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\'en 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\'en 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.

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.

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 $\gamma_b \sim$ few hundred.

Pulsars moving through ISM produce bow shocks detected in hydrogen H$\alpha$ 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$\alpha$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $\delta\rho/\rho \sim 1$, leads to significant variations of the shape of the shock seen in H$\alpha$ 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.

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 ($\sigma$ $\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.

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.

The Crab pulsar has striking radio emission properties, with the two dominant pulse components -- the main pulse and the interpulse -- consisting entirely of giant pulses. The emission is scattered in both the Crab nebula and the interstellar medium, causing multi-path propagation and thus scintillation. We study the scintillation of the Crab's giant pulses using phased Westerbork data at 1668\,MHz. We find that giant pulse spectra correlate at only $\sim 2 \%$, much lower than the $1/3$ correlation expected from a randomized signal imparted with the same impulse response function. In addition, we find that the main pulse and the interpulse appear to scintillate differently, with the scintillation of the interpulse correlating to the main pulse at a lower amplitude over a longer time and frequency range. These lines of evidence suggest that the giant pulse emission regions are extended, and that the main pulse and interpulse arise in physically distinct regions which are resolved by the scattering screen. Assuming the scattering takes place in the nebular filaments, the emission regions are of order a light cylinder radius, as projected on the sky. With further VLBI and multi-frequency data, it may be possible to measure the distance to the scattering screens, the size of giant pulse emission regions, and the physical separation between the pulse components.

We perform numerical simulations of the dynamics of magnetized and highly relativistic pulsar-like winds produced by long-lasting central engines in Gamma Ray Bursters (GRB). The very fast wind interacts with the initial relativistically expanding GRB outflow, creating a multiple-shock structure. Depending on the parameters of the model (the energy of the initial explosion, the outside density, the wind power, the delay time for the switch-on of the wind, and the magnetization of the wind) a variety of temporal behaviors can be produced. In certain regimes the dynamics of such secondary engine-driven shocks is self-similar.

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.

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, $\eta = 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 $\eta \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.

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.

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 $\omega_B$ dominates over plasma frequency $\omega_p$ and the frequency of the electromagnetic wave, $\omega_B \gg \omega_p,\, \omega$, 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 $(\omega/\omega_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 $\Delta B/B_0 = 2\pi 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.

We consider electromagnetic interaction and precursor emission of merging neutron stars. Orbital motion of the magnetized neutron stars may revive pair production within the common magnetosphere years before the merger, igniting pulsar-like magnetospheric dynamics. We identify two basic scenarios: (i) only one star is magnetized (1M-DNS scenario) and (ii) both stars are magnetized (2M-DNS scenario). Inductively created electric fields can have component along the total magnetic field (gaps) and/or the electric field may exceed the value of the local magnetic field. The key to the detection is orbital modulation of the emission. If only one star is magnetized (1M-DNS scenario) the emission is likely to be produced along the direction of the magnetic field at the location of the secondary; then, if the magnetic axis is misaligned with the orbital spin, this direction is modulated on the orbital period. For the 2M-DNS scenario, the structure of the common magnetosphere of the non-rotating neutron stars is complicated, with gaps, but no $E>B$ regions; there is strong orbital variations for the case of misaligned magnetic moments. For the same parameters of neutron stars the 2M-DNS scenario has intrinsically higher potential than the 1M-DNS one. The overall powers are not very high, $\leq 10^{45} $ erg s$^{-1}$; the best chance to detect electromagnetic precursors to the merging neutron stars is if the interaction of their magnetospheres leads to the production of pulsar-like coherent radio emission modulated at the orbital period, with luminosity of up to $\sim 1$ Jankys at the time the merger.

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{\omega} \propto - \omega$, matching both numerically and parametrically the rates observed in FBRs; a more complicated behavior of $\dot{\omega} $ 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.

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}/ \beta_0$, where $ \beta_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 ( {\beta_{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.

Mergers of white dwarfs (WDs) may lead to a variety of transient astrophysical events, SNIa being one possible outcome. Lyutikov & Toonen (2017, 2019) argued that mergers of WDs result, under various parameter regimes, in unusual central engine-powered supernova and a type of short Gamma Ray Bursts that show extended emission tails. Observations by Gvaramadze et al. (2019) of the central star and the nebula J005311 match to the details the model of Lyutikov & Toonen (2017, 2019) for the immediate product of a merger of a heavy ONeMg WD with CO WD (age, luminosity, stellar size, hydrogen deficiency and chemical composition).

We interpret recent observations of high-to-low frequency drifting features in the spectra of the repeating FRBs as evidence of sharply changing plasma properties in the emission region, presumably the neutron stars magnetospheres. The drifts are then FRBs' analogues of radius-to-frequency mapping in pulsars and Solar type-III radio burst (but not in a sense of a particular emission mechanism). The drifts rates of $\sim 100$ MHz ms$^{-1}$ at frequencies $\sim$ GHz translate to physical size of $\sim {c \omega}/{\dot{\omega}} \sim$ few $\times 10^8$ cm, matching the hypothesis of the FRB origin in the magnetospheres of neutron stars. We suggest that reconnection events result in generation of upward propagating plasma beams that produce radio emission with frequency related to the decreasing local magnetic field and plasma density.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

We outline a model of the Crab Pulsar Wind Nebula with two different populations of synchrotron emitting particles, arising from two different acceleration mechanisms: (i) Component-I due to Fermi-I acceleration at the equatorial portion of the termination shock, with particle spectral index $p_I \approx 2.2$ above the injection break corresponding to $\gamma_{wind} \sigma_{wind} \sim 10^5$, peaking in the UV ($\gamma_{wind} \sim 10^2$ is the bulk Lorentz factor of the wind, $\sigma _{wind} \sim 10^3$ is wind magnetization); (ii) Component-II due to acceleration at reconnection layers in the bulk of the turbulent Nebula, with particle index $p_{II} \approx 1.6$. The model requires relatively slow but highly magnetized wind. For both components the overall cooling break is in the infra-red at $\sim 0.01$ eV, so that the Component-I is in the fast cooling regime (cooling frequency below the peak frequency). In the optical band Component-I produces emission with the cooling spectral index of $\alpha_o \approx 0.5$, softening towards the edges due to radiative losses. Above the cooling break, in the optical, UV and X-rays, Component-I mostly overwhelms Component-II. We hypothesize that acceleration at large-scale current sheets in the turbulent nebula (Component-II) extends to the synchrotron burn-off limit of $\epsilon_s \approx 100$ MeV. Thus in our model acceleration in turbulent reconnection (Component-II) can produce both hard radio spectra and occasional gamma-ray flares. This model may be applicable to a broader class of high energy astrophysical objects, like AGNe and GRB jets, where often radio electrons form a different population from the high energy electrons.

We suggest that fast-rising blue optical transients (FBOTs) and the brightest event of the class AT2018cow result from an electron-capture collapse to a \NS\ following a merger of a massive ONeMg white dwarf (WD) with another WD. Two distinct evolutionary channels lead to the disruption of the less massive WD during the merger and the formation of a shell burning non-degenerate star incorporating the ONeMg core. During the shell burning stage a large fraction of the envelope is lost to the wind, while mass and angular momentum are added to the core. As a result, the electron-capture collapse occurs with a small envelope mass, after $\sim 10^2-10^4$ years. During the formation of a neutron star as little as $\sim 10^{-2} M_\odot $ of the material is ejected at the bounce-off with mildly relativistic velocities and total energy $\sim$ few $ 10^{50}$ ergs. This ejecta becomes optically thin on a time scale of days - this is the FBOT. During the collapse, the neutron star is spun up and magnetic field is amplified. The ensuing fast magnetically-dominated relativistic wind from the newly formed neutron star shocks against the ejecta, and later against the wind. The radiation-dominated forward shock produces the long-lasting optical afterglow, while the termination shock of the relativistic wind produces the high energy emission in a manner similar to Pulsar Wind Nebulae. If the secondary WD was of the DA type, the wind will likely have $\sim 10^{-4} M_\odot$ of hydrogen; this explains the appearance of hydrogen late in the afterglow spectrum. The model explains many of the puzzling properties of FBOTs/AT2018cow: host galaxies, a fast and light anisotropic ejecta producing a bright optical peak, afterglow high energy emission of similar luminosity to the optical, and late infra-red features.

We suggest that narrow, long radio filaments near the Galactic Center arise as kinetic jets - streams of high energy particles escaping from ram-pressure confined pulsar wind nebulae (PWNe). The reconnection between the PWN and interstellar magnetic field allows pulsar wind particles to escape, creating long narrow features. They are the low frequency analogs of kinetic jets seen around some fast-moving pulsars, such as The Guitar and The Lighthouse PWNe. The radio filaments trace a population of pulsars also responsible for the Fermi GeV excess produced by the Inverse Compton scattering by the pulsar wind particles. The magnetic flux tubes are stretched radially by the large scale Galactic winds. In addition to PWNe accelerated particles can be injected at supernovae remnants. The model predicts variations of the structure of the largest filaments on scales of $\sim$ dozens of years - smaller variations can occur on shorter time scales. We also encourage targeted observations of the brightest sections of the filaments and of the related unresolved point sources in search of the powering PWNe and pulsars.

We reconsider the structure of thermally driven rotating Parker wind. Rotation, without \Bf, changes qualitatively the structure of the subsonic region: solutions become non-monotonic and do not extend to the origin. For small angular velocities solutions have two critical points - X-point and O-points, which merge at the critical angular velocity of the central star $\Omega_{crit} = G M_\ast/(2 \sqrt{2} c_s R_{b}^2)$ (where $M_\ast$ and $R_{b}$ are mass and radius of the central star, $c_s$ is the sound speed in the wind). For larger spins there is no critical points in the solution. For disk winds (when the base of the wind rotates with Keplerian velocity) launched equatorially the coronal sound speed should be smaller than $\approx 0.22 v_K$ in order to connect to the critical curve ($v_K$ is the Keplerian velocity at a given location on the disk).

Bow-shock pulsar wind nebulae (PWNe) show a variety of morphological shapes. We attribute this diversity to the geometrical factors: relative orientations of the pulsar rotation axis, proper velocity, and the line of sight (magnetic inclination angle may also have a certain influence on the morphology). We identify three basic types of bow-shock nebulae: (i) a "Rifle Bullet" (pulsar spin axis and proper velocity are aligned); (ii) a "Frisbee" (pulsar spin axis and proper velocity are orthogonal with the spin axis lying in the plane of the sky), and (iii) a Cart Wheel" (like frisbee but the spin axis is perpendicular to the plane of the sky). Using 3D RMHD simulations, as well as analytical calculations, we reproduce the key morphological features of the bow-shock PWNe, as well as variations, are seen across different systems. magnetic stresses within the shocked pulsar wind affect the overall structure strongly, producing "whiskers", "tails", "filled-in" and "mushroom" shapes, as well as non-symmetric morphologies. On the other hand, the interstellar medium inhomogeneities and the anisotropy of the energy flux in the pulsar wind have only a mild impact of the PWN morphology. In a few cases, when we clearly identify specific morphological structures, our results do not favor alignment of the pulsar spin axis and proper velocity. Our calculations of the underlying emission processes explain the low synchrotron X-ray efficiency (in terms of the spin-down luminosity) and imply an energetical subdominant contribution of the inverse Compton process.

We obtained high temporal resolution spectroscopy of the unusual binary system AR Sco covering nearly an orbit. The H$\alpha$ emission shows a complex line structure similar to that seen in some polars during quiescence. Such emission is thought to be due to long-lived prominences originating on the red dwarf. A difference between AR Sco and these other systems is that the white dwarf in AR Sco is rapidly spinning relative to the orbital period. "Slingshot" prominences stable at 3 to 5 stellar radii require surface magnetic fields between 100 and 500 G. This is comparable to the estimated WD magnetic field strength near the surface of the secondary. Our time-resolved spectra also show emission fluxes, line equivalent widths, and continuum color varying over the orbit and the beat/spin periods of the system. During much of the orbit, the optical spectral variations are consistent with synchrotron emission with the highest energy electrons cooling between pulses. On the time-scale of the beat/spin period we detect red and blue-shifted H$\alpha$ emission flashes that reach velocities of 700 km/s. Red-shifted Balmer emission flashes are correlated with the bright phases of the continuum beat pulses while blue-shifted flashes appear to prefer the time of minimum in the beat light curve. We propose that much of the energy generated in AR Sco comes from fast magnetic reconnection events occurring near the inward face of the secondary and we show that the energy generated by magnetic reconnection can account for the observed excess luminosity from the system.

We performed calculations of the late radio and X-ray afterglow of GRB/GW170817 in the cocoon-jet paradigm, predicting the appearance of a second peak in the afterglow light curve ~ one-three years after the explosion. The model assumes that the prompt emission and early afterglows originate from a cocoon generated during break-out of the delayed magnetically powered jet. As the jet breaks out from the torus-generated wind, a nearly isotropic mildly relativistic outflow is generated; at the same time the primary jet accelerates to high Lorentz factors and avoids detection. As the fast jet slows down, it should become visible to the off-axis observer. Thus, the model has a clear prediction: the X-ray and radio afterglows should first experience a decay, as the cocoon slows down, followed by a rebrightening when the primary jet starts emitting toward an observer.

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically-dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two dimensional `ABC' structures) and zero total current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For plasma magnetization $\sigma \leq 10^2$ the spectrum power law index is $p> 2$; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\sigma \geq 10^2$, the spectra are hard, $p< 2$, yet the maximal energy $\gamma_{max}$ can still exceed the average magnetic energy per particle, $ \sim \sigma$, by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

We develop a model of gamma-ray flares of the Crab Nebula resulting from the magnetic reconnection events in highly-magnetized relativistic plasma. We first discuss physical parameters of the Crab nebula and review the theory of pulsar winds and termination shocks. We also review the principle points of particle acceleration in explosive reconnection events (Lyutikov et al. 2017a,b). It is required that particles producing flares are accelerated in highly magnetized regions of the nebula. Flares originate from the poleward regions at the base of Crab's polar outflow, where both the magnetization and the magnetic field strength are sufficiently high. The post-termination shock flow develops macroscopic (not related to the plasma properties on the skin-depth scale) kink-type instabilities. The resulting large-scales magnetic stresses drive explosive reconnection events on the light-crossing time of the reconnection region. Flares are produced at the initial stage of the current sheet development, during the X-point collapse. The model has all the ingredients needed for Crab flares: natural formation of highly magnetized regions, explosive dynamics on light travel time, development of high electric fields on macroscopic scales and acceleration of particles to energies well exceeding the average magnetic energy per particle.

The extreme properties of the gamma ray flares in the Crab Nebula present a clear challenge to our ideas on the nature of particle acceleration in relativistic astrophysical plasma. It seems highly unlikely that standard mechanisms of stochastic type are at work here and hence the attention of theorists has switched to linear acceleration in magnetic reconnection events. In this series of papers, we attempt to develop a theory of explosive magnetic reconnection in highly-magnetized relativistic plasma which can explain the extreme parameters of the Crab flares. In the first paper, we focus on the properties of the X-point collapse. Using analytical and numerical methods (fluid and particle-in-cell simulations) we extend Syrovatsky's classical model of such collapse to the relativistic regime. We find that the collapse can lead to the reconnection rate approaching the speed of light on macroscopic scales. During the collapse, the plasma particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\sigma$. For sufficiently high magnetizations and vanishing guide field, the non-thermal particle spectrum consists of two components: a low-energy population with soft spectrum, that dominates the number census; and a high-energy population with hard spectrum, that possesses all the properties needed to explain the Crab flares.

Some fast-moving pulsars, like the Guitar and the Lighthouse, exhibit asymmetric non-thermal emission features that extend well beyond their ram pressure confined pulsar wind nebulae (PWNe). Using 3D relativistic simulations we explain these features as kinetically streaming pulsar wind particles that escaped into the interstellar medium (ISM) due to reconnection between the PWN and ISM magnetic fields. The structure of the reconnecting magnetic fields at the incoming and outgoing regions produce highly asymmetric magnetic bottles, and result in asymmetric extended features. For the features to become visible, the ISM magnetic field should be sufficiently high, $B_{\rm ISM}>10$~$\mu$G. We also discuss archival observations of PWNe displaying evidence of kinetic jets: the Dragonfly PWN (PSR J2021+3651), G327.1--1.1, and MSH 11--62, the latter two of which exhibit "snail eyes" morphologies. We suggest that in those cases the pulsar is moving along the ambient magnetic field in a frisbee-type configuration.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

We consider the structure of mildly relativistic shocks in dense media, taking into account the radiation and pair loading, and diffusive radiation energy transfer within the flow. For increasing shock velocity (increasing post-shock temperature), the first important effect is the efficient energy redistribution by radiation within the shock that leads to the appearance of an isothermal jump, whereby the flow reaches the finial state through a discontinuous isothermal transition. The isothermal jump, on scales much smaller that the photon diffusion length, consists of a weak shock and a quick relaxation to the isothermal conditions. Highly radiation-dominated shocks do not form isothermal jump. Pair production can mildly increase the overall shock compression ratio to $\approx 10$ (from $4$ for matter-dominated shocks and $7$ of the radiation-dominated shocks).

We present our observations of electromagnetic transients associated with GW170817/GRB 170817A using optical telescopes of Chilescope observatory and Big Scanning Antenna (BSA) of Pushchino Radio Astronomy Observatory at 110MHz. The Chilescope observatory detected an optical transient of $\sim19^m$ on the third day in the outskirts of the galaxy NGC 4993; we continued observations following its rapid decrease. We put an upper limit of $1.5\times10^{4}$ Jy on any radio source with a duration of 10-60 s which may be associated with GW170817/GRB 170817A. The prompt gamma-ray emission consists of two distinctive components - a hard short pulse delayed by $\sim2$ seconds with respect to the LIGO signal and softer thermal pulse with $T\sim 10 $ keV lasting for another $\sim2$ seconds. The appearance of a thermal component at the end of the burst is unusual for short GRBs. Both the hard and the soft components do not satisfy the Amati relation, making GRB 170817A distinctively different from other short GRBs. Based on gamma-ray and optical observations, we develop a model for the prompt high-energy emission associated with GRB 170817A. The merger of two neutron stars creates an accretion torus of $\sim10^{-2} M_\odot$, which supplies the black hole with magnetic flux and confines the Blandford-Znajek-powered jet. We associate the hard prompt spike with the quasispherical breakout of the jet from the disk wind. As the jet plows through the wind with subrelativistic velocity, it creates a radiation-dominated shock that heats the wind material to tens of kiloelectron volts, producing the soft thermal component.

We perform spectral analyses of four bright $\gamma$-ray pulsars: PSR J0007+7303, Vela, Geminga and PSR J2021+3651 concentrating on the high-energy tails, defined as emission above 10\,GeV. The two competing models of pulsar $\gamma$-ray emission predict qualitatively different spectra well above the break energy: curvature emission predicts an exponential cut-off in the spectra, while the inverse-Compton scattering mechanism favors a power-law. We perform fits to the phase-averaged spectral energy distributions for each of the four pulsars. We find that in all cases both models fit the data equally well---the present data set does not allow any firm claim to be made about the shape of the spectra above 10\,GeV. In no case is the power-law fit or exponential cut-off fit significantly preferred over the other. The Crab pulsar remains the only known pulsar in which the power-law fit is clearly preferred over the exponential cut-off.

We point out that for sufficiently strong shocks, with Mach number $ M_1 > \sqrt{\frac{3 \gamma -1}{(3-\gamma ) \gamma }}= 1.34$ ($\gamma=5/3$), the solutions for cosmic ray-modified shocks experiences a bifurcation. As a result, for super-critical flows an isothermal jump forms (which is not a shock). The isothermal jump forms due to the energy diffusion of fast, but energetically subdominant cosmic rays. For super-critical flows the isothermal jump appears regardless of a particular feed-back mechanism from the CRs. The compression ratio at the isothermal jump is $2/ (\gamma-1)=3$, so that in the test particle regime the expected spectrum of low energy CRs experiencing first-order Fermi process is $p = 2\gamma/(3- \gamma) = 5/2$, steeper than conventional $p=2$.

Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio pulsar PSR J1119-6127 at X-ray, with XMM-Newton & NuSTAR, and at radio energies with Parkes radio telescope, during a period of magnetar-like bursts. The rotationally powered radio emission shuts off coincident with the occurrence of multiple X-ray bursts, and recovers on a time scale of ~70 seconds. These observations of related radio and X-ray phenomena further solidify the connection between radio pulsars and magnetars, and suggest that the pair plasma produced in bursts can disrupt the acceleration mechanism of radio emitting particles.

We consider structure and emission properties of a pair plasma fireball that cools due to radiation. At temperatures $T \geq 0.5 m_e c^2$ the cooling takes a form of clearly defined cooling wave, whereby the temperature and pair density experience a sharp drop within a narrow region. The surface temperature, corresponding to the location where the optical depth to infinity reaches unity, never falls much below $0.1 m_e c^2 \approx 50$ keV. The propagation velocity of the cooling wave is much smaller than the speed of light and decreases with increasing bulk temperature.

We discuss a scenario of short Gamma-Ray Bursts (GRBs) following a merger of a massive ONeMg white dwarf (WD) with a CO WD, and an ensuing accretion induced collapse (AIC). An initial system with the primary mass $M_1 \sim 6-10 \, M_\odot$ and the secondary mass $M_2 \sim 3-6\, M_\odot$ forms, via two distinct evolutionary channels, a double degenerate CO-ONeMg WD system. For sufficiently large mass ratio $q \equiv M_2/M_1 > q_{crit} \sim 0.25$ the ensuing gravitational wave-driven mass transfer is unstable, whereby the less massive CO WD is disrupted and transfers its mass to the primary ONeMg WD on a few orbital time scales. The merger product ignites shell CO burning, adding mass to the degenerate core; at the same time mass and angular momentum is lost due to powerful winds. For an ONeMg WD sufficiently close to the Chandrasekhar mass an electron-capture accretion induced collapse (AIC) follows $\sim 10^4$ years afterwards. We associate the prompt short GRB emission with a direct collapse of an ONeMg WD to a neutron star, without formation of an accretion disk. After the collapse the accretion of the unburnt part of the shell onto the newly formed NS powers the extended emission (EE). During the collapse the neutron star is spun to millisecond periods and produce long lasting relativistic winds that shock against the material lost during the shell-burning stage, and produce afterglow emission from the wind termination shock.

We consider the double-shock structure of ultra-relativistic flows produced by the interaction of magnetized wind with magnetized external medium. The contact discontinuity (CD) is a special point in the flow - density, kinetic pressure and magnetic field experience a jump or are non-analytic on the CD. To connect dynamically the outside region (the forward shock flow) with the inside region (the reverse shock flow) requires resolving flow singularities at the contact discontinuity. On the CD the pressure is communicated exclusively by the magnetic field on both sides - the CD becomes an Alfven tangential discontinuity. Thus, the dynamic amplification of the magnetic field leads to a formation of a narrow magnetosheath. We discussed a possibility that particles emitting early $X$-ray afterglows, as well as Fermi GeV photons, are accelerated via magnetic reconnection processes in the post-reverse shock region, and in the magnetosheath in particular.

We present a simple prescription for the rotation of polarization produced by a relativistically moving gravitational lens, applicable to arbitrary deflection angle and arbitrary velocity of the lens. When geometric optics is applicable, two independent components contribute to the total rotation of polarization: (i) in the frame of the lens the polarization vector experiences minimal rotation defined by the deflection angle (as measured by a set of remote observers, or no rotation if defined in terms of parallel-propagated tetrad); (ii) the effect of the motion of the lens on the polarization can be taken into account exactly using special relativistic Lorentz transformation of polarization. The effects of the gravitational lensing are thus parametrized by the deflection angle of the null geodesics (not necessarily small) and the motion of the lens (not necessarily with velocities much smaller than that of light).

In our prior papers, we considered the non-relativistic linear stability analysis of magnetized jets that do not have current sheet at the boundary. In this paper, we extend our analysis to relativistic jets. In order to find the unstable modes of current sheet-free, magnetized relativistic jets, we linearize full relativistic magnetohydrodynamics equations and solve them numerically. We find the dispersion relation of the pinch and kink mode instabilities. By comparing the dispersion relations of mildly relativistic jet (Lorentz factor 2) with moderately relativistic jet (Lorentz factor 10), we find that the jet with higher Lorentz factor is significantly more stable in both pinch and kink modes. We show that inclusion of the current sheet-free magnetic field in the jet further enhances the stability. Both pinch and kink mode instabilities become progressively more stable with increasing magnetization. We also show a scaling relation between the maximum temporal growth rate of the unstable mode and the Lorentz factor of the jet. The maximum temporal growth rates of the unstable modes are inversely proportion to the Lorentz factors for most of the modes that we study. However, for the fundamental pinch mode it is inversely proportional to the square of the Lorentz factor. This very beneficial scaling relation holds regardless of the presence of a magnetic field.

We present a model of blazar variability that can both reproduce smooth large polarization angle swings, and at the same time allow for the seemingly random behaviour of synchrotron fluxes, polarization fraction and, occasionally, $\pi$/2 polarization jumps. We associate blazar flaring activity with a jet carrying helical magnetic fields and propagating along a variable direction (and possibly with a changing bulk Lorentz factor). The model predicts that for various jet trajectories (i) EVPA can experience large smooth temporal variations while at the same time polarization fraction ($\Pi$) can be highly variable; (ii) $\Pi \sim 0$ near sudden EVPA jumps of 90$^{\circ}$, but can also remain constant for large, smoother EVPA swings; (iii) the total angle of EVPA rotation can be arbitrary large; (iv) intensity I is usually maximal at points of fastest EVPA changes, but can have a minimum. Thus, even for a regular, deterministic motion of a steadily emitting jet the observed properties can vary in a non-monotonic and/or seemingly stochastic way. Intrinsic fluctuations of the emissivity will further complicated the intensity profiles, but are expected to preserve the polarization structure.

We can probe observationally and reproduce theoretically intricate properties of the Crab Nebula nearest to the pulsar - The Inner Knot. The tiny knot is indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock that accelerates particles. It is required that the part of the wind that produces the Inner Knot has low magnetization; thus, it is not a site of gamma-ray flares. We develop a model of particle acceleration during explosive reconnection events in relativistic highly magnetized plasma and apply the model to explain the Crab gamma-ray flares. Particles are efficiently accelerated by charge-starved DC-type electric fields during initial stages of magnetic flux merges. By implication, the magnetic reconnection is an important, and possibly dominant process of particle acceleration in high energy astrophysical sources.

We develop a model of early GRB afterglows with the dominant $X$-ray contribution from the reverse shock (RS) propagating in highly relativistic (Lorentz factor $\gamma_w \sim 10^6$) magnetized wind of a long-lasting central engine. The model reproduces, in a fairly natural way, the overall trends and yet allows for variations in the temporal and spectral evolution of early optical and $X$-ray afterglows. The high energy and the optical synchrotron emission from the RS particles occurs in the fast cooling regime; the resulting synchrotron power $L_s$ is a large fraction of the wind luminosity, $L_s \approx L_w/\sqrt{1+\sigma_w}$ ($L_w$ and $\sigma_w$ are wind power and magnetization). Thus, plateaus - parts of afterglow light curves that show slowly decreasing spectral power - are a natural consequence of the RS emission. Contribution from the forward shock (FS) is negligible in the $X$-rays, but in the optical both FS and RS contribute similarly: the FS optical emission is in the slow cooling regime, producing smooth components, while the RS optical emission is in the fast cooling regime, and thus can both produce optical plateaus and account for fast optical variability correlated with the $X$-rays, e.g., due to changes in the wind properties. We discuss how the RS emission in the $X$-rays and combined FS and RS emission in the optical can explain many of puzzling properties of early GRB afterglows.

We consider fluid dynamics of relativistic double explosion - when a point explosion with energy $E_1$ is followed by a second explosion with energy $E_2$ after time $t_d$ (the second explosion could be in a form of a long lasting wind). The primary explosion creates a self-similar relativistic blast wave propagating with \Lf\ $\Gamma_1(t)$. A sufficiently strong second explosion, with total energy $E_2 \geq 10^{-2} E_1$, creates a fast second shock in the external fluid previously shocked by the primary shock. At times longer than the interval between the explosions $t_d$, yet short compared with the time when the second shock catches up the primary shock at $\sim t_d \Gamma_1^2$, the structure of the second shock is approximately self-similar. Self-similar structure of the second shock exist for the case of constant external density (in this case $\Gamma_2 \propto t^{-7/3}$), but not for the wind environment. At early times the \Lf\ of the second shock may exceed that of the primary shock and may boost the synchrotron emission of locally accelerated electrons into the Fermi LAT range.

We argue that the localization of the Repeating FRB at $\sim 1$ Gpc excludes rotationally-powered type of radio emission (e.g., analogues of Crab's giant pulses coming from very young energetic pulsars) as the origin of FRBs.

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 $\pi/2$ polarization jumps.

In this paper we describe a simple numerical approach which allows to study the structure of steady-state axisymmetric relativistic jets using one-dimensional time-dependent simulations. It is based on the fact that for narrow jets with v~c the steady-state equations of relativistic magnetohydrodynamics can be accurately approximated by the one-dimensional time-dependent equations after the substitution z=ct. Since only the time-dependent codes are now publicly available this is a valuable and efficient alternative to the development of a high-specialized code for the time-independent equations. The approach is also much cheaper and more robust compared to the relaxation method. We tested this technique against numerical and analytical solutions found in literature as well as solutions we obtained using the relaxation method and found it sufficiently accurate. In the process, we discovered the reason for the failure of the self-similar analytical model of the jet reconfinement in relatively flat atmospheres and elucidated the nature of radial oscillations of steady-state jets.

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.

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.

We investigate the dynamics of bow shock nebulae created by pulsars moving supersonically through a partially ionized interstellar medium. A fraction of interstellar neutral hydrogen atoms penetrating into the tail region of a pulsar wind will undergo photo-ionization due to the UV light emitted by the nebula, with the resulting mass loading dramatically changing the flow dynamics of the light leptonic pulsar wind. Using a quasi 1-D hydrodynamic model of relativistic flow we find that if a relatively small density of neutral hydrogen, as low as $10^{-4}$ cm$^{-3}$, penetrate inside the pulsar wind, this is sufficient to strongly affect the tail flow. Mass loading leads to the fast expansion of the pulsar wind tail, making the tail flow intrinsically non-stationary. The shapes predicted for the bow shock nebulae compare well with observations, both in H$\alpha$ and X-rays.

We study the transport of high-energy particles in pulsar wind nebulae (PWN) using three-dimensional MHD (see Porth et al. (2014) for details) and test-particle simulations, as well as a Fokker-Planck particle transport model. The latter includes radiative and adiabatic losses, diffusion, and advection on the background flow of the simulated MHD nebula. By combining the models, the spatial evolution of flux and photon index of the X-ray synchrotron emission is modelled for the three nebulae G21.5-0.9, the inner regions of Vela, and 3C 58, thereby allowing us to derive governing parameters: the magnetic field strength, average flow velocity and spatial diffusion coefficient. For comparison, the nebulae are also modelled with the semi-analytic Kennel & Coroniti (1984) model but the Porth et al. (2014) model generally yields better fits to the observational data. We find that high velocity fluctuations in the turbulent nebula (downstream of the termination shock) give rise to efficient diffusive transport of particles, with average P\'eclet number close to unity, indicating that both advection and diffusion play an important role in particle transport. We find that the diffusive transport coefficient of the order of $\sim2\times 10^{27} (L_{\rm s}/0.42\rm Ly) cm^{2}s^{-1}$ ($L_{\rm s}$ is the size of the termination shock) is independent of energy up to extreme particle Lorentz factors of $\gamma_{p}\sim10^{10}$.

We develop a model of particle acceleration in explosive reconnection events in relativistic magnetically-dominated plasmas and apply it to explain gamma-ray flares from the Crab Nebula. The model relies on development of current-driven instabilities on macroscopic scales (not related to plasma skin depths). Using analytical and numerical methods (fluid and particle-in-cell simulations), we study a number of model problems in relativistic magnetically-dominated plasma: (i) we extend Syrovatsky's classical model of explosive X-point collapse to magnetically-dominated plasmas; (ii) we consider instability of two-dimensional force-free system of magnetic flux tubes; (iii) we consider merger of two zero total poloidal current magnetic flux tubes. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization. The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced reconnection, further accelerating the particles, but proceeds at a slower reconnection rate. Crab flares result from the initial explosive stages of magnetic island mergers of magnetic flux tubes produced in the bulk of nebula at intermediate polar regions. The post-termination shock plasma flow in the wind sectors with mild magnetization naturally generates large-scale highly magnetized structures. Internal kink-like instabilities lead to the formation of macroscopic current-carrying magnetic flux tubes that merge explosively.

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}$.

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.

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.

The relativistic double neutron star binary PSR J0737-3039 shows clear evidence of orbital phase-dependent wind-companion interaction, both in radio and X-rays. In this paper we present the results of timing analysis of PSR J0737-3039 performed during 2006 and 2011 XMM-Newton Large Programs that collected ~20,000 X-ray counts from the system. We detected pulsations from PSR J0737-3039A (PSR A) through the most accurate timing measurement obtained by XMM-Newton so far, the spin period error being of 2x10^-13 s. PSR A's pulse profile in X-rays is very stable despite significant relativistic spin precession that occurred within the time span of observations. This yields a constraint on the misalignment between the spin axis and the orbital momentum axis Delta_A ~6.6^{+1.3}_{-5.4} deg, consistent with estimates based on radio data. We confirmed pulsed emission from PSR J0737-3039B (PSR B) in X-rays even after its disappearance in radio. The unusual phenomenology of PSR B's X-ray emission includes orbital pulsed flux and profile variations as well as a loss of pulsar phase coherence on time scales of years. We hypothesize that this is due to the interaction of PSR A's wind with PSR B's magnetosphere and orbital-dependent penetration of the wind plasma onto PSR B closed field lines. Finally, the analysis of the full XMM-Newton dataset provided evidences of orbital flux variability (~7%) for the first time, involving a bow-shock scenario between PSR A's wind and PSR B's magnetosphere.

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.

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, $\sigma \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.

We investigate the dynamics of bow shock nebulae created by pulsars moving supersonically through a partially ionized interstellar medium. A fraction of interstellar neutral hydrogen atoms penetrating into the tail region of a pulsar wind will undergo photo-ionization due to the UV light emitted by the nebula, with the resulting mass loading dramatically changing the flow dynamics of the light leptonic pulsar wind. Using a quasi 1-D hydrodynamic model of both non-relativistic and relativistic flow, and focusing on scales much larger than the stand-off distance, we find that if a relatively small density of neutral hydrogen, as low as $10^{-4}\,\text{cm}^{-3}$, penetrate inside the pulsar wind, this is sufficient to strongly affect the tail flow. Mass loading leads to the fast expansion of the pulsar wind tail, making the tail flow intrinsically non-stationary. The shapes predicted for the bow shock nebulae compare well with observations, both in H$\alpha$ and X-rays.

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.

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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of gamma-ray bursts. For a summary, we refer to the paper.

We present the results of 71.6 hours of observations of the Geminga pulsar (PSR J0633+1746) with the VERITAS very-high-energy gamma-ray telescope array. Data taken with VERITAS between November 2007 and February 2013 were phase-folded using a Geminga pulsar timing solution derived from data recorded by the XMM-\emph{Newton} and \emph{Fermi}-LAT space telescopes. No significant pulsed emission above 100 GeV is observed, and we report upper limits at the 95% confidence level on the integral flux above 135 GeV (spectral analysis threshold) of 4.0$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ and 1.7$\times10^{-13}$ s$^{-1}$ cm$^{-2}$ for the two principal peaks in the emission profile. These upper limits, placed in context with phase-resolved spectral energy distributions determined from five years of data from the \emph{Fermi}-LAT, constrain possible hardening of the Geminga pulsar emission spectra above $\sim$50 GeV.

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.

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 $ \epsilon_t $ to values well below the critical strain $ \epsilon_{crit}$ of a brittle fracture, $ \epsilon_t \sim 10^{-2}\epsilon_{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 $\rho \sim 10^{10} {\rm g cm}^{-3}$; this energy is detected as enhanced post-flare surface emission.

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.

In March 2013, a flaring episode from the Crab Nebula lasting ~2 weeks was detected by the Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). VERITAS provides simultaneous observations throughout this period. During the flare, the Fermi-LAT detected a 20-fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very-high-energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 * 10^-3 with 95% confidence.

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.

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.

We present solutions for Hall equilibria applicable to neutron star crusts. Such magnetic configurations satisfy a Grad-Shafranov-type equation, which is solved analytically and numerically. The solutions presented cover a variety of configurations, from purely poloidal fields connected to an external dipole to poloidal-toroidal fields connected to an external vacuum field, or fully confined within the star. We find that a dipole external field should be supported by a uniformly rotating electron fluid. The energy of the toroidal magnetic field is generally found to be a few percent of the total magnetic field energy for the fields with an external component. We discuss the evolution due to Ohmic dissipation which leads to slowing down of the electron fluid. We also find that the transition from an MHD equilibrium to a state governed by Hall effect, generates spontaneously an additional toroidal field in regions where the electron fraction changes.

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

(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.

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.

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.

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.

We present Hall equilibrium solutions for neutron stars crusts containing toroidal and poloidal magnetic field. Some simple cases are solved analytically while more complicated configurations are found numerically through a Gauss-Seidel elliptic partial differential equation solver.

We present the results of a joint observational campaign between the Green Bank radio telescope and the VERITAS gamma-ray telescope, which searched for a correlation between the emission of very-high-energy (VHE) gamma rays ($E_{\gamma} >$ 150 GeV) and Giant Radio Pulses (GRPs) from the Crab pulsar at 8.9 GHz. A total of 15366 GRPs were recorded during 11.6 hours of simultaneous observations, which were made across four nights in December 2008 and in November and December 2009. We searched for an enhancement of the pulsed gamma-ray emission within time windows placed around the arrival time of the GRP events. In total, 8 different time windows with durations ranging from 0.033 ms to 72 s were positioned at three different locations relative to the GRP to search for enhanced gamma-ray emission which lagged, led, or was concurrent with, the GRP event. Further, we performed separate searches on main pulse GRPs and interpulse GRPs and on the most energetic GRPs in our data sample. No significant enhancement of pulsed VHE emission was found in any of the preformed searches. We set upper limits of 5-10 times the average VHE flux of the Crab pulsar on the flux simultaneous with interpulse GRPs on single-rotation-period time scales. On $\sim$8-second time scales around interpulse GRPs, we set an upper limit of 2-3 times the average VHE flux. Within the framework of recent models for pulsed VHE emission from the Crab pulsar, the expected VHE-GRP emission correlations are below the derived limits.

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.

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 \Phi_\infty /(\pi c \hbar), where \Phi_\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.

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.

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 $\gamma$-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.

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.

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.

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.

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 \gamma-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 \gamma-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 \gamma-ray regions, all point to the inverse Compton origin of the high energy emission from majority of pulsars.

In the double pulsar system PSR J0737-3039A/B the strong wind produced by pulsar A distorts the magnetosphere of pulsar B. The influence of these distortions on the orbital-dependent emission properties of pulsar B can be used to determine the location of the coherent radio emission generation region in the pulsar magnetosphere. Using a model of the wind-distorted magnetosphere of pulsar B and the well defined geometrical parameters of the system, we determine the minimum emission height to be ~ 20 neutron star radii in the two bright orbital longitude regions. We can determine the maximum emission height by accounting for the amount of deflection of the polar field line with respect to the magnetic axis using the analytical magnetic reconnection model of Dungey and the semi-empirical numerical model of Tsyganenko. Both of these models estimate the maximum emission height to be ~ 2500 neutron star radii. The minimum and maximum emission heights we calculate are consistent with those estimated for normal isolated pulsars.

The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology using multi-frequency radio observations obtained with the Green Bank Telescope. We have characterized the periodic flux modulations previously discovered at 820 MHz by McLaughlin et al., and investigated the radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10^5). We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be detected. Their complex behavior is consistent with the confinement of the absorbing plasma in the dipolar magnetic field of pulsar B as suggested by Lyutikov & Thompson and such a geometric connection explains that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar A and the nominally estimated magnetic pressure of pulsar B.

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.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

The VERITAS telescope array has been operating smoothly since 2007, and has detected gamma-ray emission above 100 GeV from 40 astrophysical sources. These include blazars, pulsar wind nebulae, supernova remnants, gamma-ray binary systems, a starburst galaxy, a radio galaxy, the Crab pulsar, and gamma-ray sources whose origin remains unidentified. In 2009, the array was reconfigured, greatly improving the sensitivity. We summarize the current status of the observatory, describe some of the scientific highlights since 2009, and outline plans for the future.

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.

We present an analytical class of equilibrium solutions for the structure of relativistic sheared and rotating magnetized jets that contain no boundary current sheets. We demonstrate the overall dynamical stability of these solutions and, most importantly, a better numerical resistive stability than the commonly employed force-free structures which inevitably require the presence of dissipative surface currents. The jet is volumetrically confined by the external pressure, with no pressure gradient on the surface. We calculate the expected observed properties of such jets. Given the simplicity of these solution we suggest them as useful initial conditions for relativistic jet simulations.

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.

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 \Phi_\infty /(\pi c \hbar)$, where $\Phi_\infty \approx 2 \pi^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.

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 \gamma-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_{\gamma} scales linearly with the pulsar spin-down power \dot{E}, L_\gamma ~ \dot{E}, and not proportional to \sqrt{\dot{E}} as expected from potential-limited acceleration.

We report the detection of pulsed gamma rays from the Crab pulsar at energies above 100 Gigaelectronvolts (GeV) with the VERITAS array of atmospheric Cherenkov telescopes. The detection cannot be explained on the basis of current pulsar models. The photon spectrum of pulsed emission between 100 Megaelectronvolts (MeV) and 400 GeV is described by a broken power law that is statistically preferred over a power law with an exponential cutoff. It is unlikely that the observation can be explained by invoking curvature radiation as the origin of the observed gamma rays above 100 GeV. Our findings require that these gamma rays be produced more than 10 stellar radii from the neutron star.

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 $\Delta \theta \leq 1/\Gamma^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.

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 $\Gamma>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.

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.

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).

[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.

The oblique geometry of pulsar wind termination shock ensures that the Doppler beaming has a strong impact on the shock emission. We illustrate this using recent relativistic MHD simulations of the Crab Nebula and also show that the observed size, shape, and distance from the pulsar of the Crab Nebula inner knot are consistent with its interpretation as a Doppler-boosted emission from the termination shock. If the electrons responsible for the synchrotron gamma-rays are accelerated only at the termination shock then their short life-time ensures that these gamma-rays originate close to the shock and are also strongly effected by the Doppler beaming. As the result, bulk of the observed synchrotron gamma-rays of the Crab Nebula around 100 MeV may come from its inner knot. This hypothesis is consistent with the observed optical flux of the inner knot provided its optical-gamma spectral index is the same as the injection spectral index found in the Kennel & Coroniti model of the nebula spectrum. The observed variability of synchrotron gamma-ray emission can be caused by the instability of the termination shock discovered in recent numerical simulations. Given the small size of the knot, it is possible that the September 2010 gamma-ray flare of the Crab Nebula also came from the knot, though the actual mechanism remains unclear. The model predicts correlation of the temporal variability of the synchrotron gamma-ray flux in the Fermi and AGILE windows with the variability of the unpulsed optical flux from within 1 arcsec of the Crab pulsar.

(Abridged) The motion of a Schwarzschild black hole with velocity $v_0 = \beta_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 $\rho_{\rm ind}= B_0 \beta_0 /(2 \pi 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 $\rho_{\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 \beta_0$. The electromagnetic power of the jets is $L \approx (G M)^2 B_0^2 \beta_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.

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.

We mode Solar coronal mass ejections (CMEs) as expanding force-fee magnetic structures and find the self-similar dynamics of configurations with spatially constant \alpha, where {\bf J} =\alpha {\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.

To constrain the giant pulse (GP) emission mechanism and test the model of Lyutikov (2007) for GP emission, we have carried out a campaign of simultaneous observations of the Crab pulsar at gamma-ray (Fermi) and radio (Green Bank Telescope) wavelengths. Over 10 hours of simultaneous observations we obtained a sample of 2.1x10^4 giant pulses, observed at a radio frequency of 9 GHz, and 77 Fermi photons, with energies between 100 MeV and 5 GeV. The majority of GPs came from the interpulse (IP) phase window. We found no change in the GP generation rate within 10-120 s windows at lags of up to +-40 min of observed gamma-ray photons. The 95% upper limit for a gamma-ray flux enhancement in pulsed emission phase window around all GPs is 4 times the average pulsed gamma-ray flux from the Crab. For the subset of IP GPs, the enhancement upper limit, within the IP emission window, is 12 times the average pulsed gamma-ray flux. These results suggest that GPs, at least high-frequency IP GPs, are due to changes in coherence of radio emission rather than an overall increase in the magnetospheric particle density.

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.

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 \rho ~ r^{-\omega} steeper than \omega >4.

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.

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.

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.

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.

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.

To constrain the giant pulse (GP) emission mechanism and test the model of Lyutikov (2007) of GP emission, we are carrying out a campaign of simultaneous observations of the Crab pulsar between gamma-rays (Fermi) and radio wavelengths. The correlation between times of arrival of radio GPs and high-energy photons, whether it exists or not, will allow us to choose between different origins of GP emission and further constrain the emission physics. Our foremost goal was testing whether radio GPs are due to changes in the coherence of the radio emission mechanism, variations in the pair creation rate in the pulsar magnetosphere, or changes in the beaming direction. Accomplishing this goal requires an enormous number of simultaneous radio GPs and gamma-photons. Thus, we organized a radio observations campaign using the 42-ft telescope at the Jodrell Bank Observatory (UK), the 140-ft telescope, and the 100-m Robert C. Byrd Green Bank Telescope (GBT) at the Green Bank Observatory (WV). While the observations with the two first ones are ongoing, here we present the preliminary results of 20 hrs of observations with the GBT at the high frequency of 8.9 GHz. These particular observations were aimed to probe the model of GP emission by Lyutikov (2007) which predicts that GPs at frequencies > 4 GHz should be accompanied by gamma-ray photons of energies of 1-100 GeV.

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.

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.

(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.

The basic principles of General Theory of Relativity historically have been tested in gedanken experiments in rotating frame of references. One of the key issues, which still evokes a lot of controversy, is the centrifugal acceleration. Machabeli & Rogava (1994) argued that centrifugal acceleration reverse direction for particles moving radially with relativistic velocities within a "bead on a wire" approximation. We show that this result is frame-dependent and reflects a special relativistic dilution of time (as correctly argued by de Felice (1995)) and is analogous to freezing of motion on the black hole horizon as seen by a remote observer. It is a reversal of coordinate acceleration; there is no such effect as measured by a defined set of observers, e.g., proper and/or comoving. Frame-independent velocity of a "bead" with respect to stationary rotating observers increases and formally reaches the speed of light on the light cylinder. In general relativity, centrifugal force does reverse its direction at photon circular orbit, r=3M in Schwarzschild metric, as argued by Abramowicz (1990).

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.

We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X and soft gamma-ray data of four magnetars, including anomalous X-ray pulsars, and soft gamma repeaters. In this scenario, non-thermal magnetar spectra in the soft X-rays result from resonant cyclotron scattering of the thermal surface emission by hot magnetospheric plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while an additional component is still required to model the hard X-ray persistent magnetar emission recently discovered by INTEGRAL. The latter is an important component in terms of magnetars' luminosity, and cannot be neglected when modelling the soft X-ray part of the spectrum.

In this White Paper, we emphasize the need for and the important role of plasma astrophysics in the studies of formation, evolution of, and feedback by Active Galaxies. We make three specific recommendations: 1) We need to significantly increase the resolution of VLA, perhaps by building an EVLA-II at a modest cost. This will provide the angular resolution to study jets at kpc scales, where, for example, detailed Faraday rotation diagnosis can be done at 1GHz transverse to jets; 2) We need to build coordinated programs among NSF, NASA, and DOE to support laboratory plasma experiments (including liquid metal) that are designed to study key astrophysical processes, such as magneto-rotational instability (origin of angular momentum transport), dynamo (origin of magnetic fields), jet launching and stability. Experiments allowing access to relativistic plasma regime (perhaps by intense lasers and magnetic fields) will be very helpful for understanding the stability and dissipation physics of jets from Supermassive Black Holes; 3) Again through the coordinated support among the three Agencies, we need to invest in developing comprehensive theory and advanced simulation tools to study the accretion disks and jets in relativistic plasma physics regime, especially in connecting large scale fluid scale phenomena with relativistic kinetic dissipation physics through which multi-wavelength radiation is produced.

Multi-wavelength studies at radio, infrared, optical, X-ray, and TeV wavelengths have discovered probable counterparts to many Galactic sources of GeV emission detected by EGRET. These include pulsar wind nebulae, high mass X-ray binaries, and mixed morphology supernova remnants. Here we provide an overview of the observational properties of Galactic sources which emit across 19 orders of magnitude in energy. We also present new observations of several sources.

The double pulsar PSR J0737-3039A/B consists of two neutron stars in a highly relativistic orbit that displays a roughly 30-second eclipse when pulsar A passes behind pulsar B. Describing this eclipse of pulsar A as due to absorption occurring in the magnetosphere of pulsar B, we successfully use a simple geometric model to characterize the observed changing eclipse morphology and to measure the relativistic precession of pulsar B's spin axis around the total orbital angular momentum. This provides a test of general relativity and alternative theories of gravity in the strong-field regime. Our measured relativistic spin precession rate of 4.77 (+0.66,-0.65) degrees per year (68% confidence level) is consistent with that predicted by general relativity within an uncertainty of 13%.

(Abridged) We present a systematic fit of a model of resonant cyclotron scattering (RCS) to the X-ray data of ten magnetars, including canonical and transient anomalous X-ray pulsars (AXPs), and soft gamma repeaters (SGRs). In this scenario, non-thermal magnetar spectra in the soft X-rays (i.e. below ~10 keV) result from resonant cyclotron scattering of the thermal surface emission by hot magnetospheric plasma. We find that this model can successfully account for the soft X-ray emission of magnetars, while using the same number of free parameters than the commonly used empirical blackbody plus power-law model. However, while the RCS model can alone reproduce the soft X-ray spectra of AXPs, the much harder spectra of SGRs below ~10 keV, requires the addition of a power-law component (the latter being the same component responsible for their hard X-ray emission). Although this model in its present form does not explain the hard X-ray emission of a few of these sources, we took this further component into account in our modeling not to overlook their contribution in the ~4-10 keV band. We find that the entire class of sources is characterized by magnetospheric plasma with a density which, at resonant radius, is about 3 orders of magnitudes higher than n_{GJ}, the Goldreich-Julian electron density. The inferred values of the intervening hydrogen column densities, are also in better agreement with more recent estimates inferred from the fit of single X-ray edges. For the entire sample of observations, we find indications for a correlation between the scattering depth and the electron thermal velocity, and the field strength. Moreover, in most transient anomalous X-ray pulsars the outburst state is characterized by a relatively high surface temperature which cools down during the decay.

We consider a Hartmann layer, stationary flow of a viscose and resistive fluid between two plates with superimposed transverse magnetic field, in the limit of gyrotropic plasma, when viscosity across the field is strongly suppressed. For zero cross-field viscosity, the problem is not well posed, since viscosity then vanishes on the boundaries and in the middle of the layer, where there is no longitudinal field. An additional arbitrarily small isotropic viscosity allows one to find magnetic field and velocity profiles which are independent of this viscosity floor and different from flows with isotropic viscosity. Velocity sharply rises in a thin boundary layer, which thickness depends both on the Hartmann number and on the Lundquist number of the flow. The implication of the work is that, in simulating ICM dynamics it is imperative to use numerical schemes which take into account anisotropic viscosity. Although magnetic fields are dynamically subdominant in the ICM they do determine its the dissipative properties, stability of embedded structures and the transition to turbulence.

Relativistic outflows carrying large scale magnetic fields have large inductive potential and may accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (with jet luminosity $\geq 10^{46}$ erg s$^{-1}$). A positively charged particle carried by such a plasma is in an unstable equilibrium if ${\bf B} \cdot \nabla \times {\bf v}< 0$, so that kinetic drift along the velocity shear would lead to fast, regular energy gain. The highest rigidity particles are accelerated most efficiently implying the dominance of light nuclei for energies above the ankle in our model: from a mixed population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably.

We discuss dissipative processes in strongly gyrotropic, nearly collisionless plasma in clusters of galaxies (ICM). First, we point out that Braginsky theory, which assumes that collisions are more frequent that the system's dynamical time scale, is inapplicable to fast, sub-viscous ICM motion. Most importantly, the electron contribution to collisional magneto-viscosity dominates over that of ions for short-scale Alfvenic motions. Thus, if a turbulent cascade develops in the ICM and propagates down to scales $\leq 1$ kpc, it is damped collisionally not on ions, but on electrons. Second, in high beta plasma of ICM, small variations of the magnetic field strength, of relative value $\sim 1/\beta$, lead to development of anisotropic pressure instabilities (firehose, mirror and cyclotron). Unstable wave modes may provide additional resonant scattering of particles, effectively keeping the plasma in a state of marginal stability. We show that in this case the dissipation rate of a laminar, subsonic, incompressible flows scales as inverse of plasma beta parameter. We discuss application to the problem of ICM heating.

A number of disparate observational and theoretical pieces of evidence indicate that, contrary to the conventional wisdom, neutron stars' closed field lines are populated by dense, hot plasma and may be responsible for producing some radio and high energy emission. This conclusion is based on eclipse modeling of the binary pulsar system PSR J0737-3039A/B (Lyutikov & Thompson 2005), a quantitative theory of Crab giant pulses (Lyutikov 2007) and a number of theoretical works related to production of non-thermal spectra in magnetars through resonant scattering. In magnetars, dense pair plasma is produced by twisting magnetic field lines and associated electric fields required to lift the particles from the surface. In long period pulsars, hot particles on closed field lines can be efficiently trapped by magnetic mirroring, so that relatively low supply rate, e.g. due to a drift from open field lines, may result in high density. In short period pulsars, magnetic mirroring does not work; large densities may still be expected at the magnetic equator near the Y-point.

We propose that Crab giant pulses are generated on closed magnetic field lines near the light cylinder via anomalous cyclotron resonance on the ordinary mode. Waves are generated in a set of fine, unequally spaced, narrow emission bands at frequencies much lower than a local cyclotron frequency. Location of emission bands is fitted to spectral structures seen by Eilek et al. (2006). To reproduce the data, the required density of plasma in the giant pulses emission region is much higher, by a factor $\sim 3 \times 10^5$, than the minimal Goldreich-Julian density. Emission is generated by a population of highly energetic particles with radiation-limited Lorentz factors $\gamma \sim 7 \times 10^7$, produced during occasional reconnection close to the Y point, where the last closed field lines approach the light cylinder.

The highest energy cosmic rays could be produced by drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (i.e. FR II radiogalaxies; radio loud quasars and high power BL Lacs). We show that in such scenarios proton synchrotron radiation can give rise to detectable photon fluxes at energies ranging from hundreds of keV to tens of MeV.

Dissipation of magnetic fields in Hall plasma of neutron star crusts may power persistent high energy emission of a class of strongly magnetized neutrons stars, magnetars. We consider development of a dissipative tearing mode in Hall plasma (electron MHD) and find that its growth rate increases with the wave number of perturbations, reaching a maximum value intermediate between resistive $\tau_r$ and Hall times $\tau_H$, $\Gamma\sim 1 / \sqrt{\tau_r \tau_H}$. We argue that the tearing mode may be the principal mechanism by which strong magnetic fields are dissipated in magnetars on times scale of $\sim 10^4-10^5$ yrs powering the persistent X-ray emisison.

(Abridged)Relativistic outflows carrying large scale magnetic fields have large inductive potential and may accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (with jet luminosity $\geq 10^{46}$ erg s$^{-1}$). We point out that a positively charged particle carried by such a plasma is in an unstable equilibrium if ${\bf B} \cdot \nabla \times {\bf v}< 0$, so that kinetic drift along the velocity shear would lead to fast, {\it regular} energy gain. We show that if a seed of pre-accelerated particles with energy below the ankle $\leq 10^{18}$ eV is present, these particles can be boosted to energies above $ 10^{19}$ eV. A key feature of the mechanism is that the highest rigidity particles are accelerated most efficiently implying the dominance of light nuclei for energies above the ankle in our model: from a mixed population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably. In addition, after a particle traversed large fraction of the available potential, its Larmor radius becomes of the order of the jet thickness. In this case, {\it the maximum possible acceleration rate of inverse relativistic gyro-frequency is achieved} and a particle finally become unconfined and leave the jet.

We present the first self-consistent theoretical models of magnetar spectra that take into account the combined effects of the stellar atmosphere and its magnetosphere. We find that the proton cyclotron lines that are already weakened by atmospheric effects become indistinguishable from the continuum for moderate scattering optical depths in the magnetosphere. Furthermore, the hard excess becomes more pronounced due to resonant scattering and the resulting spectra closely resemble the observed magnetar spectra. We argue that while the absence of proton cyclotron lines in the observed spectra are inconclusive about the surface field strengths of magnetars, the continuum carries nearly unique signatures of the field strength and can thus be used to infer this quantity. We fit our theoretical spectra with a phenomenological two-blackbody model and compare our findings to source spectra with existing two-blackbody fits. The field strengths that we infer spectroscopically are in remarkable agreement with the values inferred from the period derivative of the sources assuming a dipole spindown. These detailed predictions of line energies and equivalent widths may provide an optimal opportunity for measuring directly the surface field strengths of magnetars with future X-ray telescopes such as {\em Constellation}-X.

The X-ray spectra of the magnetar candidates are customarily fitted with an empirical, two component model: an absorbed blackbody and a power-law. However, the physical interpretation of these two spectral components is rarely discussed. It has been recently proposed that the presence of a hot plasma in the magnetosphere of highly magnetized neutron stars might distort, through efficient resonant cyclotron scattering, the thermal emission from the neutron star surface, resulting in the production of non-thermal spectra. Here we discuss the Resonant Cyclotron Scattering (RCS) model, and present its XSPEC implementation, as well as preliminary results of its application to Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters.

Many astrophysical sources of high energy emission, such as black hole magnetospheres, superstrongly magnetized neutron stars (magnetars), and probably relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts involve relativistically magnetically dominated plasma. In such plasma the energy density of magnetic field greatly exceeds the thermal and the rest mass energy density of particles. Therefore the magnetic field is the main reservoir of energy and its dissipation may power the bursting emission from these sources, in close analogy to Solar flares. One of the principal dissipative instabilities that may lead to release of magnetic energy is the tearing instability. In this paper we study, both analytically and numerically, the development of tearing instability in relativistically magnetically-dominated plasma using the framework of resistive magnetodynamics. We confirm and elucidate the previously obtained result on the growth rate of the tearing mode: the shortest growth time is the same as in the case of classical non-relativistic MHD, namely $\tau =\sqrt{\tau_a \tau_d}$ where $\tau_a$ is the \Alfven crossing time and $\tau_d$ is the resistive time of a current layer.

Sharp fronts observed by Chandra satellite between dense cool cluster cores moving with near-sonic velocity through the hotter intergalactic gas, require strong suppression of thermal conductivity across the boundary. This may be due to magnetic fields tangential to the contact surface separating two plasma components. We point out that a super-Alfvenic motion of a plasma cloud (a core of a merging galaxy) through a weakly magnetized intercluster medium leads to "magnetic draping", formation of a thin, strongly magnetized boundary layer with a tangential magnetic field. For supersonic cloud motion, M_s > 1, magnetic field inside the layer reaches near-equipartition values with thermal pressure. Typical thickness of the layer is L /M_A^2 << L, where L is the size of the obstacle (plasma cloud) moving with Alfven Mach number M_A >> 1. To a various degree, magnetic draping occurs both for sub- and supersonic flows, random and ordered magnetic fields and it does not require plasma compressibility. The strongly magnetized layer will thermally isolate the two media and may contribute to the Kelvin-Helmholtz stability of the interface. Similar effects occur for radio bubbles, quasi-spherical expanding cavities blown up by AGN jets; in this case the thickness of the external magnetized layer is smaller, L /M_A^3 << L.

Resonant cyclotron scattering of the surface radiation in the magnetospheres of neutron stars may considerably modify the emergent spectra and impede efforts to constraint neutron star properties. Resonant cyclotron scattering by a non-relativistic warm plasma in an inhomogeneous magnetic field has a number of unusual characteristics: (i) in the limit of high resonant optical depth, the cyclotron resonant layer is half opaque, in sharp contrast to the case of non-resonant scattering. (ii) The transmitted flux is on average Compton up-scattered by ~ $1+ 2 beta_T$, where $\beta_T$ is the typical thermal velocity in units of the velocity of light; the reflected flux has on average the initial frequency. (iii) For both the transmitted and reflected fluxes the dispersion of intensity decreases with increasing optical depth. (iv) The emergent spectrum is appreciably non-Plankian while narrow spectral features produced at the surface may be erased. We derive semi-analytically modification of the surface Plankian emission due to multiple scattering between the resonant layers and apply the model to anomalous X-ray pulsar 1E 1048.1--5937. Our simple model fits just as well as the ``canonical'' magnetar spectra model of a blackbody plus power-law.

The Swift X-Ray Telescope often observes a rapidly decaying X-ray emission stretching to as long as $ t \sim 10^3$ seconds after a conventional prompt phase. This component is most likely due to a prompt emission viewed at large observer angles $\theta > 1/\Gamma$, where $\theta\sim 0.1$ is a typical viewing angle of the jet and$\Gamma\geq 100$ is the Lorentz factor of the flow during the prompt phase. This can be used to estimate the prompt emission radii, $r_{em} \geq 2 t c/\theta^2 \sim 6 \times 10^{15}$ cm. These radii are much larger than is assumed within a framework of a fireball model. Such large emission radii can be reconciled with a fast variability, on time scales as short as milliseconds, if the emission is beamed in the bulk outflow frame, e.g. due to a random relativistic motion of ''fundamental emitters''. This may also offer a possible explanation for X-ray flares observed during early afterglows.

I describe electromagnetic model of gamma ray bursts and contrast its main properties and predictions with hydrodynamic fireball model and its magnetohydrodynamical extension. The electromagnetic model assumes that rotational energy of a relativistic, stellar-mass central source (black-hole--accretion disk system or fast rotating neutron star) is converted into magnetic energy through unipolar dynamo mechanism, propagated to large distances in a form of relativistic, subsonic, Poynting flux-dominated wind and is dissipated directly into emitting particles through current-driven instabilities. Thus, there is no conversion back and forth between internal and bulk energies as in the case of fireball model. Collimating effects of magnetic hoop stresses lead to strongly non-spherical expansion and formation of jets. Long and short GRBs may develop in a qualitatively similar way, except that in case of long bursts ejecta expansion has a relatively short, non-relativistic, strongly dissipative stage inside the star. Electromagnetic and fireball models (as well as strongly and weakly magnetized fireballs) lead to different early afterglow dynamics, before deceleration time. Finally, I discuss the models in view of latest observational data in the Swift era.

(Abridged) We propose that giant flares on Soft Gamma-Ray Repeaters produce relativistic, strongly magnetized, weakly baryon loaded magnetic clouds, somewhat analogous to solar coronal mass ejection (CME) events. Flares are driven by unwinding of internal non-potential magnetic fields which leads to slow build-up of magnetic energy outside of the neutron star. For large magnetospheric currents, corresponding to a large twist of external magnetic field, magnetosphere becomes dynamically unstable on \Alfven crossing times scale of inner magnetosphere, $t_A \sim R_{NS}/c \sim 30 \mu$sec. Released magnetic energy results in formation of a strongly magnetized, pair-loaded, quasi-spherically expanding flux rope, topologically connected by magnetic field to the neutron star during the prompt flare emission. Magnetic stresses of the tied flux rope lead to late collimation of the expansion, on time scales longer than giant flare duration. Relativistic bulk motion of the expanding magnetic cloud, directed at an angle $\theta \sim 135^\circ$ to the line of sight (away from the observer), results in a strongly non-spherical forward shock with observed non-relativistic apparent expansion and bulk motion velocities $\beta_{app} \sim \cot \theta/2 \sim 0.4 $ at times of first radio observations approximately one week after the burst. Interaction with a shell of wind-shocked ISM and then with the unshocked ISM leads to deceleration to non-relativistic velocities approximately one month after the flare.

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.

Zhang and Kobayashi (2004) attempted to calculate early afterglow emission from a system of forward and reverse shocks in GRB outflows for the case of magnetized ejecta. We point out a fundamental error in the underlying dynamical model. According to the authors, energy and momentum carried by the magnetic field of the ejecta are not transfered to the forward shock. This is an incorrect assumption that invalidates the results.

In binary radio pulsar system J0737-3039, slow pulsar B shows orbital modulations of intensity, being especially bright at two short orbital phases. We propose 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. In our model, pulsar B is intrinsically bright at all times but its radiation beam misses the Earth at most orbital phases. We employ a simple model of distorted B magnetosphere using stretching transformations of Euler potentials of dipolar fields. To fit observations we use parameters of pulsar B derived from modeling of A eclipses (Lyutikov and Thompson 2005). The model reproduces two bright regions approximately at the observed orbital phases, explains variations of the pulse shape between them and regular timing residuals within each emission window. It also makes predictions for timing properties and secular variations of pulsar B profiles.

(Abridged) In the binary radio pulsar system J0737-3039, the faster pulsar A is eclipsed once per orbit. We construct a simple geometric model which successfully reproduces the eclipse light curves, based on the idea that the radio pulses are attenuated by synchrotron absorption on the closed magnetic field lines of pulsar B. The model explains most of the properties of the eclipse: its asymmetric form, the nearly frequency-independent duration, and the modulation of the brightness of pulsar A at both once and twice the rotation frequency of pulsar B in different parts of the eclipse. This detailed agreement confirms the dipolar structure of the star's poloidal magnetic field. The model makes clear predictions for the degree of linear polarization of the transmitted radiation. The weak frequency dependence of the eclipse duration implies that the absorbing plasma is relativistic, with a density much larger than the corotation charge density. Such hot, dense plasma can be effectively stored in the outer magnetosphere, where cyclotron cooling is slow. The gradual loss of particles inward through the cooling radius is compensated by an upward flux driven by a fluctuating component of the current, and by the pumping of magnetic helicity on the closed field lines. The trapped particles are heated to relativistic energies by the damping of magnetospheric turbulence and, at a slower rate, by the absorption of the radio emission of the companion pulsar.

We discuss the dynamical role of magnetic fields in GRB explosions in the framework of the fireball and electromagnetic models. We contrast the observational predictions of the two models and argue that only the very early afterglow observations, almost coincident with the prompt phase, can be used as decisive tests of the ejecta content.

Most of energy carried by relativistic AGN jets remains undetected until hundreds of kiloparsecs where interaction with intergalactic medium produces hot spots. The jet's hidden energy is only partially dissipated at smaller scales, from parsecs to kiloparsecs. Several media may play the role of the ``prime mover'': ions, pairs or large scale magnetic fields. Analyzing VLBI polarization structures of relativistic parsec scale jets we conclude that large-scale magnetic fields can explain the salient polarization properties of parsec-scale AGN jets. This implies that large-scale magnetic fields carry a non-negligible fraction of jet luminosity. We also discuss the possibility that relativistic AGN jets may be electromagnetically (Poynting flux) dominated. In this case, dissipation of the toroidal magnetic field (and not fluid shocks) may be responsible for particle acceleration.

The traditional view of radio pulsar wind nebulae (PWN), encouraged by the Crab nebula's X-ray and radio morphologies, is that they are a result of the integrated history of their pulsars' wind. The radio emission should therefore be largely unaffected by recent pulsar activity, and hence minimally correlated with structures in the X-ray nebulae. Observations of several PWN, both stationary (sPWN) and rapidly moving (rPWN), now show clear morphological relationships between structures in the radio and X-ray with radio intensity variations on the order of unity. We present high-resolution X-ray and radio images of several PWN of both types and discuss the morphological relationships between the two wavebands.

(Abridged) We consider the polarization properties of optically thin synchrotron radiation emitted by relativistically moving electron--positron jets carrying large-scale helical magnetic fields. In our model, the jet is cylindrical, and the emitting plasma moves parallel to the jet axis with a characteristic Lorentz factor $\Gamma$. We draw attention to the strong influence that the bulk relativistic motion of the emitting relativistic particles has on the observed polarization. We conclude that large-scale magnetic fields can explain the salient polarization properties of parsec-scale AGN jets. Since the typical degrees of polarization are $\leq 15%$, the emitting parts of the jets must have comparable rest-frame toroidal and poloidal fields. In this case, most relativistic jets are strongly dominated by the toroidal magnetic field component in the observer's frame, $B_\phi/B_z \sim \Gamma$. We also discuss the possibility that relativistic AGN jets may be electromagnetically (Poynting flux) dominated. In this case, dissipation of the toroidal magnetic field (and not fluid shocks) may be responsible for particle acceleration.

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.

At the early, ejecta dominated, stage of supernova remnant (SNR) expansion a fraction of the swept-up circumstellar magnetic field is dynamically compressed to approximate equipartition at the contact discontinuity separating the SN progenitor's wind (or the ISM) and the ejecta. We propose that the thin non-thermal X-ray filaments observed by the Chandra satellite in several young SNRs are associated with such ``pile-up'' of the magnetic field. We use a one-dimensional diffusion-convection transport equation to describe the propagation of non-thermal electrons near the contact discontinuity of a young SNR and to calculate spatially resolved emission spectra in the X-ray and TeV bands. The results suggest that the high-energy electrons are possibly accelerated at the forward shock, and emitting efficiently only when they diffuse into regions of high magnetic field near the contact discontinuity. Much more likely, though, is that they are locally accelerated at the contact discontinuity, in which case the acceleration cannot be related to Fermi-type processes and should occur due to other plasma mechanisms. As a consequence, the forward shock in young SNRs is inconspicuous and often unobservable, similar to that in the Crab nebular.

(Abridged) We interpret gamma ray bursts as relativistic, electromagnetic explosions. Specifically, we propose that they are created when a rotating, relativistic, stellar-mass progenitor loses much of its rotational energy in the form of a Poynting flux during an active period lasting $\sim 100$ s. Initially, a non-spherically symmetric, electromagnetically-dominated bubble expands non-relativistically inside the star, most rapidly along the rotational axis of the progenitor. After the bubble breaks out from the stellar surface and most of the electron-positron pairs annihilate, the bubble expansion becomes highly relativistic. After the end of the source activity most of the electromagnetic energy is concentrated in a thin shell inside the contact discontinuity between the ejecta and the shocked circumstellar material. This electromagnetic shell pushes a relativistic blast wave into the circumstellar medium. Current-driven instabilities develop in this shell at a radius $\sim3\times10^{16}$ cm and lead to dissipation of magnetic field and acceleration of pairs which are responsible for the $\gamma$-ray burst. At larger radii, the energy contained in the electromagnetic shell is mostly transferred to the preceding blast wave. Particles accelerated at the forward shock may combine with electromagnetic field from the electromagnetic shell to produce the afterglow emission.

Observations by the {\RHESSI} satellite of large polarization of the prompt $\gamma$-ray emission from the Gamma Ray Burst GRB021206 \citep{coburn03} imply that the magnetic field coherence scale is larger than the size of the visible emitting region, $\sim R/\Gamma$, where $R$ is the radius of the flow, $\Gamma$ is the associated Lorentz factor. Such fields cannot be generated in a causally disconnected, hydrodynamically dominated outflow. Electromagnetic models of GRBs \citep{lyutikov02}, in which large scale, dynamically dominant, magnetic fields are present in the outflow from the very beginning, provide a natural explanation of this large reported linear polarization. We derive Stokes parameters of synchrotron emission of a relativistically moving plasma with a given magnetic field configuration and calculate the pulse averaged polarization fraction of the emission from a relativistically expanding shell carrying global toroidal magnetic field. For viewing angles larger than $1/\Gamma$ the observed patch of the emitting shell has almost homogeneous magnetic field, producing a large fractional polarization ($56% $ for a power-law energy distribution of relativistic particles $dn/d\epsilon \propto \epsilon^{-3}$). The maximum polarization is smaller than the theoretical upper limit for a stationary plasma in uniform magnetic field due to relativistic kinematic effects.

We briefly review the electromagnetic model of Gamma Ray Bursts and then discuss how various models account for high prompt polarization. We argue that if polarization is confirmed at a level $\Pi \geq 10%$ the internal shock model is excluded.

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 $ \tau_A$ and a long resistive time scale $\tau_R$: $\Gamma \sim 1/(\tau_R \tau_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.

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.

We study the dynamics of relativistic electromagnetic explosions as a possible mechanism for the production of Gamma-Ray Bursts. We propose that a rotating relativistic stellar-mass progenitor loses much of its spin energy in the form of an electromagnetically-dominated outflow. After the flow becomes optically thin, it forms a relativistically expanding, non-spherically symmetric magnetic bubble - a ''cold fireball''. We analyze the structure and dynamics of such a cavity in the force-free approximation. During relativistic expansion, most of the magnetic energy in the bubble is concentrated in a thin shell near its surface (contact discontinuity). We suggest that either the polar current or the shell currents become unstable to electromagnetic instabilities at a radius $\sim10^{16}$ cm. This leads to acceleration of pairs and causes the $\gamma$-ray emission. At a radius $\sim10^{17}$ cm, the momentum contained in the electromagnetic shell will have been largely transferred to the surrounding blast wave propagating into the circumstellar medium. Particles accelerated at the fluid shock may combine with electromagnetic field from the electromagnetic shell to produce the afterglow emission.

Using the formulation of relativistic reconnection by Lyutikov & Uzdensky (2002) we estimate the upper possible rates of reconnection in pulsar winds using Bohm prescription for resistivity. We find that the velocity of plasma inflow into the reconnection layer may be relativistic, of the order of the speed of light in the plasma rest-frame. This in principle may allow efficient dissipation of the magnetic field energy in the wind and/or destruction of the toroidal magnetic flux. The efficiency of reconnection realized in pulsar winds remains an open question: it should depend both on the microphysical properties of plasma and on the three-dimensional structure of the reconnection flow.

The dynamics of the steady-state Sweet--Parker-type reconnection is analyzed in relativistic regime when energy density in the inflowing region is dominated by magnetic field. The structure of reconnection layer (its thickness, inflow and outflow velocities) depends on the ratio of two large dimensionless parameters of the problem - magnetization parameter $\sigma \gg 1$ (the ratio of the magnetic to particle energy-densities in the inflowing region) and the Lundquist number $S$. The inflow velocity may be relativistic (for $ S< \sigma$) or non-relativistic (for $ S > \sigma$), while the outflowing plasma is moving always relativisticly. For extremely magnetized plasmas with $\sigma \geq S^2$, the inflow four-velocity becomes of the order of the \Alfven four-velocity.

We present an X-ray and radio study of the wind nebula surrounding the central pulsar PSR J1811-1925 in the supernova remnant G11.2-0.3. Using high resolution data obtained with the Chandra X-ray observatory and with the VLA radio telescope we show the X-ray and radio emission is asymmetric around the pulsar, despite the latter's central position in the very circular shell. The new X-ray data allow us to separate the synchrotron emission of the pulsar wind nebula from the surrounding thermal emission and that from the pulsar itself. Based on X-ray data from two epochs, we observe temporal variation of the location of X-ray hot spots near the pulsar, indicating relativistic motion. We compare thermal emission observed within the shell, which may be associated with the forward shock of the pulsar wind nebula, to thermal emission from a nearby portion of the remnant shell, the temperature of which implies an expansion velocity consistent with the identification of the remnant with the historical event of 386 A.D. The measured X-ray and radio spectral indices of the nebula synchrotron emission are found to be consistent with a single synchrotron cooling break. The magnetic field implied by the break frequency is anomalously large given the apparent size and age of the nebula if a spherical morphology is assumed, but is consistent with a bipolar morphology.

We discuss properties of the expected radio emission from Soft Gamma-ray Repeaters (SGRs) during their bursting activity in the framework of the model of Thompson, Lyutikov and Kulkarni (2002), in which the high energy emission is powered by the dissipation of super-strong magnetic fields in the magnetospheres through reconnection-type events. Drawing on analogies with Solar flares we predict that coherent radio emission resembling solar type-III radio bursts may be emitted in SGRs during X-ray bursts. The radio emission should have correlated pulse profiles with X-rays, narrow-band-type radio spectrum with $\Delta \nu \leq \nu$ with the typical frequency $ \nu \geq 1 $ GHz, and, possibly, a drifting central frequency. We encourage sensitive radio observations of SGRs during the bursting activity.

(ABBREVIATED) We consider the structure of neutron star magnetospheres threaded by large-scale electrical currents, and the effect of resonant Compton scattering by the charge carriers (both electrons and ions) on the emergent X-ray spectra and pulse profiles. In the magnetar model for the SGRs and AXPs, these currents are maintained by magnetic stresses acting deep inside the star. We construct self-similar, force-free equilibria of the current-carrying magnetosphere with a power-law dependence of magnetic field on radius, B ~ r^(-2-p), and show that a large-scale twist softens the radial dependence to p < 1. The spindown torque acting on the star is thereby increased in comparison with a vacuum dipole. We comment on the strength of the surface magnetic field in the SGR and AXP sources, and the implications of this model for the narrow measured distribution of spin periods. A magnetosphere with a strong twist, B_\phi/B_\theta = O(1) at the equator, has an optical depth ~ 1 to resonant cyclotron scattering, independent of frequency (radius), surface magnetic field strength, or charge/mass ratio of the scattering charge. When electrons and ions supply the current, the stellar surface is also heated by the impacting charges at a rate comparable to the observed X-ray output of the SGR and AXP sources, if B_{dipole} ~ 10^{14} G. Redistribution of the emerging X-ray flux at the ion and electron cyclotron resonances will significantly modify the emerging pulse profile and, through the Doppler effect, generate a non-thermal tail to the X-ray spectrum. The sudden change in the pulse profile of SGR 1900+14 after the 27 August 1998 giant flare is related to an enhanced optical depth to electron cyclotron scattering, resulting from a sudden twist imparted to the external magnetic field.

The dynamics of relativistic magnetized mass loaded outflows carrying toroidal magnetic field is analyzed in the context of Pulsar Wind Nebulae (PWNs). Mass loading is very efficient in slowing down super-relativistic magnetized flows and weakening of relativistic shocks. We suggest that weakening of relativistic reverse shocks by mass loading in PWNs is responsible for the low radiative efficiencies of the majority of the PWNs. Mass loading may also result in a shock transition near the fast magnetosonic point; this is unlikely to happen in majority of PWNs. The evolution of magnetized mass loaded flows beyond the reverse shock is complicated: after initial deceleration to the minimal velocity required to transport the magnetic flux, the mass loaded flows have to {\it accelerate}. In order to be able to expand to infinity, magnetized flows should either become time dependent or destroy the toroidal magnetic flux by developing internal instabilities. Destruction of the magnetic flux initiated by mass loading may allow for the flow to slow down to sub-relativistic velocities and resolve the $\sigma$ paradox of the pulsar wind nebula.

We report the first X-ray detection of the radio pulsar PSR B1757-24 using the Chandra X-ray Observatory. The image reveals point-source emission at the pulsar position, consistent with being magnetospheric emission from the pulsar. In addition, we detect a faint tail extending nearly 20" east of the pulsar, in the same direction and with comparable morphology to the pulsar's well-studied radio tail. The X-ray tail is unlikely to be emission left behind following the passage of the pulsar, but rather is probably from synchrotron-emitting pulsar wind particles having flow velocity ~7000 km/s. Assuming the point-source X-ray emission is magnetospheric, the observed X-ray tail represents only ~0.01% of the pulsar's spin-down luminosity, significantly lower than the analogous efficiencies of most known X-ray nebulae surrounding rotation-powered pulsars.

We present high-resolution radio and X-ray studies of the composite supernova remnant G11.2-0.3. Using archival VLA data, we perform radio spectral tomography to measure for the first time the spectrum of the shell and plerion separately. We compare the radio morphology of each component to that observed in the hard and soft Chandra X-ray images. We measure the X-ray spectra of the shell and the emission in the interior and discuss the hypothesis that soft X-ray emission interior to the shell is the result of the expanding pulsar wind shocking with the supernova ejecta. We also see evidence for spatial variability in the hard X-ray emission near the pulsar, which we discuss in terms of ion mediated relativistic shocks.

We report the first X-ray detection of the radio pulsar PSR B1757-24 using the Chandra X-ray Observatory. We detect point-source emission at the pulsar position plus a faint tail extending nearly 20" east of the pulsar, in the same direction and with comparable morphology to the radio tail. Assuming the point-source X-ray emission is magnetospheric, the observed X-ray tail represents only ~0.01% of the pulsar's spin-down luminosity. This is significantly lower than the analogous efficiencies of most known X-ray nebulae surrounding rotation-powered pulsars. Assuming a non-thermal spectrum for the tail photons, we show that the tail is unlikely to be emission left behind following the passage of the pulsar, but rather is probably from synchrotron-emitting pulsar wind particles having flow velocity ~7000 km/s. We also show that there must be a significant break in the tail synchrotron spectrum between the radio and X-ray bands that is intrinsic to the particle spectrum. No emission is detected from the shell supernova remnant G5.4-1.2. The upper limits on remnant emission are unconstraining.

The dynamics of a relativistic blast wave propagating through magnetized medium is considered taking into account possible inhomogeneities of density and magnetic field and additional energy supply. Under the simplifying assumption of a spherically symmetric explosion in a medium with toroidal magnetic field self-similar solutions for the internal dynamics of the flow are derived. In the weakly magnetized case, when the bulk of the flow may be described by the unmagnetized solutions, there is a strongly magnetized sheath near the contact discontinuity (when it exists). Self-similar solutions inside the sheath are investigated. In the opposite limit of strongly magnetized upstream plasma new analytical self-similar solutions are found. Possible application to the physics of Gamma-Ray Bursts is discussed.

In the framework of the magnetar model for the Soft Gamma Repeaters and Anomalous X-ray Pulsars, we consider the structure of neutron star magnetospheres threaded by large-scale electrical currents. We construct self-similar, force-free equilibria under the assumption of axisymmetry and a power law dependence of magnetic field on radius, ${\bf B} \propto r^{-(2+p)}$. A large-scale twist of the field lines softens the radial dependence to $p < 1$, thereby accelerating the spindown torque with respect to a vacuum dipole. A magnetosphere with a strong twist ($B_\phi/B_\theta = O(1)$ at the equator) has an optical depth $\sim 1$ to resonant cyclotron scattering, independent of frequency (radius), surface magnetic field strength, or the charge/mass ratio of the scattering charge. We investigate the effects of the resonant Compton scattering by the charge carriers (both electrons and ions) on the emergent X-ray spectra and pulse profiles.

In relativistic strongly magnetized winds outflowing from the fast-rotating compact progenitors of gamma-ray bursts (GRBs) there are three regions where powerful high-frequency emission may be generated: (i) the thermal photosphere, (ii) the region of the internal wind instability and (iii) the region of the wind interaction with an ambient gas. This results in a multicomponent structure of GRBs. The emission from the thermal photosphere may be observed as a weak precursor to the main burst. The precursor should have a blackbody-like spectrum with the mean energy of photons of $\sim 1$ MeV, and its intensity should be tens to hundreds of times smaller than that of the main GRB emission. Observations of such precursors with future $\gamma$-ray missions like GLAST can clarify the nature of bursters.

Diffractive and refractive magnetospheric scintillations may allow a direct testing of the plasma inside the light cylinder. Unusual electrodynamics of the strongly magnetized electron-positron plasma allow separation of the magnetospheric and interstellar scattering. The most distinctive feature of the magnetospheric scintillations is their independence on frequency. Diffractive scattering due to small scale inhomogeneities produce a scattering angle that may be as large as 0.1 radians, and a typical decorrelation time of $10^{-8}$ seconds. Refractive scattering due to large scale inhomogeneities is also possible, with a typical angle of $10^{-3}$ radians and a correlation time of the order of $10^{-4}$ seconds. Some of the magnetospheric propagation effects may have already been observed.

We show that accretion disks with relativisticly hot electrons which are coupled to cold ions is subject to inertial-type instability in which kinematic drag between electrons and ions produces radial acceleration of electrons.

Poynting flux driven outflows from magnetized rotators are a plausible explanation for gamma-ray burst engines. We suggest a new possibility for how such outflows might transfer energy into radiating particles. We argue that the Poynting flux drives non-linearly unstable large amplitude electromagnetic waves (LAEMW) which ``break'' at radii $r_t \sim 10^{14} $ cm where the MHD approximation becomes inapplicable. In the ``foaming'' (relativisticly reconnecting) regions formed during the wave breaks the random electric fields stochastically accelerate particles to ultrarelativistic energies which then radiate in turbulent electromagnetic fields. The typical energy of the emitted photons is a fraction of the fundamental Compton energy $ \epsilon \sim f \hbar c/r_e $ with $f \sim 10^{-3}$ plus additional boosting due to the bulk motion of the medium. The emission properties are similar to synchrotron radiation, with a typical cooling time $\sim 10^{-4}$ sec. During the wave break, the plasma is also bulk accelerated in the outward radial direction and at larger radii can produce afterglows due to the interactions with external medium. The near equipartition fields required by afterglow models maybe due to magnetic field regeneration in the outflowing plasma (similarly to the field generation by LAEMW of laser-plasma interactions) and mixing with the upstream plasma.

We describe the possible electromagnetic signals expected from the magnetospheric interactions of a neutron star binary prior to merger. We find that both radio and X-ray signals of detectable strength are possible. We discuss possible links with the phenomenon of gamma-ray bursts and describe the prospects for direct detection of these signals in searches for radio and X-ray transients.

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 .

We analyze mechanisms of the excitation of Alfv\'{e}n wave in pulsar magnetospheres as a possible source of pulsar radio emission generation. We find that Cherenkov excitation of obliquely propagating Alfv\'{e}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\'{e}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\'{e}n waves into escaping radiation. Unfortunately, no decisive conclusion about the effectiveness of such conversion can be made.

We apply a theory of wave propagation through a turbulent medium to the scattering of radio waves in pulsar magnetospheres. We find that under conditions of strong density modulation the effects of magnetospheric scintillations in diffractive and refractive regimes may be observable. The most distinctive feature of the magnetospheric scintillations is their independence on frequency. Results based on diffractive scattering due to small scale inhomogeneities give a scattering angle that may be as large as 0.1 radians, and a typical decorrelation time of $10^{-8}$ seconds. Refractive scattering due to large scale inhomogeneities is also possible, with a typical angle of $10^{-3}$ radians and a correlation time of the order of $10^{-4}$ seconds. Temporal variation in the plasma density may also result in a delay time of the order of $10^{-4}$ seconds. The different scaling of the above quantities with frequency may allow one to distinguish the effects of propagation through a pulsar magnetosphere from the interstellar medium. In particular, we expect that the magnetospheric scintillations are relatively more important for nearby pulsars when observed at high frequencies.

We show that intensity dependent depolarization of single pulses (e.g., Xiluoris et al. 1994) may be due to the nonlinear decay of the "upper" ordinary mode into an unpolarized extraordinary mode and a backward propagating wave. The decay occurs in the innermost parts of the pulsar magnetosphere for obliquely propagating O waves.

Beam instabilities in the strongly magnetized electron-positron plasma of pulsar magnetospheres are considered. We analyze the resonance conditions and estimate the growth rates of the Cherenkov and cyclotron instabilities of the ordinary (O), extraordinary (X) and Alfv\'{e}n modes in two limiting regimes: kinetic and hydrodynamic. As a preliminary step, we reconsider wave dispersion and polarization properties in a one dimensional pair plasma taking into account relativistic thermal effects. We then find the location of the Cherenkov and cyclotron resonance of the X, O and Alfv\'{e}n waves with the particles from the primary beam. The importance of the different instabilities as a source of the coherent pulsar radiation generation is then estimated taking into account the angular dependence of the growth rates and the limitations on the length of the coherent wave-particle interaction imposed by the curvature of the magnetic field lines. We conclude, that in the pulsar magnetosphere the Cherenkov-type instabilities occur in the hydrodynamic regimes, while the cyclotron-type instabilities occur in the kinetic regime. We argue, that electromagnetic cyclotron-type instabilities on the X, O and probably Alfv\'{e}n waves are more likely to develop in the pulsar magnetosphere.

Electromagnetic processes associated with a charged particle moving in a strong circular magnetic field is considered in cylindrical coordinates. We investigate the relation between the vacuum curvature emission and Cherenkov emission and argue, that for the superluminal motion of a particle in the inhomogeneous magnetic field in a dielectric, the combined effects of magnetic field inhomogeneity and presence of a medium give rise to the synergetic Cherenkov-curvature emission process. We find the conditions when the operator relations between electric field and electric displacement in cylindrical coordinates may be approximated by algebraic relations. For nonresonant electromagnetic wave the interaction with particles streaming along the curved magnetic field may be described in the WKB approximation. For the resonant waves interacting with superluminal particles we find a particular case of plane-wave-like approximation. This allows a fairly simple calculation of the dielectric tensor of a plasma in a weakly inhomogenous magnetic field. We find in this approximation the polarization of normal modes in the plasma, Cherenkov-curvature and Cherenkov-drift emissivities and growth rates.

A theory of pulsar radio emission generation, in which the observed waves are produced directly by maser-type plasma instabilities operating at the anomalous cyclotron-Cherenkov resonance $\omega- k_{\parallel} v_{\parallel} + \omega_B/ \gamma_{res}=0$ and the Cherenkov-drift resonance $\omega- k_{\parallel} v_{\parallel} - k_{\perp} u_d =0$, is capable of explaining the main observational characteristics of pulsar radio emission. The instabilities are due to the interaction of the fast particles from the primary beam and the tail of the distribution with the normal modes of a strongly magnetized one-dimensional electron-positron plasma. The waves emitted at these resonances are vacuum-like, electromagnetic waves that may leave the magnetosphere directly. In this model, the cyclotron-Cherenkov instability is responsible for core emission pattern and the Cherenkov-drift instability produces conal emission. The conditions for the development of the cyclotron-Cherenkov instability are satisfied for both typical and millisecond pulsars provided that the streaming energy of the bulk plasma is not very high $\gamma_p \approx 10$. In a typical pulsar the cyclotron-Cherenkov and Cherenkov-drift resonances occur in the outer parts of magnetosphere at $r_{res} \approx 10^9 cm$. This theory can account for various aspects of pulsar phenomenology including the morphology of the pulses, their polarization properties and their spectral behavior. We propose several observational tests for the theory. The most prominent prediction are the high altitudes of the emission region and the linear polarization of conal emission in the plane orthogonal to the local osculating plane of the magnetic field.

It is shown that induced Raman scattering of electromagnetic waves in the strongly magnetized electron-positron plasma of pulsar magnetosphere may be important for wave propagation and as an effective saturation mechanism for electromagnetic instabilities. The frequencies, at which strong Raman scattering occurs in the outer parts of magnetosphere, fall into the observed radio band. The typical threshold intensities for the strong Raman scattering are of the order of the observed intensities, implying that pulsar magnetosphere may be optically thick to Raman scattering of electromagnetic waves.

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.