The nearly 100% linear polarization has been reported for a few fast radio bursts. This finding places severe limits on the emission mechanism. I argue that the totally polarized radiation could be generated in the course of relativistic magnetic reconnection in the outer magnetosphere of the magnetar. At the onset of the magnetar flare, a large scale magnetic perturbation forms a magnetic pulse, which propagates from the flare cite outwards. The pulse strongly compresses magnetospheric plasma and pushes it away. The high-frequency MHD waves are generated when the magnetic pulse reaches the current sheet separating, just beyond the light cylinder, the oppositely directed magnetic fields. Coalescence of magnetic islands in the reconnecting current sheet produces magnetosonic waves, which propagate away on the top of the magnetic pulse and escape in the far zone of the wind as radio waves polarized along the rotational axis of the magnetar. I estimate parameters of the outgoing radiation and show that they are compatible with the observed properties of FRBs.

Spectral modification of energetic magnetar flares by resonant cyclotron scattering (RCS) is considered. During energetic flares, photons emitted from the magnetically-trapped fireball near the stellar surface should resonantly interact with magnetospheric electrons or positrons. We show by a simple thought experiment that such scattering particles are expected to move at mildly relativistic speeds along closed magnetic field lines, which would slightly shift the incident photon energy due to the Doppler effect. We develop a toy model for the spectral modification by a single RCS that incorporates both a realistic seed photon spectrum from the trapped fireball and the velocity field of particles, which is unique to the flaring magnetosphere. We show that our spectral model can be effectively characterized by a single parameter; the effective temperature of the fireball, which enables us to fit observed spectra with low computational cost. We demonstrate that our single scattering model is in remarkable agreement with Swift/BAT data of intermediate flares from SGR 1900+14, corresponding to effective fireball temperatures of $T_{\rm eff}=6$-$7$ keV, whereas BeppoSAX/GRBM data of giant flares from the same source may need more elaborate models including the effect of multiple scatterings. Nevertheless, since there is no standard physically-motivated model for magnetar flare spectra, our model could be a useful tool to study magnetar bursts, shedding light on the hidden properties of the flaring magnetosphere.

The dissipation of turbulent magnetic fields is an appealing scenario to explain the origin of non-thermal particles in high-energy astrophysical sources. However, it has been suggested that the particle distribution may effectively thermalise when the radiative (synchrotron and/or Inverse Compton) losses are severe. Inspired by recent PIC simulations of relativistic turbulence, which show that electrons are impulsively heated in intermittent current sheets by a strong electric field aligned with the local magnetic field, we instead argue that in plasmas where the particle number density is dominated by the pairs (electron-positron and electron-positron-ion plasmas): (i) as an effect of fast cooling and of different injection times, the electron energy distribution is ${\rm d}n_e/{\rm d}\gamma\propto\gamma^{-2}$ for $\gamma\lesssim\gamma_{\rm heat}$ (the Lorentz factor $\gamma_{\rm heat}$ being close to the equipartition value), while the distribution steepens at higher energies; (ii) since the time scales for the turbulent fields to decay and for the photons to escape are of the same order, the magnetic and the radiation energy densities in the dissipation region are comparable; (iii) if the mass energy of the plasma is dominated by the ion component, the pairs with a Lorentz factor smaller than a critical one (of the order of the proton-to-electron mass ratio) become isotropic, while the pitch angle remains small otherwise. The outlined scenario is consistent with the typical conditions required to reproduce the Spectral Energy Distribution of blazars, and allows one to estimate the magnetisation of the emission site. Finally, we show that turbulence within the Crab Nebula may power the observed gamma-ray flares if the pulsar wind is nearly charge-separated at high latitudes.

Theoretical modelling and observations of AGN jets suggest that the non-thermal electrons emitting the observed radiation should (i) carry an amount of energy comparable to the magnetic fields ($U_e\sim U_B$), which is likely the case if the magnetic fields play a dynamically important role in the jet's acceleration process; (ii) cool efficiently in a dynamical time ($t_{\rm cool}\lesssim t_{\rm dyn}$), which is suggested by the fact that a large fraction of the jet's kinetic energy is promptly converted into radiation. These expectations are at odds with the results of the simplest one-zone Self-Synchro-Compton (SSC) model for the Spectral Energy Distribution (SED) of BL Lacs. Indeed, the model predicts $U_e\gg U_B$ and $t_{\rm cool}\gg t_{\rm dyn}$ for most of the objects. Here we closely investigate one of the key assumptions of this model, namely that the momentum distribution of the non-thermal electrons is isotropic. We find that this assumption may be an oversimplification. If the magnetic energy is dissipated via a turbulent MHD cascade, the highest energy electrons may instead retain a small pitch angle. Since the synchrotron emissivity is suppressed when the pitch angle is small, this effect may importantly affect the modelling of the SED. As an illustrative example, we present an anisotropic model for the electron momentum distribution such that $U_e\sim U_B$ and $t_{\rm cool}\lesssim t_{\rm dyn}$ at the same time. Our model manages to simultaneously solve the two problems with one only more free parameter with respect to the usual isotropic one-zone SSC model.

