Using numerical simulations we show that low-amplitude Alfv\'en waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually reconnect behind the pancake, generating a variable wind far stronger than the normal spindown wind of the magnetar. Repeating ejections drive blast waves in the amplified wind. We suggest that these ejections generate the simultaneous X-ray and radio bursts detected from SGR 1935+2154. A modest energy budget of the magnetospheric perturbation $\sim 10^{40}$ erg is sufficient to produce the observed bursts with the energy ratio $\mathcal{E}_X/\mathcal{E}_{\rm radio}\sim 10^5$. Our simulation predicts a narrow (a few ms) X-ray spike from the magnetosphere, arriving simultaneously with the radio burst emitted far outside the magnetosphere. This timing is caused by the extreme relativistic motion of the ejecta.

Recently, CHIME detected periodicity in the bursting rate of the repeating FRB 180916.J0158+65. In a popular class of models, the fast radio bursts (FRBs) are created by giant magnetic flares of a hyper-active magnetar driven by fast ambipolar diffusion in the core. We point out that in this scenario the magnetar is expected to precess freely with a period of hours to weeks. The internal magnetic field $B\sim 10^{16}$G deforms the star, and magnetic flares induce sudden changes in magnetic stresses. The resulting torques and displacements of the principal axes of inertia are capable of pumping a significant amplitude of precession. The anisotropy of the flaring FRB activity, combined with precession, implies a strong periodic modulation of the visible bursting rate. The ultra-strong field invoked in the magnetar model provides: (1) energy for the frequent giant flares, (2) the high rate of ambipolar diffusion, releasing the magnetic energy on the timescale $\sim 10^9$s, (3) the core temperature $T\approx 10^9$K, likely above the critical temperature for neutron superfluidity, (4) strong magnetospheric torques, which efficiently spin down the star, and (5) deformation with ellipticity $\epsilon> 10^{-6}$, much greater than the rotational deformation. These conditions result in a precession with negligible viscous damping, and can explain the observed 16 day period in FRB 180916.J0158+65.

Magnetars younger than one century are expected to be hyper active. Besides winds powered by rotation they generate frequent magnetic flares, which launch powerful blast waves into the wind. These internal shocks act as masers producing fast (millisecond) radio bursts (FRBs) with the following properties. (1) GHz radio emission occurs at radii $r\sim 10^{14}$ cm and lasts $\lesssim 1$ ms in observer's time. (2) Induced scattering in the surrounding wind does not suppress the radio burst. (3) The emission has linear polarization set by the magnetar rotation axis. (4) The emission drifts to lower frequencies during the burst, and its duration broadens at lower frequencies. (5) Blast waves in inhomogeneous winds may emit variable bursts; periodicity might appear on sub-ms timescales if the magnetar rotates with $\sim 1$ s period. However, the observed FRB structure is likely changed by lensing effects during propagation through the host galaxy. (6) The FRBs from magnetars are expected to repeat, with rare strong bursts (up to $\sim 10^{43}$ erg) or more frequent weak bursts. (7) When a repeating flare strikes the wind bubble in the tail of a previous flare, the FRB turns into a bright optical flash. Its luminosity may approach that of a supernova Ia and last seconds. The rate of these optical flashes in the universe is much lower than the FRB rate, however it may exceed the supernova rate. Locations of hyper-active magnetars in their host galaxies depend on how they form: magnetars created in supernovae explosions will trace star formation regions, and magnetars formed in mergers of compact objects will be offset. The merger magnetars are expected to be most energetic and particularly hyper-active.

We suggest that neutron star mergers eject an ultra-relativistic envelope of mass $m\sim 10^{-7}M_\odot$, which helps explain the gamma-ray burst from GW170817. One ejection mechanism is the ablation of the neutron star surface by the burst of neutrinos in the first $30\,\mu$s of the merger. Another, more efficient, mechanism for inflating the ultra-relativistic envelope is an internal shock in the massive ejecta from the merger. A strong shock is expected if the merger product is a magnetar, which emits a centrifugally accelerated wind. The shock propagates outward through the ejecta and accelerates in its outer layers at radii $r\sim 10^9-10^{10}$cm, launching an ultra-relativistic opaque envelope filled with $\sim 10^4$ photons per nucleon. The Lorentz factor profile of the envelope rises outward and determines its homologous expansion, which adiabatically cools the trapped photons. Once the magnetar loses its differential rotation and collapses into a black hole, a powerful jet forms. It drives a blast wave into the envelope, chasing its outer layers and eventually catching up with the envelope photosphere at $r\sim 10^{12}$cm. The ultra-relativistic photospheric breakout of the delayed blast wave emits a gamma-ray burst in a broad solid angle around the merger axis. This model explains the gamma-ray pulse from merger GW170817 with luminosity $L_\gamma\sim 10^{47}$erg/s, duration $\Delta t_{\rm obs}\sim 0.5$s, and characteristic photon energy $\sim 100$keV. The blast wave Lorentz factor at the envelope photosphere is consistent with $\Gamma\geq 5$ that we derive from the observed light curve of the burst. We suggest future tests of the model.

The remarkable null pulse coincident with the 2016 glitch in Vela rotation indicates a dynamical event involving the crust and the magnetosphere of the neutron star. We propose that a crustal quake associated with the glitch strongly disturbed the Vela magnetosphere and thus interrupted its radio emission. We present the first global numerical simulations of a neutron starquake. Our code resolves the elastodynamics of the entire crust and follows the evolution of Alfv\'en waves excited in the magnetosphere. We observe Rayleigh surface waves propagating away from the epicentre of the quake, around the circumference of the crust - an instance of the so-called whispering gallery modes. The Rayleigh waves set the initial spatial scale of the magnetospheric disturbance. Once launched, the Aflv\'en waves bounce in the closed magnetosphere, become dephased, and generate strong electric currents, capable of igniting electric discharge. Most likely, the discharge floods the magnetosphere with electron-positron plasma, quenching the radio emission. We find that the observed $\sim 0.2$ s disturbance is consistent with the damping time of the crustal waves if the crust is magnetically coupled to the superconducting core of the neutron star. The quake is expected to produce a weak X-ray burst of short duration.

We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. Magnetic reconnection under such conditions can operate in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically dominated plasmoids with a small thermal pressure. Compton drag clears cavities inside plasmoids and also affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in non-radiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two distinct components. (1) A mildly relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of ~100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives ~20% of the reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support magnetic reconnection as a mechanism powering hard X-ray emission from accreting black holes.

We present the results of our 8 year X-ray monitoring campaign on CXOU J171405.7-381031, the magnetar associated with the faint supernova remnant (SNR) CTB 37B. It is among the youngest by inferred spin-down age, and most energetic in spin-down power of magnetars, and may contribute, at least partially, to the GeV and TeV emission coincident with the SNR. We use a series of Chandra, XMM-Newton, and NuSTAR observations to characterize the timing and spectral properties of the magnetar. The spin-down rate of the pulsar almost doubled in <1 year and then decreased slowly to a more stable value. Its X-ray flux varied by approx, 50%, possibly correlated with the spin down rate. The 1-79 keV spectrum is well-characterized by an absorbed blackbody plus power-law model with an average temperature of kT=0.62+/-0.04 keV and photon index Gamma=0.92+/-0.16, or by a Comptonized blackbody with kT=0.55+/-0.04 keV and an additional hard power law with Gamma=0.70+/-0.20, In contrast with most magnetars, the pulsed signal is found to decrease with energy up to 6 keV, which is apparently caused by mixing with the hard spectral component that is pulse-phase shifted by approx. 0.43 cycles from the soft X-rays. We also analyze the spectrum of the nearby, diffuse nonthermal source XMMU J171410.8-381442, whose relation to the SNR is uncertain.

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.

Transitional millisecond pulsars provide a unique set of observational data for understanding accretion at low rates onto magnetized neutron stars. In particular, PSR~J1023+0038 exhibits a remarkable bimodality of the X-ray luminosity (low and high modes), pulsations extending from the X-ray to the optical band, GeV emission, and occasional X-ray flares. We discuss a scenario for the pulsar interaction with the accretion disk capable of explaining the observed behavior. We suggest that during the high mode the disk is truncated outside the light cylinder, allowing the pulsar wind to develop near the equatorial plane and strike the disk. The dissipative wind-disk collision energizes the disk particles and generates synchrotron emission, which peaks in the X-ray band and extends down to the optical band. The emission is modulated by the pulsar wind rotation, resulting in a pulse profile with two peaks 180\degr\ apart. This picture explains the high-mode luminosity, spectrum, and pulse profile (X-ray and optical) of PSR J1023+0038. It may also explain the X-ray flares as events of sudden increase in the effective disk cross section intercepting the wind. In contrast to previously proposed models, we suggest that the disk penetrates the light cylinder only during the low X-ray mode. This penetration suppresses the dissipation caused by the pulsar wind-disk collision, and the system enters the propeller regime. The small duty cycle of the propeller explains the low spindown rate of the pulsar.

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.

The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

Magnetar flares excite strong Alfv\'{e}n waves in the magnetosphere of the neutron star. The wave energy can (1) dissipate in the magnetosphere, (2) convert to "fast modes" and possibly escape, and (3) penetrate the neutron star crust and dissipate there. We examine and compare the three options. Particularly challenging are nonlinear interactions between strong waves, which develop a cascade to small dissipative scales. This process can be studied in the framework of force-free electrodynamics (FFE). We perform three-dimensional FFE simulations to investigate Alfv\'{e}n wave dissipation, how long it takes, and how it depends on the initial wave amplitude on the driving scale. In the simulations, we launch two large Alfv\'{e}n wave packets that keep bouncing on closed magnetic field lines and collide repeatedly until the full turbulence spectrum develops. Besides dissipation due to the turbulent cascade, we find that in some simulations spurious energy losses occur immediately in the first collisions. This effect occurs in special cases where the FFE description breaks. It is explained with a simple one-dimensional model, which we examine in both FFE and full magnetohydrodynamic settings. We find that magnetospheric dissipation through nonlinear wave interactions is relatively slow and more energy is drained into the neutron star. The wave energy deposited into the star is promptly dissipated through plastic crustal flows induced at the bottom of the liquid ocean, and a fraction of the generated heat is radiated from the stellar surface.

We examine the possible role of turbulence in feeding the emission of gamma-ray bursts (GRBs). Turbulence may develop in a GRB jet as the result of hydrodynamic or current-driven instabilities. The jet carries dense radiation and the turbulence cascade can be damped by Compton drag, passing kinetic fluid energy to photons through scattering. We identify two regimes of turbulence dissipation: (1) "Viscous" - the turbulence cascade is Compton damped on a scale $\ell_{\rm damp}$ greater than the photon mean free path $\ell_\star$. Then turbulence energy is passed to photons via bulk Comptonization by smooth shear flows on scale $\ell_\star<\ell_{\rm damp}$. (2) "Collisionless" - the cascade avoids Compton damping and extends to microscopic plasma scales much smaller than $\ell_\star$. The collisionless dissipation energizes plasma particles, which radiate the received energy; how the dissipated power is partitioned between particles needs further investigation with kinetic simulations. We show that the dissipation regime switches from viscous to collisionless during the jet expansion, at a critical value of the jet optical depth which depends on the amplitude of turbulence. Turbulent GRB jets are expected to emit nonthermal photospheric radiation. Our analysis also suggests revisions of turbulent Comptonization in black hole accretion disks discussed in previous works.

Internal shocks provide a plausible heating mechanism in the jets of gamma-ray bursts (GRBs). Shocks occurring below the jet photosphere are mediated by radiation. It was previously found that radiation mediated shocks (RMSs) inside GRB jets are inefficient photon producers, and the photons that mediate the RMS must originate from an earlier stage of the explosion. We show that this conclusion is valid only for non-magnetized jets. RMSs that propagate in moderately magnetized plasma develop a collisionless subshock which locally heats the plasma to a relativistic temperature, and the hot electrons emit copious synchrotron photons inside the RMS. We find that this mechanism is generally effective for mildly relativistic shocks and may be the main source of photons observed in GRBs. We derive a simple analytical formula for the generated photon number per proton, $Z$, which gives $Z=10^5$-$10^6$, consistent with observations. The photons are initially injected with low energies, well below the observed GRB peak. Their number is controlled by two main factors: (1) the abundance of electron-positron pairs created in the shock, which is self-consistently calculated, and (2) the upper limit on the brightness temperature of soft radiation set by induced downscattering. The injected soft photons that survive induced downscattering gain energy in the RMS via bulk Comptonization and shape its nonthermal spectrum.

On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A.

