Source author record

Maxim Lyutikov

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

ResearcherUnclaimed source record

astro-ph.HE physics.plasm-ph astro-ph astro-ph.SR gr-qc physics.flu-dyn astro-ph.CO astro-ph.EP astro-ph.IM chao-dyn nlin.CD physics.comp-ph physics.space-ph

Catalog footprint

What is connected

81works

13topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Coherent emission in pulsars, magnetars and Fast Radio Bursts: reconnection-driven free electron laser

We develop a model of the generation of coherent radio emission in the Crab pulsar, magnetars and Fast Radio Bursts (FRBs). Emission is produced by a reconnection-generated beam of particles via a variant of Free Electron Laser (FEL) mechanism, operating in a weakly-turbulent, guide-field dominated plasma. We first consider nonlinear Thomson scattering in a guide-field dominated regime, and apply to model to explain emission bands observed in Crab pulsar and in Fast Radio Bursts. We consider particle motion in a combined fields of the electromagnetic wave and thee lectromagnetic (Alfvenic) wiggler. Charge bunches, created via a ponderomotive force, Compton/Raman scatter the wiggler field coherently. The model is both robust to the underlying plasma parameters and succeeds in reproducing a number of subtle observed features: (i) emission frequencies depend mostly on the length $λ_t$ of turbulence and the Lorentz factor of the reconnection generated beam, $ω\sim γ_b^2 ( c/λ_t) $ - it is independent of the absolute value of the underlying magnetic field. (ii) The model explains both broadband emission and the presence of emission stripes, including multiple stripes observed in the High Frequency Interpulse of the Crab pulsar. (iii) The model reproduces correlated polarization properties: presence of narrow emission bands in the spectrum favors linear polarization, while broadband emission can have arbitrary polarization. (iv) The mechanism is robust to the momentum spread of the particle in the beam. We also discuss a model of wigglers as non-linear force-free Alfven solitons (light darts).

preprint2022arXiv

Faraday conversion in pair-symmetric winds of magnetars and Fast Radio Bursts

We consider propagation of polarization in the inner parts of pair-symmetric magnetar winds, close to the light cylinder. Pair plasmas in magnetic field is birefringent, a $\propto B^2$ effect. As a result, such plasmas work as phase retarders: Stokes parameters follow a circular trajectory on the Poincare sphere. In the highly magnetized regime, $ω, \, ω_p \ll ω_B$, the corresponding rotation rates are independent of the magnetic field. A plasma screen with dispersion measure DM $\sim 10^{-6}$ pc cm$^{-3}$ can induce large polarization changes, including large effective Rotation Measure (RM). The frequency scaling of the (generalized) RM, $ \propto λ^α$, mimics the conventional RM with $α=2$ for small phase shifts, but can be as small as $α=1$. In interpreting observations the frequency scaling of polarization parameters should be fitted independently. The model offers explanations for (i) large circular polarization component observed in FRBs, with right-left switching; (ii) large RM, with possible sign changes; (iii) time-depend variable polarization. Relatively dense and slow wind is needed - the corresponding effect in regular pulsars is small.

preprint2022arXiv

Faraday rotation in fast radio bursts

Fast Radio Bursts (FBRs) show highly different polarization properties: high/small RMs, high/small circular/linear fractions. We outline a complicated picture of polarization propagation in the inner parts of the magnetars' winds, at scales $\sim$ few to hundreds of light cylinder radii. The key point is the Faraday rotation of linear polarization in highly magnetized symmetric pair plasma, a $\propto B^2$ effect. Position angle (PA) rotation rate is maximal for propagation across the magnetic field and disappears only for parallel propagation. In the highly magnetized regime, $ω\ll ω_B$, it becomes independent of the magnetic field. Very specific properties of PA($λ$) (scaling of the rotation angle with the observed wavelength $λ$) can help identify/sort out the propagation effects. Two basic regimes in pair plasma predict PA $\propto λ$ and $\propto λ^3$ (depending on the magnetic dominance); both are different from the conventional plasma's PA = RM $ λ^2$. This is the main prediction of the model. A number of effects, all sensitive to the underlying parameters, contribute to the observed complicated polarization patterns: streaming of plasma along magnetic field lines near the light cylinder, Faraday depolarization, effects of limiting polarization, the associated effect of linear-circular conversion, and synchrotron absorption.

preprint2022arXiv

FRB-periodicity: mild pulsars in tight O/B-star binaries

Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $η= L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $η\leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.

preprint2022arXiv

On the nature of Fast Blue Optical Transients

Short rise times of Fast Blue Optical Transients (FBOTs) require very light ejected envelopes, $M_{ej} \leq 10^{-1} M_\odot$, much smaller than of a typical supernova. Short peak times also mean that FBOTs should be hydrodynamically, not radioactively powered. The detection by Chandra of X-ray emission in AT2020mrf of $L_X \sim 10^{42} $ erg s$^{-1}$ after 328 days implies total, overall dominant, X-ray energetics at the Gamma Ray Bursts (GRBs) level of $\sim 6 \times 10^{49}$ erg. We further develop a model of Lyutikov & Toonen (2019), whereby FBOTs are the results of a late accretion induced collapse (AIC) of the product of super-Chandrasekhar double white dwarf (WD) merger between ONeMg WD and another WD. Small ejecta mass, and the rarity of FBOTs, result from the competition between mass loss from the merger product to the wind, and ashes added to the core, on time scale of $\sim 10^3-10^4$ years. FBOTs occur only when the envelope mass before AIC is $\leq 10^{-1} M_\odot$. FBOTs proper come from central engine-powered radiation-dominated forward shock as it propagates through ejecta. FBOTs' duration is determined by the diffusion time of photons produced by the NS-driven forward shock within the expanding ejecta. All the photons produced by the central source deep inside the ejecta escape almost simultaneously, producing a short bright event. The high energy emission is generated at the highly relativistic and highly magnetized termination shock, qualitatively similar to Pulsar Wind Nebulae. The X-ray bump observed in AT2020mrf by SRG/eROSITA, predicted by Lyutikov & Toonen (2019), is coming from the break-out of the engine-powered shock from the ejecta into the preceding wind. The model requires total energetics of just few $\times 10^{50}$ ergs, slightly above the observed X-rays. We predict that the system is hydrogen poor.

preprint2022arXiv

Radius-to-frequency mapping and FRB frequency drifts

We build a model of radius-to-frequency mapping in magnetospheres of neutron stars and apply it to frequency drifts observed in Fast Radio Bursts. We assume that an emission patch propagates along the dipolar magnetic field lines producing coherent emission with frequency, direction and polarization defined by the local magnetic field. The observed temporal evolution of the frequency depends on relativistic effects of time contraction and the curvature of the magnetic field lines. The model generically produces linear scaling of the drift rate, $\dotω \propto - ω$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dotω $ is also possible. Fast rotating magnetospheres produce higher drifts rates for similar viewing parameters than the slowly rotating ones. In the case of repeaters same source may show variable drift pattens depending on the observing phase. We expect rotational of polarization position angle through a burst, though by smaller amount than in radio pulsars. All these findings compare favorably with properties of FBRs, strengthening their possible loci in the magnetospheres of neutron stars.

preprint2022arXiv

Relativistic magnetic explosions

Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using analytical, relativistic MHD as well as force-free calculations, we find that the creation of a relativistically expanding, causally disconnected flow obeys a threshold condition: it requires sufficiently high initial over-pressure and sufficiently quick decrease of the pressure in the external medium (the pre-explosion wind). In the subcritical case the magnetic cloud just "puffs-up" and quietly expands with the pre-flare wind. We also find a compact analytical solution to the Prendergast's problem - expansion of force-free plasma into vacuum.

preprint2021arXiv

Dynamics and emission of wind-powered afterglows of gamma-ray bursts: flares, plateaus and steep decays

We develop a model of early X-ray afterglows of gamma-ray bursts originating from the reverse shock (RS) propagating through ultra-relativistic, highly magnetized pulsar-like winds produced by long-lasting central engines. We first perform fluid and MHD numerical simulations of relativistic double explosions. We demonstrate that even for constant properties of the wind a variety of temporal behaviors can be produced, depending on the energy of the initial explosion and the wind power, the delay time for the switch-on of the wind, and magnetization of the wind. X-ray emission of the highly magnetized RS occurs in the fast cooling regime - this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) RS emission naturally produces light curves showing power-law temporal evolution with various temporal indices; (ii) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$ (equivalent isotropic), can reproduce the afterglows' plateau phase; (iii) termination of the wind can produce sudden steep decays; (iv) short-duration afterglow flares are due to mild variations in the wind luminosity, with small total injected energy.

preprint2021arXiv

Rotating neutron stars without light cylinders

We find a class of twisted and differentially rotating neutron star magnetospheres that do not have a light cylinder, generate no wind and thus do not spin-down. The magnetosphere is composed of embedded differentially rotating flux surfaces, with the angular velocity decreasing as $Ω\propto 1/r$ (equivalently, becoming smaller at the foot-points closer to the axis of rotation). For each given North-South self-similar twist profile there is a set of self-similar angular velocity profiles (limited from above) with a "smooth", dipolar-like magnetic field structure extending to infinity. For spin parameters larger than some critical value, the light cylinder appears, magnetosphere opens up, and the wind is generated.

preprint2020arXiv

Coherence constraints on physical parameters at bright radio sources and FRB emission mechanism

We discuss physical constrains that observations of high brightness temperature coherent radio emission, with brightness temperatures as high as $T_b \sim 10^{35}$ K, impose on the plasma parameters at relativistically moving astrophysical sources. High brightness temperatures imply a minimal plasma energy density at the source. Additional important constraints come from the fact that resonantly emitting particles lose most of their energy to non-resonant inverse Compton and synchrotron processes. We also interpret recent observations of high-to-low frequency drifting features in the spectra of repeating FRBs as analogues of type-III Solar radio bursts produced by reconnection plasma beams within magnetospheres of highly magnetized neutron stars.

preprint2020arXiv

Fast moving pulsars as probes of interstellar medium

Pulsars moving through ISM produce bow shocks detected in hydrogen H$α$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$α$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $δρ/ρ\sim 1$, leads to significant variations of the shape of the shock seen in H$α$ line emission. We successfully reproduce the morphology of the Guitar Nebula. We infer quasi-periodic density variations in the warm component of ISM with a characteristic length of $\sim0.1$~pc. Structures of this scale might be also responsible for the formation of the fine features seen at the forward shock of Tycho SNR in X-rays. Formation of such short periodic density structures in the warm component of ISM is puzzling, and bow-shock nebulae provide unique probes to study this phenomenon.