The pulsar radio emission is commonly associated with the plasma outflow in the open field line tube; then a pencil beam is emitted along the pulsar magnetic axis. Observations suggest that there is an additional radio emission mechanism specific for pulsars with high magnetic field at the light cylinder. These pulsars are known to be strong sources of non-thermal high energy radiation, which could be attributed to reconnection in the current sheet separating, just beyond the light cylinder, the oppositely directed magnetic fields. Pulsars with the highest magnetic field at the light cylinder (>100 kG) exhibit also radio pulses in phase with the high energy pulses. Moreover, giant radio pulses are observed in these pulsars. I argue that the reconnection process that produces high energy emission could also be responsible for the radio emission. Namely, coalescence of magnetic islands in the sheet produces magnetic perturbations that propagate away in the form of electro-magnetic nano-shots. I estimate the parameters of this emission and show that they are compatible with observations.

In the paradigm of magnetic acceleration of relativistic outflows, a crucial point is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma, and eventually into the observed radiation. Since the plasma is hardly accelerated beyond equipartition, MHD instabilities are often invoked to explain the dissipation of the magnetic energy. Motivated by the fast variability that is shown by the gamma-ray flares of both AGN and PWNe, different authors have proposed the Poynting flux to be dissipated in a region where the flow is converging. Here we perform a linear stability analysis of ultra-relativistic, highly magnetised outflows with such a recollimation nozzle, showing that MHD instabilities are indeed induced by the convergence of the flow. The amplitude of the perturbations increases while recollimation gets stronger, and eventually diverges when the flow is focused to a single point. Hence, depending on the geometry of the outflow, instabilities excited while the flow is converging may play an important role to dissipate the magnetic energy of the plasma.

This paper is the first in the series of papers aiming to study interaction of the electro-magnetic precursor waves generated at the front of a relativistic shock with the upstream flow. It is motivated by a simple consideration showing that the absorption of such an electro-magnetic precursor could yield an efficient transformation of the kinetic energy of the upstream flow to the energy of accelerated particles. Taking into account that the precursor is a strong wave, in which electrons oscillate with relativistic velocities, the standard plasma-radiation interaction processes should be reconsidered. In this paper, I calculate the synchrotron absorption of strong electro-magnetic waves.

We study the linear and non-linear development of the Kruskal-Schwarzchild Instability in a relativisitically expanding striped wind. This instability is the generalization of Rayleigh-Taylor instability in the presence of a magnetic field. It has been suggested to produce a self-sustained acceleration mechanism in strongly magnetized outflows found in active galactic nuclei, gamma-ray bursts, and micro-quasars. The instability leads to magnetic reconnection, but in contrast with steady-state Sweet-Parker reconnection, the dissipation rate is not limited by the current layer's small aspect ratio. We performed two-dimensional (2D) relativistic magneto-hydrodynamic (RMHD) simulations featuring two cold and highly magnetized ($1\leq\sigma\leq10^{3}$) plasma layers with an anti-parallel magnetic field separated by a thin layer of relativistically hot plasma with a local effective gravity induced by the outflow's acceleration. Our simulations show how the heavier relativistically hot plasma in the reconnecting layer drips out and allows oppositely oriented magnetic field lines to reconnect. The instability's growth rate in the linear regime matches the predictions of linear stability analysis. We find turbulence rather than an ordered bulk flow near the reconnection region, with turbulent velocities up to $\sim0.1$c, largely independent of model parameters. However, the magnetic energy dissipation rate is found to be much slower, corresponding to an effective ordered bulk velocity inflow into the reconnection region $v_{\rm in}=\beta_{\rm in}c$, of $10^{-3}\lesssim\beta_{\rm in}\lesssim 5\times10^{-3}$. This occurs due to the slow evacuation of hot plasma from the current layer, largely because of the Kelvin-Helmholtz instability experienced by the dripping plasma. 3D RMHD simulations are needed to further investigate the non-linear regime.

In the paradigm of magnetic launching of astrophysical jets, instabilities in the MHD flow are a good candidate to convert the Poynting flux into the kinetic energy of the plasma. If the magnetised plasma fills the almost entire space, the jet is unstable to helical perturbations of its body. However, the growth rate of these modes is suppressed when the poloidal component of the magnetic field has a vanishing gradient, which may be the actual case for a realistic configuration. Here we show that, if the magnetised plasma is confined into a limited region by the pressure of some external medium, the velocity shear at the contact surface excites unstable modes which can affect a significant fraction of the jet's body. We find that when the Lorentz factor of the jet is $\Gamma\sim10$ ($\Gamma\sim 100$), these perturbations typically develop after propagating along the jet for tens (hundreds) of jet's radii. Surface modes may therefore play an important role in converting the energy of the jet from the Poynting flux to the kinetic energy of the plasma, particularly in AGN. The scaling of the dispersion relation with (i) the angular velocity of the field lines and (ii) the sound speed in the confining gas is discussed.