In this work we explore the evolution of magnetic fields inside strongly magnetized neutron stars in axisymmetry. We model numerically the coupled field evolution in the core and the crust. Our code models the Hall drift and Ohmic effects in the crust, the back-reaction on the field from magnetically-induced elastic deformation of the crust, the magnetic twist exchange between the crust and the core, and the drift of superconducting flux tubes inside the core. The correct hydromagnetic equilibrium is enforced in the core. We find that i) The Hall attractor found by Gourgouliatos and Cumming in the crust exists also for configurations when the B-field penetrates into the core. However, the evolution timescale for the core-penetrating fields is dramatically different than that of the fields confined to the crust. ii) The combination of Jones' flux tube drift and Ohmic dissipation in the crust can deplete the pulsar magnetic fields on the timescale of $150$ Myr if the crust is hot ($T\sim2\times10^8$ K), but acts on much slower timescales for cold neutron stars, such as recycled pulsars ($\sim 1.8 $ Gyr, depending on impurity levels). iii) The outward motion of superfluid vortices during the rapid spin-down of a young highly magnetized pulsar, can result in a partial expulsion of flux from the core when $B\lesssim 10^{13}$ G. However for $B\gtrsim2\times 10^{13}$ G, the combination of a stronger magnetic field and a longer spin period implies that the core field cannot be expelled.

Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks.

Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $\nu F_\nu$ shock spectra. The spectrum extends up to $\epsilon_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $\epsilon \gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $v\gamma \gtrsim 1$, where $\gamma \equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $v\gamma = 3$ and a photon-to-proton ratio of $n_\gamma/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $\theta_d \equiv kT_d/m_e c^2 \approx \epsilon_d/3$ is reached at an optical depth of $\tau \gg \theta_d^{-1}$ behind the shock, where $\epsilon_d$ is the average downstream photon energy.

Accreting pulsars power relativistic jets, and display a complex spin phenomenology. These behaviours may be closely related to the large-scale configuration of the star's magnetic field, shaped by its interaction with the surrounding accretion disc. Here we present the first relativistic simulations of the interaction of a pulsar magnetosphere with an accretion flow. Our axisymmetric simulations treat the magnetospheric, or coronal, regions using a resistive extension of force-free electrodynamics. The magnetic field is also evolved inside the disc, which is a defined volume with a specified velocity field and conductivity profile, found using an alpha-disc model. We study a range of disc alpha-parameters, thicknesses, magnetic Prandtl numbers, and inner truncation radii. We find that a large fraction of the magnetic flux in the pulsar's closed zone is opened by the intrusion of the disc, leading to an enhancement of the power extracted by the pulsar wind and the spin-down torque applied to the pulsar. In our simulations, most of the spin-down contribution to the stellar torque acts on open field lines. The efficiency of field-line opening is high in the simulations' long-term quasi-steady states, which implies that a millisecond pulsar's electromagnetic wind could be strong enough to power the observed neutron-star radio jets, and may significantly affect the pulsar's spin evolution.

The persistent radio counterpart of FRB 121102 is estimated to have $N\sim 10^{52}$ particles, energy $E_N\sim 10^{48}$ erg, and size $R\sim 10^{17}$ cm. The source can be a nebula inflated and heated by an intermittent outflow from a magnetar --- a neutron star powered by its magnetic (rather than rotational) energy. The object is young and frequently liberating energy in magnetic flares driven by accelerated ambipolar diffusion in the neutron star core, feeding the nebula and producing bright millisecond bursts. The particle number in the nebula is consistent with ion ejecta from giant flares. The nebula may also contain the freeze-out of electron-positron pairs $N_\pm\sim 10^{51}$ created months after the neutron star birth; the same mechanism offers an explanation for $N_\pm$ in the Crab nebula. The persistent source around FRB 121102 is likely heated by magnetic dissipation and internal waves excited by the magnetar ejecta. The volumetric heating by waves explains the nebula's enormous efficiency in producing radio emission. The repeating radio bursts are suggested to occur much closer to the magnetar, whose flaring magnetosphere drives ultrarelativistic internal shocks into the magnetar wind. The shocks are mediated by Larmor rotation that forms a GHz maser with the observed ms duration. Furthermore, the flare ejecta can become charge-starved and then convert to electromagnetic waves.

Magnetars are young and highly magnetized neutron stars which display a wide array of X-ray activity including short bursts, large outbursts, giant flares and quasi-periodic oscillations, often coupled with interesting timing behavior including enhanced spin-down, glitches and anti-glitches. The bulk of this activity is explained by the evolution and decay of an ultrastrong magnetic field, stressing and breaking the neutron star crust, which in turn drives twists of the external magnetosphere and powerful magnetospheric currents. The population of detected magnetars has grown to about 30 objects and shows unambiguous phenomenological connection with very highly magnetized radio pulsars. Recent progress in magnetar theory includes explanation of the hard X-ray component in the magnetar spectrum and development of surface heating models, explaining the sources' remarkable radiative output.

This paper examines the mechanism of internal shocks in opaque relativistic outflows, in particular in cosmological gamma-ray bursts. The shocks produce neutrino emission and affect the observed photospheric radiation from the explosion. They develop from internal compressive waves and can be of different types depending on the composition of the outflow: (1) Shocks in "photon gas," with negligible plasma inertia, have a unique structure determined by the force-free condition---zero radiation flux in the plasma rest frame. Radiation dominance over plasma inertia suppresses formation of collisionless shocks mediated by collective electromagnetic fields. (2) If the outflow is sufficiently magnetized, a strong collisionless subshock develops, which is embedded in a thicker radiation-mediated structure. (3) Waves in outflows with a free neutron component lead to dissipation through nuclear collisions. At large optical depths, shocks have a thickness comparable to the neutron free path, with an embedded radiation-mediated and collisionless subshocks. The paper also presents first-principle simulations of magnetized flows filled with photons, demonstrating formation of shocks and their structure. Simple estimates show that magnetized sub-photospheric shocks are efficient producers of photons and have a great impact on the observed photospheric radiation. The shock structure changes as the outflow expands toward its photosphere. The dissipation is accompanied by strong $e^\pm$ pair creation, and the $e^\pm$-dressed shock carries the photosphere with it up to two decades in radius, emitting a strong pulse of nonthermal radiation.

We review the physics of GRB production by relativistic jets that start highly opaque near the central source and then expand to transparency. We discuss dissipative and radiative processes in the jet and how radiative transfer shapes the observed nonthermal spectrum released at the photosphere. A comparison of recent detailed models with observations gives estimates for important parameters of GRB jets, such as the Lorentz factor and magnetization. We also discuss predictions for GRB polarization and neutrino emission.

A radiative mechanism is proposed for magnetic flares near luminous accreting black holes. It is based on recent first-principle simulations of magnetic reconnection, which show a hierarchical chain of fast-moving plasmoids. The reconnection occurs in a compact region (comparable to the black hole radius), and the chain experiences fast Compton cooling accompanied by electron-positron pair creation. The distribution of plasmoid speeds is shaped by radiative losses, and the self-regulated chain radiates its energy in hard X-rays. The mechanism is illustrated by Monte-Carlo simulations of the transfer of seed soft photons through the reconnection layer. The emerging radiation spectrum has a cutoff near 100 keV similar to the hard-state spectra of X-ray binaries and AGN. We discuss how the chain cooling differs from previous phenomenological emission models, and suggest that it can explain the hard X-ray activity of accreting black holes from first principles. Particles accelerated at the X-points of the chain produce an additional high-energy component, explaining "hybrid Comptonization" observed in Cyg X-1.

A gamma-ray burst (GRB) jet running into an external medium is expected to generate luminous GeV-TeV emission lasting from minutes to several hours. The high-energy emission results from inverse Compton upscattering of prompt and afterglow photons by shock-heated {\it thermal} plasma. At its peak the high-energy radiation carries a significant fraction of the power dissipated at the forward shock. We discuss in detail the expected TeV luminosity, using a robust "minimal" emission model. Then, using the statistical properties of the GRB population (luminosity function, redshift distribution, afterglow energy) we simulate the expected detection rates of GRBs by current and upcoming atmospheric Cherenkov instruments. We find that GRBs exploding into a low-density interstellar medium must produce TeV emission that would have already been detected by the currently operating Cherenkov telescopes. The absence of detections is consistent with explosions into a dense wind of the GRB progenitor. If the typical environment of long GRBs is a Wolf-Rayet progenitor wind, as suggested by the recent analysis of Fermi LAT data, then the upcoming Cherenkov Telescope Array (CTA) should be able to detect a fraction of GRBs that trigger the space-borne detectors. Since absorption by the extragalactic background light limits the detectability above 0.1 TeV for all but the most nearby bursts ($z\lesssim 1$), the reduced energy threshold of CTA is the key improvement over current instruments, which should boost the detection rate by at least a factor of 3.

We review electrodynamics of rotating magnetized neutron stars, from the early vacuum model to recent numerical experiments with plasma-filled magnetospheres. Significant progress became possible due to the development of global particle-in-cell simulations which capture particle acceleration, emission of high-energy photons, and electron-positron pair creation. The numerical experiments show from first principles how and where electric gaps form, and promise to explain the observed pulsar activity from radio waves to gamma-rays.

The MeV spectral peak of gamma-ray bursts (GRBs) is best explained as photospheric emission from a dissipative relativistic jet. The observed non-blackbody spectrum shows that sub-photospheric dissipation involves both thermal plasma heating and injection of nonthermal particles, which quickly cool through inverse Compton scattering and emission of synchrotron radiation. Synchrotron photons emitted around and above the photosphere are predicted to dominate the low-energy part of the GRB spectrum, starting from roughly a decade in energy below the MeV peak. We show that this leads to a unique polarization signature: a rise in GRB polarization toward lower energies. We compute the polarization degree of GRB radiation as a function of photon energy for a generic jet model, and show the predictions for GRBs 990123, 090902B and 110721A. The expected polarization is significant in the X-ray band, in particular for bursts similar to GRB 090902B. Radiation in the MeV peak (and at higher energies) is unpolarized as long as the jet is approximately uniform on angular scales $\delta\theta \gtrsim \Gamma^{-1}$ where $\Gamma$ is the bulk Lorentz factor of the jet.

The magnetospheres of magnetars are believed to be filled with electron-positron plasma generated by electric discharge. We present a first direct numerical experiment showing how the plasma is created in an axisymmetric closed magnetosphere. The $e^\pm$ discharge occurs in response to twisting of the magnetic field lines by a shear deformation of the magnetar surface, which launches electric currents into the magnetosphere. The simulation shows the formation of an electric "gap" with unscreened electric field ($\mathbf{E}\cdot \mathbf{B}\neq 0$) that continually accelerates particles along the magnetic field lines and sustains pair creation. The accelerating voltage is self-regulated to the threshold of the $e^\pm$ discharge. It controls the rate of energy release and the lifetime of the magnetic twist. The simulation follows the global evolution of the twisted magnetosphere over a long time and demonstrates its gradual resistive untwisting. A vacuum cavity forms near the star and expands, gradually erasing magnetospheric electric currents $j$. The active j-bundle shrinks with time and its footprints form shrinking hot spots on the magnetar surface bombarded by the created particles.

We explore the interaction between Hall waves and mechanical failures inside a magnetar crust, using detailed one-dimentional models that consider temperature-sensitive plastic flow, heat transport and cooling by neutrino emission, as well as the coupling of the crustal motion to the magnetosphere. We find that the dynamics is enriched and accelerated by the fast, short-wavelength Hall waves that are emitted by each failure. The waves propagate and cause failures elsewhere, triggering avalanches. We argue that these avalanches are the likely sources of outbursts in transient magnetars.

We examine four candidate mechanisms that could explain the high surface temperatures of magnetars. (1) Heat flux from the liquid core heated by ambipolar diffusion. It could sustain the observed surface luminosity $L_s\approx 10^{35}$ erg/s if core heating offsets neutrino cooling at a temperature $T_{core}>6\times 10^8$ K. This scenario is viable if the core magnetic field exceeds $10^{16}$ G and the heat-blanketing envelope of the magnetar has a light element composition. We find however that the lifetime of such a hot core should be shorter than the typical observed lifetime of magnetars. (2) Mechanical dissipation in the solid crust. This heating can be quasi-steady, powered by gradual (or frequent) crustal yielding to magnetic stresses. We show that it obeys a strong upper limit. As long as the crustal stresses are fostered by the field evolution in the core or Hall drift in the crust, mechanical heating is insufficient to sustain persistent $L_s\approx 10^{35}$ erg/s. The surface luminosity is increased in an alternative scenario of mechanical deformations triggered by external magnetospheric flares. (3) Ohmic dissipation in the crust, in volume or current sheets. This mechanism is inefficient because of the high conductivity of the crust. Only extreme magnetic configurations with crustal fields $B>10^{16}$ G varying on a 100 meter scale could provide $L_s\approx 10^{35}$ erg/s. (4) Bombardment of the stellar surface by particles accelerated in the magnetosphere. This mechanism produces hot spots on magnetars. Observations of transient magnetars show evidence for external heating.