preprint2020arXiv

Fast Radio Bursts from reconnection events in magnetar magnetospheres

Lyutikov (2002) predicted "radio emission from soft gamma-ray repeaters (SGRs) during their bursting activity". Detection of a Mega-Jansky radio burst in temporal coincidence with high energy bursts from a Galactic magnetar SGR 1935+2154 confirms that prediction. Similarity of this radio event with Fast Radio Bursts (FRBs) suggests that FRBs are produced within magnetar magnetospheres. We demonstrate that SGR 1935+2154 satisfies the previously derived constraints on the physical parameters at the FRBs' loci. Coherent radio emission is generated in the inner parts of the magnetosphere at $r< 100 R_{\rm NS}$. The radio emission is produced by the yet unidentified plasma emission process, occurring during the initial stages of reconnection events.

preprint2020arXiv

Free electron laser in magnetars/Fast Radio Bursts

We discuss coherent free electron laser (FEL) operating during explosive reconnection events in magnetized pair plasma of magnetar magnetospheres. The model explains many salient features of Fast Radio Bursts/magnetars' radio emission: temporal coincidence of radio and high energy bursts, high efficiency of conversion of plasma kinetic energy into coherent radiation, presence of variable, narrow-band emission features drifting down in frequency, high degree of linear polarization. The model relies on magnetar-specific drifting $e^\pm$ plasma components (which generate wiggler field due to the development of the firehose instability) and the presence of reconnection-generated particle beam with mild Lorentz factor of $γ_b \sim$ few hundred.

preprint2020arXiv

Magnetospheric interaction in white dwarf binaries AR Sco and AE Aqr

We develop a model of the white dwarf (WD) - red dwarf (RD) binaries AR Sco and AE Aqr as systems in a transient propeller stage of highly asynchronous intermediate polars. The WDs are relatively weakly magnetized with magnetic field of $\sim 10^6$ G. We explain the salient observed features of the systems due to the magnetospheric interaction of two stars. Currently, the WD's spin-down is determined by the mass loading of the WD's magnetosphere from the RD's at a mild rate of $\dot{M}_{WD} \sim 10^{-11} M_\odot $/yr. Typical loading distance is determined by the ionization of the RD's wind by the WD's UV flux. The WD was previously spun up by a period of high accretion rate from the RD via Roch lobe overflow with $\dot{M} \sim 10^{-9} M_\odot $/yr, acting for as short a period as tens of thousands of years. The non-thermal X-ray and optical synchrotron emitting particles originate in reconnection events in the magnetosphere of the WD due to the interaction with the flow from the RD. In the case of AR Sco, the reconnection events produce signals at the WD's rotation and beat periods - this modulation is due to the changing relative orientation of the companions' magnetic moments and resulting variable reconnection conditions. Radio emission is produced in the magnetosphere of the RD, we hypothesize, in a way that it is physically similar to the Io-induced Jovian decametric radiation.

preprint2020arXiv

Nonlinear optics in strongly magnetized pair plasma, with applications to FRBs

Intense radiation field can modify plasma properties, the corresponding refractive index, and lead to such nonlinear propagation effects as self-focusing. We estimate the corresponding effects in pair plasma, both in unmagnetized and strongly magnetically dominated case. First, in the unmagnetized pair plasma the ponderomotive force does not lead to charge separation, but to density depletion. Second, for astrophysically relevant plasmas of pulsar magnetospheres, (and possible loci of Fast Radio Bursts), where cyclotron frequency $ω_B$ dominates over plasma frequency $ω_p$ and the frequency of the electromagnetic wave, $ω_B \gg ω_p,\, ω$, we show that (i) there is virtually no nonlinearity due to changing effective mass in the field of the wave; (ii) ponderomotive force is $F_p^{(B)} =- {m_e c^2}/({4 B_0^2}) \nabla E^2$; it is reduced by a factor $(ω/ω_B)^2$ if compared to the unmagnetized case ($B_0$ is the external magnetic field and $E$ is the electric field of the wave); (iii) for radiation beam propagating along constant magnetic field in pair plasma with density $n_\pm$, the ponderomotive force leads to appearance of circular currents that lead to the decrease of the field within the beam by a factor $ΔB/B_0 = 2πn_\pm m_e c^2 {E^2}/{B_0^4}$. Applications to the physics of FRBs are discussed; we conclude that for parameters of FRB's the dominant magnetic field completely suppresses nonlinear radiation effects.

preprint2020arXiv

Peeking Between the Pulses: The Far-UV Spectrum of the Previously Unseen White Dwarf in AR Scorpii

The compact object in the interacting binary AR Sco has widely been presumed to be a rapidly rotating, magnetized white dwarf (WD), but it has never been detected directly. Isolating its spectrum has proven difficult because the spin-down of the WD generates pulsed synchrotron radiation that far outshines the WD's photosphere. As a result, a previous study of AR Sco was unable to detect the WD in the averaged far-ultraviolet spectrum from a Hubble Space Telescope (HST) observation. In an effort to unveil the WD's spectrum, we reanalyze these HST observations by calculating the average spectrum in the troughs between synchrotron pulses. We identify weak spectral features from the previously unseen WD and estimate its surface temperature to be 11500$\pm$500K. Additionally, during the synchrotron pulses, we detect broad Lyman-$α$ absorption consistent with hot WD spectral models. We infer the presence of a pair of hotspots, with temperatures between 23000K and 28000K, near the magnetic poles of the WD. As the WD is not expected to be accreting from its companion, we describe two possible mechanisms for heating the magnetic poles. The Lyman-$α$ absorption of the hotspots appears relatively undistorted by Zeeman splitting, constraining the WD's field strength to be 100 MG, but the data are insufficient to search for the subtle Zeeman splits expected at lower field strengths.

preprint2020arXiv

Tilting Instability of Magnetically Confined Spheromaks

We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfvén time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfvén waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars' giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.

preprint2020arXiv

Turbulent model of Crab nebula radiation

We construct a turbulent model of the Crab Nebula's non-thermal emission. The present model resolves a number of long-standing problems of the Kennel-Coroniti (1984) model: (i) the sigma problem; (ii) the hard spectrum of radio electrons; (iii) the high peak energy of gamma-ray flares; (iv) and the spacial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I accelerated at the wind termination shock via Fermi-I mechanism, and Component-II accelerated in reconnecting turbulence in highly magnetized ($σ$ $\gg 1$) plasma in the central part of the Crab Nebula. The reconnecting turbulence Component-II extends from radio to gamma rays: it accelerate radio electrons with a hard spectrum, destroy the large scale magnetic flux (and thus resolves the sigma-problem), and occasionally produces gamma-ray flares (from the largest scale reconnection events). The model reproduces the broad-band spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as spatially resolved evolution of the spectral indices in IR and optical bands.

preprint2020arXiv

Wind-powered afterglows of gamma-ray bursts: flares, plateaus and steep decays

Afterglows of gamma-ray bursts often show flares, plateaus, and sudden intensity drops: these temporal features are difficult to explain as coming from the forward shock. We calculate radiative properties of early GRB afterglows with the dominant contribution from the reverse shock (RS) propagating in an ultra-relativistic (pulsar-like) wind produced by the long-lasting central engine. RS emission occurs in the fast cooling regime -- this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$, can reproduce the afterglows' plateau phase; (ii) termination of the wind can produce sudden steep decays; (iii) mild variations in the wind luminosity can produce short-duration afterglow flares.

preprint2019arXiv

Conditions for jet break-out in NS mergers

We consider conditions for jet break-out through ejecta following mergers of neutron stars and provide simple relations for the break out conditions. We demonstrate that: (i) break-out requires that the isotropic-equivalent jet energy $E_j$ exceeds the ejecta energy $E_{ej}$ by $E_j \geq E_{ej}/ β_0$, where $ β_0 = V_{ej}/c$, $V_{ej}$ is the maximum velocity of the ejecta. If the central engine terminates before the break out, the shock approaches the edge of the ejecta slowly $\propto 1/t$; late break out occurs only if at the termination moment the head of the jet was relatively close to the edge. (ii) If there is a substantial delay between the ejecta's and the jet's launching, the requirement on the jet power increases. (iii) The forward shock driven by the jet is mildly strong, with Mach number $M\approx 5/4$ (increasing with time delay $t_d$); (iii) the delay time $t_d$ between the ejecta and the jet's launching is important for $t_d > t_0= ({3 }/{16} ) {c M_{ej} V_{ej}}/{L_j} = 1.01 {\rm sec} M_{ej, -2} L_{j, 51} ^{-1} \left ( {β_{ej}} /{0.3} \right)$, where $M_{ej}$ is ejecta mass, $L_j$ is the jet luminosity (isotropic equivalent). For small delays, $t_0 $ is also an estimate of the break-out time.

preprint2019arXiv

Nonlinear force-free configurations in cylindrical geometry

We find a new family of solutions for force-free magnetic structures in cylindrical geometry. These solutions have radial power-law dependance and are periodic but non-harmonic in azimuthal direction; they generalize the vacuum $z$-independent potential fields to current-carrying configurations.

preprint2016arXiv

Emission knots and polarization swings of swinging jets

Knots (emission features in jets of active galactic nuclei) often show non-ballistic dynamics and variable emission/polarization properties. We model these features as emission pattern propagating in a jet that carries helical magnetic field and is launched along a changing direction. The model can reproduce a wide range of phenomena observed in the motion of knots: non-ballistic motion (both smooth and occasional sudden change of direction, and/or oscillatory behavior), variable brightness, confinement of knots' motion within an overlaying envelope. The model also reproduces smooth large polarization angle swings, and at the same time allows for the seemingly random behavior of synchrotron fluxes, polarization fraction and occasional $π/2$ polarization jumps.

preprint2016arXiv

Fast radio bursts as giant pulses from young rapidly rotating pulsars