In the paradigm of magnetic acceleration of relativistic jets, one of the key points is identifying a viable mechanism to convert the Poynting flux into the kinetic energy of the plasma beyond equipartition. A promising candidate is the kink instability, which deforms the body of the jet through helical perturbations. Since the detailed structure of real jets is unknown, we explore a large family of cylindrical, force-free equilibria to get robust conclusions. We find that the growth rate of the instability depends primarily on two parameters: (i) the gradient of the poloidal magnetic field; (ii) the Lorentz factor of the perturbation, which is closely related to the velocity of the plasma. We provide a simple fitting formula for the growth rate of the instability. As a tentative application, we use our results to interpret the dynamics of the jet in the nearby active galaxy M87. We show that the kink instability becomes non-linear at a distance from the central black hole comparable to where the jet stops accelerating. Hence (at least for this object), the kink instability of the jet is a good candidate to drive the transition from a Poynting-dominated to a kinetic-energy-dominated flow.

The origin of fast radio bursts remains a puzzle. Suggestions have been made that they are produced within the Earth atmosphere, in stellar coronae, in other galaxies or at cosmological distances. If they are extraterrestrial, the implied brightness temperature is very high, and therefore, the induced scattering places constraints on possible models. In this paper, constraints are obtained on flares from coronae of nearby stars. It is shown that the radio pulses with the observed power could not be generated if the plasma density within and in the nearest vicinity of the source is as high as it is necessary in order to provide the observed dispersion measure. However, one cannot exclude a possibility that the pulses are generated within a bubble with a very low density and pass through the dense plasma only in the outer corona.

In this review we describe recent observational and theoretical developments in our understanding of pulsar winds and pulsar-wind nebulae (PWNe). We put special emphasis on the results from observations of well-characterized PWNe of various types (e.g., torus-jet and bowshock-tail), the most recent MHD modeling efforts, and the status of the flaring Crab PWN puzzle.

The collapsar model explains the association of long duration gamma-Ray Bursts (GRBs) with stellar collapse. It involves a relativistic jet that forms at the core of a collapsing massive star. The jet penetrates the stellar envelope and the prompt GRB emission is produced once the jet is well outside the star. Most current models for generation of relativistic jets involve Poynting flux dominated outflows. We explore here the propagation of such a jet through a stellar envelope. The jet forms a bow shock around it. Energy dissipation at the head of this shock supplies energy to a cocoon that surrounds the jet. This cocoon exerts pressure on the jet and collimates it. While this description resembles the propagation of a hydrodynamic jets there are significant qualitative differences. Two Strong shocks, the reverse shock that slows down the hydrodynamic jet and the collimation shock that collimates it, cannot form within the Poynting flux dominated jet. As a result this jet moves much faster and dissipates much less energy while it crosses the stellar envelope. We construct here a simple analytic model that explores, self consistently, the jet-cocoon interaction and dynamics. Using this model we determine the properties of the jet, including its velocity, propagation time and shape.

Bursts of millisecond duration were recently discovered in the 1 GHz band. There is a strong evidence that they come from $\sim 1 $ Gpc distances, which implies extraordinary high brightness temperature. I propose that these bursts could be attributed to synchrotron maser emission from relativistic, magnetized shocks. At the onset of the magnetar flare, a strongly magnetized pulse is formed, which propagates away through the relativistic magnetar wind and eventually reaches the nebula inflated by the wind within the surrounding medium. I show that the observed radio bursts could be generated at shocks formed via the interaction of the magnetic pulse with the plasma within the nebula. The model predicts strong millisecond bursts in the TeV band, which could be observed even from distant galaxies.

The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to reconcile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside the Thomson photosphere. This quasi-thermal spectrum can evolve into the observed non-thermal shape by additional energy release at moderate to small Thomson optical depths, which can readily give rise to the hard spectral tail. The position of the spectral peak is determined by the temperature and Lorentz factor of the flow in the termalization zone, where the total number of photons carried by the jet is established. To reach thermalization, dissipation alone is not sufficient and photon generation requires an efficient emission/absorption process in addition to scattering. We perform a systematic study of all relevant photon production mechanisms searching for possible conditions in which thermalization can take place. We find that a significant fraction of the available energy should be dissipated at intermediate radii, $\sim 10^{10}$-- a few$\times 10^{11}$ cm and the flow there should be relatively slow: the bulk Lorentz factor could not exceed a few tens for all but the most luminous bursts with the highest $\Epk$-s. The least restrictive constraint for successful thermalization, $\Gamma\lesssim 20$, is obtained if synchrotron emission acts as the photon source. This requires, however, a non-thermal acceleration deep below the Thomson photosphere transferring a significant fraction of the flow energy to relativistic electrons with Lorentz factors between 10 and 100. Other processes require bulk flow Lorentz factors of order of a few for typical bursts. We examine the implications of these results to different GRB photospheric emission models.