We report on simultaneous NuSTAR and XMM-Newton observations of the magnetar 1E 1048.1$-$5937 obtained in July 2013, along with Rossi X-ray Timing Explorer (RXTE) data for the same source obtained from December 2002 through March 2012. The NuSTAR data provide a clear detection of this magnetar's persistent emission up to 20 keV. We detect a previously unreported small secondary peak in the average pulse profile in the 7-10 keV band, which grows to an amplitude comparable to that of the main peak in the 10--20 keV band. We show using RXTE data that this secondary peak is likely transient. We find that the pulsed fraction increases with energy from a value of $\sim$0.55 at $\sim$2 keV to a value of $\sim$0.75 near 8 keV but shows evidence for decreasing at higher energies. After filtering out multiple bright X-ray bursts during the observation, we find that the phase-averaged spectrum from combined NuSTAR and XMM data is well described by an absorbed double blackbody plus power-law model, with no evidence for the spectral turn-up near $\sim$10 keV as has been seen in some other magnetars. Our data allow us to rule out a spectral turn-up similar to those seen in magnetars 4U 0142+61 and 1E 2259+586 of $\Delta\Gamma >= 2$, where $\Delta\Gamma$ is the difference between the soft-band and hard-band photon indexes. The absence of a significant spectral turn-up is consistent with what has been observed from a particularly active subset of magnetars having high spin-inferred magnetic fields, and with previously reported trends suggesting the degree of spectral turn-up is correlated with spin-down rate and/or spin-inferred magnetic field strength.

We present global radiative transfer models for heated relativistic jets. The simulations include all relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation to and beyond the photosphere of the jet. The transfer models are compared with three gamma-ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different spectra. The models provide good fits to the observed spectra in all three cases, and we obtain estimates for the jet magnetization parameter $\varepsilon_{\rm B}$ and the Lorentz factor $\Gamma$. In the small sample of three bursts, $\varepsilon_{\rm B}$ varies between 0.01 and 0.05, and $\Gamma$ varies between 400 and 1200.

The interaction of a rotating star's magnetic field with a surrounding plasma disk lies at the heart of many questions posed by neutron stars in X-ray binaries. We consider the opening of stellar magnetic flux due to differential rotation along field lines coupling the star and disk, using a simple model for the disk-opened flux, the torques exerted on the star by the magnetosphere, and the power extracted by the electromagnetic wind. We examine the conditions under which the system enters an equilibrium spin state, in which the accretion torque is instantaneously balanced by the pulsar wind torque alone. For magnetic moments, spin frequencies, and accretion rates relevant to accreting millisecond pulsars, the spin-down torque from this enhanced pulsar wind can be substantially larger than that predicted by existing models of the disk-magnetosphere interaction, and is in principle capable of maintaining spin equilibrium at frequencies less than 1 kHz. We speculate that this mechanism may account for the non-detection of frequency increases during outbursts of SAX J1808.4-3658 and XTE J1814-338, and may be generally responsible for preventing spin-up to sub-millisecond periods. If the pulsar wind is collimated by the surrounding environment, the resulting jet can satisfy the power requirements of the highly relativistic outflows from Cir X-1 and Sco X-1. In this framework, the jet power scales relatively weakly with accretion rate, L_j ~ Mdot^{4/7}, and would be suppressed at high accretion rates only if the stellar magnetic moment is sufficiently low.

Microquasars are stellar-mass black holes accreting matter from a companion star and ejecting plasma jets at almost the speed of light. They are analogues of quasars that contain supermassive black holes of $10^6$ to $10^{10}$ solar masses. Accretion in microquasars varies on much shorter timescales than in quasars and occasionally produces exceptionally bright X-ray flares. How the flares are produced is unclear, as is the mechanism for launching the relativistic jets and their composition. An emission line near 511 kiloelectronvolts has long been sought in the emission spectrum of microquasars as evidence for the expected electron-positron plasma. Transient high-energy spectral features have been reported in two objects, but their positron interpretation remains contentious. Here we report observations of $\gamma$-ray emission from the microquasar V404 Cygni during a recent period of strong flaring activity. The emission spectrum around 511 kiloelectronvolts shows clear signatures of variable positron annihilation, which implies a high rate of positron production. This supports the earlier conjecture that microquasars may be the main sources of the electron-positron plasma responsible for the bright diffuse emission of annihilation $\gamma$-rays in the bulge region of our Galaxy. Additionally, microquasars could be the origin of the observed megaelectronvolt continuum excess in the inner Galaxy.

The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 orders of magnitude larger than in supernova remnants) render the gas far upstream of the shock neutral and shielded from ionizing radiation. The amplification of the magnetic field needed for diffusive shock acceleration requires ionized gas, thus confining the acceleration process to a narrow photo-ionized layer immediately ahead of the shock. Based on the growth rate in this layer of the hybrid non-resonant cosmic ray current-driven instability (considering also ion-neutral damping), we quantify the maximum particle energy, Emax, across the range of shock velocities and upstream densities of interest. We find values of Emax ~ 10 GeV - 10 TeV, which are broadly consistent with the inferred spectral cut-offs, but which could also in principle lead to emission extending to higher energies >100 GeV accessible to atmosphere Cherenkov telescopes, such as the planned Cherenkov Telescope Array (CTA). Detecting TeV neutrinos with IceCube in hadronic scenarios appears to be more challenging, although the prospects are improved if the shock power during the earliest, densest phases of the nova outburst is higher than is implied by the observed GeV light curves, due to downscattering of the gamma-rays by electrons within the ejecta. Novae provide ideal nearby laboratories to study magnetic field amplification and the onset of cosmic ray acceleration, because other time-dependent sources (e.g. radio supernovae) typically occur too distant to detect as gamma-ray sources.

We present temporal and spectral analysis of simultaneous 0.5-79 keV Swift-XRT and NuSTAR observations of the magnetar 4U 0142+61. The pulse profile changes significantly with photon energy between 3 and 35 keV. The pulse fraction increases with energy, reaching a value of ~20%, similar to that observed in 1E 1841-045 and much lower than the ~80% pulse fraction observed in 1E 2259+586. We do not detect the 55-ks phase modulation reported in previous Suzaku-HXD observations. The phase-averaged spectrum of 4U 0142+61 above 20 keV is dominated by a hard power law with a photon index, $\Gamma$ ~ 0.65, and the spectrum below 20 keV can be described by two blackbodies, a blackbody plus a soft power law, or by a Comptonized blackbody model. We study the full phase-resolved spectra using the electron-positron outflow model of Beloborodov (2013). Our results are consistent with the parameters of the active j-bundle derived from INTEGRAL data by Hascoet et al. (2014). We find that a significant degeneracy appears in the inferred parameters if the footprint of the j-bundle is allowed to be a thin ring instead of a polar cap. The degeneracy is reduced when the footprint is required to be the hot spot inferred from the soft X-ray data.

We report on a 350-ks NuSTAR observation of the magnetar 1E 1841-045 taken in 2013 September. During the observation, NuSTAR detected six bursts of short duration, with $T_{90}<1$ s. An elevated level of emission tail is detected after the brightest burst, persisting for $\sim$1 ks. The emission showed a power-law decay with a temporal index of 0.5 before returning to the persistent emission level. The long observation also provided detailed phase-resolved spectra of the persistent X-ray emission of the source. By comparing the persistent spectrum with that previously reported, we find that the source hard-band emission has been stable over approximately 10 years. The persistent hard X-ray emission is well fitted by a coronal outflow model, where $e^{+/-}$ pairs in the magnetosphere upscatter thermal X-rays. Our fit of phase-resolved spectra allowed us to estimate the angle between the rotational and magnetic dipole axes of the magnetar, $\alpha_{mag}=0.25$, the twisted magnetic flux, $2.5\times10^{26}\rm \ G\ cm^2$, and the power released in the twisted magnetosphere, $L_j=6\times10^{36}\rm \ erg\ s^{-1}$. Assuming this model for the hard X-ray spectrum, the soft X-ray component is well fit by a two-blackbody model, with the hotter blackbody consistent with the footprint of the twisted magnetic field lines on the star. We also report on the 3-year Swift monitoring observations obtained since 2011 July. The soft X-ray spectrum remained stable during this period, and the timing behavior was noisy, with large timing residuals.

Magnetar flares generate Alfv\'en waves bouncing in the closed magnetosphere with energy up to $\sim 10^{46}$ erg. We show that on a 10-ms timescale the waves are transmitted into the star and form a compressed packet of high energy density. This packet strongly shears the stellar crust and initiates a plastic flow, heating the crust and melting it hundreds of meters below the surface. A fraction of the deposited plastic heat is eventually conducted to the stellar surface, contributing to the surface afterglow months to years after the flare. A large fraction of heat is lost to neutrino emission or conducted into the core of the neutron star.

One possible channel for black hole formation is the collapse of a rigidly rotating massive neutron star as it loses its angular momentum or gains excessive mass through accretion. It was proposed that part of the neutron star may form a debris disk around the black hole. Such short-lived massive disks could be the sources of powerful jets emitting cosmological gamma-ray bursts. Whether the collapse creates a disk depends on the equation of state of the neutron star. We survey a wide range of equations of states allowed by observations and find that disk formation is unfeasible. We conclude that this channel of black hole formation is incapable of producing powerful jets, and discuss implications for models of gamma-ray bursts.

Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle radiative transfer simulation and can be tested against observations. Here we examine a sample of 7 bursts with best GeV+optical data and test the model. We find that the observed light curves are in agreement with the theoretical predictions and allow us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With one possible exception of GRB 090510 (which is the only short burst in the sample) the ambient medium is consistent with a wind from a Wolf-Rayet progenitor. The wind density parameter $A=\rho r^2$ varies in the sample around $10^{11}$g/cm. The initial Lorentz factor of the blast wave varies from 200 to 540 and correlates with the burst luminosity. Radiative efficiency of the prompt emission in the sample is between 0.1 and 0.8. For the two bursts with detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes.

The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface.

Bright X-ray flares are routinely detected by the Swift satellite during the early afterglow of gamma-ray bursts, when the explosion ejecta drives a blast wave into the external medium. We suggest that the flares are produced as the reverse shock propagates into the tail of the ejecta. The ejecta is expected to contain a few dense shells formed at an earlier stage of the explosion. We show an example of how such dense shells form and describe how the reverse shock interacts with them. A new reflected shock is generated in this interaction, which produces a short-lived X-ray flare. The model provides a natural explanation for the main observed features of the X-ray flares --- the fast rise, the steep power-law decline, and the characteristic peak duration \Delta t /t= (0.1-0.3).

We present the first self-consistent global simulations of pulsar magnetospheres with operating $e^\pm$ discharge. We focus on the simple configuration of an aligned or anti-aligned rotator. The star is spun up from zero (vacuum) state to a high angular velocity, and we follow the coupled evolution of its external electromagnetic field and plasma particles using the "particle-in-cell" method. A plasma magnetosphere begins to form through the extraction of particles from the star; these particles are accelerated by the rotation-induced electric field, producing curvature radiation and igniting $e^\pm$ discharge. We follow the system evolution for several revolution periods, longer than required to reach a quasi-steady state. Our numerical experiment puts to test previous ideas for the plasma flow and gaps in the pulsar magnetosphere. We first consider rotators capable of producing pairs out to the light cylinder through photon-photon collisions. We find that their magnetospheres are similar to the previously obtained force-free solutions with a current sheet and the Y-point near the light cylinder. The magnetosphere continually ejects $e^\pm$ pairs and ions. Pair creation is sustained by a strong electric field along the current sheet. We observe powerful curvature and synchrotron emission from the current sheet, consistent with Fermi observations of gamma-ray pulsars. We then study pulsars that can only create pairs in the strong-field region near the neutron star, well inside the light cylinder. We find that both aligned and anti-aligned rotators relax to the "dead" state with suppressed pair creation and electric currents, regardless of the discharge voltage.

Ultraluminous X-ray sources (ULX) are off-nuclear point sources in nearby galaxies whose X-ray luminosity exceeds the theoretical maximum for spherical infall (the Eddington limit) onto stellar-mass black holes. Their luminosity ranges from $10^{40}$ erg s$^{-1} < L_X$(0.5 - 10 keV) $<10^{40}$ erg s$^{-1}$. Since higher masses imply less extreme ratios of the luminosity to the isotropic Eddington limit theoretical models have focused on black hole rather than neutron star systems. The most challenging sources to explain are those at the luminous end ($L_X$ > $10^{40}$ erg s$^{-1}$), which require black hole masses MBH >50 solar masses and/or significant departures from the standard thin disk accretion that powers bright Galactic X-ray binaries. Here we report broadband X-ray observations of the nuclear region of the galaxy M82, which contains two bright ULXs. The observations reveal pulsations of average period 1.37 s with a 2.5-day sinusoidal modulation. The pulsations result from the rotation of a magnetized neutron star, and the modulation arises from its binary orbit. The pulsed flux alone corresponds to $L_X$(3 - 30 keV) = $4.9 \times 10^{39}$ erg s$^{-1}$. The pulsating source is spatially coincident with a variable ULX which can reach $L_X$ (0.3 - 10 keV) = $1.8 \times 10^{40}$ erg s$^{-1}$. This association implies a luminosity ~100 times the Eddington limit for a 1.4 solar mass object, or more than ten times brighter than any known accreting pulsar. This finding implies that neutron stars may not be rare in the ULX population, and it challenges physical models for the accretion of matter onto magnetized compact objects.

We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on timescales ~ 10-10^3 r_ g/c, where r_g is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.