We discuss possible association of fast radio bursts (FRBs) with supergiant pulses emitted by young pulsars (ages $\sim$ tens to hundreds of years) born with regular magnetic field but very short -- few milliseconds -- spin periods. FRBs are extra-Galactic events coming from distances $d \lesssim 100$ Mpc. Most of the dispersion measure (DM) comes from the material in the freshly ejected SNR shell; for a given burst the DM should decrease with time. FRBs are not expected to be seen below $\sim 300 $ MHz due to free-free absorption in the expanding ejecta. A supernova might have been detected years before the burst; FRBs are mostly associated with star forming galaxies. The model requires that some pulsars are born with very fast spins, of the order of few milliseconds. The observed distribution of spin-down powers $\dot{E}$ in young energetic pulsars is consistent with equal birth rate per decade of $\dot{E}$. Accepting this injection spectrum and scaling the intrinsic brightness of FRBs with $\dot{E}$, we predict the following properties of a large sample of FRBs: (i) the brightest observed events come from a broad distribution in distances; (ii) for repeating bursts brightness either remains nearly constant (if the spin-down time is longer than the age of the pulsar) or decreases with time otherwise; in the latter case DM $\propto \dot{E}$.

preprint2016arXiv

Fermi GBM signal contemporaneous with GW150914 - an unlikely association

We argue that the physical constraints required by the association of the Fermi GBM signal contemporaneous with GW150914 - radiative power of $10^{49} $ erg s$^{-1}$, and corresponding magnetic fields on the black hole of the order of $10^{12}$ Gauss - are astrophysical highly implausible. Combined with the relatively high random probability of coincidence of 0.22 percents, we conclude that the electromagnetic signal is likely unrelated to the BH merger.

preprint2016arXiv

How else can we detect Fast Radio Bursts?

We discuss possible electromagnetic signals accompanying Fast Radio Bursts (FRBs) that are expected in the scenario where FRBs originate in neutron star magnetospheres. For models involving Crab-like giant pulses, no appreciable contemporaneous emission is expected at other wavelengths. Magnetar giant flares, driven by the reconfiguration of the magnetosphere, however, can produce both contemporaneous bursts at other wavelengths as well as afterglow-like emission. We conclude that the best chances are: (i) prompt short GRB-like emission; (ii) a contemporaneous optical flash that can reach naked eye peak luminosity (but only for a few milliseconds); (iii) a high energy afterglow emission. Case (i) could be tested by coordinated radio and high-energy experiments. Case (ii) could be seen in a coordinated radio-optical surveys, \eg\ by the Palomar Transient Factory in a 60-second frame as a transient object of $m=15-20$ magnitude with an expected optical detection rate of about 0.1~hr$^{-1}$, an order of magnitude higher than in radio. Shallow, but large-area sky surveys such as ASAS-SN and EVRYSCOPE could also detect prompt optical flashes from the more powerful Lorimer-burst clones. The best constraints on the optical-to-radio power for this kind of emission could be provided by future observations with facilities like LSST. Case (iii) might be seen in relatively rare cases that the relativistically ejected magnetic blob is moving along the line of sight.

preprint2016arXiv

On the Linear Stability of Sheared and Magnetized Jets Without Current Sheets - Non-Relativistic Case

In a prior paper (Kim et al. 2015) we considered the linear stability of magnetized jets that carry no net electric current and do not have current sheets. In this paper, in addition to physically well-motivated magnetic field structures, we also include the effects of jet shear. The jets we study have finite thermal pressure in addition to having realistic magnetic field structures and velocity shear. We find that shear has a strongly stabilizing effect on various modes of jet instability. Increasing shear stabilizes the fundamental pinch modes at long wavelengths and short wavelengths. Increasing shear also stabilizes the first reflection pinch modes at short wavelengths. Increasing shear has only a very modest stabilizing effect on the fundamental kink modes at long wavelengths; however, increasing shear does have a strong stabilizing effect on the fundamental kink modes at short wavelengths. The first reflection kink modes are strongly stabilized by increasing shear at shorter wavelengths. Overall, we find that the combined effect of magnetic field and shear stabilizes jets more than shear alone. In addition to the results from a formal linear stability analysis, we present a novel way of visualizing and understanding jet stability. This gives us a deeper understanding of the enhanced stability of sheared, magnetized jets. We also emphasize the value of our numerical approach in understanding the linear stability of jets with realistic structure.

preprint2016arXiv

Relativistic Rotating Vector Model

The direction of polarization produced by a moving source rotates with the respect to the rest frame. We show that this effect, induced by pulsar rotation, leads to an important correction to polarization swings within the framework of rotating vector model (RVM); this effect has been missed by previous works. We construct relativistic RVM taking into account finite heights of the emission region that lead to aberration, time-of-travel effects and relativistic rotation of polarization. Polarizations swings at different frequencies can be used, within the assumption of the radius-to-frequency mapping, to infer emission radii and geometry of pulsars.

preprint2015arXiv

Magnetar Activity via the Density-Shear Instability in Hall-MHD

We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing electron number density and magnetic field, a plane-parallel analogue of an azimuthal or meridional magnetic field in the crust of a magnetar. The growth rate of the instability depends on the Hall properties of the field (magnetic field intensity, electron number density and the corresponding scale-heights), while being insensitive to weak resistivity. Since the Hall effect is the driving process for the evolution of the crustal magnetic field of magnetars, we argue that this instability is critical for systems containing strong meridional or azimuthal fields. We find that this process mediates the formation of localised structures with much stronger magnetic field than the average, which can lead to magnetar activity and accelerate the dissipation of the field and consequently the production of Ohmic heating. Assuming a $5\times10^{14}$G magnetic field at the base of crust, we anticipate that magnetic field as strong as $10^{15}$G will easily develop in regions of typical size of a few $10^{2}$ meters, containing magnetic energy of $10^{43}$erg, sufficient to power magnetar bursts. These active regions are more likely to appear in the magnetic equator where the tangential magnetic field is stronger.

preprint2015arXiv

On the Linear Stability of Magnetized Jets Without Current Sheets: Non-Relativistic Case

In this paper we consider stability of magnetized jets that carry no net electric current and do not have current sheets (Gourgouliatos et al. 2012). The non-relativistic MHD equations are linearized around the background velocity and the magnetic field structure of the jet. The resulting linear equations are solved numerically inside the jet. We find that introduction of current-sheet-free magnetic field significantly improves jet stability relative to unmagnetized jets or magnetized jets with current sheets at their surface. This particularly applies to the fundamental pinch and kink modes - they become completely suppressed in a wide range of long wavelengths that are known to become most pernicious to jet stability when the evolution enters the non-linear regime. The reflection modes, both for the pinch and kink instability, also become progressively more stable with increased magnetization.

preprint2015arXiv

The inner knot of the Crab nebula

We model the inner knot of the Crab Nebula as a synchrotron emission coming from the non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock appearance. The model can reproduce the knot location, size, elongation, brightness distribution, luminosity and polarization provided the effective magnetization of the section of the pulsar wind producing the knot is low, $σ\leq 1$. In the striped wind model, this implies that the striped zone is rather wide, with the magnetic inclination angle of the Crab pulsar $\ge 45^\circ$; this agrees with the previous model-dependent estimate based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock and that this shock is a site of efficient particle acceleration. On the other hand, the deduced low magnetization of the knot plasma implies that this is an unlikely site for the Crab's gamma-ray flares, if they are related to the fast relativistic magnetic reconnection events.

preprint2015arXiv

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

preprint2014arXiv

Density-shear instability in electron MHD

We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growth-rate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is "at least as stable" as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD.

preprint2014arXiv

Magnetar activity mediated by plastic deformations of neutron star crust

We advance a "Solar flare" model of magnetar activity, whereas a slow evolution of the magnetic field in the upper crust, driven by electron MHD (EMHD) flows, twists the external magnetic flux tubes, producing persistent emission, bursts and flares. At the same time the neutron star crust plastically relieves the imposed magnetic field stress, limiting the strain $ ε_t $ to values well below the critical strain $ ε_{crit}$ of a brittle fracture, $ ε_t \sim 10^{-2}ε_{crit} $. Magnetar-like behavior, occurring near the magnetic equator, takes place in all neutron stars, but to a different extent. The persistent luminosity is proportional to cubic power of the magnetic field (at a given age), and hence is hardly observable in most rotationally powered neutron stars. Giant flares can occur only if the magnetic field exceeds some threshold value, while smaller bursts and flares may take place in relatively small magnetic fields. Bursts and flares are magnetospheric reconnection events that launch Alfven shocks which convert into high frequency whistlers upon hitting the neutron star surface. The resulting whistler pulse induces a strain that increases with depth both due to the increasing electron density (and the resulting slowing of the waves), and due to the increasing coherence of a whistler pulse with depth. The whistler pulse is dissipated on a time scale of approximately a day at shallow depths corresponding to $ρ\sim 10^{10} {\rm g cm}^{-3}$; this energy is detected as enhanced post-flare surface emission.

preprint2014arXiv

Super-Eddington accretion in ultra-luminous neutron star binary

We discuss properties of the ultra-luminous $X$-ray source in the galaxy M82, NuSTAR J095551+6940.8, containing an accreting neutron star. The neutron star has surface magnetic field $ B_{NS} \approx 1.4 \times 10^{13 } \, {\rm G}$ and experiences accretion rate of $9 \times 10^{-7} M_\odot {\rm \, yr}^{-1} $. The magnetospheric radius, close to the corotation radius, is $\sim 2 \times 10^8$ cm. The accretion torque on the neutron star is reduce well below what is expected in a simple magnetospheric accretion due to effective penetration of the stellar magnetic field into the disk beyond the corotation radius. As a result, the radiative force of the surface emission does not lead to strong coronal wind, but pushes plasma along magnetic field lines towards the equatorial disk. The neutron star is nearly an orthogonal rotator, with the angle between the rotation axis and the magnetic moment $\geq 80$ degrees. Accretion occurs through optically thick -- geometrically thin and flat accretion "curtain", which cuts across the polar cap. High radiation pressure from the neutron star surface is nevertheless smaller than that the ram pressure of the accreting material flowing through the curtain, and thus fails to stop the accretion. At distances below few stellar radii the magnetic suppression of the scattering cross-section becomes important. The $X$-ray luminosity (pulsed and persistent components) comes both from the neutron star surface as a hard $X$-ray component and as a soft component from reprocessing by the accretion disk.

preprint2013arXiv

A Giant Sample of Giant Pulses from the Crab Pulsar

