Search

Collisionless tearing instability with a power-law distribution function in a relativistic pair plasma with a guide field is studied. When the current sheet is supported by plasma pressure, the tearing mode is suppressed as the particle spectrum hardens. In the force-free limit, the instability growth rate becomes independent of the particle spectrum. We apply these results to relativistic MHD turbulence, where magnetic energy greatly exceeds plasma rest energy, and derive an expression for the transverse size of turbulent eddies unstable to tearing mode. We also establish the critical plasma magnetization parameter above which charge starvation prevents the tearing instability. These results might be useful for developing more accurate models of particle acceleration in relativistic astrophysical sources.

We study the tearing instability of a current sheet in a relativistic pair plasma with a power law distribution function. We first estimate the growth rate analytically and then confirm the analytical results by solving numerically the dispersion equation, taking into account all exact particle trajectories within the reconnecting layer. We found that the instability is suppressed when the particle spectrum becomes harder., Comment: 20 pages, 6 figures

We study the nonlinear decay of the fast magnetosonic into the Alfv\'en waves in relativistic force-free magnetohydrodynamics. The work has been motivated by models of pulsar radio emission and fast radio bursts (FRBs), in which the emission is generated in neutron star magnetospheres at conditions when not only the Larmor but also the plasma frequencies significantly exceed the radiation frequency. The decay process places limits on the source luminosity in these models. We estimated the decay rate and showed that the phase volume of Alfv\'en waves available for the decay of an fms wave is infinite. Therefore the energy of fms waves could be completely transferred to the small-scale Alfv\'en waves not via a cascade, as in the Kolmogorov turbulence, but directly. Our results explain the anomalously low radio efficiency of the Crab pulsar and show that FRBs could not be produced well within magnetar magnetospheres., Comment: ApJ, in press

We present a common unifying macroscopic framework for precursors in relativistic shock waves. These precursors transfer energy and momentum from the hot downstream to the cold upstream, modifying the shock structure. Derishev & Piran (2016} have shown that in a steady state, there is a maximal fraction of the downstream energy flux that the precursor can carry. We show that at this critical value, the shock disappears, and the flow passes through a sonic point. This behavior resembles the classical Newtonian Rayleigh flow problem. At the critical value, the transition is unstable as perturbations in the upstream accumulate at the sonic point. Thus, if such a point is reached, the shock structure is drastically modified, and the flow becomes turbulent., Comment: Revised. Several minor points clarified. Accepted for publication in MNRAS

The gravitational waves from the binary neutron star merger GW170817 were accompanied by a multi-wavelength electromagnetic counterpart, which confirms the association of the merger with a short gamma-ray burst (sGRB). The afterglow observations implied that the event was accompanied by a narrow, $\sim 5~$deg, and powerful, $\sim 10^{50}$ erg, jet. We study the propagation of a Poynting flux-dominated jet within the merger ejecta (kinematic, neutrino-driven and MRI turbulence-driven) of a neutrino-radiation-GR-MHD simulation of two coalescing neutron stars. We find that the presence of a post-merger low-density/low-pressure polar cavity, that arose due to angular momentum conservation, is crucial to let the jet break out. At the same time the ejecta collimates the jet to a narrow opening angle. The collimated jet has a narrow opening angle of $\sim 4$-$7$ deg and an energy of $10^{49}$-$10^{50}~$erg, in line with the observations of GW170817 and other sGRBs., Comment: 10 pages, 5 figures, 2 tables, accepted to ApJ Letters

The luminosity of the Soft Gamma Repeater (SGR) flares significantly exceeds the Eddington luminosity. This is because they emit mainly in the E-mode, for which the radiative cross-sections are strongly suppressed. The energy is released in the magnetosphere forming a magnetically trapped pair fireball, and the surface of the star is illuminated by the powerful radiation from the fireball. We study the ablation of the matter from the surface by this radiation. The E-mode photons are scattered within the surface layer, partly being converted into O-photons, whose scattering cross-section is of the order of the Thomson cross-section. The high radiation pressure of the O-mode radiation expels the plasma upwards. The uplifted matter forms a thick baryon sheath around the fireball. If an illuminated fraction of the star's surface includes the polar cap, a heavy, mildly relativistic baryonic wind is formed., Comment: 11 pages, 5 figures

Nonlinear effects are crucial for the propagation of Fast Radio Bursts (FRBs) near the source. We study the filamentation of FRBs in the relativistic winds of magnetars, which are commonly invoked as the most natural FRB progenitors. As a result of filamentation, the particle number density and the radiation intensity develop strong gradients along the direction of the wind magnetic field. A steady state is reached when the plasma pressure balances the ponderomotive force. In such a steady state, particles are confined into periodically spaced thin sheets, and electromagnetic waves propagate between them as in a waveguide. We show that: (i) The dispersion relation resembles that in the initial homogeneous plasma, but the effective plasma frequency is determined by the separation of the sheets, not directly by the mean particle density. (ii) The contribution of relativistic magnetar winds to the dispersion measure of FRBs could be several orders of magnitude larger than previously thought. The dispersion measure of the wind depends on the properties of individual bursts (e.g. the luminosity), and therefore can change significantly among different bursts from repeating FRBs. (iii) Induced Compton scattering is suppressed because most of the radiation propagates in near vacuum regions., Comment: 9 pages, accepted for publication in ApJL

Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and discuss implications for FRBs. As an effect of the instability, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. When the FRB front expands outside the radius where the instability ends, the radiation sheets are scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $\tau_{\rm sc}\sim{\rm\; \mu s-ms}$ at the frequency $\nu\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $\tau_{\rm sc}\propto\nu^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $\tau_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $\Delta\nu\sim 1/\tau_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind., Comment: Section 3 is modified. Accepted for publication in MNRAS