Magnetar activity is generated by shear motions of the neutron star surface, which relieve internal magnetic stresses. An analogy with earthquakes and faults is problematic, as the crust is permeated by strong magnetic fields, which greatly constrain crustal displacements. We describe a new deformation mechanism that is specific to strongly magnetized neutron stars. The magnetically stressed crust begins to move because of a thermoplastic instability, which launches a wave that shears the crust and burns its magnetic energy. The propagating wave front resembles the deflagration front in combustion physics. We describe the conditions for the instability, the front structure and velocity, and discuss implications for observed magnetar activity.

We report on the first hard X-ray detection of the Geminga pulsar above 10 keV using a 150 ks observation with the NuSTAR observatory. The double-peaked pulse profile of non-thermal emission seen in the soft X-ray band persists at higher energies. Broadband phase-integrated spectra over the 0.2-20 keV band with NuSTAR and archival XMM-Newton data do not fit to a conventional two-component model of a blackbody plus power-law, but instead exhibit spectral hardening above ~5 keV. We find two spectral models fit the data well: (1) a blackbody (kT1 ~ 42 eV) with a broken power-law (Gamma1 ~ 2.0, Gamma2 ~ 1.4 and Ebreak ~ 3.4 keV), and (2) two blackbody components (kT1 ~ 44 eV and kT2 ~ 195 eV) with a power-law component (Gamma ~ 1.7). In both cases, the extrapolation of the Rayleigh-Jeans tail of the thermal component is consistent with the UV data, while the non-thermal component overpredicts the near-infrared data, requiring a spectral flattening at E ~ 0.01 - 1 keV. While strong phase variation of the power-law index is present below ~5 keV, our phase-resolved spectroscopy with NuSTAR indicates that another hard non-thermal component with Gamma ~ 1.3 emerges above ~5 keV. The spectral hardening in non-thermal X-ray emission as well as spectral flattening between the optical and X-ray bands argue against the conjecture that a single power-law may account for multi-wavelength non-thermal spectra of middle-aged pulsars.

We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double-blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggests that an additional component, such as a power-law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap.

Prompt or early optical emission in gamma-ray bursts is notoriously difficult to measure, and observations of the dozen cases show a large variety of properties. Yet, such early emission promises to help us achieve a better understanding of the GRB emission process(es). We performed dedicated observations of the ultra-long duration (T90 about 7000 s) GRB 130925A in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope. We detect an optical/NIR flare with an amplitude of nearly 2 mag which is delayed with respect to the keV--MeV prompt emission by about 300--400 s. The decay time of this flare is shorter than the duration of the flare (500 s) or its delay. While we cannot offer a straightforward explanation, we discuss the implications of the flare properties and suggest ways toward understanding it.

We report the detection of eight bright X-ray bursts from the 6.5-s magnetar 1E 1048.1-5937, during a 2013 July observation campaign with the Nuclear Spectroscopic Telescope Array (NuSTAR). We study the morphological and spectral properties of these bursts and their evolution with time. The bursts resulted in count rate increases by orders of magnitude, sometimes limited by the detector dead time, and showed blackbody spectra with kT=6-8 keV in the T90 duration of 1-4 s, similar to earlier bursts detected from the source. We find that the spectra during the tail of the bursts can be modeled with an absorbed blackbody with temperature decreasing with flux. The bursts flux decays followed a power-law of index 0.8-0.9. In the burst tail spectra, we detect a ~13 keV emission feature, similar to those reported in previous bursts from this source as well as from other magnetars observed with the Rossi X-ray Timing Explorer (RXTE). We explore possible origins of the spectral feature such as proton cyclotron emission, which implies a magnetic field strength of B~2X10^15 G in the emission region. However, the consistency of the energy of the feature in different objects requires further explanation.

We show that the light curve of the double GeV+optical flash in GRB 130427A is consistent with radiation from the blast wave in a wind-type medium with density parameter $A=\rho r^2\sim 5\times 10^{10}$ g cm$^{-1}$. The peak of the flash is emitted by copious $e^\pm$ pairs created and heated in the blast wave; our first-principle calculation determines the pair-loading factor and temperature of the shocked plasma. Using detailed radiative transfer simulations we reconstruct the observed double flash. The optical flash is dominated by synchrotron emission from the thermal plasma behind the forward shock, and the GeV flash is produced via inverse Compton (IC) scattering by the same plasma. The seed photons for IC scattering are dominated by the prompt MeV radiation during the first tens of seconds, and by the optical to X-ray afterglow thereafter. IC cooling of the thermal plasma behind the forward shock reproduces all GeV data from a few seconds to $\sim 1$ day. We find that the blast wave Lorentz factor at the peak of the flash is $\Gamma\approx 200$, and the forward shock magnetization is $\epsilon_B\sim 2\times 10^{-4}$. An additional source is required by the data in the optical and X-ray bands at times $>10^2$ s; we speculate that this additional source may be a long-lived reverse shock in the explosion ejecta.