We observed the Crab pulsar with the 43-m telescope in Green Bank, WV over a timespan of 15 months. In total we obtained 100 hours of data at 1.2 GHz and seven hours at 330 MHz, resulting in a sample of about 95000 giant pulses (GPs). This is the largest sample, to date, of GPs from the Crab pulsar taken with the same telescope and backend and analyzed as one data set. We calculated power-law fits to amplitude distributions for main pulse (MP) and interpulse (IP) GPs, resulting in indices in the range of 2.1-3.1 for MP GPs at 1.2 GHz and in the range of 2.5-3.0 and 2.4-3.1 for MP and IP GPs at 330 MHz. We also correlated the GPs at 1.2 GHz with GPs from the Robert C. Byrd Green Bank Telescope (GBT), which were obtained simultaneously at a higher frequency (8.9 GHz) over a span of 26 hours. In total, 7933 GPs from the 43-m telescope at 1.2 GHz and 39900 GPs from the GBT were recorded during these contemporaneous observations. At 1.2 GHz, 236 (3%) MP GPs and 23 (5%) IP GPs were detected at 8.9 GHz, both with zero chance probability. Another 15 (4%) low-frequency IP GPs were detected within one spin period of high-frequency IP GPs, with a chance probability of 9%. This indicates that the emission processes at high and low radio frequencies are related, despite significant pulse profile shape differences. The 43-m GPs were also correlated with Fermi gamma-ray photons to see if increased pair production in the magnetosphere is the mechanism responsible for GP emission. A total of 92022 GPs and 393 gamma-ray photons were used in this correlation analysis. No significant correlations were found between GPs and gamma-ray photons. This indicates that increased pair production in the magnetosphere is likely not the dominant cause of GPs. Possible methods of GP production may be increased coherence of synchrotron emission or changes in beaming direction.

preprint2013arXiv

Afterglow jet breaks, GeV photons and the electromagnetic model of short GRBs

We discuss three topics: (i) the dynamics of afterglow jet breaks; (ii) the origin of Fermi-LAT photons; (iii) the electromagnetic model of short GRBs

preprint2013arXiv

Electron MHD: dynamics and turbulence

(Abridged) We consider dynamics and turbulent interaction of whistler modes within the framework of inertialess electron MHD (EMHD). We argue there is no energy principle in EMHD: any stationary closed configuration is neutrally stable. We consider the turbulent cascade of whistler modes. We show that (i) harmonic whistlers are exact non-linear solutions; (ii) co-linear whistlers do not interact (including counter-propagating); (iii) waves with the same value of the wave vector, $k_1=k_2$, do not interact; (iv) whistler modes have a dispersion that allows a three-wave decay, including into a zero frequency mode; (v) the three-wave interaction effectively couples modes with highly different wave numbers and propagation angles. In addition, linear interaction of a whistler with a single zero-mode can lead to spatially divergent structures via parametric instability. All these properties are drastically different from MHD, so that the qualitative properties of the Alfven turbulence cannot be transferred to the EMHD turbulence. We derive the Hamiltonian formulation of EMHD, and using Bogolyubov transformation reduce it to the canonical form; we calculate the matrix elements for the three-wave interaction of whistlers. We solve numerically the kinetic equation and show that, generally, the EMHD cascade develops within a broad range of angles, while transiently it may show anisotropic, nearly two dimensional structures. Development of a cascade depends on the forcing (non-universal) and often fails to reach a steady state. Analytical estimates predict the spectrum of magnetic fluctuations for the quasi-isotropic cascade $\propto k^{-2}$. The cascade remains weak (not critically-balanced). The cascade is UV-local, while the infrared locality is weakly (logarithmically) violated.

preprint2013arXiv

I. Inverse Compton origin of pulsar γ-ray emission. II. Reconnection model of Crab flares

I. There is growing evidence that pulsars' high energy emission is generated via Inverse Compton mechanism. II. The particles producing Crab flares, and possibly most of the Crab Nebula's high energy emission, are accelerated via reconnection events, and not at shock via Fermi mechanisms.

preprint2013arXiv

Magnetospheric "anti-glitches" in magnetars

We attribute the rapid spindown of magnetar 1E 2259+586 observed by Archibald et al. (2013), termed the "anti-glitch", to partial opening of the magnetosphere during the X-ray burst, followed by changes of the structure of the closed field line region. To account for the observed spin decrease during the X- ray flare all that is needed is the transient opening, for just one period, of a relatively small fraction of the magnetosphere, of the order of only few percent. More generally, we argue that in magnetars all timing irregularities have magnetospheric origin and are induced either by (i) the fluctuations in the current structure of the magnetosphere (similar to the long term torque variations in the rotationally powered pulsars); or, specifically to magnetars, by (ii) opening of a fraction of the magnetosphere during bursts and flares - the latter events are always accompanied by rapid spindown, an "anti-glitch". Slow rotational motion of the neutron star crust, driven by crustal magnetic fields, leads beyond some twist limit to explosive instability of the external magnetic fields and transient opening of a large magnetic flux in a CME-type event, the post-flares increase of magnetospheric currents accompanied by enhanced X-luminosity and spindown rate, changing profiles, as well as spectral hardening, all in agreement with the magnetospheric model of torque fluctuations.

preprint2013arXiv

Radio Emission Region Exposed: Courtesy of the Double Pulsar

The double pulsar system PSR J0737-3039A/B offers exceptional possibilities for detailed probes of the structure of the pulsar magnetosphere, pulsar winds and relativistic reconnection. We numerically model the distortions of the magnetosphere of pulsar B by the magnetized wind from pulsar A, including effects of magnetic reconnection and of the geodetic precession. Geodetic precession leads to secular evolution of the geometric parameters and effectively allows a 3D view of the magnetosphere. Using the two complimentary models of pulsar B's magnetosphere, adapted from the Earth's magnetosphere models by Tsyganenko (ideal pressure confinement) and Dungey (highly resistive limit), we determine the precise location and shape of the coherent radio emission generation region within pulsar B's magnetosphere. We successfully reproduce orbital variations and secular evolution of the profile of B, as well as subpulse drift (due to reconnection between the magnetospheric and wind magnetic fields), and determine the location and the shape of the emission region. The emission region is located at about 3750 stellar radii and has a horseshoe-like shape, which is centered on the polar magnetic field lines. The best fit angular parameters of the emission region indicate that radio emission is generated on the field lines which, according to the theoretical models, originate close to the poles and carry the maximum current. When considered together, the results of the two models converge and can explain why the modulation of B's radio emission at A's period is observed only within a certain orbital phase region. Our results imply that the wind of pulsar A has a striped structure only 1000 light cylinder radii away. We discuss the implications of these results for pulsar magnetospheric models, mechanisms of coherent radio emission generation, and reconnection rates in relativistic plasma.

preprint2013arXiv

Topics in microphysics of relativistic plasmas

Astrophysical plasmas can have parameters vastly different from the more studied laboratory and space plasmas. In particular, the magnetic fields can be the dominant component of the plasma, with energy-density exceeding the particle rest-mass energy density. Magnetic fields then determine the plasma dynamical evolution, energy dissipation and acceleration of non-thermal particles. Recent data coming from astrophysical high energy missions, like magnetar bursts and Crab nebula flares, point to the importance of magnetic reconnection in these objects. In this review we outline a broad spectrum of problems related to the astrophysical relevant processes in magnetically dominated relativistic plasmas. We discuss the problems of large scale dynamics of relativistic plasmas, relativistic reconnection and particle acceleration at reconnecting layers, turbulent cascade in force-fee plasmas. A number of astrophysical applications are also discussed.

preprint2012arXiv

Crab Nebula gamma-ray flares as relativistic reconnection minijets

The unusually short durations, high luminosities, and high photon energies of the Crab Nebula gamma-ray flares require relativistic bulk motion of the emitting plasma. We explain the Crab flares as the result of randomly oriented relativistic "minijets" originating from reconnection events in a magnetically dominated plasma. We develop a statistical model of the emission from Doppler boosted reconnection minijets and find analytical expressions for the moments of the resulting nebula light curve (e.g. time average, variance, skewness). The light curve has a flat power spectrum that transitions at short timescales to a decreasing power-law of index 2. The flux distribution from minijets follows a decreasing power-law of index ~ 1, implying the average flux from flares is dominated by bright rare events. The predictions for the flares' statistics can be tested against forthcoming observations. We find the observed flare spectral energy distributions (SEDs) have several notable features: A hard power-law index of p ~< 1 for accelerated particles that is expected in various reconnection models, including some evidence of a pile-up near the radiation reaction limit. Also, the photon energy at which the SED peaks is higher than that implied by the synchrotron radiation reaction limit, indicating the flare emission regions' Doppler factors are ~> few. We conclude that magnetic reconnection can be an important, if not dominant, mechanism of particle acceleration within the nebula.

preprint2012arXiv

Hair of astrophysical black holes

The "no hair" theorem is not applicable to black holes formed from collapse of a rotating neutron star. Rotating neutron stars can self-produce particles via vacuum breakdown forming a highly conducting plasma magnetosphere such that magnetic field lines are effectively "frozen-in" the star both before and during collapse. In the limit of no resistivity, this introduces a topological constraint which prohibits the magnetic field from sliding off the newly-formed event horizon. As a result, during collapse of a neutron star into a black hole, the latter conserves the number of magnetic flux tubes N_B = e Φ_\infty /(πc \hbar), where Φ_\infty is the initial magnetic flux through the hemispheres of the progenitor and out to infinity. The black hole's magnetosphere subsequently relaxes to the split monopole magnetic field geometry with self-generated currents outside the event horizon. The dissipation of the resulting equatorial current sheet leads to a slow loss of the anchored flux tubes, a process that balds the black hole on long resistive time scales rather than the short light-crossing time scales expected from the vacuum "no-hair" theorem.

preprint2012arXiv

IC model of pulsar high energy emission

We discuss growing evidence that pulsar high energy is emission is generated via Inverse Compton mechanism. We reproduce the broadband spectrum of Crab pulsar, from UV to very high energy gamma-rays - nearly ten decades in energy, within the framework of the cyclotron-self-Compton model. Emission is produced by two counter-streaming beams within the outer gaps, at distances above ~ 20 NS radii. The outward moving beam produces UV-X-ray photons via Doppler-booster cyclotron emission, and GeV photons by Compton scattering the cyclotron photons produced by the inward going beam. The scattering occurs in the deep Klein-Nishina regime, whereby the IC component provides a direct measurement of particle distribution within the magnetosphere. The required plasma multiplicity is high, ~ 10^6-10^7, but is consistent with the average particle flux injected into the pulsar wind nebula.

