Source author record

Yuri Lyubarsky

Yuri Lyubarsky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2025arXiv

Cavitation instability in unmagnetized relativistic pair shocks

We investigate the formation of plasma cavities in unmagnetized relativistic pair shocks and demonstrate that these cavities emerge naturally as a nonlinear consequence of asymmetric Weibel instability. We provide an intuitive physical picture and a minimal fluid model that captures the essential features of this cavitation instability and compare them with PIC results. This mechanism may provide the missing link between kinetic Weibel turbulence and macroscopic magnetic fields in astrophysical shocks.

preprint2022arXiv

Filamentation of Fast Radio Bursts in magnetar winds

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.

preprint2021arXiv

Emission Mechanisms of Fast Radio Bursts

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.

preprint2021arXiv

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

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.

preprint2020arXiv

Fast radio bursts from reconnection in magnetar magnetosphere

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.

preprint2020arXiv

Spectral Modification of Magnetar Flares by Resonant Cyclotron Scattering

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.

preprint2015arXiv

Induced scattering limits on fast radio bursts from stellar coronae

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.

preprint2014arXiv

A model for fast extragalactic radio bursts

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.

preprint2014arXiv

The dynamics of a highly magnetized jet propagating inside a star

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

Stream instabilities in relativistically hot plasma

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.

preprint2010arXiv

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

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.

preprint2010arXiv

Asymptotic theory of relativistic, magnetized jets

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.

preprint2010arXiv

Curvature-drift instability fails to generate pulsar radio emission

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.

preprint2010arXiv

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

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.

preprint2009arXiv

Adjustment of the electric charge and current in pulsar magnetospheres

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.

preprint2009arXiv

Asymptotic structure of Poynting dominated jets

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.