We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 for the first ~4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR)} and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.60 +/- 0.07 over our data span. We further show that the change in spin evolution was likely coincident with a bright X-ray burst observed in 2013 June by Swift, and if so, there was no accompanying discontinuity in the frequency. We find that the source 3-10 keV flux has declined monotonically by a factor of ~2 over an 80-day period post-outburst accompanied by a ~20% decrease in the source's blackbody temperature, although there is evidence for both flux and kT having levelled off. We argue that the torque variations are likely to be magnetospheric in nature and will dominate over any dynamical signatures of orbital motion around Sgr A*.

Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated thermal X-ray, optical, and radio emission. The forward shock is radiative initially when the density of shocked gas is highest, at which times radio emission originates from the dense cooling layer immediately downstream of the shock. The radio light curve is characterized by sharper rises to maximum and later peak times at progressively lower frequencies, with a peak brightness temperature that is approximately independent of frequency. We apply our model to the recent gamma-ray classical nova V1324 Sco, obtaining an adequate fit to the early radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a velocity ~1e3 km/s and mass ~ 2e-4 M_sun; the former is consistent with the velocities of narrow line absorption systems observed previously in nova spectra, while the total ejecta mass of the DES and fast outflow is consistent with that inferred independently by modeling the late radio peak. Rapid evolution of the early radio light curves require the DES possess a steep outer density profile, which may indicate that the onset of mass loss from the white dwarf was rapid, providing indirect evidence that the DES was expelled by the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times may contribute significantly to the optical/UV emission, which we speculate is responsible for the previously unexplained `plateaus' and secondary maxima in nova optical light curves.

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available.

The most common progenitors of gamma-ray bursts (GRBs) are massive stars with strong stellar winds. We show that the GRB blast wave in the wind should emit a bright GeV flash. It is produced by inverse Compton cooling of the thermal plasma behind the forward shock wave. The main part of the flash is shaped by scattering of the prompt MeV radiation (emitted at smaller radii) which streams through the external blast wave. The inverse-Compton flash is bright due to the huge e+- enrichment of the external medium. At late times, the blast wave switches to normal synchrotron-self-Compton cooling. The mechanism is demonstrated by a detailed transfer simulation. The observed prompt MeV radiation is taken as an input of the simulation; we use GRB 080916C as an example. The result reproduces the GeV flash observed by the Fermi telescope. It explains the delayed onset, the steep rise, the peak flux, the time of the peak, the long smooth decline, and the spectral slope of GeV emission. The wind density required to reproduce all these features is typical of Wolf-Rayet stars. Our simulation predicts strong TeV emission 1 min after the burst trigger; then a cutoff in the observed high-energy spectrum is expected from absorption by extragalactic background light. In addition, a bright optical counterpart of the GeV flash is predicted for plausible values of the magnetic field; such a double (optical+GeV) flash has been observed in GRB 130427A.

We test a model recently proposed for the persistent hard X-ray emission from magnetars. In the model, hard X-rays are produced by a decelerating electron-positron flow in the closed magnetosphere. The flow decelerates as it radiates its energy away via resonant scattering of soft X-rays, then it reaches the top of the magnetic loop and annihilates there. We test the model against observations of three magnetars: 4U 0142+61, 1RXS J1708-4009, and 1E 1841-045. We find that the model successfully fits the observed phase-resolved spectra. We derive constraints on the angle between the rotational and magnetic axes of the neutron star, the object inclination to the line of sight, and the size of the active twisted region filled with the plasma flow. Using the fit of the hard X-ray component of the magnetar spectrum, we revisit the remaining soft X-ray component. We find that it can be explained by a modified two-temperature blackbody model. The hotter blackbody is consistent with a hot spot covering 1-10% of the neutron star surface. Such a hot spot is expected at the base of the magnetospheric e+- outflow, as some particles created in the e+- discharge flow back and bombard the stellar surface.

The peak time of optical afterglow may be used as a proxy to constrain the Lorentz factor Gamma of the gamma-ray burst (GRB) ejecta. We revisit this method by including bursts with optical observations that started when the afterglow flux was already decaying; these bursts can provide useful lower limits on Gamma. Combining all analyzed bursts in our sample, we find that the previously reported correlation between Gamma and the burst luminosity L_gamma does not hold. However, the data clearly shows a lower bound Gamma_min which increases with L_gamma. We suggest an explanation for this feature: explosions with large jet luminosities and Gamma < Gamma_min suffer strong adiabatic cooling before their radiation is released at the photosphere; they produce weak bursts, barely detectable with present instruments. To test this explanation we examine the effect of adiabatic cooling on the GRB location in the L_gamma - Gamma plane using a Monte Carlo simulation of the GRB population. Our results predict detectable on-axis "orphan" afterglows. We also derive upper limits on the density of the ambient medium that decelerates the explosion ejecta. We find that the density in many cases is smaller than expected for stellar winds from normal Wolf-Rayet progenitors. The burst progenitors may be peculiar massive stars with weaker winds or there might exist a mechanism that reduces the stellar wind a few years before the explosion.

AE Aquarii is a cataclysmic variable with the fastest known rotating magnetized white dwarf (P_spin = 33.08 s). Compared to many intermediate polars, AE Aquarii shows a soft X-ray spectrum with a very low luminosity (L_X ~ 10^{31} erg/s). We have analyzed overlapping observations of this system with the NuSTAR and the Swift X-ray observatories in September of 2012. We find the 0.5-30 keV spectra to be well fitted by either an optically thin thermal plasma model with three temperatures of 0.75 +0.18 -0.45, 2.29 +0.96 -0.82, and 9.33 +6.07 -2.18 keV, or an optically thin thermal plasma model with two temperatures of 1.00 +0.34 -0.23 and 4.64 +1.58 -0.84 keV plus a power-law component with photon index of 2.50 +0.17 -0.23. The pulse profile in the 3-20 keV band is broad and approximately sinusoidal, with a pulsed fraction of 16.6 +/- 2.3%. We do not find any evidence for a previously reported sharp feature in the pulse profile.

Magnetar magnetospheres are believed to be strongly twisted, due to shearing of the stellar crust by internal magnetic stresses. We present time-dependent axisymmetric simulations showing in detail the evolution of relativistic force-free magnetospheres subjected to slow twisting through large angles. When the twist amplitude is small, the magnetosphere moves quasi-statically through a sequence of equilibria of increasing free energy. At some twist amplitude the magnetosphere becomes tearing-mode unstable to forming a resistive current sheet, initiating large-scale magnetic reconnection in which a significant fraction of the magnetic free energy can be dissipated. This "critical" twist angle is insensitive to the resistive length scale. Rapid shearing temporarily stabilizes the magnetosphere beyond the critical angle, allowing the magnetosphere of a rapidly differentially rotating star to store, and dissipate, more free energy. In addition to these effects, shearing the surface of a rotating star increases the spindown torque applied to the star. If shearing is much slower than rotation, the resulting spikes in spindown rate can occur on timescales anywhere from the long twisting timescale to the stellar spin period or shorter, depending both on the stellar shear distribution and the existing distribution of magnetospheric twists. A model in which energy is stored in the magnetosphere, and released by a magnetospheric instability, therefore predicts large changes in the measured spindown rate before SGR giant flares.

We report new spectral and temporal observations of the magnetar 1E 1841-045 in the Kes 73 supernova remnant obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR). Combined with new Swift and archival XMM-Newton and Chandra observations, the phase-averaged spectrum is well characterized by a blackbody plus double power-law model, in agreement with previous, multi-mission X-ray results. However, we are unable to reproduce the spectral results reported using Suzaku observations. The pulsed fraction of the source is found to increase with photon energy. The measured rms pulsed fraction is ~12% and ~17% at ~20 keV and ~50 keV, respectively. We detect a new feature in the 24--35 keV band pulse profile that is uniquely double-peaked. This feature may be associated with a possible absorption or emission feature in the phase-resolved spectrum. We fit the X-ray data using the recently developed electron-positron outflow model of Beloborodov (2013) for the hard X-ray emission from magnetars. This produces a satisfactory fit allowing a constraint on the angle between the rotation and magnetic axes of the neutron star of ~20 degrees and on the angle between the rotation axis and line-of-sight of ~50 degrees. In this model, the soft X-ray component is inconsistent with a single blackbody; adding a second blackbody or a power-law component fits the data. The two-blackbody interpretation suggests a hot spot of temperature kT~0.9 keV occupying ~1% of the stellar surface.

Merger of a white dwarf binary creates a differentially rotating object which is expected to generate strong magnetic fields. Kinetic energy stored in differential rotation is partially dissipated in the magnetically dominated corona, which forms a hot variable outflow with ejection velocity comparable to $10^9$ cm s$^{-1}$. The outflow should carry significant mass and energy for hours to days, creating an expanding fireball with the following features. (i) The fireball is initially opaque and its internal energy is dominated by the trapped thermal radiation. The stored heat is partially converted to kinetic energy of the flow (through adiabatic cooling) and partially radiated away. (ii) Internal shocks develop in the fireball and increase its radiative output. (iii) A significant fraction of the emitted energy is in the optical band. As a result, a bright optical transient with luminosity $L\sim 10^{41}-10^{42}$ erg s$^{-1}$ and a characteristic peak duration comparable to 1 day may be expected from the merger. In contrast to classical novae or supernovae, the transient does not involve nuclear energy. The decay after its peak reflects the damping of differential rotation in the merger remnant. Such outbursts may be detected in the local Universe with current and upcoming optical surveys.

The twisted magnetospheres of magnetars must sustain a persistent flow of electron-positron plasma. The flow dynamics is controlled by the radiation field around the hot neutron star. The problem of plasma motion in the self-consistent radiation field is solved using the method of virtual beams. The plasma and radiation exchange momentum via resonant scattering and self-organize into the "radiatively locked" outflow with a well-defined, decreasing Lorentz factor. There is an extended zone around the magnetar where the plasma flow is ultra-relativistic; its Lorentz factor is self-regulated so that it can marginally scatter thermal photons. The flow becomes slow and opaque in an outer equatorial zone, where the decelerated plasma accumulates and annihilates; this region serves as a reflector for the thermal photons emitted by the neutron star. The e+- flow carries electric current, which is sustained by a moderate induced electric field. The electric field maintains a separation between the electron and positron velocities, against the will of the radiation field. The two-stream instability is then inevitable, and the induced turbulence can generate low-frequency emission. In particular, radio emission may escape around the magnetic dipole axis of the star. Most of the flow energy is converted to hard X-ray emission, which is examined in the accompanying paper.

Magnetic spin-down of a millisecond neutron star has been proposed as the power source of hydrogen-poor "superluminous" supernovae (SLSNe-I). However, producing an unambiguous test that can distinguish this model from alternatives, such as circumstellar interaction, has proven challenging. After the supernova explosion, the pulsar wind inflates a hot cavity behind the expanding stellar ejecta: the nascent millisecond pulsar wind nebula. Electron/positron pairs injected by the wind cool through inverse Compton scattering and synchrotron emission, producing a pair cascade and hard X-ray spectrum inside the nebula. These X-rays ionize the inner exposed side of the ejecta, driving an ionization front that propagates outwards with time. Under some conditions this front can breach the ejecta surface within months after the optical supernova peak, allowing ~0.1-1 keV photons to escape the nebula unattenuated with a characteristic luminosity L_X ~ 1e43-1e45 erg/s. This "ionization break-out" may explain the luminous X-ray emission observed from the transient SCP 06F, providing direct evidence that this SLSN was indeed engine-powered. Luminous break-out requires a low ejecta mass and that the spin-down time of the pulsar be comparable to the photon diffusion timescale at optical maximum, the latter condition similar to that required for a supernova with a high optical fluence. These relatively special requirements may explain why most SLSNe-I are not accompanied by detectable X-ray emission. Global asymmetry of the supernova ejecta increases the likelihood of an early break-out along the direction of lowest density. Atomic states with lower threshold energies are more readily ionized at earlier times near optical maximum, allowing UV break-out across a wider range of pulsar and ejecta properties than X-ray break-out, possibly contributing to the blue/UV colors of SLSNe-I.

Jet heating via nuclear collisions may be the main mechanism for gamma-ray burst (GRB) emission. Besides producing the observed gamma-rays, collisional heating must generate 10-100 GeV neutrinos, implying a close relation between the neutrino and gamma-ray luminosities. We exploit this theoretical relation to make predictions for possible GRB detections by IceCube+DeepCore. To estimate the expected neutrino signal, we use the largest sample of bursts observed by BATSE in 1991-2000. A GRB neutrino could have been detected if IceCube+DeepCore operated at that time. Detection of 10-100 GeV neutrinos would have significant implications, shedding light on the composition of GRB jets and their Lorentz factors. This could be an important target in designing future upgrades of the IceCube+DeepCore observatory.

We report the discovery of 3.76-s pulsations from a new burst source near Sgr A* observed by the NuSTAR Observatory. The strong signal from SGR J1745-29 presents a complex pulse profile modulated with pulsed fraction 27+/-3 % in the 3-10 keV band. Two observations spaced 9 days apart yield a spin-down rate of Pdot = (6.5+/-1.4)x10^{-12}. This implies a magnetic field B = 1.6x10^14 G, spin-down power Edot = 5x10^33 erg/s, and characteristic age P/2Pdot = 9x10^3 yr, for the rotating dipole model. However, the current Pdot may be erratic, especially during outburst. The flux and modulation remained steady during the observations and the 3-79 keV spectrum is well fitted by a combined blackbody plus power-law model with temperature kT_BB = 0.96+/-0.02 keV and photon index = 1.5+/-0.4, respectively. The neutral hydrogen column density (nH ~ 1.4x10^23 cm^{-2}) measured by NuSTAR and Swift suggests that SGR J1745-29 is located at or near the Galactic Center. The lack of an X-ray counterpart in the published Chandra survey catalog sets a quiescent 2-8 keV luminosity limit of Lx ~< 10^32 erg/s. The bursting, timing, and spectral properties indicate a transient magnetar undergoing an outburst with 2-79 keV luminosity up to 3.5x10^35 erg/s for a distance of 8 kpc. SGR J1745-29 joins a growing subclass of transient magnetars, indicating that many magnetars in quiescence remain undetected in the X-ray band or have been detected as high-B radio pulsars. The peculiar location of SGR J1745-29 has important implications for the formation and dynamics of neutron stars in the Galactic Center region.

We study plasma flows above pulsar polar caps using time-dependent simulations of plasma particles in the self-consistent electric field. The flow behavior is controlled by the dimensionless parameter alpha=(j/c rho_GJ) where j is the electric current density and rho_GJ is the Goldreich-Julian charge density. The region of the polar cap where 0<alpha<1 is a "dead zone" --- in this zone particle acceleration is inefficient and pair creation is not expected even for young, rapidly rotating pulsars. Pulsars with polar caps near the rotation axis are predicted to have a hollow-cone structure of radio emission, as the dead zone occupies the central part of the polar cap. Our results apply to charge-separated flows of electrons (j<0) or ions (j>0). In the latter case, we consider the possibility of a mixed flow consisting of different ion species, and observe the development of two-stream instability. The dead zone at the polar cap is essential for the development of an outer gap near the null surface rho_GJ=0.

Observations indicate that the peak of gamma-ray burst spectrum forms in the opaque region of an ultra-relativistic jet. Recent radiative transfer calculations support this picture and show that the spectral peak is inherited from initially thermal radiation, which is changed by heating into a broad photon distribution with a high-energy tail. We discuss the processes that regulate the observed position of the spectral peak E_pk. The opaque jet has three radial zones: (1) Planck zone r<R_P where a blackbody spectrum is enforced; this zone ends where Thomson optical depth decreases to tau~10^5. (2) Wien zone R_P<r<R_W with Kompaneets parameter y>>1 where radiation has a Bose-Einstein spectrum, and (3) Comptonization zone r>R_W where the radiation spectrum develops the high-energy tail. Besides the initial jet temperature, an important factor regulating E_pk is internal dissipation (of bulk motions and magnetic energy) at large distances from the central engine. Dissipation in the Planck zone reduces E_pk, and dissipation in the Wien zone increases E_pk. In jets with sub-dominant magnetic fields, the predicted E_pk varies around 1 MeV up to a maximum value of about 10 MeV. If the jet carries an energetically important magnetic field, E_pk can be additionally increased by dissipation of magnetic energy. This increase is hinted by observations, which show E_pk up to about 20 MeV. We also consider magnetically dominated jets; then a simple model of magnetic dissipation gives E_pk~30 Gamma_W keV where Gamma_W is the jet Lorentz factor at the Wien radius R_W.

We present the first simulations of evolving, strongly twisted magnetar magnetospheres. Slow shearing of the magnetar crust is seen to lead to a series of magnetospheric expansion and reconnection events, corresponding to X-ray flares and bursts. The axisymmetric simulations include rotation of the neutron star and the magnetic wind through the light cylinder. We study how the increasing twist affects the spindown rate of the star, finding that a dramatic increase in spindown occurs. Particularly spectacular are explosive events caused by the sudden opening of large amounts of overtwisted magnetic flux, which may be associated with the observed giant flares. These events are accompanied by a short period of ultra-strong spindown, resulting in an abrupt increase in spin period, such as was observed in the giant flare of SGR 1900+14.

Persistent activity of magnetars is associated with electric discharge that continually injects relativistic particles into the magnetosphere. Large active magnetic loops around magnetars must be filled with outflowing particles that interact with radiation via resonant scattering and spawn electron-positron pairs. The outflow energy is processed into copious e+- until the plasma enters outer parts of the loop where the magnetic field is reduced below 10^13 G. In the outer zone, photons scattered by the outflow do not convert to e+- pairs and the outflow radiates its energy away. The escaping radiation forms a distinct hard X-ray peak in the magnetar spectrum. It has the following features: (1) Its luminosity L=10^35-10^36 erg/s can easily exceed the thermal luminosity from the magnetar surface. (2) Its spectrum extends from 10 keV to the MeV band with a hard spectral slope, which depends on the object inclination to the line of sight. (3) The anisotropic hard X-ray emission exhibits strong pulsations as the magnetar spins. (4) The emission spectrum typically peaks around 1 MeV, but the peak position significantly oscillates with the spin period. (5) The emission is dominated by the extraordinary polarization mode at photon energies below 1 MeV. (6) The decelerated pairs accumulate and annihilate at the top of the magnetic loop, and emit the 511-keV line with luminosity L_ann\sim0.1L. Features (1)-(3) agree with available data; (4)-(6) can be tested by future observations.

We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial derivatives, allowing the use of an unstaggered mesh and the complete force-free current density. The method has low numerical dissipation and diffusion outside of singular current sheets. We present a range of one- and two-dimensional tests, and demonstrate convergence to both smooth and discontinuous analytic solutions. As a first application, we revisit the aligned rotator problem, obtaining a steady solution with resistivity localised in the equatorial current sheet outside the light cylinder.

Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim that collisional heating generates the observed Band-type spectrum of GRBs. We extend the model to study the effects of magnetic fields on the emitted spectrum. We find that the spectrum peak remains near 1 MeV for the entire range of the magnetization parameter $0<\epsilon_B<2$ that is explored in our simulations. The low-energy part of the spectrum softens with increasing $\epsilon_B$, and a visible soft excess appears in the keV band. The high-energy part of the spectrum extends well above the GeV range and can contribute to the prompt emission observed by Fermi/LAT. Overall, the radiation spectrum created by the collisional mechanism appears to agree with observations, with no fine-tuning of parameters.

Analytical and numerical solutions are obtained for the equation of radiative transfer in ultra-relativistic opaque jets. The solution describes the initial trapping of radiation, its adiabatic cooling, and the transition to transparency. Two opposite regimes are examined: (1) Matter-dominated outflow. Surprisingly, radiation develops enormous anisotropy in the fluid frame before decoupling from the fluid. The radiation is strongly polarized. (2) Radiation-dominated outflow. The transfer occurs as if radiation propagated in vacuum, preserving the angular distribution and the blackbody shape of the spectrum. The escaping radiation has a blackbody spectrum if (and only if) the outflow energy is dominated by radiation up to the photospheric radius.

The curvature of a relativistic blast wave implies that its emission arrives to observers with a spread in time. This effect is believed to wash out fast variability in the lightcurves of GRB afterglows. We note that the spreading effect is reduced if emission is anisotropic in the rest-frame of the blast wave (i.e. if emission is limb-brightened or limb-darkened). In particular, synchrotron emission is almost certainly anisotropic, and may be strongly anisotropic, depending on details of electron acceleration in the blast wave. Anisotropic afterglows can display fast and strong variability at high frequencies (above the 'fast-cooling' frequency). This may explain the existence of bizarre features in the X-ray afterglows of GRBs, such as sudden drops and flares. We also note that a moderate anisotropy can significantly delay the 'jet break' in the lightcurve, which makes it harder to detect.

Hyper-accretion discs around black holes emit copious neutrinos and anti-neutrinos. A fraction of the emitted neutrinos convert to electron-positron plasma above the disc through the annihilation reaction $\nu\bar\nu\to e^+e^-$. This process may drive relativistic jets associated with GRB explosions. We calculate the efficiency of energy deposition by neutrinos. Our calculation is fully relativistic and based on a geodesic-tracing method. We find that the efficiency of neutrino heating is a well-defined function of (i) accretion rate and (ii) spin of the black hole. It is practically independent of the details of neutrino transport in the opaque zone of the disc. The results help identify accretion discs whose neutrino emission can power GRBs.

Like the solar corona, the external magnetic field of magnetars is twisted by surface motions of the star. The twist energy is dissipated over time. We discuss the theory of this activity and its observational status. (1) Theory predicts that the magnetosphere tends to untwist in a peculiar way: a bundle of electric currents (the "j-bundle") is formed with a sharp boundary, which shrinks toward the magnetic dipole axis. Recent observations of shrinking hot spots on magnetars are consistent with this behavior. (2) Continual discharge fills the j-bundle with electron-positron plasma, maintaining a nonthermal corona around the neutron star. The corona outside a few stellar radii strongly interacts with the stellar radiation and forms a "radiatively locked" outflow with a high e+- multiplicity. The locked plasma annihilates near the apexes of the closed magnetic field lines. (3) New radiative-transfer simulations suggest a simple mechanism that shapes the observed X-ray spectrum from 0.1 keV to 1 MeV: part of the thermal X-rays emitted by the neutron star are reflected from the outer corona and then upscattered by the inner relativistic outflow in the j-bundle, producing a beam of hard X-rays.

White dwarfs inspiraling into black holes of mass $\MBH\simgt 10^5M_\odot$ are detectable sources of gravitational waves in the LISA band. In many of these events, the white dwarf begins to lose mass during the main observational phase of the inspiral. The mass loss starts gently and can last for thousands of orbits. The white dwarf matter overflows the Roche lobe through the $L_1$ point at each pericenter passage and the mass loss repeats periodically. The process occurs very close to the black hole and the released gas can accrete, creating a bright source of radiation with luminosity close to the Eddington limit, $L\sim 10^{43}$~erg~s$^{-1}$. This class of inspirals offers a promising scenario for dual detections of gravitational waves and electromagnetic radiation.

Nuclear and Coulomb collisions in GRB jets create a hot electron-positron plasma. This collisional heating starts when the jet is still opaque and extends to the transparent region. The e+- plasma radiates its energy. As a result, a large fraction of the jet energy is converted to escaping radiation with a well-defined spectrum. The process is simulated in detail using the known rates of collisions and accurate calculations of radiative transfer in the expanding jet. The result reproduces the spectra of observed GRBs that typically peak near 1 MeV and extend to much higher energies with a photon index \beta ~ -2.5. This suggests that collisional heating may be the main mechanism for GRB emission.

Magnetospheres of neutron stars are anchored in the rigid crust and can be twisted by sudden crustal motions ("starquakes"). The twisted magnetosphere does not remain static and gradually untwists, dissipating magnetic energy and producing radiation. The equation describing this evolution is derived, and its solutions are presented. Two distinct regions coexist in untwisting magnetospheres: a potential region where curl(B)=0 ("cavity") and a current-carrying bundle of field lines ("j-bundle"). The cavity has a sharp boundary, which expands with time and eventually erases all of the twist. In this process, the electric current of the j-bundle is sucked into the star. Observational appearance of the untwisting process is discussed. A hot spot forms at the footprints of the j-bundle. The spot shrinks with time toward the magnetic dipole axis, and its luminosity and temperature gradually decrease. As the j-bundle shrinks, the amplitude of its twist can grow to the maximum possible value ~ 1. The strong twist near the dipole axis increases the spindown rate of the star and can generate a broad beam of radio emission. The model explains the puzzling behavior of magnetar XTE J1810-197 -- a canonical example of magnetospheric evolution following a starquake. We also discuss implications for other magnetars. The untwisting theory suggests that the nonthermal radiation of magnetars is preferentially generated on a bundle of extended closed field lines near the dipole axis.

Accretion onto black holes in wind-fed binaries and in collapsars forms small rotating discs with peculiar properties. Such ``mini-discs'' accrete on the free-fall time without help of viscosity and nevertheless can have a high radiative efficiency. The inviscid mini-disc model was previously constructed for a non-rotating black hole. We extend the model to the case of a spinning black hole, calculate the structure and radiative efficiency of the disc and find their dependence on the black hole spin. If the angular momenta of the disc and the black hole are anti-aligned, a hydrodynamic analog of Penrose process takes place.

We highlight recent theoretical and observational progress in several areas of neutron star astrophysics, and discuss the prospect for advances in the next decade.

A fraction of neutrino emission from GRB accretion disks annihilates above the disk, creating e+- plasma that can drive GRB explosions. We calculate the efficiency of this annihilation using the recent detailed model of hyper-accretion disks around Kerr black holes. Our calculation is fully relativistic and based on a geodesic-tracing method. We find that the efficiency is a well-defined function of (1) accretion rate and (2) spin of the black hole. It is practically independent of the details of neutrino transport in the opaque zone of the disk. The results help identify the accretion disks whose neutrino emission can power GRBs.

Hyper-accretion disks are short-lived, powerful sources of neutrinos and magnetized jets. Such disks are plausible sources of gamma-ray bursts. This review describes the disk structure, the neutrino conversion to electron-positron plasma around the disk, and the post-burst evolution.

Electric currents j flow along the open magnetic field lines from the polar caps of neutron stars. Activity of a polar cap depends on the ratio \alpha=j/c\rho_GJ, where \rho_GJ is the corotation charge density. The customary assumption \alpha\approx 1 is not supported by recent simulations of pulsar magnetospheres and we study polar caps with arbitrary \alpha. We argue that no significant activity is generated on field lines with 0<\alpha<1. Charges are extracted from the star and flow along such field lines with low energies. By contrast, if \alpha>1 or \alpha<0, a high voltage is generated, leading to unsteady e^+- discharge on a scale-height smaller than the size of the polar cap. The discharge can power observed pulsars. Voltage fluctuations in the discharge imply unsteady twisting of the open flux tube and generation of Alfven waves. These waves are ducted along the tube and converted to electromagnetic waves, providing a new mechanism for pulsar radiation.

The standard model of afterglow production by the forward shock wave is not supported by recent observations. We propose a model in which the forward shock is invisible and afterglow is emitted by a long-lived reverse shock in the burst ejecta. It explains observed optical and X-ray light curves, including the plateau at 10^3-10^4 s with a peculiar chromatic break, and the second break that was previously associated with a beaming angle of the explosion. The plateau forms following a temporary drop of the reverse-shock pressure much below the forward-shock pressure. A simplest formalism that can describe such blast waves is the ``mechanical'' model (Beloborodov, Uhm 2006); we use it in our calculations.

Approximate analytical formulae are derived for the pulse profile produced by small hot spots on a rapidly rotating neutron star. Its Fourier amplitudes and phases are calculated. The proposed formalism takes into account gravitational bending of light, Doppler effect, anisotropy of emission, and time delays. Its accuracy is checked with exact numerical calculations.

Relativistic blast waves can be described by a mechanical model. In this model, the "blast" -- the compressed gas between the forward and reverse shocks -- is viewed as one hot body. Equations governing its dynamics are derived from conservation of mass, energy, and momentum. Simple analytical solutions are obtained in the two limiting cases of ultra-relativistic and non-relativistic reverse shock. Equations are derived for the general explosion problem.

Persistent high-energy emission of magnetars is produced by a plasma corona around the neutron star, with total energy output of ~10^{36}erg/s. The corona forms as a result of occasional starquakes that twist the external magnetic field of the star and induce electric currents in the closed magnetosphere. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. The voltage along the magnetic field lines is maintained near threshold for ignition of pair production, in the regime of self-organized criticality. The voltage is found to be about 1 GeV which is in agreement with the observed dissipation rate \~10^{36}erg/s. The coronal particles impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona also emits curvature radiation up to 10^{14}Hz and can supply the observed IR-optical luminosity.

We develop a theoretical model that explains the formation of hot coronae around strongly magnetized neutron stars -- magnetars. The starquakes of a magnetar shear its external magnetic field, which becomes non-potential and is threaded by an electric current. Once twisted, the magnetosphere cannot untwist immediately because of its self-induction. The self-induction electric field lifts particles from the stellar surface, accelerates them, and initiates avalanches of pair creation in the magnetosphere. The created plasma corona maintains the electric current demanded by curl(B) and regulates the self-induction e.m.f. by screening. This corona persists in dynamic equilibrium: it is continually lost to the stellar surface on the light-crossing time of 10^{-4} s and replenished with new particles. In essence, the twisted magnetosphere acts as an accelerator that converts the toroidal field energy to particle kinetic energy. Using a direct numerical experiment, we show that the corona self-organizes quickly (on a millisecond timescale) into a quasi-steady state, with voltage ~1 GeV along the magnetic lines. The heating rate of the corona is ~10^{36} erg/s, in agreement with the observed persistent, high-energy output of magnetars. We deduce that a static twist that is suddenly implanted into the magnetosphere will decay on a timescale of 1-10 yrs. The particles accelerated in the corona impact the solid crust, knock out protons, and regulate the column density of the hydrostatic atmosphere of the star. The transition layer between the atmosphere and the corona is the likely source of the observed 100-keV emission from magnetars. The corona emits curvature radiation and can supply the observed IR-optical luminosity. (Abridged)

We calculate the structure of accretion disks around Kerr black holes for accretion rates 0.001 - 10 M_sun/s. Such disks are plausible candidates for the central engine of gamma-ray bursts. Our disk model is fully relativistic and treats accurately microphysics of the accreting matter: neutrino emissivity, opacity, electron degeneracy, and nuclear composition. The neutrino-cooled disk forms above a critical accretion rate that depends on the black hole spin. The disk has the ``ignition'' radius r_ign where neutrino flux rises dramatically, cooling becomes efficient, and the proton-to-nucleon ratio Y_e drops. Other characteristic radii are r_alpha where most of alpha-particles are disintegrated, r_nu where the disk becomes neutrino-opaque, and r_trap where neutrinos get trapped and advected into the black hole. We find r_alpha, r_ign, r_nu, r_trap and show their dependence on the accretion rate. We discuss the qualitative picture of accretion and present sample numerical models of the disk structure. All neutrino-cooled disks regulate themselves to a characteristic state such that: (1) electrons are mildly degenerate, (2) Y_e ~ 0.1, and (3) neutrons dominate the pressure in the disk.

Two disks of young stars have recently been discovered in the Galactic Center. The disks are rotating in the gravitational field of the central black hole at radii r=0.1-0.3 pc and thus open a new opportunity to measure the central mass. We find that the observed motion of stars in the clockwise disk implies M=4.3+/-0.5 million solar masses for the fiducial distance to the Galactic Center R_0=8 kpc and derive the scaling of M with R_0. As a tool for our estimate we use orbital roulette, a recently developed method. The method reconstructs the three-dimensional orbits of the disk stars and checks the randomness of their orbital phases. We also estimate the three-dimensional positions and orbital eccentricities of the clockwise-disk stars.

We report the definite spectroscopic identification of 41 OB supergiants, giants and main sequence stars in the central parsec of the Galaxy. Detection of their absorption lines have become possible with the high spatial and spectral resolution and sensitivity of the adaptive optics integral field spectrometer SPIFFI/SINFONI on the ESO VLT. Several of these OB stars appear to be helium and nitrogen rich. Almost all of the ~80 massive stars now known in the central parsec (central arcsecond excluded) reside in one of two somewhat thick (<|h|/R>~0.14) rotating disks. These stellar disks have fairly sharp inner edges (R~1") and surface density profiles that scale as R^{-2}. We do not detect any OB stars outside the central 0.5 pc. The majority of the stars in the clockwise system appear to be on almost circular orbits, whereas most of those in the `counter-clockwise' disk appear to be on eccentric orbits. Based on its stellar surface density distribution and dynamics we propose that IRS 13E is an extremely dense cluster (core density > 3x10^8 sunmass/pc^3), which has formed in the counter-clockwise disk. The stellar contents of both systems are remarkably similar, indicating a common age of ~6+/-2 Myr. The K-band luminosity function of the massive stars suggests a top-heavy mass function and limits the total stellar mass contained in both disks to ~1.5x10^4 sunmass. Our data strongly favor in situ star formation from dense gas accretion disks for the two stellar disks. This conclusion is very clear for the clockwise disk and highly plausible for the counter-clockwise system.

