Romain Hascoet

ON THE ORIGIN OF GeV EMISSION IN GAMMA-RAY BURSTS

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. 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 by the prompt radiation ahead of the blast wave. 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 minute 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.

Prompt thermal emission in gamma-ray bursts

Context. Gamma-ray burst (GRB) spectra globally appear non-thermal, but recent observations of a few bursts with Fermi GBM have confirmed previous indications from BATSE of the presence of an underlying thermal component. Photospheric emission is indeed expected when the relativistic outflow emerging from the central engine becomes transparent to its own radiation, with a quasi- blackbody spectrum in absence of additional sub-photospheric dissipation. However, its intensity strongly depends on the acceleration mechanism - thermal or magnetic - of the flow. Aims. We aim to compute the thermal and non-thermal emissions (light curves and spectra) produced by an outflow with a variable Lorentz factor, where the power u Eiso injected at the origin is partially thermal (fractionth ≤ 1) and partially magnetic (fraction 1− � th). The thermal emission is produced at the photosphere, and the non-thermal emission in the optically thin regime. Apart from the value ofth, we want to test how the other model parameters affect the observed ratio of the thermal to non-thermal emission. Methods. We followed the adiabatic cooling of the flow from the origin to the photosphere and computed the emitted radiation, which is a sum of modified black bodies at different temperatures (as the temperature strongly depends on the Lorentz factor of each shell at transparency). If the non-thermal emission comes from internal shocks, it is obtained from a multi-shell model where a fraction of the energy dissipated in shell collision is transferred to electrons and radiated via the synchrotron mechanism. If, conversely, the non-thermal emission originates in magnetic reconnection, the lack of any detailed theory for this process forced us to use a very simple parametrisation to estimate the emitted spectrum. Results. If the non-thermal emission is made by internal shocks, we self-consistently obtained the light curves and spectra of the thermal and non-thermal components for any distribution of the Lorentz factor in the flow. If the non-thermal emission results from magnetic reconnection we were unable to produce a light curve and could only compare the respective non-thermal and thermal spectra. In the different considered cases, we varied the model parameters to see when the thermal component in the light curve and/or spectrum is likely to show up or, on the contrary, to be hidden. We finally compared our results to the proposed evidence for the presence of a thermal component in GRB spectra. Focussing on GRB 090902B and GRB 10072B, we showed how these observations can be used to constrain the nature and acceleration mechanism of GRB outflows.

Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

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.

Shocks in nova outflows – I. Thermal emission

Growing evidence for shocks in nova outflows include (1) mult iple velocity components in the optical spectra; (2) hard X-ray emission starting wee ks to months after the outburst; (3) an early radio flare on timescales of months, in excess of t hat predicted from the freely expanding photo-ionized gas; and, perhaps most dramatically, (4)∼ GeV gamma-ray emission. We present a one dimensional model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated th ermal X-ray, optical, and radio emission. The lower velocity DES could represent an earlier stage of mass loss from the white dwarf or ambient material not directly related to the thermonuclear runaway. 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. Our predicted 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 producing classical nova V1324 Sco, obtaining an adequate fit to the ear ly radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a characteristic velocity∼ 10 3 km s −1 and mass∼ few 10 −4 M⊙; 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 inferr ed independently by modeling the late radio peak as uniformly expanding photo-ionized gas. Rapid evolution of the early radio light curves require the DES to 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 as the result of the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times are found to contribute significantly to the optical/UV emission, which we speculate may be responsible for the previously unexplained ‘plateaus’ and secondary maxima in nova optical light curve s.

NuSTAR OBSERVATIONS OF MAGNETAR 1E 1841−045

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. 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, in agreement with previous multimission X-ray results. However, we are unable to reproduce the spectral results reported based on Suzaku observations. The pulsed fraction of the source is found to increase with photon energy. The measured rms pulsed fractions are ~12% and ~17% at ~20 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 by Beloborodov 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° and on the angle between the rotation axis and line of sight of ~50°. 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.

Phase-resolved X-ray spectra of magnetars and the coronal outflow model

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.

NuSTAR Observations of the Magnetar 1E 2259+586

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, suggest 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.

Accounting for the XRT early steep decay in models of the prompt gamma-ray burst emission

Context. The Swift-XRT observations of the early X-ray afterglow of GRBs show that it usually begins with a steep decay phase. Aims. A possible origin of this early steep decay is the high latitude emission that subsists when the on-axis emission of the last dissipating regions in the relativistic outflow has been switched-off. We wish to establish which of various models of the prompt emission are compatible with this interpretation. Methods. We successively consider internal shocks, photospheric emission, and magnetic reconnection and obtain the typical decay timescale at the end of the prompt phase in each case. Results. Only internal shocks naturally predict a decay timescale comparable to the burst duration, as required to explain XRT observations in terms of high latitude emission. The decay timescale of the high latitude emission is much too short in photospheric models and the observed decay must then correspond to an effective and generic behavior of the central engine. Reconnection models require some ad hoc assumptions to agree with the data, which will have to be validated when a better description of the reconnection process becomes available.

Measuring Ambient Densities and Lorentz Factors of Gamma-Ray Bursts from GeV and Optical Observations

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 ± 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 = ρr 2 varies in the sample around 10 11 g cm 1 . 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.

DEEP NUSTAR AND SWIFT MONITORING OBSERVATIONS OF THE MAGNETAR 1E 1841−045

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 T90 ≤ 1 s. An elevated level of emission tail is detected after the brightest burst, persisting for ∼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 phaseresolved 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 for over approximately 10 yr. 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, amag = 0.25, the twisted magnetic flux, 2.5 × 10^(26) G cm^2, and the power released in the twisted magnetosphere, L_j = 6 × 10^(36) 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 yr 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.