preprint2012arXiv

Inverse Compton model of pulsar high energy emission

We reproduce the broadband spectrum of Crab pulsar, from UV to very high energy gamma-rays - nearly ten decades in energy, within the framework of the cyclotron-self-Compton model. Emission is produced by two counter-streaming beams within the outer gaps, at distances above $\sim$ 20 NS radii. The outward moving beam produces UV-$X$-ray photons via Doppler-booster cyclotron emission, and GeV photons by Compton scattering the cyclotron photons produced by the inward going beam. The scattering occurs in the deep Klein-Nishina regime, whereby the IC component provides a direct measurement of particle distribution within the magnetosphere. The required plasma multiplicity is high, $\sim 10^6-10^7$, but is consistent with the average particle flux injected into the pulsar wind nebula. The importance of Compton scattering in the Klein-Nishina regime also implies the importance of pair production in the outer gaps. We suggest that outer gaps are important sources of pairs in pulsar magnetospheres. Cyclotron motion of particles in the pulsar magnetosphere may be excited due to coherent emission of radio waves by streaming particles at the anomalous cyclotron resonance. Thus, a whole range of Crab non-thermal emission, from coherent radio waves to very high energy $γ$-rays - nearly eighteen decades in energy - may be a manifestation of inter-dependent radiation processes. The present model, together with the observational evidence in favor of the IC scattering (Lyutikov et al. 2012; Lyutikov 2012), demonstrates that the inverse Compton scattering can be the dominant high energy emission mechanism in majority of pulsars.

preprint2012arXiv

Magnetic structure of Coronal Mass Ejections

We present several models of the magnetic structure of solar coronal mass ejections (CMEs). First, we model CMEs as expanding force-free magnetic structures. While keeping the internal magnetic field structure of the stationary solutions, expansion leads to complicated internal velocities and rotation, while the field structures remain force-free. Second, expansion of a CME can drive resistive dissipation within the CME changing the ionization states of different ions. We fit in situ measurements of ion charge states to the resistive spheromak solutions. Finally, we consider magnetic field structures of fully confined stable magnetic clouds containing both toroidal and poloidal magnetic fields and having no surface current sheets. Expansion of such clouds may lead to sudden onset of reconnection events.

preprint2012arXiv

On the dynamics of mechanical failures in magnetized neutron-star crusts

We consider the dynamics of a mechanical failure induced by a shear stress in a strongly magnetized neutron-star crust. We show that even if the elastic properties of the crust allow the creation of a shear crack, the strongly sheared magnetic field around the crack leads to a back-reaction from the Lorentz force which does not allow large relative displacement of the crack surfaces. Instead, the global evolution of the crack proceeds on a slow resistive time scale, and is unable to release any substantial mechanical energy. Our calculations demostrate that for {\it some} magnetic-field configurations, the magnetic forces cause, effectively, a plastic deformation of the crust when the resulting elastic shear stress exceeds the critical value for mechanical failure.

preprint2012arXiv

The electromagnetic model of short GRBs, the nature of prompt tails, supernova-less long GRBs and highly efficient episodic accretion

Many short GRBs show prompt tails lasting up to hundreds of seconds that can be energetically dominant over the initial sub-second spike. In this paper we develop an electromagnetic model of short GRBs that explains the two stages of the energy release, the prompt spike and the prompt tail. The key ingredient of the model is the recent discovery that an isolated black hole can keep its open magnetic flux for times much longer than the collapse time and, thus, can spin-down electromagnetically, driving the relativistic wind. First, the merger is preceded by an electromagnetic precursor wind. If a fraction of the wind power is converted into pulsar-like coherent radio emission, this may produce an observable radio burst of few milliseconds. At the active stage of the merger, two neutron stars produces a black hole surrounded by an accretion torus in which the amplified magnetic field extracts the rotational energy of the black hole and drives an axially-collimated electromagnetic wind. For observers nearly aligned with the orbital normal this is seen as a classical short GRB. After the accretion of the torus, the isolated black hole keeps the open magnetic flux and drives the equatorially (not axially) collimated outflow, which is seen by an observer at intermediate polar angles as a prompt tail. The tail carries more energy than the prompt spike, but its emission is de-boosted for observers along the orbital normal. Observers in the equatorial plane miss the prompt spike and interpret the prompt tail as an energetic supernova-less long GRB. We also demonstrate that episodic accretion onto the BH of magnetized clouds that carry non-zero magnetic flux can be highly efficient in extracting the spin energy of the BH, producing the outflows with the power exceeding the average accretion power.

preprint2012arXiv

The gamma-ray spectrum of Geminga and the inverse Compton model of pulsar high energy emission

We reanalyze the Fermi spectra of the Geminga and Vela pulsars. We find that the spectrum of Geminga above the break is exceptionally well approximated by a simple power law without the exponential cut-off, making Geminga's spectrum similar to that of Crab. Vela's broadband gamma-ray spectrum is equally well fit with both the exponential cut-off and the double power law shapes. In the broadband double power-law fits, for a typical Fermi spectrum of a bright γ-ray pulsar, most of the errors accumulate due to the arbitrary parametrization of the spectral roll-off. In addition, a power law with an exponential cut-off gives an acceptable fit for the underlying double power-law spectrum for a very broad range of parameters, making such fitting procedures insensitive to the underlying Fermi photon spectrum. Our results have important implications for the mechanism of pulsar high energy emission. A number of observed properties of γ-ray pulsars, i.e., the broken power law spectra without exponential cut-offs and stretching in case of Crab beyond the maximal curvature limit, spectral breaks close to or exceeding the maximal breaks due to curvature emission, a Crab patterns of relative intensities of the leading and trailing pulses repeated in the X-ray and γ-ray regions, all point to the inverse Compton origin of the high energy emission from majority of pulsars.

preprint2011arXiv

Coronal mass ejections as expanding force-free structures

We mode Solar coronal mass ejections (CMEs) as expanding force-fee magnetic structures and find the self-similar dynamics of configurations with spatially constant α, where {\bf J} =α{\bf B}, in spherical and cylindrical geometries, expanding spheromaks and expanding Lundquist fields correspondingly. The field structures remain force-free, under the conventional non-relativistic assumption that the dynamical effects of the inductive electric fields can be neglected. While keeping the internal magnetic field structure of the stationary solutions, expansion leads to complicated internal velocities and rotation, induced by inductive electric field. The structures depends only on overall radius R(t) and rate of expansion \dot{R}(t) measured at a given moment, and thus are applicable to arbitrary expansion laws. In case of cylindrical Lundquist fields, the flux conservation requires that both axial and radial expansion proceed with equal rates. In accordance with observations, the model predicts that the maximum magnetic field is reached before the spacecraft reaches the geometric center of a CME.

preprint2011arXiv

Crab GeV flares from corrugated termination shock

Very high energy gamma-ray flares from the Crab nebular detected by AGILE and Fermi satellites challenge our understanding of the pulsar wind nebulae. The short duration of the flares, only few days, is particularly puzzling since it is much shorter than the dynamical times scale of the nebular. In this work we investigate analytically and via numerical simulations the electromagnetic signatures expected from the large amplitude low frequency magnetosonic waves generated within the Crab nebular which induce the corrugation perturbations of the termination shock. As a result, the oblique termination shock produces time-dependent, mildly relativistic post-shock flow. Using the relativistic MHD version of the RIEMANN code, we simulate the interaction of the termination shock with downstream perturbations. We demonstrate that mild Doppler boosting of the synchrotron emission in the post-shock flow can produce bright, short time scale flares.

preprint2011arXiv

Dynamics of rising magnetized cavities and UHECR acceleration in clusters of galaxies

We study the expansion of low density cavities produced by Active Galactic Nuclei jets in clusters of galaxies. The long term stability of these cavities requires the presence of linked magnetic fields. We find solutions describing the self-similar expansion of structures containing large-scale electromagnetic fields. Unlike the force-free spheromak-like configurations, these solutions have no surface currents and, thus, are less susceptible to resistive decay. The cavities are internally confined by external pressure, with zero gradient at the surface. If the adiabatic index of the plasma within the cavity is $Γ>4/3$, the expansion ultimately leads to the formation of large-scale current sheets. The resulting dissipation of the magnetic field can only partially offset the adiabatic and radiative losses of radio emitting electrons. We demonstrate that if the formation of large-scale current sheets is accompanied by explosive reconnection of the magnetic field, the resulting reconnection layer can accelerate cosmic rays to ultra high energies. We speculate that the enhanced flux of UHECRs towards Centaurus A originates at the cavities due to magnetic reconnection.

preprint2011arXiv

Electromagnetic power of merging and collapsing compact objects

[Abridged] Electromagnetic emission can be produced as a precursor to the merger, as a prompt emission during the collapse of a NS and at the spin-down stage of the resulting BH. We demonstrate that the time evolution of the axisymmetric force-free magnetic fields can be expressed in terms of the hyperbolic Grad-Shafranov equation. We find exact non-linear time-dependent split-monopole structure of magnetosphere driven by spinning and collapsing NS in Schwarzschild geometry. Based on this solution, we argue that the collapse of a NS into the BH happens smoothly, without natural formation of current sheets or other dissipative structures on the open field lines and, thus, does not allow the magnetic field to become disconnected from the star and escape to infinity. Thus, as long as an isolated Kerr BH can produce plasma and currents, it does not lose its open magnetic field lines, its magnetospheric structure evolved towards a split monopole and the BH spins down electromagnetically. The "no hair theorem", which assumes that the outside medium is a vacuum, is not applicable in this case: highly conducting plasma introduces a topological constraint forbidding the disconnection of the magnetic field lines from the BH. Eventually, a single random large scale spontaneous reconnection event will lead to magnetic field release, shutting down the electromagnetic BH engine forever. We also discuss the nature of short Gamma Ray Bursts and suggest that the similarity of the early afterglows properties of long and short GRBs can be related to the fact that in both cases a spinning BH can retains magnetic field for sufficiently long time to extract a large fraction of its rotation energy and produce high energy emission via the internal dissipation in the wind.

preprint2011arXiv

In Situ Heating of the 2007 May 19 CME Ejecta Detected by STEREO/PLASTIC and ACE