We propose a simple mechanical model for relativistic explosions with both forward and reverse shocks, which allows one to do fast calculations of GRB afterglow. The blast wave in the model is governed by pressures P_F and P_R at the forward and reverse shocks. We show that the simplest assumption P_F=P_R is in general inconsistent with energy conservation law. The model is applied to GRBs with non-uniform ejecta. Such "stratified fireballs" are likely to emerge with a monotonic velocity profile after an internal-shock stage. We calculate the early afterglow emission expected from stratified fireballs.

We calculate the structure of accretion disk around a spinning black hole for accretion rates 0.01 - 10 M_sun/s. The model is fully relativistic and treats accurately the disk microphysics including neutrino emissivity, opacity, electron degeneracy, and nuclear composition. We find that the accretion flow always regulates itself to a mildly degenerate state with the proton-to-nucleon ratio Y_e ~ 0.1 and becomes very neutron-rich. The disk has a well defined "ignition" radius where neutrino flux raises dramatically, cooling becomes efficient, and Y_e suddenly drops. We also calculate other characteristic radii of the disk, including the neutrino-opaque and neutrino-trapping radii, and show their dependence on the accretion rate. Accretion disks around fast-rotating black holes produce intense neutrino fluxes which may deposit enough energy above the disk to generate a GRB jet.