The recently discovered gamma-ray flares from the Crab nebula are generally attributed to the magnetic energy release in a highly magnetized region within the nebula. I argue that such a region naturally arises in the polar region of the inner nebula. In pulsar winds, efficient dissipation of the Poynting flux into the plasma energy occur only in the equatorial belt where the energy is predominantly transferred by alternating fields. At high latitudes, the pulsar wind remains highly magnetized therefore the termination shock in the polar region is weak and the postshock flow remains relativistic. I study the structure of this flow and show that the flow at first expands and decelerates and then it converges and accelerates. In the converging part of the flow, the kink instability triggers the magnetic dissipation. The energy release zone occurs at the base of the observed jet. A specific turbulence of relativistically shrinking magnetic loops efficiently accelerates particles so that the synchrotron emission in the hundreds MeV band, both persistent and flaring, comes from this site.

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

The subpulse modulation of pulsar radio emission goes to prove that the plasma flow in the open field line tube breaks into isolated narrow streams. I propose a model which attributes formation of streams to the process of the electric current adjustment in the magnetosphere. A mismatch between the magnetospheric current distribution and the current injected by the polar cap accelerator gives rise to reverse plasma flows in the magnetosphere. The reverse flow shields the electric field in the polar gap and thus shuts up the plasma production process. I assume that a circulating system of streams is formed such that the upward streams are produced in narrow gaps separated by downward streams. The electric drift is small in this model because the potential drop in narrow gaps is small. The gaps have to drift because by the time a downward stream reaches the star surface and shields the electric field, the corresponding gap has to shift. The transverse size of the streams is determined by the condition that the potential drop in the gaps is sufficient for the pair production. This yields the radius of the stream roughly 10% of the polar cap radius, which makes it possible to fit in the observed morphological features such as the "carousel" with 10-20 subbeams and the system of the core - two nested cone beams.

The instabilities of relativistic ion beams in a relativistically hot electron background are derived for general propagation angles. It is shown that the Weibel instability in the direction perpendicular to the streaming direction is the fastest growing mode, and probably the first to appear, consistent with the aligned filaments that are seen in PIC simulations. Oblique, quasiperpendicular modes grow almost as fast, as the growth rate varies only moderately with angle, and they may distort or corrugate the filaments after the perpendicular mode saturates.

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter \sigma, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.

The structure of a relativistically hot, strongly magnetized jet is investigated at large distances from the source. Asymptotic equations are derived describing collimation and acceleration of the externally confined jet. Conditions are found for the transformation of the thermal energy into the fluid kinetic energy or into the Poynting flux. Simple scalings are presented for the jet collimation angle and Lorentz factors.

Reconnection of alternating magnetic fields is an important energy transformation mechanism in Poynting dominated outflows. We show that the reconnection is facilitated by the Kruskal-Schwarzschild instability of current sheets separating the oppositely directed fields. This instability, which is a magnetic counterpart of the Rayleigh-Taylor instability, develops if the flow is accelerated. Then the plasma drips out of the current sheet providing conditions for rapid reconnection. Since the magnetic dissipation leads to the flow acceleration, the process is self-sustaining. In pulsar winds, this process could barely compete with the earlier proposed dissipation mechanisms. However, the novel mechanism turns out to be very efficient at AGN and GRB conditions.

The curvature drift instability has long been considered as a viable mechanism for pulsar radio emission. We reconsidered this mechanism by finding an explicit solution describing propagation of short-wave electro-magnetic waves in a plasma flow along curved magnetic field lines. We show that even though the waves could be amplified, the amplification factor remains very close to unity therefore this mechanism is unable to generate high brightness temperature emission from initial weak fluctuations.

The acceleration of relativistic jets from the Poynting to the matter dominated stage is considered. The are generally two collimation regimes, which we call equilibrium and non-equilibrium, correspondingly. In the first regime, the jet is efficiently accelerated till the equipartition between the kinetic and electro-magnetic energy. We show that after the equilibrium jet ceases to be Poynting dominated, the ratio of the electro-magnetic to the kinetic energy decreases only logarithmically so that such jets become truly matter dominated only at extremely large distances. Non-equilibrium jets remain generally Poynting dominated till the logarithmically large distances. In the only case when a non-equilibrium jet is accelerated till the equipartition level, we found that the flow is not continued to the infinity but is focused towards the axis at a finite distance from the origin.

We present a simple numerical model of the plasma flow within the open field line tube in the pulsar magnetosphere. We study how the plasma screens the rotationally induced electric field and maintains the electric current demanded by the global structure of the magnetosphere. We show that even though bulk of the plasma moves outwards with relativistic velocities, a small fraction of particles is continuously redirected back forming reverse plasma flows. The density and composition (positrons or electrons, or both) of these reverse flows are determined by the distribution of the Goldreich-Julian charge density along the tube and by the global magnetospheric current. These reverse flows could significantly affect the process of the pair plasma production in the polar cap accelerator. Our simulations also show that formation of the reverse flows is accompanied by the generation of long wavelength plasma oscillations, which could be converted, via the induced scattering on the bulk plasma flow, into the observed radio emission.