This paper is the first in a series devoted to the numerical study of nonlinear interactions of electromagnetic waves with plasma. We start with non-magnetized pair plasmas, where the primary processes are induced (Compton) scattering and the filamentation instability. In this paper, we consider waves in which electron oscillations are non-relativistic. Here, the numerical results can be compared to analytical theory, facilitating the development of appropriate numerical tools and framework. We distill the analytic theory, reconciling plasma and radiative transfer pictures of induced scattering and developing in detail the kinetic theory of modulation/filamentation instability. We carry out homogeneous numerical simulations using the particle-in-cell codes EPOCH and Tristan-MP, for both monochromatic waves and wave packets. We show that simulations of both processes are consistent with theoretical predictions, setting the stage for analyzing the highly nonlinear regime., Comment: 18 pages, 15 figures, Comments welcome

Fast radio bursts (FRBs) are recently discovered mysterious single pulses of radio emission, mostly coming from cosmological distances ($\sim 1$ Gpc). Their short duration, $\sim 1$ ms, and large luminosity evidence coherent emission. I review the basic physics of coherent emission mechanisms proposed for FRBs. In particular, I discuss the curvature emission of bunches, the synchrotron maser, and the emission of radio waves by variable currents in the course of magnetic reconnection. Special attention is paid to magnetar flares as the most promising sources of FRBs. Non-linear effects are outlined that could place bounds on the power of the outgoing radiation., Comment: To be published in Universe

The fate of relativistic pair beams produced in the intergalactic medium by very high energy emission from blazars remains controversial in the literature. The possible role of resonance beam plasma instability has been studied both analytically and numerically but no consensus has been reached. In this paper, we thoroughly analyze the development of this type of instability. This analysis takes into account that a highly relativistic beam loses energy only due to interactions with the plasma waves propagating within the opening angle of the beam (we call them parallel waves), whereas excitation of oblique waves results merely in an angular spreading of the beam, which reduces the instability growth rate. For parallel waves, the growth rate is a few times larger than for oblique ones, so they grow faster than oblique waves and drain energy from the beam before it expands. However, the specific property of extragalactic beams is that they are extraordinarily narrow; the opening angle is only $\Delta\theta\sim 10^{-6}-10^{-5}$. In this case, the width of the resonance for parallel waves, $\propto\Delta\theta^2$, is too small for them to grow in realistic conditions. We perform both analytical estimates and numerical simulations in the quasilinear regime. These show that for extragalactic beams, the growth of the waves is incapable of taking a significant portion of the beam's energy. This type of instability could at best lead to an expansion of the beam by some factor but the beam's energy remains nearly intact., Comment: To be published in MNRAS

Fast Radio Bursts (FRBs) are extreme astrophysical phenomena entering the realm of non-linear optics, a field developed in laser physics. A classical non-linear effect is self-modulation. We examine the propagation of FRBs through the circumburst environment using the idealised setup of a monochromatic linearly-polarised GHz wave propagating through a uniform plasma slab of density $N$ at distance $R$ from the source. We find that self-modulation occurs if the slab is located within a critical radius $R_{\rm crit}\sim 10^{17}(N/10^2{\rm\; cm}^{-3})(L/10^{42}{\rm\; erg\; s}^{-1}){\rm\; cm}$, where $L$ is the isotropic equivalent of the FRB luminosity. Self-modulation breaks the burst into pancakes transverse to the radial direction. When $R\lesssim R_{\rm crit}$, the transverse size of the pancakes is smaller than the Fresnel scale. The pancakes are strongly diffracted as the burst exits the slab, and interference between the pancakes produces a frequency modulation of the observed intensity with a sub-GHz bandwidth. When $R\sim R_{\rm crit}$, the transverse size of the pancakes becomes comparable with the Fresnel scale, and the effect of diffraction is weaker. The observed intensity is modulated on a timescale of ten microseconds, which corresponds to the radial width of the pancakes. Our results suggest that self-modulation may cause the temporal and frequency structure observed in FRBs., Comment: accepted for publication in MNRAS

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., Comment: 12 pages, 6 figures, accepted to MNRAS; Appendix B is added

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., Comment: Sect. 3 is modified and Appendix B is added

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., Comment: accepted for publication by MNRAS (minor changes)

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., Comment: MNRAS, in press

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., Comment: 7 pages, accepted by MNRAS

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., Comment: 8 pages, accepted for publication in MNRAS

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., Comment: 12 pages, 12 figures, Accepted for Publication in MNRAS

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., Comment: Accepted ApJ

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., Comment: 60 pages, 20 figures; slightly modified version of the article published in Space Science Reviews; citations added in the revised version

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., Comment: 31 pages, 6 figures

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., Comment: 5 pages

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., Comment: 18 pages, 10 figures, accepted to ApJ

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., Comment: 30 pages, 11 figures, and 1 table, accepted to ApJ

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., Comment: 8 pages, 1 figure

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., Comment: 10 pages, 6 figures

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., Comment: 14 pages, 5 figures, accepted for publication in ApJ

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., Comment: To appear in Phys.RevE

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., Comment: to appear in ApJ

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., Comment: 24 pages, 11 figures

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., Comment: 24 pages, 10 figures, and 1 table, submitted to ApJ

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., Comment: Submitted to ApJ; the second version: corrected typos, minor changes in the text

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., Comment: 15 pages, 9 figures, to appear in the ApJ August 10, 2008, v683n 1 issue

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., Comment: submitted to ApJ

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., Comment: Accepted by A&A

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., Comment: Invited talk at conference on "Solitary Neutron Stars", London, 2006

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., Comment: 22 pages, 10 figures. To appear in ApJ vol. 653

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., Comment: 13 pages, 3 figures

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.