We outline a model of the Crab Pulsar Wind Nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index $p_I \approx 2.2$ above the injection break corresponding to $\gamma_{wind} \sigma_{wind} \sim 10^5$, peaking in the UV ($\gamma_{wind} \sim 10^2$ is the bulk Lorentz factor of the wind, $\sigma _{wind} \sim 10^3$ is wind magnetization); (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index $p_{II} \approx 1.6$. The model requires relatively slow but highly magnetized wind. For both components the overall cooling break is in the infra-red at $\sim 0.01$ eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band Component-I produces emission with the cooling spectral index of $\alpha_o \approx 0.5$, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of $\epsilon_s \approx 100$ MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high energy astrophysical objects, like AGNe and GRB jets, where often radio electrons form a different population from the high energy electrons.

Recently it was shown that the centrifugal instability may be important in the dynamics of astrophysical jets undergoing reconfinement by external pressure. However, these studies were limited to the case of unmagnetised flows. Here we explore the role of the magnetic field within both the Newtonian and relativistic frameworks. Since the jet problem is rather complicated, we focus instead on the simpler problem of cylindrical rotation and axial magnetic field, which shares significant similarity with the jet problem, and consider only axisymmetric perturbations. The studied equilibrium configurations involve a cylindrical interface and they are stable to non-magnetic centrifugal and magneto-rotational instabilities everywhere except this interface. We use a heuristic approach to derive the local stability criterion for the interface in the magnetic case and numerical simulations to verify the role of the magnetic field. The theory and simulations agree quite well for Newtonian models but indicate a potential discrepancy for the relativistic models in the limit of high Lorentz factor of the rotational motion at the interface. In general, the magnetic field sets a critical wavelength below which the centrifugal modes are stabilised. We discuss the implication of our findings for the astrophysical jets, which suggest that the centrifugal instability develops only in jets with relatively low magnetisation. Namely, the magnetic pressure has to be below the thermal one and for the relativistic case the jets have to be kinetic-energy dominated.

Jets are a ubiquitous part of the accretion process, seen in a wide variety of objects ranging from active galaxies (AGN) to X-ray binary stars and even newly formed stars. AGN jets are accelerated by the supermassive black hole of their host galaxy by a coupling between the magnetic field and inflowing material. They are the source for many exciting phenomena and can profoundly influence the larger galaxy and surrounding cluster. This White Paper points out what advances can be achieved in the field by new technologies, concentrating on the zone where jets are accelerated to relativistic speeds and collimated. The ngVLA and new space VLBI missions will give higher angular resolution, sensitivity and fidelity in the radio, penetrating this zone for additional objects and allowing us to resolve fundamental questions over the physics of jet acceleration and collimation. Interferometry in other bands would allow us to probe directly flaring components. We also emphasize the need for polarimetry, which is essential to revealing the role and configuration of magnetic fields.

In this paper we study the reconfinement of initially freely-expanding unmagnetised relativistic jets by the pressure of non-uniform external gas using numerical approach. The results are compared with the simple semi-analytic model proposed by Komissarov & Falle (1997). In particular, we explore the reconfinement in power-law atmospheres and in the King atmosphere, which describes the X-ray coronas of giant elliptic galaxies. The results show significant deviations from the KF model, which systematically underestimates the reconfinement scale. For the power-law atmospheres the disagreement increases with the power index k, exceeding two orders of magnitude for k=1.5. For the King model, strong deviations are found on the outskirts of the atmosphere, where the distribution approaches a power law. However for jets reconfined inside the core, the reconfinement scale is increased only by the factor of two. When the King model is modified by adding a central cusp, this has little impact on the jets which are reconfined outside of the cusp region but inside the cusp the reconfinement scale significantly reduces.

Jets powered by active galactic nuclei appear impressively stable compared with their terrestrial and laboratory counterparts-they can be traced from their origin to distances exceeding their injection radius by up to a billion times. However, some less energetic jets get disrupted and lose their coherence on the scale of their host galaxy. Quite remarkably, on the same scale, these jets are expected to become confined by the thermal pressure of the intra-galactic gas. Motivated by these observations, we have started a systematic study of active galactic nuclei jets undergoing reconfinement via computer simulations. Here, we show that in the case of unmagnetized relativistic jets, the reconfinement is accompanied by the development of an instability and transition to a turbulent state. During their initial growth, the perturbations have a highly organized streamwise-oriented structure, indicating that it is not the Kelvin-Helmholtz instability, the instability which has been the main focus of the jet stability studies so far. Instead, it is closely related to the centrifugal instability. This instability is likely to be behind the division of active galactic nuclei jets into two morphological types in the Fanaroff-Riley classification.

We develop a model of gamma-ray flares of the Crab Nebula resulting from the magnetic reconnection events in highly-magnetized relativistic plasma. We first discuss physical parameters of the Crab nebula and review the theory of pulsar winds and termination shocks. We also review the principle points of particle acceleration in explosive reconnection events (Lyutikov et al. 2017a,b). It is required that particles producing flares are accelerated in highly magnetized regions of the nebula. Flares originate from the poleward regions at the base of Crab's polar outflow, where both the magnetization and the magnetic field strength are sufficiently high. The post-termination shock flow develops macroscopic (not related to the plasma properties on the skin-depth scale) kink-type instabilities. The resulting large-scales magnetic stresses drive explosive reconnection events on the light-crossing time of the reconnection region. Flares are produced at the initial stage of the current sheet development, during the X-point collapse. The model has all the ingredients needed for Crab flares: natural formation of highly magnetized regions, explosive dynamics on light travel time, development of high electric fields on macroscopic scales and acceleration of particles to energies well exceeding the average magnetic energy per particle.

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically-dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two dimensional `ABC' structures) and zero total current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For plasma magnetization $\sigma \leq 10^2$ the spectrum power law index is $p> 2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\sigma \geq 10^2$, the spectra are hard, $p< 2$, yet the maximal energy $\gamma_{max}$ can still exceed the average magnetic energy per particle, $ \sim \sigma$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

The extreme properties of the gamma ray flares in the Crab Nebula present a clear challenge to our ideas on the nature of particle acceleration in relativistic astrophysical plasma. It seems highly unlikely that standard mechanisms of stochastic type are at work here and hence the attention of theorists has switched to linear acceleration in magnetic reconnection events. In this series of papers, we attempt to develop a theory of explosive magnetic reconnection in highly-magnetized relativistic plasma which can explain the extreme parameters of the Crab flares. In the first paper, we focus on the properties of the X-point collapse. Using analytical and numerical methods (fluid and particle-in-cell simulations) we extend Syrovatsky's classical model of such collapse to the relativistic regime. We find that the collapse can lead to the reconnection rate approaching the speed of light on macroscopic scales. During the collapse, the plasma particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For sufficiently high magnetizations and vanishing guide field, the non-thermal particle spectrum consists of two components: a low-energy population with soft spectrum, that dominates the number census; and a high-energy population with hard spectrum, that possesses all the properties needed to explain the Crab flares.

Near the central engine, many astrophysical jets are expected to rotate about their axis. Further out they are expected to go through the processes of reconfinement and recollimation. In both these cases, the flow streams along a concave surface and hence, it is subject to the centrifugal force. It is well known that such flows may experience the Centrifugal Instability (CFI), to which there are many laboratory examples. The recent computer simulations of relativistic jets from Active Galactic Nuclei undergoing the process of reconfinement show that in such jets CFI may dominate over the Kelvin-Helmholtz instability associated with velocity shear (Gourgouliatos and Komissarov 2018). In this letter, we generalise the Rayleigh criterion for CFI in rotating fluids to relativistic flows using a heuristic analysis. We also present the results of computer simulations which support our analytic criterion for the case of an interface separating two uniformly-rotating cylindrical flows. We discuss the difference between CFI and the Rayleigh-Taylor instability in flows with curved streamlines.

Relativistic jets associated with active galactic nuclei and gamma-ray bursts propagate over huge distances without significant loss of momentum. At the same time they are bright emitters, which is indicative of strong energy dissipation. This points towards a mechanism of internal dissipation which does not result in a global disruption of the flow. One possibility is internal shocks and another one is turbulence driven by local instabilities. Such instabilities can be triggered when a freely expanding jet is reconfined by either the cocoon or external gas pressure. In this paper we study the dynamics of two-component spine-sheath hydrodynamic jets coming into pressure equilibrium with external gas using 2D computer simulations. We find that the jet oscillations lead to a rapid onset of Rayleigh-Taylor-type instabilities, which results in additional internal dissipation and mixing of the jet components. Although slightly different in details, this outcome holds both for the heavy-spine-light-sheath and light-spine-heavy-sheath configurations. The results may provide an explanation to the spatial flaring observed in some AGN jets on kpc-scales.

In this contribution we review the recent progress in the modeling of Pulsar Wind Nebulae (PWN). We start with a brief overview of the relevant physical processes in the magnetosphere, the wind-zone and the inflated nebula bubble. Radiative signatures and particle transport processes obtained from 3D simulations of PWN are discussed in the context of optical and X-ray observations. We then proceed to consider particle acceleration in PWN and elaborate on what can be learned about the particle acceleration from the dynamical structures called "wisps" observed in the Crab nebula. We also discuss recent observational and theoretical results of gamma-ray flares and the inner knot of the Crab nebula, which had been proposed as the emission site of the flares. We extend the discussion to GeV flares from binary systems in which the pulsar wind interacts with the stellar wind from a companion star. The chapter concludes with a discussion of solved and unsolved problems posed by PWN.

The jets from active galactic nuclei exhibit stability which seems to be far superior compared to that of terrestrial and laboratory jets. They manage to propagate over distances up to a billion of initial jet radii. Yet this may not be an indication of some exotic physics but mainly a reflection of the specific environment these jets propagate through. The key property of this environment is a rapid decline of density and pressure along the jet, which promotes its rapid expansion. Such an expansion can suppress global instabilities, which require communication across the jet, and hence ensure its survival over huge distances. At kpc scales, some AGN jets do show signs of strong instabilities and even turn into plumes. This could be a result of the flattening of the external pressure distribution in their host galaxies or inside the radio lobes. In this regard, we discuss the possible connection between the stability issue and the Fanaroff-Riley classification of extragalactic radio sources. The observations of AGN jets on sub-kpc scale do not seem to support their supposed lack of causal connectivity. When interpreted using simple kinematic models, they reveal a rather perplexing picture with more questions than answers on the jets dynamics.

In our prior papers, we considered the non-relativistic linear stability analysis of magnetized jets that do not have current sheet at the boundary. In this paper, we extend our analysis to relativistic jets. In order to find the unstable modes of current sheet-free, magnetized relativistic jets, we linearize full relativistic magnetohydrodynamics equations and solve them numerically. We find the dispersion relation of the pinch and kink mode instabilities. By comparing the dispersion relations of mildly relativistic jet (Lorentz factor 2) with moderately relativistic jet (Lorentz factor 10), we find that the jet with higher Lorentz factor is significantly more stable in both pinch and kink modes. We show that inclusion of the current sheet-free magnetic field in the jet further enhances the stability. Both pinch and kink mode instabilities become progressively more stable with increasing magnetization. We also show a scaling relation between the maximum temporal growth rate of the unstable mode and the Lorentz factor of the jet. The maximum temporal growth rates of the unstable modes are inversely proportion to the Lorentz factors for most of the modes that we study. However, for the fundamental pinch mode it is inversely proportional to the square of the Lorentz factor. This very beneficial scaling relation holds regardless of the presence of a magnetic field.

In this paper we describe a simple numerical approach which allows to study the structure of steady-state axisymmetric relativistic jets using one-dimensional time-dependent simulations. It is based on the fact that for narrow jets with v~c the steady-state equations of relativistic magnetohydrodynamics can be accurately approximated by the one-dimensional time-dependent equations after the substitution z=ct. Since only the time-dependent codes are now publicly available this is a valuable and efficient alternative to the development of a high-specialized code for the time-independent equations. The approach is also much cheaper and more robust compared to the relaxation method. We tested this technique against numerical and analytical solutions found in literature as well as solutions we obtained using the relaxation method and found it sufficiently accurate. In the process, we discovered the reason for the failure of the self-similar analytical model of the jet reconfinement in relatively flat atmospheres and elucidated the nature of radial oscillations of steady-state jets.

We develop a model of particle acceleration in explosive reconnection events in relativistic magnetically-dominated plasmas and apply it to explain gamma-ray flares from the Crab Nebula. The model relies on development of current-driven instabilities on macroscopic scales (not related to plasma skin depths). Using analytical and numerical methods (fluid and particle-in-cell simulations), we study a number of model problems in relativistic magnetically-dominated plasma: (i) we extend Syrovatsky's classical model of explosive X-point collapse to magnetically-dominated plasmas; (ii) we consider instability of two-dimensional force-free system of magnetic flux tubes; (iii) we consider merger of two zero total poloidal current magnetic flux tubes. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization. The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced reconnection, further accelerating the particles, but proceeds at a slower reconnection rate. Crab flares result from the initial explosive stages of magnetic island mergers of magnetic flux tubes produced in the bulk of nebula at intermediate polar regions. The post-termination shock plasma flow in the wind sectors with mild magnetization naturally generates large-scale highly magnetized structures. Internal kink-like instabilities lead to the formation of macroscopic current-carrying magnetic flux tubes that merge explosively.

In a prior paper (Kim et al. 2015) we considered the linear stability of magnetized jets that carry no net electric current and do not have current sheets. In this paper, in addition to physically well-motivated magnetic field structures, we also include the effects of jet shear. The jets we study have finite thermal pressure in addition to having realistic magnetic field structures and velocity shear. We find that shear has a strongly stabilizing effect on various modes of jet instability. Increasing shear stabilizes the fundamental pinch modes at long wavelengths and short wavelengths. Increasing shear also stabilizes the first reflection pinch modes at short wavelengths. Increasing shear has only a very modest stabilizing effect on the fundamental kink modes at long wavelengths; however, increasing shear does have a strong stabilizing effect on the fundamental kink modes at short wavelengths. The first reflection kink modes are strongly stabilized by increasing shear at shorter wavelengths. Overall, we find that the combined effect of magnetic field and shear stabilizes jets more than shear alone. In addition to the results from a formal linear stability analysis, we present a novel way of visualizing and understanding jet stability. This gives us a deeper understanding of the enhanced stability of sheared, magnetized jets. We also emphasize the value of our numerical approach in understanding the linear stability of jets with realistic structure.

The stability of current sheets in collisionless relativistic pair plasma was studied via two-dimensional two-fluid relativistic magnetohydrodynamic simulations with vanishing internal friction between fluids. In particular, we investigated the linear growth of the tearing and drift-kink modes in the current sheets both with and without the guide field and obtained the growth rates which are very similar to what has been found in the corresponding PIC simulations. This suggests that the two-fluid simulations can be useful in studying the large-scale dynamics of astrophysical relativistic plasmas in problems involving magnetic reconnection.

We model the inner knot of the Crab Nebula as a synchrotron emission coming from the non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock appearance. The model can reproduce the knot location, size, elongation, brightness distribution, luminosity and polarization provided the effective magnetization of the section of the pulsar wind producing the knot is low, $\sigma \leq 1$. In the striped wind model, this implies that the striped zone is rather wide, with the magnetic inclination angle of the Crab pulsar $\ge 45^\circ$; this agrees with the previous model-dependent estimate based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock and that this shock is a site of efficient particle acceleration. On the other hand, the deduced low magnetization of the knot plasma implies that this is an unlikely site for the Crab's gamma-ray flares, if they are related to the fast relativistic magnetic reconnection events.

In stark contrast to their laboratory and terrestrial counterparts, the cosmic jets appear to be very stable. We propose that the reason behind this remarkable property is the loss of causal connectivity across these jets, caused by their rapid expansion in response to fast decline of external pressure with the distance from the "jet engine". In atmospheres with power-law pressure distribution, the total loss of causal connectivity occurs, when the power index k>2 - the steepness which is expected to be quite common for many astrophysical environments. This conclusion does not seem to depend on the physical nature of jets - it applies both to relativistic and non-relativistic flows, both magnetically-dominated and unmagnetized jets. In order to verify it, we have carried out numerical simulations of moderately magnetized and moderately relativistic jets. Their results give strong support to our hypothesis and provide with valuable insights. In particular, we find that the z-pinched inner cores of magnetic jets expand slower than their envelopes and become susceptible to instabilities even when the whole jet is stable. This may result in local dissipation and emission without global disintegration of the flow. Cosmic jets may become globally unstable when they enter flat sections of external atmospheres. We propose that the Fanaroff-Riley morphological division of extragalactic radio sources into two classes is related to this issue. In particular, we argue that the low power FR-I jets become re-confined, causally connected and globally unstable on the scale of galactic X-ray coronas, whereas more powerful FR-II jets re-confine much further out and remain largely intact until they terminate at hot spots. Using this idea, we derived the relationship between the critical jet power and the optical luminosity of the host galaxy, which is in a very good agreement with the observations.

In this paper we consider stability of magnetized jets that carry no net electric current and do not have current sheets (Gourgouliatos et al. 2012). The non-relativistic MHD equations are linearized around the background velocity and the magnetic field structure of the jet. The resulting linear equations are solved numerically inside the jet. We find that introduction of current-sheet-free magnetic field significantly improves jet stability relative to unmagnetized jets or magnetized jets with current sheets at their surface. This particularly applies to the fundamental pinch and kink modes - they become completely suppressed in a wide range of long wavelengths that are known to become most pernicious to jet stability when the evolution enters the non-linear regime. The reflection modes, both for the pinch and kink instability, also become progressively more stable with increased magnetization.

In this paper we discuss the development of Rayleigh-Taylor filaments in axisymmetric simulations of Pulsar wind nebulae (PWN). High-resolution adaptive mesh refinement magnetohydrodynamic (MHD) simulations are used to resolve the non-linear evolution of the instability. The typical separation of filaments is mediated by the turbulent flow in the nebula and hierarchical growth of the filaments. The strong magnetic dissipation and field-randomization found in recent global three-dimensional simulations of PWN suggests that magnetic tension is not strong enough to suppress the growth of RT filaments, in agreement with the observations of prominent filaments in the Crab nebula. The long-term axisymmetric results presented here confirm this finding.

The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy.

In this paper we give a detailed account of the first 3D relativistic magnetohydrodynamic (MHD) simulations of Pulsar Wind Nebulae (PWN), with parameters most suitable for the Crab Nebula. In order to clarify the new features specific to 3D models, reference 2D simulations have been carried out as well. Compared to the previous 2D simulations, we considered pulsar winds with much stronger magnetisation, up to \sigma=3, and accounted more accurately for the anticipated magnetic dissipation in the striped zone of these winds. While the 3D models preserve the separation of the post termination shock flow into the equatorial and polar components, their relative strength and significance differ. Whereas the highly magnetised 2D models produce highly coherent and well collimated polar jets capable of efficient "drilling" through the supernova shell, in the corresponding 3D models the jets are disrupted by the kink mode current driven instability and "dissolve" into the main body of PWN after propagation of several shock radii. Our results show that contrary to the expectations based on 1D analytical and semi-analytical models, our numerical solutions with highly magnetized pulsar winds still produce termination shocks with radii comparable to those deduced from the observations. We present polarized synchrotron maps constructed from the 3D simulations, showing that the wealth of morphological features observed with Hubble and Chandra can well be reproduced with 3D-MHD models.

We present first results of three dimensional relativistic magnetohydrodynamical simulations of Pulsar Wind Nebulae. They show that the kink instability and magnetic dissipation inside these nebulae may be the key processes allowing to reconcile their observations with the theory of pulsar winds. In particular, the size of the termination shock, obtained in the simulations, agrees very well with the observations even for Poynting-dominated pulsar winds. Due to magnetic dissipation the total pressure in the simulated nebulae is particle-dominated and more or less uniform. While in the main body of the simulated nebulae the magnetic field becomes rather randomized, close to the termination shock, it is dominated by the regular toroidal field freshly injected by the pulsar wind. This field is responsible for driving polar outflows and may explain the high polarization observed in pulsar wind nebulae.

Magnetic dissipation is frequently invoked as a way of powering the observed emission of relativistic flows in Gamma Ray Bursts and Active Galactic Nuclei. Pulsar Wind Nebulae provide closer to home cosmic laboratories which can be used to test the hypothesis. To this end, we analyze the observational data on the spindown power of the Crab pulsar, energetics of the Crab nebula, and its magnetic field. We show that unless the magnetic inclination angle of the Crab pulsar is very close to 90 degrees the overall magnetization of the striped wind after total dissipation of its stripes is significantly higher than that deduced in the Kennel-Coroniti model and recent axisymmetric simulations of Pulsar Wind Nebulae. On the other hand, higher wind magnetization is in conflict with the observed low magnetic field of the Crab nebula, unless it is subject to efficient dissipation inside the nebula as well. For the likely inclination angle of 45 degrees the data require magnetic dissipation on the timescale about 80 years, which is short compared to the life-time of the nebula but long compared to the time scale of Crab's gamma-ray flares.

We have revisited the issue of shock dissipation and emission and its implications for the internal shock model of the prompt GRB emission and studied it in the context of impulsive Poynting-dominated flows. Our results show that unless the magnetization of GRB jets is extremely high, \sigma > 100 in the prompt emission zone, the magnetic model may still be compatible with the observations. The main effect of reduced dissipation efficiency is merely an increase in the size of the dissipation zone and even for highly magnetised GRB jets this size may remain below the external shock radius, provided the central engine can emit magnetic shells on the time scale well below the typical observed variability scale of one second. Our analytical and numerical results suggest that magnetic shells begin strongly interact with each other well before they reach the coasting radius. As the result, the impulsive jet in the dissipation zone is best described not as a collection of shells but as a continuous highly magnetised flow with a high amplitude magnetosonic wave component. How exactly the dissipated wave energy is distributed between the radiation and the bulk kinetic energy of radial jets depends on the relative rates of radiative and adiabatic cooling. In the fast radiative cooling regime, the corresponding radiative efficiency can be as high as the wave contribution to their energy budget, independently of the magnetization. Moreover, after leaving the zone of prompt emission the jet may still remain Poynting-dominated, leading to weaker emission from the reverse shock compared to non-magnetic models.

The Magnetodynamics, or Force-Free Degenerate Electrodynamics, is recognized as a very useful approximation in studies of magnetospheres of relativistic stars. In this paper we discuss various forms of Magnetodynamic equations which can be used to study magnetospheres of black holes. In particular, we focus on the 3+1 equations which allow for curved and dynamic spacetime.

In this paper we propose a new plausible mechanism of supernova explosions specific to close binary systems. The starting point is the common envelope phase in the evolution of a binary consisting of a red super giant and a neutron star. As the neutron star spirals towards the center of its companion it spins up via disk accretion. Depending on the specific angular momentum of gas captured by the neutron star via the Bondi-Hoyle mechanism, it may reach millisecond periods either when it is still inside the common envelope or after it has merged with the companion core. The high accretion rate may result in strong differential rotation of the neutron star and generation of the magnetar-strength magnetic field. The magnetar wind can blow away the common envelope if its magnetic field is as strong as $10^{15}\,$G, and can destroy the entire companion if it is as strong as $10^{16}\,$G. The total explosion energy can be comparable to the rotational energy of a millisecond pulsar and reach $10^{52}\,$erg. However, only a small amount of $^{56}$Ni is expected to be produced this way. The result is an unusual type-II supernova with very high luminosity during the plateau phase, followed by a sharp drop in brightness and a steep light-curve tail. The remnant is either a solitary magnetar or a close binary involving a Wolf-Rayet star and a magnetar. When this Wolf-Rayet star explodes this will be a third supernovae explosion in the same binary.

The oblique geometry of pulsar wind termination shock ensures that the Doppler beaming has a strong impact on the shock emission. We illustrate this using recent relativistic MHD simulations of the Crab Nebula and also show that the observed size, shape, and distance from the pulsar of the Crab Nebula inner knot are consistent with its interpretation as a Doppler-boosted emission from the termination shock. If the electrons responsible for the synchrotron gamma-rays are accelerated only at the termination shock then their short life-time ensures that these gamma-rays originate close to the shock and are also strongly effected by the Doppler beaming. As the result, bulk of the observed synchrotron gamma-rays of the Crab Nebula around 100 MeV may come from its inner knot. This hypothesis is consistent with the observed optical flux of the inner knot provided its optical-gamma spectral index is the same as the injection spectral index found in the Kennel & Coroniti model of the nebula spectrum. The observed variability of synchrotron gamma-ray emission can be caused by the instability of the termination shock discovered in recent numerical simulations. Given the small size of the knot, it is possible that the September 2010 gamma-ray flare of the Crab Nebula also came from the knot, though the actual mechanism remains unclear. The model predicts correlation of the temporal variability of the synchrotron gamma-ray flux in the Fermi and AGILE windows with the variability of the unpulsed optical flux from within 1 arcsec of the Crab pulsar.

The strong variability of magnetic central engines of AGN and GRBs may result in highly intermittent strongly magnetized relativistic outflows. We find a new magnetic acceleration mechanism for such impulsive flows that can be much more effective than the acceleration of steady-state flows. This impulsive acceleration results in kinetic-energy-dominated flows at astrophysically relevant distances from the central source. For a spherical flow, a discrete shell ejected from the source over a time t_0 with Lorentz factor Gamma~1 and initial magnetization sigma_0 = B_0^2/(4 pi rho_0 c^2) >> 1 quickly reaches a typical Lorentz factor Gamma ~ sigma_0^{1/3} and magnetization sigma ~ sigma_0^{2/3} at the distance R_0 ~ ct_0. At this point the magnetized shell of width Delta ~ R_0 in the lab frame loses causal contact with the source and continues to accelerate by spreading significantly in its own rest frame. The expansion is driven by the magnetic pressure gradient and leads to relativistic relative velocities between the front and back of the shell. While the expansion is roughly symmetric in the center of momentum frame, in the lab frame most of the energy and momentum remain in a region (or shell) of width Delta ~ R_0 at the head of the flow. This acceleration proceeds as Gamma ~ (sigma_0 R/R_0)^{1/3} and sigma ~ sigma_0^{2/3} (R/R_0)^{-1/3} until reaching a coasting radius R_c ~ R_0 sigma_0^2 where the kinetic energy becomes dominant: Gamma ~ sigma_0 and sigma ~ 1 at R_c. Then the shell starts coasting and spreading (radially) causing its magnetization to drop as sigma ~ R_c/R at R>R_c. Given the typical variability time-scales of AGN and GRBs, the magnetic acceleration in these sources is a combination of the quasi-steady-state collimation acceleration close to the source and the impulsive acceleration further out. The interaction with the external medium is briefly discussed.

This is a brief review of the recent developments in the theory of magnetic acceleration of relativistic jets. We attempt to explain the key results of this complex theory using basic physical arguments and simple calculations. The main focus is on the standard model, which describes steady-state axisymmetric ideal MHD flows. We argue that this model is over-restrictive and discuss various alternatives.

When a magnetically-dominated super-fast magnetosonic GRB jet leaves the progenitor star the external pressure support may drop and the jet may enter the regime of ballistic expansion during which its magnetic acceleration becomes highly ineffective. However, recent numerical simulations suggested that the transition to this regime is accompanied by a sudden "burst" of acceleration. We confirm this finding and attribute the acceleration to the sideways expansion of the jet - the magnetic energy is converted into the kinetic one in the strong magnetosonic rarefaction wave, which is launched when the jet loses its external support. This type of acceleration, the rarefaction acceleration, is specific to relativistic jets because their energy budget can still be dominated by magnetic energy even in highly super-fast magnetosonic regime. Just like the collimation acceleration of externally confined magnetized jets, it is connected with the geometry of magnetic flux sufaces. In both cases, in the acceleration zone the poloidal field lines diverge faster than in the monopolar configuration. On the other hand, whereas the collimation acceleration keeps the product of jet opening angle and Lorentz factor somewhat below unity, the rarefaction acceleration allows to make it significantly larger, in agreement with the standard model of jet breaks in afterglow light curves.

The very first stars in the Universe can be very massive, up to $10^3M_\odot$. If born in large numbers, such massive stars can have strong impact on the subsequent star formation, producing strong ionising radiation and contaminating the primordial gas with heavy elements. They would leave behind massive black holes that could act as seeds for growing supermassive black holes of active galactic nuclei. Given the anticipated fast rotation, such stars would end their live as supermassive collapsars and drive powerful magnetically-dominated jets. In this letter we investigate the possibility of observing the bursts of high-energy emission similar to the Long Gamma Ray Bursts associated with normal collapsars. We show that during the collapse of supercollapsars, the Blandford-Znajek mechanism can produce jets as powerful as few$\times10^{52}$erg/s and release up to $10^{54}$erg of the black hole rotational energy. Due to the higher intrinsic time scale and higher redshift the initial bright phase of the burst can last for about $10^4$ seconds, whereas the central engine would remain active for about one day. Due to the high redshift the burst spectrum is expected to be soft, with the spectral energy distribution peaking at around 20-30keV. The peak total flux density is relatively low, $10^{-7}{erg} {cm}^{-2} {s}^{-1}$, but not prohibitive. If one supercollapsar is produced per every mini-halo of dark matter arising from the 3-$\sigma$ cosmological fluctuations then the whole sky frequency of such bursts could reach several tens per year.

In this letter, we describe results of new high-resolution axisymmetric relativistic MHD simulations of Pulsar Wind Nebulae. The simulations reveal strong breakdown of the equatorial symmetry and highly variable structure of the pulsar wind termination shock. The synthetic synchrotron maps, constructed using a new more accurate approach, show striking similarity with the well known images of the Crab Nebula obtained by Chandra, and the Hubble Space Telescope. In addition to the jet-torus structure, these maps reproduce the Crab's famous moving wisps whose speed and rateof production agree with the observations. The variability is then analyzed using various statistical methods, including the method of structure function and wavelet transform. The results point towards the quasi-periodic behaviour with the periods of 1.5-3yr and MHD turbulence on scales below 1yr. The full account of this study will be presented in a follow up paper.

The strong dependence of the neutrino annihilation mechanism on the mass accretion rate makes it difficult to explain the LGRBs with duration in excess of 100 seconds as well as the precursors separated from the main gamma-ray pulse by few hundreds of seconds. Even more difficult is to explain the Swift observations of the shallow decay phase and X-ray flares, if they indeed indicate activity of the central engine for as long as 10,000 seconds. These data suggest that some other, most likely magnetic mechanisms have to be considered. The magnetic models do not require the development of accretion disk within the first few seconds of the stellar collapse and hence do not require very rapidly rotating stellar cores at the pre-supernova state. This widens the range of potential LGRB progenitors. In this paper, we re-examine the close binary scenario allowing for the possibility of late development of accretion disks in the collapsar model and investigate the available range of mass accretion rates, black hole masses, and spins. A particularly interesting version of the binary progenitor involves merger of a WR star with an ultra-compact companion, neutron star or black hole. In this case we expect the formation of very long-lived accretion disks, that may explain the phase of shallow decay and X-ray flares observed by Swift. Similarly long-lived magnetic central engines are expected in the current single star models of LGRB progenitors due to their assumed exceptionally fast rotation.

Collapse of massive stars may result in formation of accreting black holes in their interior. The accreting stellar matter may advect substantial magnetic flux onto the black hole and promote 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 particularly, we show that the Blandford-Znajek mechanism is activated when the rest mass-energy density of matter drops below the energy density of magnetic field in the very 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 stellar core required in the collapsar model and suggest that the conflict can be avoided if the progenitor star is a component of 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 disk 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 additional power provided by the disk.

We present a relativistic-MHD numerical study of axisymmetric, magnetically driven jets with parameters applicable to gamma-ray burst (GRB) flows. We also present analytic expressions for the asymptotic jet shape and other flow parameters that agree very well with the numerical results. All current-carrying outflows exhibit self-collimation and consequent acceleration near the rotation axis, but unconfined outflows lose causal connectivity across the jet and therefore do not collimate or accelerate efficiently in their outer regions. Magnetically accelerated jets confined by an external pressure that varies with distance with a power-law index < 2 assume a paraboloidal shape and have an acceleration efficiency > 50%. They attain Lorentz factors > 30 on scales 10^9-3x10^10 cm, consistent with the possibility that short/hard GRB jets are accelerated on scales where they can be confined by moderately relativistic winds from accretion discs, and > 100 on scales 10^10-10^12 cm, consistent with the possibility that long/soft GRB jets are accelerated within the envelopes of collapsing massive stars. We also find that the Lorentz factor of a magnetically accelerated jet is approximately inversely proportional to the opening half-angle of the poloidal streamlines. This implies that the gamma-ray emitting components of GRB outflows are very narrow, with a half-angle < 1 degree in regions where the Lorentz factor exceeds 100, and that the afterglow light curves of these components would either exhibit a very early jet break or show no jet break at all.

In this study we explore the magnetic mechanism of hypernovae and relativistic jets of long duration gamma ray bursts within the collapsar scenario. This is an extension of our earlier work [1]. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state and takes into account the cooling associated with emission of neutrinos and the energy losses due to dissociation of nuclei. The neutrino heating is not included. We describe solutions with different black hole rotation, mass accretion rate, and strength of progenitor's magnetic field. Some of them exhibits strong explosions driven by Poynting-dominated jets with power up to $12\times10^{51} {erg s}^{-1}$. These jets originate from the black hole and powered via the Blandford-Znajek mechanism. A provisional criterion for explosion is derived. A number of simulation movies can be downloaded from http://www.maths.leeds.ac.uk/~serguei/research/movies/anim.html

We present numerical simulations of cold, axisymmetric, magnetically driven relativistic outflows. The outflows are initially sub-Alfv\'enic and Poynting flux-dominated, with total--to--rest-mass energy flux ratio up to $\mu \sim 620$. To study the magnetic acceleration of jets we simulate flows confined within a funnel with rigid wall of prescribed shape, which we take to be $z\propto r^a$ (in cylindrical coordinates, with $a$ ranging from 1 to 2). This allows us to eliminate the numerical dissipative effects induced by a free boundary with an ambient medium. We find that in all cases they converge to a steady state characterized by a spatially extended acceleration region. For the jet solutions the acceleration process is very efficient - on the outermost scale of the simulation more than half of the Poynting flux has been converted into kinetic energy flux, and the terminal Lorentz factor approached its maximum possible value ($\Gamma_\infty \simeq \mu$). The acceleration is accompanied by the collimation of magnetic field lines in excess of that dictated by the funnel shape. The numerical solutions are generally consistent with the semi-analytic self-similar jets solutions and the spatially extended acceleration observed in some astrophysical relativistic jets. In agreement with previous studies we also find that the acceleration is significantly less effective for wind solutions suggesting that pulsar winds may remain Poynting dominated when they reach the termination shock.

During the three decades since its theoretical discovery the Blandford-Znajek process of extracting the rotational energy of black holes has become one of the foundation stones in the building of modern relativistic astrophysics. However, it is also true that for a long time its physics was not well understood, as evidenced by the controversy that surrounded it since 1990s. Thanks to the efforts of many theorists during the last decade the state of affairs is gradually improving. In this lecture I attempt to explain the key ingredients of this process in more or less systematic, rigorous, and at the same time relatively simple fashion. A particular attention is paid to the similarities and differences between the Blandford-Znajek and Penrose processes. To this purpose I formulate the notion of energy counter flow. The concept of horizon membrane is replaced with the concept of vacuum as an electromagnetically active medium. The effect of negative phase velocity of electromagnetic waves in the black hole ergosphere is also discussed.

In this letter we briefly describe the first results of our numerical study on the possibility of magnetic origin of relativistic jets of long duration gamma ray bursters within the collapsar scenario. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state of stellar matter, takes into account the cooling associated with emission of neutrinos, and the energy losses due to dissociation of nuclei. The neutrino heating is not included. We describe the solution for one particular model where the progenitor star has magnetic field $B=10^{10}$G. The solution exhibits strong explosion driven by the Poynting-dominated jets whose power exceeds $2\times10^{51} {erg/s}$. The jets originate mainly from the black hole and they are powered via the Blandford-Znajek mechanism. The full details of the simulations together with the results of parameter study will be presented elsewhere. A number of simulation movies can be downloaded from http://www.maths.leeds.ac.uk/~serguei/research/movies/anim.html

In this paper we report on the early evolution of core-collapse supernova explosion following the birth of a magnetar with the dipolar magnetic field of B=10^{15}G and the rotational period of 2ms, which was studied by means of axisymmetric general relativistic MHD simulations. The numerical models exhibit highly collimated magnetically-driven jets very early on. The jets are super-Alfvenic but remain sub-fast until the end of the simulations (t=0.2s). The power released in the jets is about 3x10^{50}erg/s which implies the spin-down time of ~37s. The total rotational energy of the magnetar, E~10^{52}erg, is sufficient to drive hypernova but it is not clear as to how large a fraction of this energy can be transfered to the stellar envelope. Given the observed propagation speed of the jets, v_p~0.17c, they are expected to traverse the progenitor in few seconds and after this most of the released rotational energy would be simply carried away by these jets into the surrounding space. Our results provide the first more or less self-consistent numerical model of a central engine capable of producing, in the supernova setting and on a long-term basis, collimated jets with sufficient power to explain long duration GRBs and their afterglows. Although the flow speed of our jets is relatively low, v_j~0.5c$, the cooling of proto-neutron star will eventually result in much higher magnetization of its magnetospheres and ultra-relativistic asymptotic speeds of the jets. Given the relatively long cooling time-scale we still expect the jets to be only weakly relativistic by the time of break out. This leads to a model of GRB jets with systematic longitudinal variation of Lorentz factor which may have specific observational signatures both in the prompt and the afterglow emission.

The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric field consistent with the electric charge distribution via the method of Generalized Lagrange Multiplier. The hyperbolic fluxes are computed using the HLL prescription and the source terms are accounted via the time-splitting technique. The results of test simulations show that the scheme can handle equally well both resistive current sheets and shock waves and thus can be a useful tool for studying phenomena of relativistic astrophysics that involve both colliding supersonic flows and magnetic reconnection.

We present numerical simulations of axisymmetric, magnetically driven relativistic jets. To eliminate the dissipative effects induced by a free boundary with an ambient medium we assume that the flow is confined by a rigid wall of a prescribed shape, which we take to be $z\propto r^a$ (in cylindrical coordinates, with $a$ ranging from 1 to 3). The outflows are initially cold, sub-Alfv\'enic and Poynting flux-dominated, with a total--to--rest-mass energy flux ratio $\mu \sim 15$. We find that in all cases they converge to a steady state characterized by a spatially extended acceleration region. The acceleration process is very efficient: on the outermost scale of the simulation as much as $\sim 77%$ of the Poynting flux has been converted into kinetic energy flux, and the terminal Lorentz factor approaches its maximum possible value ($\Gamma_\infty \simeq \mu$). We also find a high collimation efficiency: all our simulated jets develop a cylindrical core. We argue that this could be the rule for current-carrying outflows that start with a low initial Lorentz factor ($\Gamma_0 \sim 1$). Our conclusions on the high acceleration and collimation efficiencies are not sensitive to the particular shape of the confining boundary or to the details of the injected current distribution, and they are qualitatively consistent with the semi-analytic self-similar solutions derived by Vlahakis & K\"onigl. We apply our results to the interpretation of relativistic jets in AGNs: we argue that they naturally account for the spatially extended accelerations inferred in these sources ($\Gamma_\infty \ga 10$ attained on radial scales $R\ga 10^{17} {\rm cm}$) and are consistent with the transition to the matter-dominated regime occurring already at $R\ga 10^{16} {\rm cm}$.

The expulsion of axisymmetric magnetic field from the event horizons of rapidly rotating black holes has been seen as an astrophysically important effect that may significantly reduce or even nullify the efficiency of the Blandford-Znajek mechanism of powering the relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts. However, this Meissner-like effect is seen in vacuum solutions of black hole electrodynamics whereas the Blandford-Znajek mechanism is concerned with plasma-filled magnetospheres. In this paper we argue that conductivity dramatically changes the properties of axisymmetric electromagnetic solutions -- even for a maximally rotating Kerr black hole the magnetic field is pulled inside the event horizon. Moreover, the conditions resulting in outgoing Poynting flux in the Blandford-Znajek mechanism exist not on the event horizon but everywhere within the black hole ergosphere. Thus, the Meissner effect is unlikely to be of interest in astrophysics of black holes, at least not in the way this has been suggested so far. These conclusions are supported by the results of time-dependent numerical simulations with three different computer codes. The test problems involve black holes with the rotation parameter ranging from $a=0.999$ to $a=1$. The pure electrodynamic simulations deal with the structure of conductive magnetospheres of black holes placed in a uniform-at-infinity magnetic field (Wald's problem) and the magnetohydrodynamic simulations are used to study the magnetospheres arising in the problem of disk accretion.

Many astrophysical sources of high energy emission, such as black hole magnetospheres, superstrongly magnetized neutron stars (magnetars), and probably relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts involve relativistically magnetically dominated plasma. In such plasma the energy density of magnetic field greatly exceeds the thermal and the rest mass energy density of particles. Therefore the magnetic field is the main reservoir of energy and its dissipation may power the bursting emission from these sources, in close analogy to Solar flares. One of the principal dissipative instabilities that may lead to release of magnetic energy is the tearing instability. In this paper we study, both analytically and numerically, the development of tearing instability in relativistically magnetically-dominated plasma using the framework of resistive magnetodynamics. We confirm and elucidate the previously obtained result on the growth rate of the tearing mode: the shortest growth time is the same as in the case of classical non-relativistic MHD, namely $\tau =\sqrt{\tau_a \tau_d}$ where $\tau_a$ is the \Alfven crossing time and $\tau_d$ is the resistive time of a current layer.

The dynamics of accretion discs around galactic and extragalactic black holes may be influenced by their magnetic field. In this paper we generalise the fully relativistic theory of stationary axisymmetric tori in Kerr metric of Abramowicz et al.(1978) by including strong toroidal magnetic field and construct analytic solutions for barotropic tori with constant angular momentum. This development is particularly important for the general relativistic computational magnetohydrodynamics that suffers from the lack of exact analytic solutions that are needed to test computer codes.

In this paper we present the results of time-dependent simulations of dipolar axisymmetric magnetospheres of neutron stars carried out both within the framework of relativistic magnetohydrodynamics and within the framework of 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 over the robustness of this approach. On the other hand, the MHD approach allows to make a significant progress. We start with a nonrotating magnetically dominated dipolar magnetospheres and follow its evolution as the stellar rotation is switched on. We find that the time-dependent solution gradually approaches the steady state that is very close to the stationary solution of the Pulsar Equation found by Contopoulos et al.(1999). This result suggests that other stationary solutions that have the y-point located well inside the light cylinder are unstable. The role of the particle inertia and pressure on the structure and dynamics of MHD magnetospheres is studied in details, as well as the potential implications of the 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 sub-pulse phenomena.

In this paper we report the results of axisymmetric relativistic MHD simulations for the problem of Kerr black hole immersed into a rarefied plasma with ''uniform'' magnetic field. The long term solution shows properties which are significantly different from those of the initial transient phase studied recently by Koide(2003). The topology of magnetic field lines within the ergosphere is similar to that of the split-monopole model with a strong current sheet in the equatorial plane. Closer inspection reveals a system of isolated magnetic islands inside the sheet and ongoing magnetic reconnection. No regions of negative hydrodynamic ''energy at infinity'' are seen inside the ergosphere and the so-called MHD Penrose process does not operate. Yet, the rotational energy of the black hole continues to be extracted via purely electromagnetic mechanism of Blandford and Znajek(1977). However, this is not followed by development of strong relativistic outflows from the black hole. Combined with results of other recent simulations this signals a potential problem for the standard MHD model of relativistic astrophysical jets should they still be observed at distances as small as few tens of gravitational radii from the central black hole.

Numerical simulations combined with simple analytical arguments are used to reach a number of important conclusions on the nature of the Blandford-Znajek mechanism. We show that, just like in the Penrose mechanism and in the MHD models of Punsly and Coroniti, the key role in this mechanism is played by the black hole ergosphere. The poloidal currents are driven by the gravitationally induced electric field which cannot be screened within the ergosphere by any static distribution of the electric charge of locally created pair plasma. Contrary to what is expected in the Membrane paradigm, the energy and angular momentum are extracted not only along the magnetic field lines penetrating the event horizon but along all field lines penetrating the ergosphere. In dipolar magnetic configurations symmetric relative to the equatorial plane the force-free approximation breaks down within the ergosphere where a strong current sheet develops along the equatorial plane. This current sheet supplies energy and angular momentum at infinity to the surrounding force-free magnetosphere. The Blandford-Znajek monopole solution is found to be asymptotically stable and causal. The so-called horizon boundary condition of Znajek is shown to be a regularity condition at fast critical surface.

In this paper we report the results of the first ever time-dependent general relativistic magnetohydrodynamic simulations of the magnetically dominated monopole magnetospheres of black holes. It is found that the numerical solution evolves towards a stable steady-state solution which is very close to the corresponding force-free solution found by Blandford and Znajek. Contrary to the recent claims, the particle inertia does not become dynamically important near the event horizon and the force-free approximation provides a proper framework for magnetically dominated magnetospheres of black holes. For the first time, our numerical simulations show the development of an ultra-relativistic particle wind from a rotating black hole. However, the flow remains Poynting dominated all the way up to the fast critical point. This suggests that the details of the so-called ``astrophysical load'', where the electromagnetic energy is transferred to particles, may have no effect on the efficiency of the Blandford-Znajek mechanism.

Recent discoveries of the intriguing ``jet-torus'' structure in the Crab Nebula and other pulsar nebulae prompted calls for re-examining of their theory. The most radical proposals involve abolishing of the MHD approximation altogether and developing of purely electromagnetic models. However, the classical MHD models of the Crab Nebula were hampered by the assumption of spherical symmetry made in order to render the flow equations easily integrable. The impressive progress in computational relativistic magnetohydrodynamics in recent years has made it possible to study the Crab nebula via numerical simulations without making such a drastic simplification of the problem. In this letter we present the results of the first study of such kind. They show that the jet-torus pattern can be explained within MHD approximation when anisotropy of pulsar winds is taken into account. They also indicate that the flow in the nebula is likely to be much more intricate than it has been widely believed.

It is widely accepted in the astrophysical community that the event horizon plays crucial role in the Blandford-Znajek mechanism of extraction of rotational energy of black holes. In fact, this view is a quintessence of the Membrane Paradigm of black holes which suggests that the event horizon, or rather the so-called stretched horizon, is similar to a rotating conducting sphere of finite resistivity. In this paper we argue that this interpretation is rather misleading and the proper explanation of the Blandford-Znajek mechanism has to be sought in the properties of the ergospheric region of black holes.

In this paper we examine various issues closely related to the ongoing discussion on the nature of the Blandford-Znajek mechanism of extraction of rotational energy of black holes. In particular, we show that switch-on and switch-off shocks are allowed by the shock equations of relativistic MHD and have similar properties to their Newtonian counterparts. Just like in Newtonian MHD they are limits of fast and slow shock solutions and as such they may be classified as weakly evolutionary shocks. The analysis of Punsly's MHD waveguide problem shows that its solution cannot have the form of a traveling step wave and that both fast and Alfv\'en waves are essential for generating the flow in the guide. Causality considerations are used to argue that the Blandford-Znajek perturbative solution is in conflict with the membrane paradigm. An alternative interpretation is presented according to which the role of an effective unipolar inductor in the Blandford-Znajek mechanism is played by the ergospheric region of a rotating black hole. Various implications of this are discussed.

This paper formulates time-dependent, force-free, degenerate electrodynamics as a hyperbolic system of conservation laws. It is shown that this system has four characteristic modes, a pair of fast waves propagating with the speed of light and a pair of Alfven waves. All these modes are linearly degenerate. The results of this analytic study can be used in developing upwind numerical schemes for the electrodynamics of black holes and pulsar magnetospheres. As an example, this paper describes a simple one-dimensional numerical scheme based on linear and exact Riemann solvers.

Relativistic shocks accelerate particles by the first order Fermi mechanism. These particles then emit synchrotron emission in the post shock gas. We have developed a numerical code which integrates the relativistic Euler equations for fluid dynamics with a general equation of state, together with the Liouville equation for the accelerated particles. We present tests of this code and, in addition, we use it to study the gamma ray burst afterglow predicted by the fireball model, along with the hydrodynamics of a relativistic blastwave. We find that, while, broadly speaking, the behaviour of the emission is similar to that already predicted with semi-analytic approaches, the detailed behaviour is somewhat different. The ``breaks'' in the synchrotron spectrum behave differently with time, and the spectrum above the final break is harder than previously expected. These effects are due to the incorporation of the geometry of the (spherical) blastwave, along with relativistic beaming and adiabatic cooling of the energetic particles leading to a mix, in the observed spectrum, between recently injected "uncooled" particles and the older "cooled" population in different parts of the evolving, inhomogeneous flow.