In relativistic, Poynting dominated outflows, acceleration and collimation are intimately connected. An important point is that the Lorentz force is nearly compensated by the electric force therefore the acceleration zone spans a large range of scales. We derived the asymptotic equations describing relativistic, axisymmetric MHD flows far beyond the light cylinder. These equations do not contain either intrinsic small scales (like the light cylinder radius) or terms that nearly cancel each other (like the electric and magnetic forces) therefore they could be easily solved numerically. They also suit well for qualitative analysis of the flow and in many cases, they could even be solved analytically or semi-analytically. We show that there are generally two collimation regimes. In the first regime, the residual of the hoop stress and the electric force is counterbalanced by the pressure of the poloidal magnetic field so that at any distance from the source, the structure of the flow is the same as the structure of an appropriate cylindrical equilibrium configuration. In the second regime, the pressure of the poloidal magnetic field is negligible small so that the flow could be conceived as composed from coaxial shrinking magnetic loops. In the two collimation regimes, the flow is accelerated in different ways. We study in detail the structure of jets confined by the external pressure with a power law profile. In particular, we obtained simple scalings for the extent of the acceleration zone, for the terminal Lorentz factor and for the collimation angle.

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

We study the compression driven magnetic reconnection in the relativistic electron-positron plasma. Making use of a 2.5D particle-in-cell code, we simulated compression of a magnetized plasma layer containing a current sheet within it. We found that the particle spectrum within the reconnecting sheet becomes non-thermal; it could be approximated by a power-law distribution with an index of -1 and an exponential cutoff.

Effect of the induced Compton and Raman scattering on short, bright radio pulses is investigated. It is shown that when a single pulse propagates through the scattering medium, the effective optical depth is determined by the duration of the pulse but not by the scale of the medium. The induced scattering could hinder propagation of the radio pulse only if close enough to the source a dense enough plasma is presented. The induced scattering within the relativistically moving source places lower limits on the Lorentz factor of the source. The results are applied to the recently discovered short extragalactic radio pulse.

We review current theoretical ideas on pulsar winds and their surrounding nebulae. Relativistic MHD models of the wind of the aligned rotator, and of the striped wind, together with models of magnetic dissipation are discussed. It is shown that the observational signature of this dissipation is likely to be point-like, rather than extended, and that pulsed emission may be produced. The possible pulse shapes and polarisation properties are described. Particle acceleration at the termination shock of the wind is discussed, and it is argued that two distinct mechanisms must be operating, with the first-order Fermi mechanism producing the high-energy electrons (above 1 TeV) and either magnetic annihilation or resonant absorption of ion cyclotron waves responsible for the 100 MeV to 1 TeV electrons. Finally, MHD models of the morphology of the nebula are discussed and compared with observation.

Most of the rotational luminosity of a pulsar is carried away by a relativistic magnetised wind in which the matter energy flux is negligible compared to the Poynting flux. Near the equatorial plane of an obliquely rotating pulsar magnetosphere, the magnetic field reverses polarity with the pulsar period, forming a wind with oppositely directed field lines. This structure is called a striped wind; dissipation of alternating fields in the striped wind is the object of our study. The aim of this paper is to study the conditions required for magnetic energy release at the termination shock of the striped pulsar wind. Magnetic reconnection is considered via analytical methods and 1D relativistic PIC simulations. An analytical condition on the upstream parameters for partial and full magnetic reconnection is derived from the conservation laws of energy, momentum and particle number density across the relativistic shock. Furthermore, by using a 1D relativistic PIC code, we study in detail the reconnection process at the termination shock. We found a very simple criterion for dissipation of alternating fields at the termination shock, depending on the upstream parameters of the flow. 1D relativistic PIC simulations are in agreement with our criterion. Thus, alternating magnetic fields annihilate easily at relativistic highly magnetised shocks.

Slow dissipation of non-potential magnetic fields in the magnetosphere of the magnetar is assumed to accelerate particles to hundreds MeV along the magnetic field lines. We consider interaction of fast particles with the surface of the magnetar. We argue that the collisionless dissipation does not work in the atmosphere of the neutron star because the two-stream instability is stabilized by the inhomogeneity of the atmosphere. Rather, the dominant dissipation mechanism is collisional Landau levelexcitations followed by pair production via the deexcitation gamma-rays ultimately leading to electrons with the energy below the Landau energy. We show that, because of the effects of the superstrong magnetic field, these electrons could emerge from the surface carrying most of the original energy so that a hot corona arises with the temperature of 1 - 2 MeV. This extended corona is better suited than a thin atmosphere to convert most of the primary beam energy to non-thermal radiation and, as we show, most of the coronal energy release is radiated away in the hard X-ray and the soft gamma-ray bands by Comptonization and bremsstrahlung. The radiation spectrum is a power-law with the photon index $1<\alpha<2$. The model may account for the persistent hard X-ray emission discovered recently from the soft gamma-ray repeaters and anomalous X-ray pulsars and predicts that the radiation spectrum is extended into the MeV band.