In situ measurements of ion charge states can provide unique insight into the heating and evolution of coronal mass ejections when tested against realistic non-equilibrium ionization modeling. In this work we investigate the representation of the CME magnetic field as an expanding spheromak configuration, where the plasma heating is prescribed by the choice of anomalous resistivity and the spheromak dynamics. We chose as a test case, the 19 May 2007 CME observed by STEREO and ACE. The spheromak is an appealing physical model, because the location and degree of heating is fixed by the choice of anomalous resistivity and the spheromak expansion rate which we constrain with observations. This model can provide the heating required between 1.1$R_{\sun}$ and earth orbit to produce charge states observed in the CME flux rope. However this source of heating in the spheromak alone has difficulty accounting for the rapid heating to Fe$^{8 - 11+}$ at lower heights, as observed in STEREO EUVI due to the rapid radiative cooling that occurs at the high densities involved. Episodes of heating and cooling clearly unrelated to spheromak expansion are observed prior to the eruption, and presumably still play a role during the eruption itself. Spheromak heating is also not capable of reproducing the high Fe charge states (Fe$^{16+}$ and higher) seen in situ exterior to the flux rope in this CME. Thus while the spheromak configuration may be a valid model for the magnetic topology, other means of energization are still required to provide much of the rapid heating observed.

preprint2011arXiv

On the sideways expansion of relativistic non-spherical shocks and GRB afterglows

Expansion of non-spherical relativistic blast waves is considered in the Kompaneets (the thin shell) approximation. We find that the relativistic motion effectively "freezes out" the lateral dynamics of the shock front: only extremely strongly collimated shocks, with the opening angles $Δθ\leq 1/Γ^2$, show appreciable modification of profiles due to sideways expansion. For less collimated profiles the propagation is nearly ballistic; the sideways expansion of relativistic shock becomes important only when they become mildly relativistic.

preprint2011arXiv

Relativistic magnetohydrodynamics in one dimension

We derive a number of solution for one-dimensional dynamics of relativistic magnetized plasma that can be used as benchmark estimates in relativistic hydrodynamic and magnetohydrodynamic numerical codes. First, we analyze the properties of simple waves of fast modes propagating orthogonally to the magnetic field in relativistically hot plasma. The magnetic and kinetic pressures obey different equations of state, so that the system behaves as a mixture of gases with different polytropic indices. We find the self-similar solutions for the expansion of hot strongly magnetized plasma into vacuum. Second, we derive linear hodograph and Darboux equations for the relativistic Khalatnikov potential, which describe arbitrary one-dimensional isentropic relativistic motion of cold magnetized plasma and find their general and particular solutions. The obtained hodograph and Darboux equations are very powerful: system of highly non-linear, relativistic, time dependent equations describing arbitrary (not necessarily self-similar) dynamics of highly magnetized plasma reduces to a single linear differential equation.

preprint2011arXiv

Schwarzschild black holes as unipolar inductors: expected electromagnetic power of a merger

(Abridged) The motion of a Schwarzschild black hole with velocity $v_0 = β_0 c$ through a constant magnetic field $B_0$ in vacuum induces a component of the electric field along the magnetic field, generating a non-zero second Poincare electromagnetic invariant $ ^* F \cdot F \neq 0$. This will produce (e.g., via radiative effects and vacuum breakdown) an electric charge density of the order of $ρ_{\rm ind}= B_0 β_0 /(2 πe R_G)$, where $R_G = 2 G M/c^2$ is the Schwarzschild radius and $M$ is the mass of the black hole; the charge density $ρ_{\rm ind}$ is similar to the Goldreich-Julian density. The magnetospheres of moving black holes resemble in many respects the magnetospheres of rotationally-powered pulsars, with pair formation fronts and outer gaps, where the sign of the induced charge changes. As a result, the black hole will generate bipolar electromagnetic jets each consisting of two counter-aligned current flows (four current flows total), each carrying an electric current of the order $I \approx e B_0 R_G β_0$. The electromagnetic power of the jets is $L \approx (G M)^2 B_0^2 β_0^2/c^3$; for a particular case of merging black holes the resulting Poynting power is $ L \approx {(G M)^3 B_0^2 /(c^5 R)}$, where $R$ is the radius of the orbit. In addition, in limited regions near the horizon the first electromagnetic invariant changes sign, so that the induced electric field becomes larger than the magnetic field, $E>B$. The total energy loss from a system of merging BHs is a sum of two components with similar powers, one due to the rotation of space-time within the orbit, driven by the non-zero angular momentum in the system, and the other due to the linear motion of the BHs through the magnetic field.

preprint2011arXiv

Signatures of large-scale magnetic fields in AGN jets: transverse asymmetries

We investigate the emission properties that a large-scale helical magnetic field imprints on AGN jet synchrotron radiation. A cylindrically symmetric relativistic jet and large-scale helical magnetic field produce significant asymmetrical features in transverse profiles of fractional linear polarization, intensity, Faraday rotation, and spectral index. The asymmetrical features of these transverse profiles correlate with one another in ways specified by the handedness of the helical field, the jet viewing angle (theta_ob), and the bulk Lorentz factor of the flow (Gamma). Thus, these correlations may be used to determine the structure of the magnetic field in the jet. In the case of radio galaxies (theta_ob~1 radian) and a subclass of blazars with particularly small viewing angles (theta_ob << 1/Gamma), we find an edge-brightened intensity profile that is similar to that observed in the radio galaxy M87. We present observations of the AGNs 3C 78 and NRAO 140 that display the type of transverse asymmetries that may be produced by large-scale helical magnetic fields.

preprint2011arXiv

Simple waves in relativistic fluids

We consider the Riemann problem for relativistic flows of polytropic fluids and find relations for the flow characteristics. Evolution of physical quantities take especially simple form for the case of cold magnetized plasmas. We find exact, explicit analytical solutions for one dimensional expansion of magnetized plasma into vacuum, valid for arbitrary magnetization. We also consider expansion into cold unmagnetized external medium both for stationary initial conditions and for initially moving plasma, as well as reflection of rarefaction wave from a wall. We also find self-similar structure of three-dimensional magnetized outflows into vacuum, valid close to the plasma-vacuum interface. The key results of this work, the self-similar solutions, were incorporated post-initial submission into appendices of the published version of Granot et al. (2010).

preprint2011arXiv

Slowly balding black holes

The "no hair" theorem, a key result in General Relativity, states that an isolated black hole is defined by only three parameters: mass, angular momentum, and electric charge; this asymptotic state is reached on a light-crossing time scale. We find that the "no hair" theorem is not formally applicable for black holes formed from collapse of a rotating neutron star. Rotating neutron stars can self-produce particles via vacuum breakdown forming a highly conducting plasma magnetosphere such that magnetic field lines are effectively "frozen-in" the star both before and during collapse. In the limit of no resistivity, this introduces a topological constraint which prohibits the magnetic field from sliding off the newly-formed event horizon. As a result, during collapse of a neutron star into a black hole, the latter conserves the number of magnetic flux tubes $N_B = e Φ_\infty /(πc \hbar)$, where $Φ_\infty \approx 2 π^2 B_{NS} R_{NS}^3 /(P_{\rm NS} c)$ is the initial magnetic flux through the hemispheres of the progenitor and out to infinity. We test this theoretical result via three-dimensional general relativistic plasma simulations of rotating black holes that start with a neutron star dipole magnetic field with no currents initially present outside the event horizon. The black hole's magnetosphere subsequently relaxes to the split monopole magnetic field geometry with self-generated currents outside the event horizon. The dissipation of the resulting equatorial current sheet leads to a slow loss of the anchored flux tubes, a process that balds the black hole on long resistive time scales rather than the short light-crossing time scales expected from the vacuum "no-hair" theorem.

preprint2011arXiv

Spin Tilts in the Double Pulsar Reveal Supernova Spin Angular-Momentum Production

The system PSR J0737-3039 is the only binary pulsar known to consist of two radio pulsars (PSR J0737-3039 A and PSR J0737-3039 B). This unique configuration allows measurements of spin orientation for both pulsars: pulsar A's spin is tilted from the orbital angular momentum by no more than 14 degrees at 95% confidence; pulsar B's by 130 +/- 1 degrees at 99.7% confidence. This spin-spin misalignment requires that the origin of most of B's present-day spin is connected to the supernova that formed pulsar B. Under the simplified assumption of a single, instantaneous kick during the supernova, the spin could be thought of as originating from the off-center nature of the kick, causing pulsar B to tumble to its misaligned state. With this assumption, and using current constraints on the kick magnitude, we find that pulsar B's instantaneous kick must have been displaced from the center of mass of the exploding star by at least 1 km and probably 5--10 km. Regardless of the details of the kick mechanism and the process that produced pulsar B's current spin, the measured spin-spin misalignment in the double pulsar system provides an empirical, direct constraint on the angular momentum production in this supernova. This constraint can be used to guide core-collapse simulations and the quest for understanding the spins and kicks of compact objects.

preprint2011arXiv

The very-high energy emission from pulsars: a case for inverse Compton scattering

The observations of gamma-ray emission from pulsars with the Fermi-LAT detector and the detection of the Crab pulsar with the VERITAS array of Cherenkov telescopes at energies above 100 GeV make it unlikely that curvature radiation is the main source of photons above GeV energies in the Crab and many other pulsars. We outline a model in which the broad UV-X-ray component and the very high energy γ-ray emission of pulsars are explained within the Synchrotron-Self-Compton (SSC) framework. We argue that the bulk of the observed radiation is generated by the secondary plasma, which is produced in cascades in the outer gaps of the magnetosphere. We find that the inverse-Compton (IC) scattering occurs in the Klein-Nishina regime, which favors synchrotron photons in the UV band as target field for the scattering process. The primary beam is accelerated in a modest electric field, with a field strength that is of the order of a few percent of the magnetic field near the light cylinder. Overall, in the Klein-Nishina regime of the IC scattering the particle distribution in the gap does not evolve towards a stationary distribution and thus is intrinsically time-dependent. We point out that in a radiation reaction-limited regime of particle acceleration the gamma-ray luminosity L_γ scales linearly with the pulsar spin-down power \dot{E}, L_γ~ \dot{E}, and not proportional to \sqrt{\dot{E}} as expected from potential-limited acceleration.

preprint2010arXiv

Double explosions and jet formation in GRB-supernova progenitors

