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. [ABSTRACT FROM AUTHOR]

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. It was shown that in a steady state, there is a maximal fraction of the downstream energy flux that the precursor can carry. We show here that at this critical value, the shock disappears, and the flow passes through a sonic point. This behaviour 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. [ABSTRACT FROM AUTHOR]

We study the nonlinear decay of the fast magnetosonic (fms) into the Alfvén 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én 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én 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. [ABSTRACT FROM AUTHOR]

The luminosity of the Soft Gamma Repeater 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. [ABSTRACT FROM AUTHOR]

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 ≲ 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 time-scale depends on the properties of the magnetar wind. In a cold wind, the typical scattering time-scale is τsc ∼   |$\mu$| s–ms at the frequency |$\nu \sim 1\, {\rm GHz}$|⁠. The scattering time-scale increases at low frequencies, with the scaling τsc ∝  ν −2. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering time-scale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering time-scale 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}$|⁠. Broad-band frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind. [ABSTRACT FROM AUTHOR]

High brightness-temperature radiation is observed in various astrophysical sources: active galactic nuclei, pulsars, interstellar masers, and flaring stars; the discovery of fast radio bursts renewed interest in the nonlinear interaction of intense radiation with plasma. In astronomical systems, the radiation frequency is typically well above the plasma frequency and its spectrum is broad, so nonlinear processes differ considerably from those typically studied in laboratory plasma. This paper is the first in a series devoted to the numerical study of nonlinear interactions of electromagnetic waves with plasma. We start with nonmagnetized pair plasmas, where the primary processes are induced (Compton) scattering and filamentation instability. In this paper, we consider waves in which electron oscillations are nonrelativistic. 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 the 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. [ABSTRACT FROM AUTHOR]

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 analyse 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 Δθ ∼ 10−6 to 10−5. In this case, the width of the resonance for parallel waves, ∝Δθ2, is too small for them to grow in realistic conditions. We perform both analytical estimates and numerical simulations in the quasi-linear 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. [ABSTRACT FROM AUTHOR]

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 idealized setup of a monochromatic linearly polarized 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 crit ∼ 1017(N /102 cm−3)(L /1042 erg s−1) 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 ≲ R 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 ∼ R 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 time-scale of 10 µm, 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. [ABSTRACT FROM AUTHOR]

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 distortion 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 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. [ABSTRACT FROM AUTHOR]

The nearly 100% linear polarization has been reported for a few fast radio bursts (FRBs). This finding places severe limits on the emission mechanism. I argue that the completely 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 outward. 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. [ABSTRACT FROM AUTHOR]

This is the second in the series of papers aiming to study interaction of the electromagnetic precursor waves from relativistic shocks with the upstream flow. Here, I consider the induced scattering of strong waves. In such a wave, the electrons oscillate with relativistic velocities therefore, the scattering generally occurs in harmonics of the incident wave. I show that the induced scattering occurs predominantly in the first harmonics. I also show that even though in the weak case regime, the induced scattering rate is proportional to the intensity of the incident wave, in the strong wave case, the rate decreases as the wave amplitude grows. [ABSTRACT FROM AUTHOR]

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 electromagnetic nano-shots. I estimate the parameters of this emission and show that they are compatible with observations. [ABSTRACT FROM AUTHOR]

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 magnetohydrodynamic (RMHD) simulations featuring two cold and highlymagnetized (1 ≤ σ ≤ 103) 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 ~0.1c, 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 vin =βinc of 10-3 ≲βin ≲5×10-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. [ABSTRACT FROM AUTHOR]

This paper is the first in the series of papers aiming to study interaction of the electromagnetic 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 electromagnetic 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 electromagnetic waves. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 in front of it. Energy dissipation at the shock generates an energetic cocoon that surrounds the jet. This cocoon exerts pressure on the jet and collimates it. While this description resembles the propagation of a hydrodynamic jet 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. [ABSTRACT FROM PUBLISHER]

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 longitudinal 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. [ABSTRACT FROM AUTHOR]

The growth rates for Weibel and Buneman instabilities of relativistic ion beams in a relativistically hot electron background are derived analytically for general propagation angles. The Weibel instability perpendicular to the streaming direction is found to be the fastest growing mode and probably the first to appear. 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. The growth rate of the purely longitudinal (Buneman) mode is significantly smaller, contrary to the non-relativistic case. The results are consistent with simulations, which display aligned magnetic filaments and their subsequent disruption. [ABSTRACT FROM AUTHOR]

We have investigated the development of current-driven (CD) kink instability in relativistic jets via 3D RMHD simulations. In this investigation 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 CD kink instability. The linear growth and nonlinear evolution of the CD kink instability depend moderately on the radial density profile and strongly on the magnetic pitch profile. Kink amplitude growth in the nonlinear 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 nonlinear regime nearly ceases for increasing magnetic pitch. [ABSTRACT FROM AUTHOR]

Most of the rotational luminosity of a pulsar is carried away by a relativistic magnetized wind in which the matter energy flux is negligible compared to the Poynting flux. However, observations indicate that most of the Poynting flux is eventually converted into ultrarelativistic particles. The mechanism responsible for this transformation remains poorly understood. 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, called a striped wind. We 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 for different upstream Lorentz factors and magnetizations, and find good agreement with our analytical criterion. Thus, alternating magnetic fields annihilate easily at relativistic highly magnetized shocks. We apply our criterion for dissipation to deduce bounds on the pair multiplicity, κ, in the pulsar wind. [ABSTRACT FROM AUTHOR]

The Crab nebula is regarded as one of the most important “cosmic laboratories” in astrophysics, which has made a bigger impact on the development of astronomy than any other single object beyond the solar system. The most intriguing recent result is the completely unexpected discovery of a peculiar “jet-torus” structure in the inner part of the nebula. Similar structures were found later in other Crab-like nebulae. This discovery clearly indicates significant anisotropy of the wind from the Crab pulsar which has been ignored so far in simplified theoretical models of the nebula. Fortunately, the impressive progress in computational relativistic magnetohydrodynamics in recent years has made possible to study the Crab nebula without making such a drastic simplification of the problem. In this paper we present the results of the first study of such kind. They provide a likely explanation of the jet-torus pattern and show that the flow in the nebula is much more complex than it has been widely believed. [ABSTRACT FROM AUTHOR]

We have investigated the development of current-driven (CD) kink instability in relativistic jets via 3D RMHD simulations. In this investigation a static force-free equilibrium helical magnetic field 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 CD kink instability. The linear growth and nonlinear evolution of the CD kink instability depends moderately on the radial density profile and strongly on the magnetic pitch profile. Kink amplitude growth in the nonlinear 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 nonlinear regime nearly ceases for increasing magnetic pitch. [ABSTRACT FROM PUBLISHER]

Bursts of millisecond duration were recently discovered in the 1 GHz band. There is a strong evidence that they come from ∼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. [ABSTRACT FROM AUTHOR]

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 accordance with the linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outward from the axis, then the initial small perturbations grow strongly, and if multiple wavelengths are excited, then nonlinear 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. [ABSTRACT FROM AUTHOR]