A model is presented for the quasiperiodic component of magnetar emission during the tail phase of giant flares. The model invokes modulation of the particle number density in the magnetosphere. The magnetospheric currents are modulated by torsional motion of the surface and we calculate that the amplitude of neutron star surface oscillation should be ~1% of the NS radius in order to produce the observed features in the power spectrum. Using an axisymmetric analytical model for structure of the magnetosphere of an oscillating NS, we calculate the angular distribution of the optical depth to the resonant Compton scattering. The anisotropy of the optical depth may be why QPO are observed only at particular rotational phases.

X-ray emission following giant flares of magnetars can be categorized into three categories of time scales (a) short term afterglow (b) medium term afterglow and (c) long term afterglow. Short term afterglow, which declines over several hours, seems to correspond to gravitational resettling of uplifted material. Medium term afterglow, which declines over several weeks or months, appears to be the cooling of the heated outer crust, and long term afterglow, which declines over a period of many years, can be understood to be the cooling of the inner crust. The long term afterglow profile may be a very sensitive indicator of neutron star mass.

It is argued and demonstrated by particle-in-cell simulations that the synchrotron maser instability could develop at the front of a relativistic, magnetized shock. The instability generates strong low-frequency electromagnetic waves propagating both upstream and downstream of the shock. Upstream of the shock, these waves make electrons lag behind ions so that a longitudinal electric field arises and the electrons are accelerated up to the ion kinetic energy. Then thermalization at the shock front results in a plasma with equal temperatures of electrons and ions. Downstream of the shock, the amplitude of the maser-generated wave may exceed the strength of the shock-compressed background magnetic field. In this case the shock-accelerated particles radiate via nonlinear Compton scattering rather than via a synchrotron mechanism. The spectrum of the radiation differs, in the low-frequency band, from that of the synchrotron radiation, providing possible observational tests of the model.

It is estimated that the Weibel instability is not generally an effective mechanism for generating ultrarelativistic astrophysical shocks. Even if the upstream magnetic field is as low as in the interstellar medium, the shock is mediated not by the Weibel instability but by the Larmor rotation of protons in the background magnetic field. Future simulations should be able to verify or falsify our conclusion.

The 2004 Dec. 27 giant Gamma-ray flare detected from the magnetar SGR 1806-20 created an expanding radio nebula which we have monitored with the Australia Telescope Compact Array and the Very Large Array. These data indicate that there was an increase in the observed flux ~25 days after the initial flare that lasted for ~8 days, which we believe is the result of ambient material swept-up and shocked by this radio nebula. For a distance to SGR 1806-20 of 15 kpc, using the properties of this rebrightening we infer that the initial blast wave was dominated by baryonic material of mass M>10^{24.5} g. For an initial expansion velocity v~0.7c (as derived in an accompanying paper), we infer this material had an initial kinetic energy E>10^{44.5} ergs. If this material originated from the magnetar itself, it may have emitted a burst of ultra-high energy (E > 1 TeV) neutrinos far brighter than that expected from other astrophysical sources.

On 2004 Dec. 27, the soft gamma repeater (SGR) 1806-20 emitted the brightest giant flare (GF) ever detected from an SGR, with an (isotropic) energy release $\sim 100$ times greater than the only two other known SGR GFs. It was followed by a very bright, fading radio afterglow. Extensive follow-up radio observations provided a wealth of information with unprecedented astrometric precision, revealing the temporal evolution of the source size, along with densely sampled light curves and spectra. Here we expand on our previous work on this source, by explaining these observations within one self-consistent dynamical model. In this scenario, the early radio emission is due to the outflow ejected during the GF energizing a thin shell surrounding a pre-existing cavity, where the observed steep temporal decay of the radio emission seen beginning on day 9 is attributed to the adiabatic cooling of the shocked shell. The shocked ejecta and external shell move outward together, driving a forward shock into the ambient medium, and are eventually decelerated by a reverse shock. As we show in Gelfand et al. (2005), the radio emission from the shocked external medium naturally peaks when significant deceleration occurs, and then decays relatively slowly. The dynamical modeling of the collision between the ejecta and the external shell together with the observed evolution of the source size (which is nicely reproduced in our model) suggest that most of the energy in the outflow was in mildly relativistic material, with an initial expansion velocity $v/c \lesssim 0.7d_{15}$, for a distance of $15d_{15}$ kpc to SGR 1806-20. An initially highly relativistic outflow would not have produced a long coasting phase at a mildly relativistic expansion velocity, as was observed.