Progenitors of long GRBs, and core-collapse supernovae in general, may have two separate mechanisms driving the outflows: quasi-isotropic neutrino-driven supernova explosions followed by a highly collimated relativistic outflow driven by the GRB central engine, a black hole or a magnetar. We consider the dynamics of the second GRB-driven explosion propagating through expanding envelope generated by the passage of the primary supernova shock. Beyond the central core, in the region of steep density gradient created by the SN shock breakout, the accelerating secondary quasi-spherical GRB shock become unstable to corrugation and under certain conditions may form a highly collimated jet, a "chimney", when a flow expands almost exclusively along a nearly cylindrically collimated channel. Thus, weakly non-spherical driving and/or non-spherical initial conditions of the wind cavity may produce highly non-spherical, jetted outflows. For a constant luminosity GRB central engine, this occurs for density gradient in the envelope ρ~ r^{-ω} steeper than ω>4.

preprint2010arXiv

Dynamics of strongly magnetized ejecta in Gamma Ray Bursts

We consider dynamical scales in magnetized GRB outflows, using the solutions to the Riemann problem of expanding arbitrarily magnetized outflows (Lyutikov 2010). For high ejecta magnetization, the behavior of the forward shock closely resembles the so-called thick shell regime of the hydrodynamical expansion. The exception is at small radii, where the motion of the forward shock is determined by the dynamics of subsonic relativistic outflows. The behaviors of the reverse shock is different in fluid and magnetized cases: in the latter case, even for medium magnetization, sigma ~ 1, the reverse shock forms at fairly large distances, and may never form in a wind-type external density profile.

preprint2010arXiv

Electromagnetic extraction of energy from merging black holes

We calculate the evolution of the prompt intrinsic Poynting power generated by merging black holes. Orbiting black holes induce rotation of the space-time. In a presence of magnetic field supported by an accretion disk outside of the orbit, this results in a generation of an electromagnetic outflow via the Blandford-Znajek-type process with luminosity $L_{EM} \sim G^3 M^3 B^2/(c^5 R_{orb}) $ and reaching a fairly low maximum values of $L_{EM} = 10^{37}-10^{39} m_6 $ erg s$^{-1}$ ($m_6$ is the masses of black holes in millions of Solar mass) at the time of the merger. Dissipation of the wind power may produce two types of observed signatures: a highly variable collimated emission coming from the internal dissipation within the jets and a broad-band near-isotropic emission generated at the termination shocks.

preprint2010arXiv

Gamma Ray Bursts: back to the blackboard

Exceptional observational breakthroughs in the field of Gamma Ray Burst research are not paralleled by theoretical advances. In this review, based on the introductory talk given at the "The Shocking Universe" meeting, I argue that any present day model of GRBs, especially of Short type, is grossly incomplete. I will highlight various contradictions with observations that many models face and briefly mention a number of ideas that might or might not work. In particular, I will discuss (i) a possibility that early X-ray afterglows are coming from internal dissipation, and not from the forwards shock; (ii) that prompt radiation is beamed in the outflow frame.

preprint2010arXiv

Probing Shock Properties with Non-thermal X-ray Filaments in Cas A

Thin non-thermal X-ray filaments are often seen in young supernova remnants. We used data from the 1 Ms Chandra observation of Cassiopeia A to study spectral properties of some of the filaments in this remnant. For all the cases that we examined, the X-ray spectrum across the filaments hardens, at about 10% level, going outward, while observed filament widths depend only weakly on the photon energy. Using a model that includes radiative cooling, advection and diffusion of accelerated particles behind the shock, we estimated the magnetic field, turbulence level and shock obliquity.

preprint2010arXiv

Resolving Doppler-factor crisis in AGNs: non-steady magnetized outflows

Magnetically-driven non-stationary acceleration of jets in AGNs results in the leading parts of the flow been accelerated to much higher Lorentz factors than in the case of steady state acceleration with the same parameters. The higher Doppler-boosted parts of the flow may dominate the high energy emission of blazar jets. We suggest that highly variable GeV and TeV emission in blazars is produced by the faster moving leading edges of highly magnetized non-stationary ejection blobs, while the radio data trace the slower-moving bulk flow. Model predictions compare favorably with the latest Fermi gamma-ray and MOJAVE radio VLBI results.

preprint2010arXiv

Structure of magnetic fields in intracluster cavities

Observations of clusters of galaxies show ubiquitous presence of X-ray cavities, presumably blown by the AGN jets. We consider magnetic field structures of these cavities. Stability requires that they contain both toroidal and poloidal magnetic fields, while realistic configurations should have vanishing magnetic field on the boundary. For axisymmetric configurations embedded in unmagnetized plasma, the continuity of poloidal and toroidal magnetic field components on the surface of the bubble then requires solving the elliptical Grad-Shafranov equation with both Dirichlet and Neumann boundary conditions. This leads to a double eigenvalue problem, relating the pressure gradients and the toroidal magnetic field to the radius of the bubble. We have found fully analytical stable solutions. This result is confirmed by numerical simulation. We present synthetic X-ray images and synchrotron emission profiles and evaluate the rotation measure for radiation traversing the bubble.

preprint2009arXiv

High sigma model of pulsar wind nebulae

Pulsars and central engines of long gamma ray burst -- collapsars -- may produce highly magnetized (Poynting flux dominated) outflows expanding in a dense surrounding (interstellar medium or stellar material). For certain injection conditions, the magnetic flux of the wind cannot be accommodated within the cavity. In this case, ideal (non-dissipative) MHD models, similar to the Kennel and Coroniti (1984) model of the Crab nebular, break down (the so-called sigma problem). This is typically taken to imply that the wind should become particle-dominated on scales much smaller than the size of the cavity. The wind is then slowed down by a fluid-type (low magnetization) reverse shock. Recent Fermi results, indicating that synchrotron spectrum of the Crab nebula extends well beyond the upper limit of the most efficient radiation reaction-limited acceleration, contradict the presence of a low sigma reverse shock. We propose an alternative possibility, that the excessive magnetic flux is destroyed in a reconnection-like process in two regions: near the rotational axis and near the equator. We construct an example of such highly magnetized wind having two distinct reconnection regions and suggest that these reconnection cites are observed as tori and jets in pulsar wind nebulae. The model reproduces, qualitatively, the observed morphology of the Crab nebula. In parts of the nebular the dissipation occurs in a relativistically moving wind, alleviating the requirements on the acceleration rate.

preprint2009arXiv

Magnetocentrifugal launching of jets from disks around Kerr black holes

Strong magnetic fields modify particle motion in the curved space-time of spinning black holes and change the stability conditions of circular orbits. We study conditions for magnetocentrifugal jet launching from accretion disks around black holes, whereby large scale black hole lines anchored in the disk may fling tenuous coronal gas outward. For a Schwarzschild black hole, magnetocentrifugal launching requires that the poloidal component of magnetic fields makes an angle less than $60^\circ$ to the outward direction at the disk surface, similar to the Newtonian case. For prograde rotating disks around Kerr black holes, this angle increases and becomes $90^\circ$ for footpoints anchored to the disk near the horizon of a critically spinning $a=M$ black hole. Thus, a disk around a critically spinning black hole may centrifugally launch a jet even along the rotation axis.

preprint2009arXiv

Structure of Cosmic Ray-modified Perpendicular Shocks

Kinetic diffusion of cosmic rays ahead of perpendicular shocks induces charge non-neutrality, which is mostly, yet not completely, screened by the bulk plasma via polarization drift current. Hydrodynamic shear instabilities as well as modified Buneman instability of the polarization current generate the turbulence necessary for a Fermi-type acceleration. Thus, similar to the case of parallel shocks, in perpendicular shocks the diffusing cosmic rays generate unstable plasma currents that in turn excite turbulence. This allows a self-consistent evolution of a shock-cosmic rays system. In the kinetic regime of the modified Buneman instability, electrons may be heated in the cosmic ray precursor.

preprint2009arXiv

Structure of magnetic fields in non-convective stars

(Abridged) We develop a theoretical framework to construct axisymmetric magnetic equilibria in stars, consisting of both poloidal and toroidal magnetic field components. In a stationary axisymmetric configuration, the poloidal current is a function of the poloidal magnetic flux only, and thus should vanish on field lines extending outside of the star. Non-zero poloidal current is limited to a set of toroid-shape flux surfaces fully enclosed inside the star. If we demand that there are no current sheets, then on the separatrix delineating the regions of zero and finite toroidal magnetic field both the poloidal flux function and its derivative should match. Thus, for a given magnetic field in the bulk of the star, the elliptical Grad-Shafranov equation that describes magnetic field structure inside the toroid is an ill-posed problem, with both Dirichlet and Newman boundary conditions and {\it a priori} unknown distribution of toroidal and poloidal electric currents. We discuss a procedure which allows to solve this ill-posed problem by adjusting the unknown current functions. We find a poloidal current-carrying solution that leaves the shape of the flux function and, correspondingly, the toroidal component of the electric current the same as in the case of no poloidal current. The equilibria discussed in this paper may have arbitrary large toroidal magnetic field, and may include a set of stable equilibria. The method developed here can also be applied to magnetic structure of differentially rotating stars, as well as to calculate velocity field in incompressible isolated fluid vortex with a swirl.

preprint2005arXiv

Eclipses and orbital modulations in binary pulsar PSR J0737-3039

Binary radio pulsar system J0737-3039 provides an exceptional opportunity to study innermost structure of pulsar magnetospheres due to very tight orbit, favorable directions of pulsars' rotation and magnetic axes and extremely fortuitous orientation of the orbit. In this system the millisecond pulsar A is eclipsed once per orbit. During eclipse a clear modulation at the 2.77 s period of pulsar B is seen, pointing unambiguously to magnetospheric origin of eclipses. A simple geometric model, based on the idea that the radio pulses are attenuated by synchrotron absorption on the closed magnetic field lines of pulsar B, can successfully reproduces the eclipse light curves down to intricate details. This detailed agreement confirms the dipolar structure of the neutron star's magnetic field. The model gives clear predictions for temporal evolution of eclipse profile due to geodetic precession of pulsar B. In addition, pulsar B shows orbital modulations of intensity, being especially bright at two short orbital phases. We showed that these modulations are due to distortion of pulsar B magnetosphere by pulsar A wind which produces orbital phase-dependent changes of the direction along which radio waves are emitted. Thus, pulsar B is intrinsically bright at all times but its radiation beam misses the Earth at most orbital phases.

