Benoit Cerutti

The relativistic jet of Cygnus X−3 in gamma‐rays

High-energy gamma-rays have been detected from Cyg X−3, a system composed of a Wolf–Rayet star and a black hole or neutron star. The gamma-ray emission is linked to the radio emission from the jet launched in the system. The flux is modulated with the 4.8-h orbital period, as expected if high-energy electrons are upscattering photons emitted by the Wolf–Rayet star to gamma-ray energies. This modulation is computed assuming that high-energy electrons are located at some distance along a relativistic jet of arbitrary orientation. Modelling shows that the jet must be inclined and that the gamma-ray emitting electrons cannot be located within the system. This is consistent with the idea that the electrons gain energy, where the jet is recollimated by the stellar wind pressure and forms a shock. Jet precession should strongly affect the gamma-ray modulation shape at different epochs. The power in non-thermal electrons represents a small fraction of the Eddington luminosity only if the inclination is low, i.e. if the compact object is a black hole.

Particle acceleration in axisymmetric pulsar current sheets

The equatorial current sheet in pulsar magnetospheres is often regarded as an ideal site for particle acceleration via relativistic reconnection. Using 2D spherical particle-in-cell simulations, we investigate particle acceleration in the axisymmetric pulsar magnetosphere as a function of the injected plasma multiplicity and magnetization. We observe a clear transition from a highly charge-separated magnetosphere for low plasma injection with little current and spin-down power, to a nearly force-free solution for high plasma multiplicity characterized by a prominent equatorial current sheet and high spin-down power. We find significant magnetic dissipation in the current sheet, up to 30% within 5 light-cylinder radii in the high-multiplicity regime. The simulations unambiguously demonstrate that the dissipated Poynting flux is efficiently channeled to the particles in the sheet, close to the Y-point within about 1-2 light cylinder radii from the star. The mean particle energy in the sheet is given by the upstream plasma magnetization at the light cylinder. The study of particle orbits shows that all energetic particles originate from the boundary layer between the open and the closed field lines. Energetic positrons always stream outward, while high-energy electrons precipitate back towards the star through the sheet and along the separatrices, which may result in auroral-like emission. Our results suggest that the current sheet and the separatrices may be the main source of high-energy radiation in young pulsars.

Relativistic Doppler-boosted emission in gamma-ray binaries

Context. Gamma-ray binaries could be compact pulsar wind nebulae formed when a young pulsar orbits a massive star. The pulsar wind is contained by the stellar wind of the O or Be companion, creating a relativistic comet-like structure accompanying the pulsar along its orbit. Aims. The X-ray and the very high energy (>100 GeV, VHE) gamma-ray emission from the binary LS 5039 are modulated on the orbital period of the system. Maximum and minimum flux occur at the conjunctions of the orbit, suggesting that the explanation is linked to the orbital geometry. The VHE modulation has been proposed to be due to the combined effect of Compton scattering and pair production on stellar photons, both of which depend on orbital phase. The X-ray modulation could be due to relativistic Doppler boosting in the comet tail where both the X-ray and VHE photons would be emitted. Methods. Relativistic aberrations change the seed stellar photon flux in the comoving frame so Doppler boosting affects synchrotron and inverse Compton emission differently. The dependence with orbital phase of relativistic Doppler-boosted (isotropic) synchrotron and (anisotropic) inverse Compton emission is calculated, assuming that the flow is oriented radially away from the star (LS 5039) or tangentially to the orbit (LSI +61°303, PSR B1259-63). Results. Doppler boosting of the synchrotron emission in LS 5039 produces a lightcurve whose shape corresponds to the X-ray modulation. The observations imply an outflow velocity of 0.15-0.33c consistent with the expected flow speed at the pulsar wind termination shock. In LS I +61°303, the calculated Doppler boosted emission peaks in phase with the observed VHE and X-ray maximum. Conclusions. Doppler boosting is not negligible in gamma-ray binaries, even for mildly relativistic speeds. The boosted modulation reproduces the X-ray modulation in LS 5039 and could also provide an explanation for the puzzling phasing of the VHE peak in LS I +61 °303.