The extraordinary giant flare (GF) of 2004 December 27 from the soft gamma repeater (SGR) 1806-20 was followed by a bright radio afterglow. We present an analysis of VLA observations of this radio afterglow from SGR 1806-20, consisting of previously reported 8.5 GHz data covering days 7 to 20 after the GF, plus new observations at 8.5 and 22 GHz from day 24 to 81. We detect motion in the flux centroid of the afterglow, at an average velocity of 0.26 +/- 0.03 c (assuming a distance of 15 kpc) at a position angle of -45 degrees. This motion, in combination with the growth and polarization measurements, suggests an asymmetric outflow, mainly from one side of the magnetar. We find a deceleration in the expansion, from ~9 mas/day to <5 mas/day. The time of deceleration is roughly coincident with the rebrightening in the radio light curve, as expected to result when the ejecta from the GF sweeps up enough of the external medium, and transitions from a coasting phase to the Sedov-Taylor regime. The radio afterglow is elongated and maintains a 2:1 axis ratio with an average position angle of -40 degrees (north through east), oriented perpendicular to the average intrinsic linear polarization angle.

Soft gamma repeaters (SGRs) are "magnetars", a small class of slowly spinning neutron stars with extreme surface magnetic fields, ~10^15 gauss. On 2004 December 27, a giant flare was detected from the magnetar SGR 1806-20, the third such event ever recorded. This burst of energy was detected by a variety of instruments and even caused an ionospheric disturbance in the Earth's upper atmosphere recorded around the globe. Here we report the detection of a fading radio afterglow produced by this outburst, with a luminosity 500 times larger than the only other detection of a similar source. From day 6 to day 19 after the flare from SGR 1806-20, a resolved, linearly polarized, radio nebula was seen, expanding at approximately a quarter the speed of light. To create this nebula, at least 4x10^43 ergs of energy must have been emitted by the giant flare in the form of magnetic fields and relativistic particles. The combination of spatially resolved structure and rapid time evolution allows a study in unprecedented detail of a nearby analog to supernovae and gamma-ray bursts.

Magnetars comprise two classes of rotating neutron stars (Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars), whose X-ray emission is powered by an ultrastrong magnetic field, B ~ 10^15 G. Occasionally SGRs enter into active episodes producing many short X-ray bursts; extremely rarely (about once per 50 years per source), SGRs emit a giant flare, an event with total energy at least 1000 times higher than their typical bursts. Here we report that, on 2004 December 27, SGR 1806-20 emitted the brightest extra-solar transient event ever recorded, even surpassing the full moon brightness for 0.2 seconds. The total (isotropic) flare energy is 2x10^46 erg, 100 times higher than the only two previous events, making this flare a once in a century event. This colossal energy release likely occurred during a catastrophic reconfiguration of the magnetar's magnetic field. Such an event would have resembled a short, hard Gamma Ray Burst (GRB) if it had occurred within 40 Mpc, suggesting that extragalactic SGR flares may indeed form a subclass of GRBs.

Reconnection of the magnetic lines of force is considered in case the magnetic energy exceeds the rest energy of the matter. It is shown that the classical Sweet-Parker and Petschek models are generalized straightforwardly to this case and the reconnection rate may be estimated by substituting the Alfven velocity in the classical formulas by the speed of light. The outflow velocity in the Sweet-Parker configuration is mildly relativistic. In the Petschek configuration, the outflow velocity is ultrarelativistic whereas the angle between the slow shocks is very small. Due to the strong compression, the plasma outflow in the Petschek configuration may become strongly magnetized if the reconnecting fields are not exactly antiparallel.

SGR 1627-41 was discovered in 1998 after a single active episode which lasted \~6 weeks. We report here our monitoring results of the decay trend of the persistent X-ray luminosity of the source during the last 5 years. We find an initial temporal power law decay with index 0.47, reaching a plateau which is followed by a sharp (factor of ten) flux decline ~800 days after the source activation. The source spectrum is best described during the entire period by a single power law with high absorption (N_H=9.0(7)x10^(22) cm^(-2)); the spectral index, however, varies dramatically between 2.2-3.8 spanning the entire range for all known SGR sources. We discuss the cooling behavior of the neutron star assuming a deep crustal heating initiated by the burst activity of the source during 1998.

We consider the curvature radiation of the point-like charge moving relativistically along curved magnetic field lines through a pulsar magnetospheric electron-positron plasma. We demonstrate that the radiation power is largely suppressed as compared with the vacuum case, but still at a considerable level, high enough to explain the observed pulsar luminosities. The emitted radiation is polarized perpendicularly to the plane of the curved magnetic filed lines coincides with $ which can freely escape from the magnetospheric plasma. Our results strongly support the coherent curvature radiation by the spark-associated solitons as a plausible mechanism of pulsar radio emission.

