Showing 1–50 of 64 results for author: Lyubarsky, Y

On the escape of low-frequency waves from magnetospheres of neutron stars

Authors: Ephim Golbraikh, Yuri Lyubarsky

Abstract: We study the nonlinear decay of the fast magnetosonic 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.… ▽ More We study the nonlinear decay of the fast magnetosonic 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. △ Less

Submitted 17 September, 2023; originally announced September 2023.

Comments: ApJ, in press

The precursor structure in relativistic shocks

Authors: Barel Skuratovsky, Yuri Lyubarsky, Tsvi Piran

Abstract: 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… ▽ More 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. △ Less

Submitted 6 October, 2023; v1 submitted 30 July, 2023; originally announced July 2023.

Comments: Revised. Several minor points clarified. Accepted for publication in MNRAS

The collimation of relativistic jets in post-neutron star binary merger simulations

Authors: Matteo Pais, Tsvi Piran, Yuri Lyubarsky, Kenta Kiuchi, Masaru Shibata

Abstract: 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-domin… ▽ More 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. △ Less

Submitted 13 March, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 10 pages, 5 figures, 2 tables, accepted to ApJ Letters

Radiatively driven evaporation from magnetar's surface

Authors: Ivan Demidov, Yuri Lyubarsky

Abstract: 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 stud… ▽ More 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. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: 11 pages, 5 figures

Saturation of the filamentation instability and dispersion measure of Fast Radio Bursts

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei M. Beloborodov, Lorenzo Sironi, Masanori Iwamoto

Abstract: 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 stead… ▽ More 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. △ Less

Submitted 12 January, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 9 pages, accepted for publication in ApJL

Filamentation of Fast Radio Bursts in magnetar winds

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei Beloborodov, Lorenzo Sironi