Relativistic neutron-loaded outflows in gamma-ray bursts are studied at their early stages, before deceleration by a surrounding medium. The outflow has four components: radiation, electrons, protons and neutrons. The components interact with each other and exchange energy as the outflow expands. The presence of neutrons significantly changes the outflow evolution. Before neutrons decouple from protons, friction between the two components increases their temperatures by many orders of magnitude. After the decoupling, the gradual neutron decay inside the outflow has a drag effect on the protons and reduces their final Lorentz factor.

The recently discovered soft gamma-ray emission from the anomalous X-ray pulsar 1E 1841-045 has a luminosity L_g ~ 10^{36} ergs/s. This luminosity exceeds the spindown power by three orders of magnitude and must be fed by an alternative source of energy such as an ultrastrong magnetic field. A gradual release of energy in the stellar magnetosphere is expected if it is twisted and a strong electric current is induced on the closed field lines. We examine two mechanisms of gamma-ray emission associated with the gradual dissipation of this current. (1) A thin surface layer of the star is heated by the downward beam of current-carrying charges, which excite Langmuir turbulence in the layer. As a result, it can reach a temperature kT ~ 100 keV and emit bremsstrahlung photons up to this characteristic energy. (2) The magnetosphere is also a source of soft gamma rays at a distance of ~100 km from the star, where the electron cyclotron energy is in the keV range. A large electric field develops in this region in response to the outward drag force felt by the current-carrying electrons from the flux of keV photons leaving the star. A seed positron injected in this region undergoes a runaway acceleration and upscatters keV photons above the threshold for pair creation. The created pairs emit a synchrotron spectrum consistent with the observed 20-100 keV emission. This spectrum is predicted to extend to higher energies and reach a peak at ~1 MeV.

The MeV radiation front of gamma-ray bursts creates copious e+- pairs as it propagates through an ambient medium. The created pairs enrich the leptonic component of the medium by a large factor at distances R < R_load ~ 10^{16} cm from the burst center. The following blast wave sweeps up the pair-rich medium and then emits the observed afterglow radiation. We find that the afterglow has a "memory" of e+- loading outside R_load. The e+- pairs remain in the swept-up material and slowly cool down by emitting synchrotron radiation. They are likely to dominate the blast-wave emission in IR, optical, and UV bands during the first minutes of the observed afterglow. The expected e+- radiation is described by a simple formula, which is derived analytically and checked by numerical integration of synchrotron emission over the blast material; a suitable Lagrangian formalism is developed for such calculations. The main signature of e+- radiation is its flat ("white") spectrum in a broad range of frequencies from IR to UV and possibly soft X-rays. This radiation can be detected by Swift satellite, which would enable new observational tests for the explosion physics.

The synchrotron optical flash caught in GRB 990123 overlaps with the MeV radiation front, and the optical-emitting electrons must also produce GeV-TeV emission by inverse Compton scattering of MeV photons. The ultra-high-energy flash can be much stronger than its optical counterpart. We also note that Compton cooling by MeV photons immediately terminates the optical emission unless the fireball Lorentz factor exceeds 10^3. Severe Compton losses may explain the non-detections of optical flashes in several long GRBs. Such failed optical flashes should be especially efficient GeV producers and likely to develop e+- cascades. This probably happened in GRB 941017, and its mysterious high-energy component is well explained by Compton upscattering of GRB photons at the fireball deceleration radius. The proposed mechanism of GeV emission should not work for short gamma-ray bursts that early decouple from the fireball and avoid interaction with the electrons in the deceleration flash. Observations by Swift and GLAST will provide an opportunity to test these expectations. The existing data for GRB 990123 already impose interesting constraints on the explosion.

The traditional way of estimating the gravitational field from observed motions of test objects is based on the virial relation between their kinetic and potential energy. We find a more efficient method. It is based on the natural presumption that the objects are observed at a random moment of time and therefore have random orbital time phases. The proposed estimator, which we call "orbital roulette", checks the randomness of the phases. The method has the following advantages: (1) It estimates accurately Keplerian (point-mass) potentials as well as non-Keplerian potentials where the unknown gravitating mass is distributed in space. (2) It is a complete statistical estimator: it checks a trial potential and accepts it or rules it out with a certain significance level; the best-fit measurement is thus supplemented with error bars at any confidence level. (3) It needs no a priori assumptions about the distribution of orbital parameters of the test bodies. We test our estimator with Monte-Carlo-generated motions and demonstrate its efficiency. Useful applications include the Galactic Center, dark-matter halo of the Galaxy, and clusters of stars or galaxies.

Fireballs with high energy per baryon rest mass (> 400) contain a relatively slow neutron component. We show here that in this situation the thermal history of fireballs is very different from the standard adiabatic cooling.

Physics of GRB blast waves is discussed with a focus on two effects: (1) pair creation in the external medium by the gamma-ray front and (2) decay of neutrons ahead of the decelerating blast wave. Both effects impact the afterglow mechanism at radii up to 10^{17}cm.

Cosmological GRBs are discussed with an emphasis on their plausible connection with black holes. GRBs can be triggered by collapse of stellar-mass objects that leads to formation of a black hole and a transient debris disk with a huge accretion rate. The disk is believed to produce a relativistic jet (``fireball'') that expands and emits to infinity the observed burst of gamma-rays. This accretion-jet picture is similar to quasars and X-ray binaries, however, there are important differences: the physical conditions and the cooling mechanism in the disk are very different. The observed radiation is emitted when the expanding fireball becomes transparent, at distances much larger than the Schwarzschild radius. The burst is then observed as a powerful relativistic explosion and the transient accretion disk in its center serves as a brief source of energy that drives the explosion. The explosion picture depends on the fireball nuclear composition which is shaped close to the black hole. A large amount of free neutrons survive till the emission phase and link the physics of the central engine to observed radiation.

The customary picture of the fireball deceleration by an external medium neglects two physical agents in the problem that are important at radii R<10^{17} cm: a radiative precursor (the prompt GRB) and a relativistic front of free neutrons. The radiative precursor preaccelerates the medium and loads it with e+- pairs. This provides a new explanation for the early optical flashes in GRBs. The front of free neutrons must form in a standard baryonic fireball and change the mechanism of the fireball deceleration. The neutron effect, however, disappears if the fireball is strongly dominated by the Poynting flux. The neutrons thus provide a unique link between the progenitor physics and the observed external blast wave. The explosion mechanism at R<10^{17} cm is testable with upcoming observations of early afterglows.

Observations of the galactic center revealed a population of young massive stars within 0.4 pc from Sgr A* -- the presumed location of a supermassive black hole. The origin of these stars is a puzzle as their formation in citu should be suppressed by the black hole's tidal field. We find that out of 13 stars whose 3-dimensional velocities have been measured by Genzel et. al. (2000), 10 lie in a thin disk. The half-opening angle of the disk is consistent with zero within the measurement errors, and does not exceed 10 degrees. We propose that a recent burst of star formation has occurred in a dense gaseous disk around Sgr A*. Such a disk is no longer present because, most likely, it has been accreted by the central black hole. The three-dimensional orbit of S2, the young star closest to Sgr A*, has been recently mapped out with high precision. It is inclined to the stellar disk by 75 degrees. We find that the orbit should undergo Lense-Thirring precession with the period of (5/a) Myr, where a<1 is the dimensionless spin of the black hole. Therefore it is possible that originally S2 orbit lay in the disk plane. If so, we can constrain the black hole spin $a$ be greater than 0.2(t_{S2}/5 Myr), where t_{S2} is the age of S2.

Standard GRB fireballs must carry free neutrons. This crucially changes the mechanism of fireball deceleration by an external medium. As the ion fireball decelerates, the coasting neutrons form a leading front. They gradually decay, leaving behind a relativistic trail of decay products mixed with the ambient medium. The ion fireball sweeps up the trail and drives a shock wave in it. Thus, observed afterglow emission is produced in the neutron trail. The impact of neutrons turns off at 10^{17}cm from the explosion center, and here a spectral transition is expected in GRB afterglows. Absence of neutron signatures would point to absence of baryons and a dominant Poynting flux in the fireballs.

Fireballs of gamma-ray bursts are partially made of free neutrons. Their presence crucially changes the mechanism of the fireball deceleration by an external medium. At late stages of the explosion, neutrons fully decouple from ions and coast with a constant Lorentz factor \Gamma_n. As the ion fireball decelerates, the neutrons form a leading front. This front gradually decays, leaving behind a trail of decay products mixed with the ambient medium. The kinetic energy of the decay products far exceeds the medium rest-mass energy, and the trail has a high Lorentz factor at radii R < R_{trail} = 10R_\beta \approx 10^{17} cm where R_\beta = 10^{16}(\Gamma_n/300) cm is the mean radius of neutron decay. The ion fireball sweeps up the trail material and drives a shock wave in it. Thus, observed afterglow emission is produced in the neutron trail. It can naturally re-brighten or display a spectral transition at R \approx R_{trail} where the impact of neutrons turns off. Absence of any neutron signatures would point to an extremely low baryon loading of the fireballs and a strong dominance of a Poynting flux.

We study three processes that shape the nuclear composition of GRB fireballs: (1) neutronization in the central engine, (2) nucleosynthesis in the fireball as it expands and cools, and (3) spallation of nuclei in subsequent internal shocks. The fireballs are found to have a neutron excess and a marginally successful nucleosynthesis. They are composed of free nucleons, alpha-particles, and deuterium. A robust result is the survival of a significant neutron component, which has important implications. First, as shown in previous works, neutrons can lead to observable multi-GeV neutrino emission. Second, as we show in an accompanying paper, neutrons impact the explosion dynamics at radii up to 10^{17} cm and change the mechanism of the GRB afterglow emission.

The standard afterglow model neglects the presence of the GRB prompt radiation ahead of the blast wave. In fact, the blast wave is dramatically influenced by the leading gamma-ray front which preaccelerates the ambient medium and loads it with electron-positron pairs. The front sweeps the medium outward with a high Lorentz factor and results in a spectacular effect: the GRB ejecta moves freely in a cavity behind its own radiation front. When the front expands sufficiently and gets diluted, a blast wave develops, and it does it differently from the standard model used before. The afterglow should initially appear as a steep rise of soft emission (from infrared to soft X-rays) at a time comparable to the prompt GRB duration and then the emission should quickly evolve to a normal X-ray afterglow. This may explain the prompt optical flash observed in GRB 990123 and allows one to infer the ejecta Lorentz factor in this burst: \Gamma=200. The effect of the gamma-ray front on the afterglow emission is especially pronounced if the GRB has a massive progenitor and the blast wave propagates in the progenitor wind. We emphasize the importance of early afterglow observations in soft bands, as they will allow one to test different progenitor models.

Three classical modes of accretion are briefly discussed: wind-fed, spherical, and disk. The three modes are illustrated with the mass transfer onto black holes in high-mass X-ray binaries. Then a new regime of mini-disk accretion is described and it is argued that observed wind-fed X-ray sources are likely to accrete in this regime. Switching from one accretion mode to another can cause the observed spectral state transitio ns.

Observations and models of GRBs are reviewed with an emphasis on the early phase of the burst. Recent revision of the early afterglow model is described.

A photon emitted near a compact object at an angle $\alpha$ with respect to the radial direction escapes to infinity at a different angle $\psi>\alpha$. This bending of light is caused by a strong gravitational field. We show that, in a Schwarzschild metric, the effect is described by $1-\cos\alpha=(1-\cos\psi)(1-r_g/R)$ where $R/r_g$ is the emission radius in Schwarzschild units. The formula is approximate and it applies at $R\geq 2r_g$ only, however at these radii it has amazing accuracy, fully sufficient in many applications. As one application we develop a new formulation for the light bending effects in pulsars. It reveals the simple character of these effects and gives their quantitative description with practically no losses of accuracy (for the typical radius of a neutron star $R=3r_g$ the error is 1%). The visible fraction of a star surface is shown to be $S_v/4\pi R^2=[2(1-r_g/R)]^{-1}$ which is 3/4 for $R=3r_g$. The instantaneous flux of a pulsar comes from one or two antipodal polar caps that rotate in the visible zone. The pulse produced by one blackbody cap is found to be sinusoidal (light bending impacts the pulse amplitude but not its shape). When both caps are visible, the pulse shows a plateau: the variable parts of the antipodal emissions precisely cancel each other. The pulsed fraction of blackbody emission with antipodal symmetry has an upper limit $A_{max}=(R-2r_g)/(R+2r_g)$. Pulsars with $A>A_{max}$ must be asymmetric.