The origin of radio emission from plerions is considered. Recent observations suggest that radio emitting electrons are presently accelerated rather than having been injected at early stages of the plerion evolution. The observed flat spectra without a low frequency cutoff imply an acceleration mechanism that raises the average particle energy by few orders of magnitude but leaves most of the particles at the energy less than about few hundred MeV. It is suggested that annihilation of the alternating magnetic field at the pulsar wind termination shock provides the necessary mechanism. Toroidal stripes of opposite magnetic polarity are formed in the wind emanated from an obliquely rotating pulsar magnetosphere (the striped wind). At the termination shock, the flow is compressed and the magnetic field annihilates by driven reconnection. Jump conditions are obtained for the shock in a striped wind. It is shown that postshock MHD parameters of the flow are the same as if the energy of alternating field has already been converted into the plasma energy upstream the shock. Therefore the available estimates of the ratio of the Poynting flux to the matter energy flux, $\sigma$, should be attributed not to the total upstream Poynting flux but only to that associated with the average magnetic field. A simple model for the particle acceleration in the shocked striped wind is presented.

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.

The X-ray afterglow that is observed following large flares on magnetars can be accurately fit by simple and quantitative theoretical models: The long term afterglow, lasting of order weeks, can be understood as thermal radiation of a heated neutron star crust that is re-scattered in the magnetosphere. Short term afterglow is well fit by the cooling of a non-degenerate, pair-rich layer, which gradually shrinks and releases heat to a pair-free zone above it. Measurements of persistent optical and infrared emission directly probe long lived currents in the magnetosphere which are a likely source of collective plasma emission. The superstrong magnetic field plays an important role in generating these various emissions, and previous inference of its strength in magnetars is supported by the good fits with observation.

Transient X-ray emission, with an approximate t^{-0.7} decay, was observed from SGR 1900+14 over 40 days following the the giant flare of 27 Aug 1998. We calculate in detail the diffusion of heat to the surface of a neutron star through an intense 10^{14}-10^{15} G magnetic field, following the release of magnetic energy in its outer layers. We show that the power law index, the fraction of burst energy in the afterglow, and the return to persistent emission can all be understood if the star is composed of normal baryonic material.

Fast magnetosonic waves in a magnetically-dominated plasma are investigated. In the pulsar wind, these waves may transport a significant fraction of the energy flux. It is shown that the nonlinear steepening and subsequent formation of multiple shocks is a viable mechanism for the wave dissipation in the pulsar wind. The wave dissipation both in the free pulsar wind and beyond the wind termination shock is considered.

Generation of the pulsar radio emission from plasma waves excited by the two-stream instability is considered. Special attention is given to propagation effects.

It is proposed that magnetospheric currents above the surfaces of magnetars radiate coherent emission in analogy to pulsars. Scaling the magnetospheric parameters suggests that the coherent emission from magnetars would emerge in the infra-red or optical.

The recent Chandra image of the Crab nebula shows a striking, axisymmetric polar jet. It is shown that jets are formed in axisymmetric, magnetized pulsar winds and that the jet luminosity scales relative to the total as (\gamma_0\sigma_{eq})^{-4/3}, where \sigma_{eq} is the ratio of Poynting flux to particle kinetic energy output at the equator at the base of the flow and \gamma_0 the initial Lorentz factor of the flow. The results are applied to the image of the Crab nebula, and the limit is set for the Crab pulsar of \sigma_{eq} \leq 100. It is argued that conventional pulsar theory needs to be reexamined in light of these limits.

The spin-down power of a pulsar is thought to be carried away in an MHD wind in which, at least close to the star, the energy transport is dominated by Poynting flux. The pulsar drives a low-frequency wave in this wind, consisting of stripes of toroidal magnetic field of alternating polarity, propagating in a region around the equatorial plane. The current implied by this configuration falls off more slowly with radius than the number of charged particles available to carry it, so that the MHD picture must, at some point, fail. Recently, magnetic reconnection in such a structure has been shown to accelerate the wind significantly. This reduces the magnetic field in the comoving frame and, consequently, the required current, enabling the solution to extend to much larger radius. This scenario is discussed and, for the Crab Nebula, the range of validity of the MHD solution is compared with the radius at which the flow appears to terminate. For sufficiently high particle densities, it is shown that a low frequency entropy wave can propagate out to the termination point. In this case, the "termination shock" itself must be responsible for dissipating the wave.

It is generally thought that most of the spin-down power of a pulsar is carried away in an MHD wind dominated by Poynting flux. In the case of an oblique rotator, a significant part of this energy can be considered to be in a low-frequency wave, consisting of stripes of toroidal magnetic field of alternating polarity, propagating in a region around the equatorial plane. Magnetic reconnection in such a structure has been proposed as a mechanism for transforming the Poynting flux into particle energy in the pulsar wind. We have re-examined this process and conclude that the wind accelerates significantly in the course of reconnection. This dilates the timescale over which the reconnection process operates, so that the wind requires a much larger distance than was previously thought in order to convert the Poynting flux to particle flux. In the case of the Crab, the wind is still Poynting-dominated at the radius at which a standing shock is inferred from observation. An estimate of the radius of the termination shock for other pulsars implies that all except the milli-second pulsars have Poynting-flux dominated winds all the way out to the shock front.