preprint2004arXiv

On the nature of eclipses in binary pulsar J0737-3039

We consider magnetohydrodynamical interaction between relativistic pulsar wind and static magnetosphere in binary pulsar system PSR J0737-3039. We construct semi-analytical model describing the form of the interface separating the two pulsars. An assumption of vacuum dipole spin down for Pulsar B leads to eclipse duration ten times longer than observed. We discuss possible Pulsar B torque modification and magnetic field estimates due to the interaction with Pulsar A wind. Unless the orbital inclination is $\leq 86 ^\circ$, the duration of eclipses is typically shorter than the one implied by the size of the eclipsing region. We propose that eclipses occur due to synchrotron absorption by mildly relativistic particles in the shocked Pulsar A wind. The corresponding optical depth may be high enough if Pulsar A wind density is at the upper allowed limits. We derive jump conditions at oblique, relativistic, magnetohydrodynamical shocks and discuss the structure of the shocked Pulsar A wind. Finally, we speculate on a possible mechanism of orbital modulation of Pulsar B radio emission.

preprint2003arXiv

Explosive reconnection in magnetars

X-ray activity of Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters may result from the heating of their magnetic corona by direct currents dissipated by magnetic reconnection. We investigate the possibility that X-ray flares and bursts observed from AXPs and SGRs result from magnetospheric reconnection events initiated by development of tearing mode in magnetically-dominated relativistic plasma. We formulate equations of resistive force-free electrodynamics, discuss its relation to ideal electrodynamics, and give examples of both ideal and resistive equilibria. Resistive force-free current layers are unstable toward the development of small-scale current sheets where resistive effects become important. Thin current sheets are found to be unstable due to the development of resistive force-free tearing mode. The growth rate of tearing mode is intermediate between the short \Alfven time scale $ τ_A$ and a long resistive time scale $τ_R$: $Γ\sim 1/(τ_R τ_A)^{1/2}$, similar to the case of non-relativistic non-force-free plasma. We propose that growth of tearing mode is related to the typical rise time of flares, $\sim 10$ msec. Finally, we discuss how reconnection may explain other magnetar phenomena and ways to test the model.

preprint2002arXiv

Role of Reconnection in AGN Jets

We discuss the possible role of reconnection in electro-magnetically dominated cores of relativistic AGN jets. We suggest that reconnection may proceed in a two-fold fashion: initial explosive collapse on the Alfven time-scale of a current-carrying jet (which is of the order of the light crossing time) and subsequent slow quasi-steady reconnection. Sites of explosive collapse are associated with bright knots, while steady-state reconnection re-energizes particles in the ``bridges'' between the knots. Ohmic dissipation in reconnection layers leads to particle acceleration either by inductive electric fields or by stochastic particle acceleration in the ensuing electro-magnetic turbulence.

preprint1999arXiv

Coherent Mechanisms of Pulsar Radio Emission

Relativistic plasma masers operating on the anomalous cyclotron-Cherenkov resonance and the Cherenkov-drift resonance are capable of explaining the main observational characteristics of pulsar radio emission. Both electromagnetic instabilities are due to the interaction of the fast particles from the primary beam and from the tail of the secondary pairs distribution with the normal modes of a strongly magnetized one-dimensional electron-positron plasma. In a typical pulsar both resonances occur in the outer parts of magnetosphere .

preprint1999arXiv

Excitation of Alfven Waves and Pulsar Radio Emission

We analyze mechanisms of the excitation of Alfvén wave in pulsar magnetospheres as a possible source of pulsar radio emission generation. We find that Cherenkov excitation of obliquely propagating Alfvén waves is inefficient, while excitation at the anomalous cyclotron resonance by the particles from the primary beam and from the tail of the bulk distribution function may have a considerable growth rate. The cyclotron instability on Alfvén waves occurs in the kinetic regime still not very closed to the star: $r \geq 50 R_{NS}$. We also discuss various mechanisms of conversion of Alfvén waves into escaping radiation. Unfortunately, no decisive conclusion about the effectiveness of such conversion can be made.

preprint1999arXiv

Turbulent 4-wave Interaction of Two Type of Waves

We consider turbulent 4-wave interaction of two types of waves: acoustic waves (dispersion $ω= k$) and electromagnetic-type waves (dispersion $Ω^2 = m^2 + p^2$). For large wave vectors ($ k \gg m$), when the dispersion of EM-type waves becomes quasiacoustic, the stationary spectra are obtained following a standard Zakharov approach. In the small wave number region ($ k, p \ll m$) we derive nonlinear differential Kompaneets-type kinetic equation for arbitrary distributions of the interacting fields and find when this equation has Kolmogorov-type power law solutions.

preprint1998arXiv

Cherenkov-drift emission mechanism

Emission of a charged particle propagating in a medium with a curved magnetic field is reconsidered stressing the analogy between this emission mechanism and collective Cherenkov-type plasma emission. It is explained how this mechanism differs from conventional Cherenkov, cyclotron or curvature emission and how it includes, to some extent, the features of each of these mechanisms. Presence of a medium supporting subluminous waves is essential for the possibility of wave amplification by particles streaming along the curved magnetic field with a finite curvature drift. We suggest an analogy between the curvature drift emission and the anomalous cyclotron-Cherenkov emission. Treating the emission in cylindrical coordinates in the plane-wave-like approximation allows one to compute the single particle emissivity and growth rate of the Cherenkov-drift instability. We compare the growth rates calculated using the single particle emissivity and using the dielectric tensor of one dimensional plasma streaming along the curved field. In calculating the single particle emissivity it is essential to know the normal modes of the medium and their polarization which can be found from the dielectric tensor of the medium. This emission mechanism may be important for the problem of pulsar radio emission generation.

Maxim Lyutikov

What is connected

Connect this record

See the researcher in context

Building this map preview

81 published item(s)

Coherent emission in pulsars, magnetars and Fast Radio Bursts: reconnection-driven free electron laser

Faraday conversion in pair-symmetric winds of magnetars and Fast Radio Bursts

Faraday rotation in fast radio bursts

FRB-periodicity: mild pulsars in tight O/B-star binaries

On the nature of Fast Blue Optical Transients

Radius-to-frequency mapping and FRB frequency drifts

Relativistic magnetic explosions

Dynamics and emission of wind-powered afterglows of gamma-ray bursts: flares, plateaus and steep decays

Rotating neutron stars without light cylinders

Coherence constraints on physical parameters at bright radio sources and FRB emission mechanism

Fast moving pulsars as probes of interstellar medium

Fast Radio Bursts from reconnection events in magnetar magnetospheres

Free electron laser in magnetars/Fast Radio Bursts

Magnetospheric interaction in white dwarf binaries AR Sco and AE Aqr

Nonlinear optics in strongly magnetized pair plasma, with applications to FRBs

Peeking Between the Pulses: The Far-UV Spectrum of the Previously Unseen White Dwarf in AR Scorpii

Tilting Instability of Magnetically Confined Spheromaks

Turbulent model of Crab nebula radiation

Wind-powered afterglows of gamma-ray bursts: flares, plateaus and steep decays

Conditions for jet break-out in NS mergers

Nonlinear force-free configurations in cylindrical geometry

Emission knots and polarization swings of swinging jets

Fast radio bursts as giant pulses from young rapidly rotating pulsars

Fermi GBM signal contemporaneous with GW150914 - an unlikely association

How else can we detect Fast Radio Bursts?

On the Linear Stability of Sheared and Magnetized Jets Without Current Sheets - Non-Relativistic Case

Relativistic Rotating Vector Model

Magnetar Activity via the Density-Shear Instability in Hall-MHD

On the Linear Stability of Magnetized Jets Without Current Sheets: Non-Relativistic Case

The inner knot of the Crab nebula

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

Density-shear instability in electron MHD

Magnetar activity mediated by plastic deformations of neutron star crust

Super-Eddington accretion in ultra-luminous neutron star binary

A Giant Sample of Giant Pulses from the Crab Pulsar

Afterglow jet breaks, GeV photons and the electromagnetic model of short GRBs

Electron MHD: dynamics and turbulence

I. Inverse Compton origin of pulsar γ-ray emission. II. Reconnection model of Crab flares

Magnetospheric "anti-glitches" in magnetars

Radio Emission Region Exposed: Courtesy of the Double Pulsar

Topics in microphysics of relativistic plasmas

Crab Nebula gamma-ray flares as relativistic reconnection minijets

Hair of astrophysical black holes

IC model of pulsar high energy emission

Inverse Compton model of pulsar high energy emission

Magnetic structure of Coronal Mass Ejections

On the dynamics of mechanical failures in magnetized neutron-star crusts

The electromagnetic model of short GRBs, the nature of prompt tails, supernova-less long GRBs and highly efficient episodic accretion

The gamma-ray spectrum of Geminga and the inverse Compton model of pulsar high energy emission

Coronal mass ejections as expanding force-free structures

Crab GeV flares from corrugated termination shock

Dynamics of rising magnetized cavities and UHECR acceleration in clusters of galaxies

Electromagnetic power of merging and collapsing compact objects

In Situ Heating of the 2007 May 19 CME Ejecta Detected by STEREO/PLASTIC and ACE

On the sideways expansion of relativistic non-spherical shocks and GRB afterglows

Relativistic magnetohydrodynamics in one dimension

Schwarzschild black holes as unipolar inductors: expected electromagnetic power of a merger

Signatures of large-scale magnetic fields in AGN jets: transverse asymmetries

Simple waves in relativistic fluids

Slowly balding black holes

Spin Tilts in the Double Pulsar Reveal Supernova Spin Angular-Momentum Production

The very-high energy emission from pulsars: a case for inverse Compton scattering

Double explosions and jet formation in GRB-supernova progenitors

Dynamics of strongly magnetized ejecta in Gamma Ray Bursts

Electromagnetic extraction of energy from merging black holes

Gamma Ray Bursts: back to the blackboard

Probing Shock Properties with Non-thermal X-ray Filaments in Cas A

Resolving Doppler-factor crisis in AGNs: non-steady magnetized outflows

Structure of magnetic fields in intracluster cavities

High sigma model of pulsar wind nebulae

Magnetocentrifugal launching of jets from disks around Kerr black holes

Structure of Cosmic Ray-modified Perpendicular Shocks

Structure of magnetic fields in non-convective stars

Eclipses and orbital modulations in binary pulsar PSR J0737-3039