Abstract: 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… ▽ More 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 $τ_{\rm sc}\sim{\rm\; μs-ms}$ at the frequency $ν\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $τ_{\rm sc}\proptoν^{-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 $τ_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $Δν\sim 1/τ_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind. △ Less

Submitted 25 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Section 3 is modified. Accepted for publication in MNRAS

Nonlinear electromagnetic-wave interactions in pair plasma: (I) Non-relativistic regime

Authors: Arka Ghosh, Daniel Kagan, Uri Keshet, Yuri Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2021; originally announced November 2021.

Comments: 18 pages, 15 figures, Comments welcome

Journal ref: The Astrophysical Journal, Volume 930, Issue 2, id.106, 18 pp. (2022)

Emission Mechanisms of Fast Radio Bursts

Authors: Yuri Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 28 February, 2021; originally announced March 2021.

Comments: To be published in Universe

The role of resonant plasma instabilities in the evolution of blazar induced pair beams

Authors: Roy Perry, Yuri Lyubarsky

Abstract: 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 i… ▽ More 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 $Δθ\sim 10^{-6}-10^{-5}$. In this case, the width of the resonance for parallel waves, $\proptoΔθ^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. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: To be published in MNRAS

Self-modulation of Fast Radio Bursts

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei M. Beloborodov, Lorenzo Sironi

Abstract: 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… ▽ More 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. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: accepted for publication in MNRAS

Spectral Modification of Magnetar Flares by Resonant Cyclotron Scattering

Authors: Shotaro Yamasaki, Yuri Lyubarsky, Jonathan Granot, Ersin Gogus

Abstract: 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 speed… ▽ More 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. △ Less

Submitted 24 July, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: 12 pages, 6 figures, accepted to MNRAS; Appendix B is added

Fast radio bursts from reconnection in magnetar magnetosphere

Authors: Yuri Lyubarsky

Abstract: 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,… ▽ More 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. △ Less

Submitted 27 May, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: Sect. 3 is modified and Appendix B is added

Magnetic energy dissipation and origin of non-thermal spectra in radiatively efficient relativistic sources

Authors: Emanuele Sobacchi, Yuri E. Lyubarsky

Abstract: 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 el… ▽ More 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}γ\proptoγ^{-2}$ for $γ\lesssimγ_{\rm heat}$ (the Lorentz factor $γ_{\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. △ Less

Submitted 8 January, 2020; v1 submitted 26 November, 2019; originally announced November 2019.

Comments: accepted for publication by MNRAS (minor changes)

On the magnetisation and the radiative efficiency of BL Lac jets

Authors: E. Sobacchi, Y. E. Lyubarsky

Abstract: 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 (… ▽ More 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. △ Less

Submitted 29 December, 2018; originally announced December 2018.

Comments: Submitted to MNRAS. Revised after the first referee report

Radio emission from the Crab and Crab-like pulsars

Authors: Yuri Lyubarsky

Abstract: 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 coul… ▽ More 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. △ Less

Submitted 27 November, 2018; originally announced November 2018.

Comments: MNRAS, in press

Instability induced by recollimation in highly magnetised outflows

Authors: Emanuele Sobacchi, Yuri E. Lyubarsky

Abstract: 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 v… ▽ More 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. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: 7 pages, accepted by MNRAS

Interaction of the electro-magnetic precursor from a relativistic shock with the upstream flow. I. Synchrotron absorption of strong electromagnetic waves

Authors: Yuri Lyubarsky

Abstract: 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 accel… ▽ More 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. △ Less

Submitted 27 November, 2017; originally announced November 2017.

Comments: 8 pages, accepted for publication in MNRAS

External Confinement and Surface Modes in Magnetised Force-Free Jets

Authors: E. Sobacchi, Y. E. Lyubarsky

Abstract: 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 ha… ▽ More 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 $Γ\sim10$ ($Γ\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. △ Less

Submitted 3 October, 2017; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: 9 pages, accepted by MNRAS

2D Relativistic MHD Simulations of the Kruskal-Schwarzschild Instability in a Relativistic Striped Wind

Authors: Ramandeep Gill, Jonathan Granot, Yuri Lyubarsky

Abstract: 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… ▽ More 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σ\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}=β_{\rm in}c$, of $10^{-3}\lesssimβ_{\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. △ Less

Submitted 13 November, 2017; v1 submitted 10 August, 2017; originally announced August 2017.

Comments: 12 pages, 12 figures, Accepted for Publication in MNRAS

Kink Instability of Force-Free Jets: a Parameter Space Study

Authors: E. Sobacchi, Y. E. Lyubarsky, M. C. Sormani

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2017; v1 submitted 10 February, 2017; originally announced February 2017.

Comments: 7 pages, 5 figures, accepted by MNRAS. Minor revision after referee comments

Induced scattering limits on fast radio bursts from stellar coronae

Authors: Yuri Lyubarsky, Sofiya Ostrovska

Abstract: 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 flar… ▽ More 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. △ Less

Submitted 23 December, 2015; originally announced December 2015.

Comments: Accepted ApJ

Pulsar-Wind Nebulae: Recent Progress in Observations and Theory

Authors: Oleg Kargaltsev, Benoit Cerutti, Yuri Lyubarsky, Edoardo Striani

Abstract: 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. 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. △ Less

Submitted 30 July, 2015; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: 60 pages, 20 figures; slightly modified version of the article published in Space Science Reviews; citations added in the revised version

The dynamics of a highly magnetized jet propagating inside a star

Authors: Omer Bromberg, Jonathan Granot, Yuri Lyubarsky, Tsvi Piran

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2014; originally announced February 2014.

Comments: 31 pages, 6 figures

A model for fast extragalactic radio bursts

Authors: Yuri Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 26 January, 2014; originally announced January 2014.

Comments: 5 pages

Highly magnetized region in pulsar wind nebulae and origin of the Crab gamma-ray flares

Authors: Y. E. Lyubarsky

Abstract: 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 predominant… ▽ More 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. △ Less

Submitted 7 September, 2012; originally announced September 2012.

Comments: Submitted to MNRAS

On thermalization in gamma-ray burst jets and the peak energies of photospheric spectra

Authors: Indrek Vurm, Yuri Lyubarsky, Tsvi Piran

Abstract: 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 t… ▽ More 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, $Γ\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. △ Less

Submitted 21 December, 2012; v1 submitted 4 September, 2012; originally announced September 2012.

Comments: 18 pages, 10 figures, accepted to ApJ

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. III. Rotating Relativistic Jets

Authors: Yosuke Mizuno, Yuri Lyubarsky, Ken-Ichi Nishikawa, Philip E. Hardee

Abstract: 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… ▽ More 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. △ Less

Submitted 20 July, 2012; originally announced July 2012.

Comments: 30 pages, 11 figures, and 1 table, accepted to ApJ

Adjustment of the electric current in pulsar magnetospheres and origin of subpulse modulation

Authors: Yuri Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 7 June, 2012; originally announced June 2012.

Comments: 8 pages, 1 figure

Stream instabilities in relativistically hot plasma

Authors: Rashid Shaisultanov, Yuri Lyubarsky, David Eichler

Abstract: 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… ▽ More 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. △ Less

Submitted 4 April, 2011; originally announced April 2011.

Comments: 10 pages, 6 figures

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. II. Relaxation of Pulsar Wind Nebula

Authors: Yosuke Mizuno, Yuri Lyubarsky, Ken-Ichi Nishikawa, Philip E. Hardee

Abstract: 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… ▽ More 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 σ, 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. △ Less

Submitted 13 December, 2010; originally announced December 2010.

Comments: 14 pages, 5 figures, accepted for publication in ApJ

Asymptotic theory of relativistic, magnetized jets

Authors: Yuri Lyubarsky

Abstract: 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 a… ▽ More 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. △ Less

Submitted 12 December, 2010; originally announced December 2010.

Comments: To appear in Phys.RevE

A new mechanism for dissipation of alternating fields in Poynting dominated outflows

Authors: Yuri Lyubarsky

Abstract: 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 t… ▽ More 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. △ Less

Submitted 7 December, 2010; originally announced December 2010.

Journal ref: ApJL, 725, L234 (2010)

Curvature-drift instability fails to generate pulsar radio emission

Authors: Alexander Kaganovich, Yuri Lyubarsky

Abstract: 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 t… ▽ More 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. △ Less

Submitted 29 August, 2010; originally announced August 2010.

Comments: to appear in ApJ

Transformation of the Poynting flux into the kinetic energy in relativistic jets

Authors: Y. E. Lyubarsky

Abstract: 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 d… ▽ More 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. △ Less

Submitted 30 September, 2009; v1 submitted 25 September, 2009; originally announced September 2009.

Comments: Submitted to MNRAS Minor changes in the Conclusions

Adjustment of the electric charge and current in pulsar magnetospheres

Authors: Yuri Lyubarsky

Abstract: 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 i… ▽ More 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. △ Less

Submitted 16 April, 2009; originally announced April 2009.

Comments: 24 pages, 11 figures

Journal ref: Astrophys.J.696:320-327, 2009

3D Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. I. Instability of a static column

Authors: Yosuke Mizuno, Yuri Lyubarsky, Ken-Ichi Nishikawa, Philip E. Hardee

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2009; originally announced March 2009.

Comments: 24 pages, 10 figures, and 1 table, submitted to ApJ

Journal ref: Astrophys.J.700:684-693,2009

Asymptotic structure of Poynting dominated jets

Authors: Yuri Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2009; v1 submitted 19 February, 2009; originally announced February 2009.

Comments: Submitted to ApJ; the second version: corrected typos, minor changes in the text

Journal ref: Astrophys.J.698:1570-1589,2009

Particle acceleration in the driven relativistic reconnection

Authors: Yuri Lyubarsky, Michael Liverts

Abstract: 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 an… ▽ More 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. △ Less

Submitted 7 May, 2008; v1 submitted 1 May, 2008; originally announced May 2008.

Comments: 15 pages, 9 figures, to appear in the ApJ August 10, 2008, v683n 1 issue

Induced scattering of short radio pulses

Authors: Yuri Lyubarsky

Abstract: 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 enou… ▽ More 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. △ Less

Submitted 13 April, 2008; originally announced April 2008.

Comments: submitted to ApJ

Magnetic reconnection at the termination shock in a striped pulsar wind

Authors: Jerome Petri, Yuri Lyubarsky

Abstract: 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 win… ▽ More 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. △ Less

Submitted 12 July, 2007; originally announced July 2007.

Comments: Accepted by A&A

Journal ref: Int.J.Mod.Phys.D17:1961-1967,2008

On the nature of QPO in the tail of SGR giant flares

Authors: A. N. Timokhin, D. Eichler, Yu. Lyubarsky

Abstract: 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… ▽ More 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. △ Less

Submitted 25 June, 2007; originally announced June 2007.

Comments: 7 pages, 4 figures, emulateapj style

Journal ref: Astrophys.J.680:1398-1404,2008

On the corona of magnetars

Authors: Yury Lyubarsky, David Eichler

Abstract: 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… ▽ More 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<α<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. △ Less

Submitted 25 June, 2007; originally announced June 2007.

Comments: 23 pages, 2 fig. Submitted to ApJ

The theory of pulsar winds and nebulae

Authors: J. G. Kirk, Y. Lyubarsky, J. Petri

Abstract: 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… ▽ More 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. △ Less

Submitted 18 July, 2007; v1 submitted 6 March, 2007; originally announced March 2007.

Comments: 33 pages, to appear in Springer Lecture Notes on "Neutron stars and pulsars, 40 years after the discovery", ed W.Becker

Long Term Cooling of Magnetar Crusts

Authors: David Eichler, Yuri Lyubarsky, Chryssa Kouveliotou, Colleen A. Wilson

Abstract: 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 th… ▽ More 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. △ Less

Submitted 23 November, 2006; originally announced November 2006.

Comments: Invited talk at conference on "Solitary Neutron Stars", London, 2006

Electron-ion coupling upstream of relativistic collisionless shocks

Authors: Yuri Lyubarsky

Abstract: 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 arise… ▽ More 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. △ Less

Submitted 1 November, 2006; originally announced November 2006.

Comments: 22 pages, 10 figures. To appear in ApJ vol. 653

Journal ref: Astrophys.J.652:1297-1305,2006

Are GRB shocks mediated by the Weibel instability?

Authors: Yuri Lyubarsky, David Eichler

Abstract: 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… ▽ More 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. △ Less

Submitted 22 December, 2005; originally announced December 2005.

Comments: 13 pages, 3 figures

Journal ref: Astrophys.J.647:1250-1254,2006

The Growth, Polarization, and Motion of the Radio Afterglow from the Giant Flare from SGR 1806-20

Authors: G. B. Taylor, J. D. Gelfand, B. M. Gaensler, J. Granot, C. Kouveliotou, R. P. Fender, E. Ramirez-Ruiz, D. Eichler, Y. E. Lyubarsky, M. Garrett, R. A. M. J. Wijers

Abstract: 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 ce… ▽ More 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. △ Less

Submitted 10 October, 2005; v1 submitted 16 April, 2005; originally announced April 2005.

Comments: 5 pages, 3 figures, as accepted to ApJ Letters

Journal ref: Astrophys.J. 634 (2005) L93-L96

A Re-brightening of the Radio Nebula associated with the 2004 December 27 giant flare from SGR 1806--20

Authors: J. D. Gelfand, Y. E. Lyubarsky, D. Eichler, B. M. Gaensler, G. B. Taylor, J. Granot, K. J. Newton-McGee, E. Ramirez-Ruiz, C. Kouveliotou, R. A. M. J. Wijers

Abstract: 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 sho… ▽ More 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. △ Less

Submitted 31 October, 2005; v1 submitted 11 March, 2005; originally announced March 2005.

Comments: 7 pages, 2 figures. Accepted by ApJL, minor changes due to referee comments

Journal ref: Astrophys.J. 634 (2005) L89-L92

Diagnosing the outflow from the SGR 1806-20 Giant Flare with Radio Observations

Authors: J. Granot, E. Ramirez-Ruiz, G. B. Taylor, D. Eichler, Y. E. Lyubarsky, R. A. M. J. Wijers, B. M. Gaensler, J. D. Gelfand, C. Kouveliotou

Abstract: 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, reveal… ▽ More 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. △ Less

Submitted 11 October, 2005; v1 submitted 13 March, 2005; originally announced March 2005.

Comments: 15 pages, 2 figures, clearer and more elaborate version; ApJ, in press

Journal ref: Astrophys.J. 638 (2006) 391-396

Gamma-Ray Observations of a Giant Flare from The Magnetar SGR 1806-20

Authors: D. M. Palmer, S. Barthelmy, N. Gehrels, R. M. Kippen, T. Cayton, C. Kouveliotou, D. Eichler, R. A. M. J. Wijers, P. M. Woods, J. Granot, Y. E. Lyubarsky, E. Ramirez-Ruiz, L. Barbier, M. Chester, J. Cummings, E. E. Fenimore, M. H. Finger, B. M. Gaensler, D. Hullinger, H. Krimm, C. B. Markwardt, J. A. Nousek, A. Parsons, S. Patel, T. Sakamoto , et al. (3 additional authors not shown)

Abstract: 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 ti… ▽ More 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. △ Less

Submitted 1 March, 2005; originally announced March 2005.

Comments: Submitted to Nature 2005-02-02, revised 2005-03-01. 21 pp, incl. 6 figures

Report number: LA-UR-05-0629

Journal ref: Nature 434 (2005) 1107-1109