Modelling high-energy pulsar light curves from first principles

Current models of gamma-ray lightcurves in pulsars suffer from large uncertainties on the precise location of particle acceleration and radiation. Here, we present an attempt to alleviate these difficulties by solving for the electromagnetic structure of the oblique magnetosphere, particle acceleration, and the emission of radiation self-consistently, using 3D spherical particle-in-cell simulations. We find that the low-energy radiation is synchro-curvature radiation from the polar-cap regions within the light cylinder. In contrast, the high-energy emission is synchrotron radiation that originates exclusively from the Y-point and the equatorial current sheet where relativistic magnetic reconnection accelerates particles. In most cases, synthetic high-energy lightcurves contain two peaks that form when the current sheet sweeps across the observers line of sight. We find clear evidence of caustics in the emission pattern from the current sheet. High-obliquity solutions can present up to two additional secondary peaks from energetic particles in the wind region accelerated by the reconnection-induced flow near the current sheet. The high-energy radiative efficiency depends sensitively on the viewing angle, and decreases with increasing pulsar inclination. The high-energy emission is concentrated in the equatorial regions where most of the pulsar spindown is released and dissipated. These results have important implications for the interpretation of gamma-ray pulsar data.

The gamma‐ray emitting region of the jet in Cyg X‐3

We study models of the ?-ray emission of Cyg X-3 observed by Fermi. We calculate the average X-ray spectrum during the ?-ray active periods. Then, we calculate spectra from Compton scattering of a photon beam into a given direction by isotropic relativistic electrons with a power-law distribution, both based on the KleinNishina cross-section and in the Thomson limit. Applying the results to scattering of stellar blackbody radiation in the inner jet of Cyg X-3, we find that a low-energy break in the electron distribution at a Lorentz factor of similar to 300103 is required by the shape of the observed X-ray/?-ray spectrum in order to avoid overproducing the observed X-ray flux. The electrons giving rise to the observed ?-rays are efficiently cooled by Compton scattering, and the power-law index of the acceleration process is ?2.53. The bulk Lorentz factor of the jet and the kinetic power before the dissipation region depend on the fraction of the dissipation power supplied to the electrons; if it is ?1/2, the Lorentz factor is similar to 2.5, and the kinetic power is similar to 1038 erg s-1, which represents a firm lower limit on the jet power, and is comparable to the bolometric luminosity of Cyg X-3. Most of the power supplied to the electrons is radiated. The broad-band spectrum constrains the synchrotron and self-Compton emission from the ?-ray emitting electrons, which requires the magnetic field to be relatively weak, with the magnetic energy density ? a few times 10-3 of that in the electrons. The actual value of the magnetic field strength can be inferred from a future simultaneous measurement of the infrared and ?-ray fluxes.

Pulsar-Wind Nebulae

In this review we describe recent observational and theoretical developments in our understanding of pulsar winds and pulsar-wind nebulae (PWNe). We put special emphasis on the results from observations of well-characterized PWNe of various types (e.g., torus-jet and bowshock-tail), the most recent MHD modeling efforts, and the status of the flaring Crab PWN puzzle.

Gamma-ray flares in the Crab Nebula: A case of relativistic reconnection?

The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by shocks. Here, it is argued that the flares may be the smoking gun of magnetic dissipation in the Nebula. Using 2D and 3D particle-in-cell simulations, it is shown that the observations are consistent with relativistic magnetic reconnection, where pairs are subject to strong radiative cooling. The Crab flares may highlight the importance of relativistic magnetic reconnection in astrophysical sources.

Spectral signature of a free pulsar wind in the gamma-ray binaries LS 5039 and LSI +61 ◦ 303