Gamma-ray bursts (GRBs) are emitted by relativistic ejecta from powerful cosmic explosions. Their light curves suggest that the gamma-ray emission occurs at early stages of the ejecta expansion, well before it decelerates in the ambient medium. If so, the launched gamma-ray front must overtake the ejecta and sweep the ambient medium outward. As a result a gap is opened between the ejecta and the medium that surfs the radiation front ahead. Effectively, the ejecta moves in a cavity until it reaches a radius R_{gap}=10^{16}E_{54}^{1/2} cm where E is the isotropic energy of the GRB. At R=R_{gap} the gap is closed, a blast wave forms and collects the medium swept by radiation. Further development of the blast wave is strongly affected by the leading radiation front: the front plays the role of a precursor where the medium is loaded with e+- pairs and preaccelerated just ahead of the blast. It impacts the emission from the blast at R < R_{load}=5R_{gap} (the early afterglow). A spectacular observational effect results: GRB afterglows should start in optical/UV and evolve fast (< min) to a normal X-ray afterglow. The early optical emission observed in GRB 990123 may be explained in this way. The impact of the front is especially strong if the ambient medium is a wind from a massive progenitor of the GRB. In this case three phenomena are predicted: (1) The ejecta decelerates at R<R_{load} producing a lot of soft radiation. (2) The light curve of soft emission peaks at t_{peak}=40(1+z)E_{54}^{1/2}(Gamma_{ej}/100)^{-2} s where Gamma_{ej} is the Lorentz factor of the ejecta. Given measured redshift z and t_{peak}, one finds Gamma_{ej}. (3) The GRB acquires a spectral break at 5 - 50 MeV because harder photons are absorbed by radiation scattered in the wind.

Models of X-ray production in black-hole sources are reviewed and compared with recent observations. Possible diagnostics of the hot-disk and disk-corona models are discussed in the light of spectral and temporal data. A new model of a small-scale inviscid accretion disk is presented.

A very high activity of Cygnus X-1 on 1999 April 19-21 was recorded by BATSE Large Area Detectors onboard the Compton Gamma-Ray Observatory. The peak luminosity was one order of magnitude higher than the normal luminosity of Cyg X-1. This fact can be critical for models of the hard state of Cyg X-1. The longest outburst lasted ~1000 s and demonstrated very unusual temporal and spectral behavior which indicates the presence of two emission components. One component is relatively soft (with a cutoff below 100 keV) and highly variable, and the other one is hard (extending above 100 keV), with much slower variability.

We compute the X-ray spectra produced by non-static coronae atop accretion discs around black holes and neutron stars. The hot corona is radiatively coupled to the underlying disc (the reflector) and generates an X-ray spectrum which is sensitive to the bulk velocity of the coronal plasma, beta=v/c. We show that an outflowing corona reproduces the hard-state spectrum of Cyg X-1 and similar objects. The dynamic model predicts a correlation between the observed amplitude of reflection R and the X-ray spectrum slope Gamma since both strongly depend on beta. A similar correlation was observed and its shape is well fitted by the dynamic model. The scattering of soft radiation in an outflowing corona can also account for the observed optical-UV polarization pattern in active galactic nuclei.

Gas falling quasi-spherically onto a black hole forms an inner accretion disc if its specific angular momentum $l$ exceeds $\lmin\sim r_gc$ where $r_g$ is the Schwarzschild radius. The standard disc model assumes $l\gg\lmin$. We argue that, in many black-hole sources, accretion flows have angular momenta just above the threshold for disc formation, $l\simgt\lmin$, and assess the accretion mechanism in this regime. In a range $\lmin<l<\lcr$, a small-scale disc forms in which gas spirals fast into the black hole without any help of horizontal viscous stresses. Such an `inviscid' disc, however, interacts inelastically with the feeding infall. The disc-infall interaction determines the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as $14 r_g$, and the energy release peaks at $2r_g$. The disc emits a Comptonized X-ray spectrum with a break at $\sim 100$ keV. This accretion regime is likely to take place in wind-fed X-ray binaries and is also possible in active galactic nuclei.

In wind-fed X-ray binaries the accreting matter is Compton cooled and falls freely onto the compact object. The matter has a modest angular momentum $l$ and accretion is quasi-spherical at large distances from the compact object. Initially small non-radial velocities grow in the converging supersonic flow and become substantial in the vicinity of the accretor. The streamlines with $l>(GMR_*)^{1/2}$ (where $M$ and $R_*$ are the mass and radius of the compact object) intersect outside $R_*$ and form a two-dimensional caustic which emits X-rays. The streamlines with low angular momentum, $l<(GMR_*)^{1/2}$, run into the accretor. If the accretor is a neutron star, a large X-ray luminosity results. We show that the distribution of accretion rate/luminosity over the star surface is sensitive to the angular momentum distribution of the accreting matter. The apparent luminosity depends on the side from which the star is observed and can change periodically with the orbital phase of the binary. The accretor then appears as a `Moon-like' X-ray source.

The fraction of a fireball kinetic kinetic energy that can be radiated away by internal shocks is sensitive to the amplitude of initial fluctuations in the fireball. We give a simple analytical description for dissipation of modest-amplitude fluctuations and confirm it with direct numerical simulations. At high amplitudes, the dissipation occurs in a non-linear regime with efficiency approaching 100 %. Most of the explosion energy is then radiated away by the prompt GRB and only a fraction remains to be radiated by the afterglow.

Power density spectra (PDSs) of long gamma-ray bursts (GRBs) provide useful information on GRBs, indicating their self-similar temporal structure. The best power-law PDSs are displayed by the longest bursts (T_90>100 s) in which the range of self-similar time scales covers more than 2 decades. Shorter bursts have apparent PDS slopes more strongly affected by statistical fluctuations. The underlying power law can then be reproduced with high accuracy by averaging the PDSs for a large sample of bursts. This power law has a slope approximately equal to -5/3 and a sharp break at about 1 Hz. The power-law PDS provides a new sensitive tool for studies of GRBs. In particular, we calculate the PDSs of bright bursts in separate energy channels. The PDS flattens in the hard channel (h\nu>300 keV) and steepens in the soft channel (h\nu<50 keV), while the PDS of bolometric light curves approximately follows the -5/3 law. We then study dim bursts and compare them to the bright ones. We find a strong correlation between the burst brightness and the PDS slope. This correlation shows that the bursts are far from being standard candles and dim bursts should be intrinsically weak. The time dilation of dim bursts is probably related to physical processes occurring in the burst rather than to a cosmological redshift.

Power density spectra (PDSs) of long gamma-ray bursts provide useful information on GRBs, indicating their self-similar temporal structure. The best power-law PDSs are displayed by the longest bursts (T_90 > 100 s) in which the range of self-similar time scales covers more than 2 decades. Shorter bursts have apparent PDS slopes more strongly affected by statistical fluctuations. The underlying power law can then be reproduced with high accuracy by averaging the PDSs for a large sample of bursts. This power-law has a slope alpha\approx -5/3 and a sharp break at 1 Hz. The power-law PDS provides a new sensitive tool for studies of gamma-ray bursts. In particular, we calculate the PDSs of bright bursts in separate energy channels. The PDS flattens in the hard channel (h\nu > 300 keV) and steepens in the soft channel (h\nu < 50 keV), while the PDS of bolometric light curves approximately follows the -5/3 law. We then study dim bursts and compare them to the bright ones. We find a strong correlation between the burst brightness and the PDS slope. This correlation shows that the bursts are far from being standard candles and dim bursts should be intrinsically weak. The time dilation of dim bursts is probably related to physical processes occurring in the burst rather than to a cosmological redshift. Finally, we test the internal shock model against the observed PDS. We demonstrate how the model can reproduce the -5/3 power law.

An electron-positron atmosphere is inevitably created around a black hole accretion disc, the spectrum of which extends to MeV energies. Pairs created in photon-photon collisions outside the disc are blown away by soft radiation (which dominates the bolometric luminosity of the disc) and form a semirelativistic outflow. We simulate numerically the conversion of the MeV radiation into a vertical e+- outflow above a disc-like source. The outflowing e+- plasma becomes optically thick to Thomson scattering if the compactness of the gamma-ray source exceeds 30. The scattering by e+- then collimates the bulk of soft radiation along the disc axis, and the apparent bolometric luminosity of the disc depends strongly on its inclination to the line of sight. The anisotropic central emission may account for the lack of Fe K-alpha lines in the X-ray spectra of bright radio-quiet quasars. The scattering in e+- outflows may also explain the orientation of optical polarization in non-blazar active galactic nuclei.

We apply Fourier analysis to 214 light curves of long gamma-ray bursts and study the statistical properties of their power density spectra (PDSs). The averaged PDS is found to follow a power law of index -5/3 over almost 2 decades of frequency with a break at about 2 Hz. Individual PDSs are exponentially distributed around the power law. It provides evidence that the diversity of the bursts is due to random realizations of the same process which is self-similar over the full range of time-scales. The -5/3 slope of the average spectrum may indicate that gamma-ray bursts are related to a phenomenon well studied in hydrodynamics - fully developed turbulence.

Optical/UV radiation from accretion disks in quasars is likely to be partly scattered by a hot plasma enveloping the disk. We investigate whether the scattering may produce the steep rises in polarization observed blueward of the Lyman limit in some quasars. We suggest and assess two models. In the first model, primary disk radiation with a Lyman edge in absorption passes through a static ionized "skin" covering the disk, which has a temperature about 3 keV and a Thomson optical depth about unity. Electron scattering in the skin smears out the edge and produces a steep rise in polarization at lambda < 912 A. In the second model, the scattering occurs in a hot coronal plasma outflowing from the disk with a mildly relativistic velocity. We find that the second model better explains the data. The ability of the models to fit the observed rises in polarization is illustrated with the quasar PG 1630+377.

Models of black hole accretion disks are reviewed, with an emphasis on the theory of hard X-ray production. The following models are considered: i) standard, ii) super-critical, iii) two-temperature, and iv) disk+corona. New developments have recently been made in hydrodynamical models of accretion and in phenomenological radiative models fitting the observed X-ray spectra. Attempts to unify the two approaches are only partly successful.

The hard X-rays in Cyg X-1 and similar black hole sources are possibly produced in an active corona atop an accretion disk. We suggest that the observed weakness of X-ray reflection from the disk is due to bulk motion of the emitting hot plasma away from the reflector. A mildly relativistic motion causes aberration reducing X-ray emission towards the disk. This in turn reduces the reprocessed radiation from the disk and leads to a hard spectrum of the X-ray source. The resulting spectral index is Gamma=1.9B^{1/2} where B=gamma(1+beta) is the aberration factor for a bulk velocity beta=v/c. The observed Gamma=1.6 and the amount of reflection, R=0.3, in Cyg X-1 in the hard state can both be explained assuming a bulk velocity beta=0.3. We discuss one possible scenario: the compact magnetic flares are dominated by e+- pairs which are ejected away from the reflector by the pressure of the reflected radiation. We also discuss physical constraints on the disk-corona model and argue that the magnetic flares are related to magneto-rotational instabilities in the accretion disk.

We calculate the structure of accretion disc around a rapidly rotating black hole with a super-Eddington accretion rate. The luminosity and height of the disc are reduced by the advection effect. In the case of a large viscosity parameter, alpha > 0.03, the accretion flow strongly deviates from thermodynamic equilibrium and overheats in the central region. With increasing accretion rate, the flow temperature steeply increases, reaches a maximum, and then falls off. The maximum is achieved in the advection dominated regime of accretion. The maximum temperature in the disc around a massive black hole, M=10^8 M_sun, with alpha=0.3 is of order 3 x 10^8 K. Discs with large accretion rates can emit X-rays in quasars as well as in galactic black hole candidates.

A fraction of the radiation produced by an accretion disk may be Thomson scattered by a wind flowing away from the disk. Employing a simple plane-parallel model of the wind, we calculate the polarization of the scattered radiation and find that its sign depends on the wind velocity, beta=v/c. In the case, 0.12 < beta < 0.78, the polarization is parallel to the disk normal, i.e., it is orthogonal to the standard Chandrasekhar's polarization expected from accretion disks. The velocity of an electron-positron wind is likely to saturate near the equilibrium value beta_* of order 0.5 for which the accelerating radiation pressure is balanced by the Compton drag. Then the change of polarization by the wind is most pronounced. This may help to reconcile the standard accretion disk model with the optical polarimetric observations of non-blazar AGNs.

In the innermost region of the advective accretion disk orbiting a black hole of high spin, the inertia of heat stored in the accreting gas is comparable to that of the gas rest mass itself. Accounting for this effect, we derive additional terms in the disk structure equations, and show that the heat inertia plays a significant role in the global energy conservation and dynamics of accretion in the relativistic advective disks.

A simple re-scaling of velocities calculated from the Paczy{\'n}ski \& Wiita (1980) pseudo-Newtonian potential makes them consistent with special relativity and greatly improves the agreement with the exact relativistic calculations. The improvement is relevant in calculations of the Doppler effect, spectra and radiation transfer which all involve effects of special relativity.