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.

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 reanalyze the prompt emission of two of the brightest Fermi GRBs (080916C and 090926A) with a new model composed of 3 components: (i) a thermal-like component--approximated with a black body (BB)--interpreted as the jet photosphere emission of a magnetized relativistic outflow, (ii) a non-thermal component--approximated with a Band function--interpreted as synchrotron radiation in an optically thin region above the photosphere either from internal shocks or magnetic field dissipation, and (iii) an extra power law (PL) extending from low to high energies likely of inverse Compton origin, even though it remains challenging. Through fine-time spectroscopy down to the 100 ms time scale, we follow the smooth evolution of the various components. From this analysis the Band function is globally the most intense component, although the additional PL can overpower the others in sharp time structures. The Band function and the BB component are the most intense at early times and globally fade across the burst duration. The additional PL is the most intense component at late time and may be correlated with the extended high-energy emission observed thousands of seconds after the burst with the Fermi/Large Area Telescope (LAT). Unexpectedly, this analysis also shows that the additional PL may be present from the very beginning of the burst. We investigate the effect of the three components on the new time-resolved luminosity-hardness relation in both the observer and rest frames and show that a strong correlation exists between the flux of the non-thermal component and its E$_{peak}$ only when the three components are fitted simultaneously to the data (i.e., F$_i^{NT}$-E$_{peak,i}^{NT}$ relation). In addition, this result points toward a universal relation between those two quantities for all GRBs when transposed to the central engine rest frame (i.e., L$_i^{NT}$-E$_{peak,i}^{rest,NT}$ relation).

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.

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.

The early X-ray afterglow of gamma-ray bursts revealed by Swift carried many surprises. We focus in this paper on the plateau phase whose origin remains highly debated. We confront several newly discovered correlations between prompt and afterglow quantities (isotropic emitted energy in gamma-rays, luminosity and duration of the plateau) to several models proposed for the origin of plateaus in order to check if they can account for these observed correlations. We first show that the scenario of plateau formation by energy injection into the forward shock leads to an efficiency crisis for the prompt phase and therefore study two possible alternatives: the first one still takes place within the framework of the standard forward shock model but allows for a variation of the microphysics parameters to reduce the radiative efficiency at early times; in the second scenario the early afterglow results from a long-lived reverse shock. Its shape then depends on the distribution of energy as a function of Lorentz factor in the ejecta. In both cases, we first present simple analytical estimates of the plateau luminosity and duration and then compute detailed light curves. In the two considered scenarios we find that plateaus following the observed correlations can be obtained under the condition that specific additional ingredients are included. In the forward shock scenario, the preferred model supposes a wind external medium and a microphysics parameter epsilon_e that first varies as n^{-\nu} (n being the external density), with \nu~1 to get a flat plateau, before staying constant below a critical density n_0. To produce a plateau in the reverse shock scenario the ejecta must contain a tail of low Lorentz factor with a peak of energy deposition at \Gamma >~ 10.

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.

We study the observed correlations between the duration and luminosity of the early afterglow plateau and the isotropic gamma-ray energy release during the prompt phase. We discuss these correlations in the context of two scenarios for the origin of the plateaus. In the first one the afterglow is made by the forward shock and the plateau results from variations of the microphysics parameters while in the second one the early afterglow is made by a long-lived reverse shock propagating in a low Lorentz factor tail of the ejecta.

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.

We perform a detailed study on the dynamics of a relativistic blast wave with the presence of a long-lived reverse shock (RS). Although a short-lived RS has been widely considered, the RS is believed to be long-lived as a consequence of a stratification expected on the ejecta Lorentz factors. The existence of a long-lived RS makes the forward shock (FS) dynamics to deviate from a self-similar Blandford-McKee solution. Employing the "mechanical model" that correctly incorporates the energy conservation, we present an accurate solution for both the FS and RS dynamics. We conduct a sophisticated calculation of the afterglow emission. Adopting a Lagrangian description of the blast wave, we keep track of an adiabatic evolution of numerous shells between the FS and RS. An evolution of the electron spectrum is also followed individually for every shell. We then find the FS and RS light curves by integrating over the entire FS and RS shocked regions, respectively. In particular, we make use of an analytic expression for observed spectral flux, which we derive here in terms of an observed frequency and observer time. Exploring a total of 20 different ejecta stratifications, we explain in detail how a stratified ejecta affects its blast wave dynamics and afterglow light curves. We show that, while the FS light curves are not sensitive to the ejecta stratifications, the RS light curves exhibit much richer features, including steep declines, plateaus, bumps, re-brightenings, and a variety of temporal decay indices. These distinctive RS features may be observable if the RS has higher values of the micophysics parameters than the FS. We discuss possible applications of our results in understanding the GRB afterglow data.

GRB 041219A is one of the longest and brightest gamma-ray bursts (GRBs) ever observed. It was discovered by the INTEGRAL satellite, and thanks to a precursor happening about 300 s before the bulk of the burst, ground based telescopes were able to catch the rarely-observed prompt emission in the optical and in the near infrared bands. Here we present the detailed analysis of its prompt gamma-ray emission, as observed with IBIS on board INTEGRAL, and of the available X-ray afterglow data collected by XRT on board Swift. We then present the late-time multi-band near infrared imaging data, collected at the TNG, and the CFHT, that allowed us to identify the host galaxy of the GRB as an under-luminous, irregular galaxy of about 5x10^9 M_Sun at best fit redshift of z=0.31 -0.26 +0.54. We model the broad-band prompt optical to gamma-ray emission of GRB 041219A within the internal shock model. We were able to reproduce the spectra and light curve invoking the synchrotron emission of relativistic electrons accelerated by a series of propagating shock waves inside a relativistic outflow. On the other hand, it is less easy to simultaneously reproduce the temporal and spectral properties of the infrared data.