Context. LS 5039 and LSI +61 ◦ 303 are two binaries that have been detected in the TeV energy domain. These binaries are composed of a massive star and a compact object, possibly a young pulsar. The gamma-ray emission would be due to particle acceleration at the collision site between the relativistic pulsar wind and the stellar wind of the massive star. Part of the emission may also originate from inverse Compton scattering of stellar photons on the unshocked (free) pulsa r wind. Aims. The purpose of this work is to constrain the bulk Lorentz factor of the pulsar wind and the shock geometry in the compact pulsar wind nebula scenario for LS 5039 and LSI +61 ◦ 303 by computing the unshocked wind emission and comparing it to observations. Methods. Anisotropic inverse Compton losses equations are derived and applied to the free pulsar wind in binaries. The unshocked wind spectra seen by the observer are calculated taking into account theγ−γ absorption and the shock geometry. Results. A pulsar wind composed of monoenergetic pairs produces a typical sharp peak at an energy which depends on the bulk Lorentz factor and whose amplitude depends on the size of the emitting region. This emission from the free pulsar wind is found to be strong and diffi cult to avoid in LS 5039 and LSI+61 ◦ 303. Conclusions. If the particles in the pulsar are monoenergetic then the obs ervations constrain their energy to roughly 10-100 GeV. For more complex particle distributions, the free pulsar wind emiss ion will be diffi cult to distinguish from the shocked pulsar wind emission.

What caused the GeV flare of PSR B1259-63?

PSR B1259-63 is a gamma-ray binary system composed of a high spindown pulsar and a massive star. Non-thermal emission up to TeV energies is observed near periastron passage, attributed to emission from high energy e+e- pairs accelerated at the shock with the circumstellar material from the companion star, resulting in a small-scale pulsar wind nebula. Weak gamma-ray emission was detected by the Fermi/LAT at the last periastron passage, unexpectedly followed 30 days later by a strong flare, limited to the GeV band, during which the luminosity nearly reached the spindown power of the pulsar. The origin of this GeV flare remains mysterious. We investigate whether the flare could have been caused by pairs, located in the vicinity of the pulsar, up-scattering X-ray photons from the surrounding pulsar wind nebula rather than UV stellar photons, as usually assumed. Such a model is suggested by the geometry of the interaction region at the time of the flare. We compute the gamma-ray lightcurve for this scenario, based on a simplified description of the interaction region, and compare it to the observations. The GeV lightcurve peaks well after periastron with this geometry. The pairs are inferred to have a Lorentz factor ~500. They also produce an MeV flare with a luminosity ~1e34 erg/s prior to periastron passage. A significant drawback is the very high energy density of target photons required for efficient GeV emission. We propose to associate the GeV-emitting pairs with the Maxwellian expected at shock locations corresponding to high pulsar latitudes, while the rest of the non-thermal emission arises from pairs accelerated in the equatorial region of the pulsar wind termination shock.

Absorption of high-energy gamma rays in Cygnus X-3

Context. The microquasar Cygnus X-3 was detected at high energies by the gamma-ray space telescopes AGILE and Fermi .T he gamma-ray emission is transient, modulated with the orbital period and seems related to major radio flares, i.e. to the relativistic jet. The GeV gamma-ray flux can be substantially attenuated by internal absorption with the ambient X-rays. Aims. We examine quantitatively the effect of pair production in Cygnus X-3 and put constraints on the location of the gamma-ray source. Methods. Cygnus X-3 exhibits complex temporal and spectral patterns in X-rays. During gamma-ray flares, the X-ray emission can be approximated by a bright disk black-body component and a non-thermal tail extending in hard X-rays, which is possibly related to a corona above the disk. We calculate numerically the exact optical depth for gamma rays above a standard accretion disk. Emission and absorption in the corona are also investigated. Results. GeV gamma rays are significantly absorbed by soft X-rays emitted from the inner parts of the accretion disk. The absorption pattern is complex and anisotropic. Isotropization of X-rays caused by Thomson scattering in the companion-star wind tends to increase the gamma-ray opacity. Gamma rays from the corona suffer from strong absorption by photons from the disk and cannot explain the observed high-energy emission, unless the corona is unrealistically extended. Conclusions. The lack of an absorption feature in the GeV emission indicates that high-energy gamma rays should be located at a minimum distance ∼10 8 −10 10 cm from the compact object. The gamma-ray emission is unlikely to have a coronal origin.