Pulsars have been identified as good candidates for the acceleration of cosmic rays, up to ultra-high energies. However, a precise description of the acceleration processes at play is still to be established. Using 2D particle-in-cell simulations, we study proton acceleration in axisymmetric pulsar magnetospheres. Protons and electrons are extracted from the neutron star surface by the strong electric field induced by the rotation of the star, and electrons and positrons are produced in the magnetosphere through pair production process. As pair production has a crucial impact on electromagnetic fields, on gaps and thus on particle acceleration, we study its influence on the maximum energy and luminosity of protons escaping the magnetosphere. Protons are accelerated and escape in all our simulations. However, the acceleration sites are different for the protons and the pairs. As shown in previous studies, pairs are accelerated to their highest energies at the Y-point and in the equatorial current sheet, where magnetic reconnection plays and important role. In contrast, protons gain most of their kinetic energy below the light-cylinder radius within the separatrix current layers, but they are not confined within the equatorial current sheet. Their maximum Lorentz factors can reach $15\%$ to $75\%$ of the maximum Lorentz factor obtained by acceleration through the full vacuum potential drop from pole to equator, and increase with decreasing pair production. Their luminosity can reach $0.2\%$ to $4\%$ of the theoretical spin down luminosity of an aligned pulsar, and the minimum luminosity is obtained at the transition between the force-free and electrosphere regimes. These estimates support that millisecond pulsars could accelerate cosmic rays up to PeV energies and that new born millisecond pulsars could accelerate cosmic rays up to ultra-high energies.

Black holes are known to launch powerful relativistic jets and emit highly variable gamma radiation. How these jets are loaded with plasma remains poorly understood. Spark gaps are thought to drive particle acceleration and pair creation in the black-hole magnetosphere. In this paper, we perform 2D axisymmetric general-relativistic particle-in-cell simulations of a monopole black-hole magnetosphere with a realistic treatment of inverse Compton scattering and pair production. We find that the magnetosphere can self-consistently fill itself with plasma and activate the Blandford-Znajek mechanism. A highly time-dependent spark gap opens near the inner light surface which injects pair plasma into the magnetosphere. These results may account for the high-energy activity observed in active galactic nuclei and explain the origin of plasma at the base of the jet.

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.

In this Letter we propose that coherent radio emission of Crab, other young energetic pulsars, and millisecond pulsars is produced in the magnetospheric current sheet beyond the light cylinder. We carry out global and local two-dimensional kinetic plasma simulations of reconnection to illustrate the coherent emission mechanism. Reconnection in the current sheet beyond the light cylinder proceeds in the very efficient plasmoid-dominated regime, and current layer gets fragmented into a dynamic chain of plasmoids which undergo successive coalescence. Mergers of sufficiently large plasmoids produce secondary perpendicular current sheets, which are also plasmoid-unstable. Collisions of plasmoids with each other and with the upstream magnetic field eject fast-magnetosonic waves, which propagate upstream across the background field and successfully escape from the plasma as electromagnetic waves that fall in the radio band. This model successfully explains many important features of the observed radio emission from Crab and other pulsars with high magnetic field at the light cylinder: phase coincidence with the high-energy emission, nano-second duration (nanoshots), and extreme instantaneous brightness of individual pulses.

The recent detection of gravitational waves associated with a binary neutron star merger revives interest in interacting pulsar magnetospheres. Current models predict that a significant amount of magnetic energy should be released prior to the merger, leading to electromagnetic precursor emission. In this paper, we revisit this problem in the light of the recent progress in kinetic modeling of pulsar magnetospheres. We limit our work to the case of aligned magnetic moments and rotation axes, and thus neglect the orbital motion. We perform global two-dimensional axisymmetric particle-in-cell simulations of two pulsar magnetospheres merging at a rate consistent with the emission of gravitational waves. Both symmetric and asymmetric systems are investigated. Simulations show a significant enhancement of magnetic dissipation within the magnetospheres as both stars get closer. Even though the magnetospheric configuration depends on the relative orientations of the pulsar spins and magnetic axes, all configurations present nearly the same radiative signature, indicating that a common dissipation mechanism is at work. The relative motion of both pulsars drives magnetic reconnection at the boundary between the two magnetospheres, leading to efficient particle acceleration and high-energy synchrotron emission. Polar-cap discharge is also strongly enhanced in asymmetric configurations, resulting in vigorous pair production and potentially additional high-energy radiation. We observe an increase in the pulsar radiative efficiency by two orders of magnitude over the last orbit before the merger exceeding the spindown power of an isolated pulsar. The expected signal is too weak to be detected at high energies even in the nearby universe. However, if a small fraction of this energy is channeled into radio waves, it could be observed as a non-repeating fast radio burst.

The origin of the pulsed gamma-ray emission in pulsars remains an open issue. The combination of sensitive observations in the GeV domain by AGILE and {\em Fermi}-LAT and increasingly sophisticated numerical simulations have recently brought new insights into our understanding of the pulsed emission and particle acceleration processes in pulsars. Particle-in-cell simulations of pulsar magnetospheres show that the equatorial current sheet forming beyond the light cylinder is the main culprit for magnetic dissipation, particle acceleration and bright high-energy synchrotron radiation all together. The shinning current sheet naturally results in a pulse of light each time the sheet crosses our line of sight, which happens twice in most cases. Synthetic lightcurves present robust features reminiscent of observed gamma-ray pulsars by the {\em Fermi}-LAT and AGILE, opening up new perspectives for direct comparison between simulations and observations.

Black holes drive powerful plasma jets to relativistic velocities. This plasma should be collisionless, and self-consistently supplied by pair creation near the horizon. We present general-relativistic collisionless plasma simulations of Kerr-black-hole magnetospheres which begin from vacuum, inject electron-positron pairs based on local unscreened electric fields, and reach steady states with electromagnetically powered Blandford-Znajek jets and persistent current sheets. Particles with negative energy-at-infinity are a general feature, and can contribute significantly to black-hole rotational-energy extraction in a variant of the Penrose process. The generated plasma distribution depends on the pair-creation environment, and we describe two distinct realizations of the force-free electrodynamic solution. This sensitivity suggests that plasma kinetics will be useful in interpreting future horizon-resolving submillimeter and infrared observations.

We investigate the dynamics and emission of a starved magnetospheric region (gap) formed in the vicinity of a Kerr black hole horizon, using a new, fully general relativistic particle-in-cell code that implements Monte Carlo methods to compute gamma-ray emission and pair production through the interaction of pairs and gamma rays with soft photons emitted by the accretion flow. It is found that when the Thomson length for collision with disk photons exceeds the gap width, screening of the gap occurs through low-amplitude, rapid plasma oscillations that produce self-sustained pair cascades, with quasi-stationary pair and gamma-ray spectra, and with a pair multiplicity that increases in proportion to the pair production opacity. The gamma-ray spectrum emitted from the gap peaks in the TeV band, with a total luminosity that constitutes a fraction of about $10^{-5}$ of the corresponding Blandford-Znajek power. This stage is preceded by a prompt discharge phase of duration $\sim r_g/c$, during which the potential energy initially stored in the gap is released as a flare of curvature TeV photons. We speculate that the TeV emission observed in M87 may be produced by pair discharges in a spark gap.

Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron-ion plasma, spanning a wide range of ambient ion magnetizations $\sigma_i$, from the semirelativistic regime (ultrarelativistic electrons but nonrelativistic ions, 0.001<<$\sigma_i$<<1) to the fully relativistic regime (both species are ultrarelativistic, $\sigma_i$>>1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and nonthermal particle acceleration depend on $\sigma_i$. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvenic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the semirelativistic, small-$\sigma_i$ regime; (3) while the released magnetic energy is partitioned equally between electrons and ions in the ultrarelativistic limit, the electron energy fraction declines gradually with decreased $\sigma_i$ and asymptotes to about 0.25 in the semirelativistic regime; (4) reconnection leads to efficient nonthermal electron acceleration with a $\sigma_i$-dependent power-law index, $p(\sigma_i) \simeq $const$+0.7 {\sigma_i}^{-1/2}$. These findings are important for understanding black hole systems and lend support to semirelativistic reconnection models for powering nonthermal emission in blazar jets, offering a natural explanation for the spectral indices observed in these systems.

Rapidly rotating neutron stars blow a relativistic, magnetized wind mainly composed of electron-positron pairs. The free expansion of the wind terminates far from the neutron star where a weakly magnetized pulsar wind nebula forms, implying efficient magnetic dissipation somewhere upstream. The wind current sheet that separates the two magnetic polarities is usually considered as the most natural place for magnetic dissipation via relativistic reconnection, but its efficiency remains an open question. Here, the goal of this work is to revisit this issue in light of the most recent progress in the understanding of reconnection and pulsar electrodynamics. We perform large two-dimensional particle-in-cell simulations of the oblique rotator to capture the multi-scale evolution of the wind. We find that the current sheet breaks up into a dynamical chain of magnetic islands separated by secondary thin current sheets. The sheet thickness increases linearly with radius while the Poynting flux decreases monotonically as reconnection proceeds. The radius of complete annihilation of the stripes is given by the plasma multiplicity parameter at the light cylinder. Current starvation within the sheets does not occur before complete dissipation as long as there is enough charges where the sheets form. Particles are efficiently heated up to a characteristic energy set by the magnetization parameter at the light cylinder. Energetic pulsed synchrotron emission peaks close to the light cylinder, and presents sub-pulse variability associated with the formation of plasmoids in the sheet. This study suggests that the striped component of the wind dissipates far before reaching the termination shock in isolated pulsars, even in very-high-multiplicity systems such as the Crab pulsar. Pulsars in binary systems may provide the best environments to study magnetic dissipation in the wind.

We discuss the role of particle-in-cell (PIC) simulations in unveiling the origin of the emitting particles in PWNe. After describing the basics of the PIC technique, we summarize its implications for the quiescent and the flaring emission of the Crab Nebula, as a prototype of PWNe. A consensus seems to be emerging that, in addition to the standard scenario of particle acceleration via the Fermi process at the termination shock of the pulsar wind, magnetic reconnection in the wind, at the termination shock and in the Nebula plays a major role in powering the multi-wavelength signatures of PWNe.

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.

Polarization is a powerful diagnostic tool to constrain the site of the high-energy pulsed emission and particle acceleration in gamma-ray pulsars. Recent particle-in-cell simulations of pulsar magnetosphere suggest that high-energy emission results from particles accelerated in the equatorial current sheet emitting synchrotron radiation. In this study, we re-examine the simulation data to compute the phase-resolved polarization properties. We find that the emission is mildly polarized and that there is an anticorrelation between the flux and the degree of linear polarization (on-pulse: ~15%, off-pulse: ~30%). The decrease of polarization during pulses is mainly attributed to the formation of caustics in the current sheet. Each pulse of light is systematically accompanied by a rapid swing of the polarization angle due to the change of the magnetic polarity when the line of sight passes through the current sheet. The optical polarization pattern observed in the Crab can be well-reproduced for a pulsar inclination angle ~60 degrees and an observer viewing angle ~130 degrees. The predicted high-energy polarization is a robust feature of the current sheet emitting scenario which can be tested by future X-ray and gamma-ray polarimetry instruments.

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 observer's 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.

Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $\sigma$ and system sizes $L$. The particle spectra are well-represented by a power law $\gamma^{-\alpha}$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to $\sigma$ and $L$, respectively. For large $L$ and $\sigma$, the power-law index $\alpha$ approaches about 1.2.

We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.

It has recently been demonstrated that self-consistent particle-in-cell simulations of low-obliquity pulsar magnetospheres in flat spacetime show weak particle acceleration and no pair production near the poles. We investigate the validity of this conclusion in a more realistic spacetime geometry via general-relativistic particle-in-cell simulations of the aligned pulsar magnetospheres with pair formation. We find that the addition of frame-dragging effect makes local current density along the magnetic field larger than the Goldreich-Julian value, which leads to unscreened parallel electric fields and the ignition of a pair cascade. When pair production is active, we observe field oscillations in the open field bundle which could be related to pulsar radio emission. We conclude that general relativistic effects are essential for the existence of pulsar mechanism in low obliquity rotators.

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.

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.

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

We present first-principles relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron-positron pairs. We inject secondary pairs at locations of primary energetic particles, whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere, with the Y-point and current sheet. Solutions with obliquities $\lt 40^{\circ}$ do not show pair production in the open field line region, because the local current density along magnetic field is below the Goldreich-Julian value. The bulk outflow in these solutions is charge separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved space-time or multipolar surface fields.

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.

The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and 2D particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in 3D, using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight.

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.

It is generally accepted that astrophysical sources cannot emit synchrotron radiation above 160 MeV in their rest frame. This limit is given by the balance between the accelerating electric force and the radiation reaction force acting on the electrons. The discovery of synchrotron gamma-ray flares in the Crab Nebula, well above this limit, challenges this classical picture of particle acceleration. To overcome this limit, particles must accelerate in a region of high electric field and low magnetic field. This is possible only with a non-ideal magnetohydrodynamic process, like magnetic reconnection. We present the first numerical evidence of particle acceleration beyond the synchrotron burnoff limit, using a set of 2D particle-in-cell simulations of ultra-relativistic pair plasma reconnection. We use a new code, Zeltron, that includes self-consistently the radiation reaction force in the equation of motion of the particles. We demonstrate that the most energetic particles move back and forth across the reconnection layer, following relativistic Speiser orbits. These particles then radiate >160 MeV synchrotron radiation rapidly, within a fraction of a full gyration, after they exit the layer. Our analysis shows that the high-energy synchrotron flux is highly variable in time because of the strong anisotropy and inhomogeneity of the energetic particles. We discover a robust positive correlation between the flux and the cut-off energy of the emitted radiation, mimicking the effect of relativistic Doppler amplification. A strong guide field quenches the emission of >160 MeV synchrotron radiation. Our results are consistent with the observed properties of the Crab flares, supporting the reconnection scenario.

We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides estimates of the total jet power and the jet magnetic field strength. Energetic constraints for the VHE flare are extremely tight: for an isotropic particle distribution they require a huge co-moving energy density in the emitting region and a very efficient radiative process. We disfavor hadronic processes due to their low radiative efficiency, as well as the synchrotron scenario recently proposed for the case of HE flares in the Crab Nebula, since the parameters needed to overcome the radiative losses are quite extreme. The VHE emission can be explained by the Synchrotron Self-Compton (SSC) mechanism for D_{VHE} ~ 20 or by the External Radiation Compton (ERC) mechanism involving the infrared radiation of the dusty torus for D_{VHE} ~ 50. After discussing several alternative scenarios, we propose that the extreme energy density constraint can be satisfied when the emission comes from highly anisotropic short-lived bunches of particles formed by the kinetic beaming mechanism in magnetic reconnection sites. By focusing the emitting particles into very narrow beams, this mechanism allows one to relax the causality constraint on the source size, decreasing the required energy density by 4 orders of magnitude.

We report on the first study of the angular distribution of energetic particles and radiation generated in relativistic collisionless electron-positron pair plasma reconnection, using two-dimensional particle-in-cell simulations. We discover a strong anisotropy of the particles accelerated by reconnection and the associated strong beaming of their radiation. The focusing of particles and radiation increases with their energy; in this sense, this "kinetic beaming" effect differs fundamentally from the relativistic Doppler beaming usually invoked in high-energy astrophysics, in which all photons are focused and boosted achromatically. We also present, for the first time, the modeling of the synchrotron emission as seen by an external observer during the reconnection process. The expected lightcurves comprise several bright symmetric sub-flares emitted by the energetic beam of particles sweeping across the line of sight intermittently, and exhibit super-fast time variability as short as about one tenth of the system light-crossing time. The concentration of the energetic particles into compact regions inside magnetic islands and particle anisotropy explain the rapid variability. This radiative signature of reconnection can account for the brightness and variability of the gamma-ray flares in the Crab Nebula and in blazars.

We study models of the gamma-ray emission of Cyg X-3 observed by Fermi. We calculate the average X-ray spectrum during the gamma-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 Klein-Nishina 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 ~ 300--1000 is required by the shape of the observed X-ray/gamma-ray spectrum in order to avoid overproducing the observed X-ray flux. The electrons giving rise to the observed \g-rays are efficiently cooled by Compton scattering, and the power-law index of the acceleration process is ~ 2.5--3. 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 ~ 2.5, and the kinetic power is ~ 10^38 erg/s, 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 gamma-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 IR and gamma-ray fluxes.

The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with classical models of particle acceleration and pulsar wind nebulae. We investigate whether particle acceleration in a magnetic reconnection layer can account for the puzzling properties of the flares. We numerically integrate relativistic test-particle orbits in the vicinity of the layer, including the radiation reaction force, and using analytical expressions for the large-scale electromagnetic fields. As they get accelerated by the reconnection electric field, the particles are focused deep inside the current layer where the magnetic field is small. The electrons suffer less from synchrotron losses and are accelerated to extremely high energies. Population studies show that, at the end of the layer, the particle distribution piles up at the maximum energy given by the electric potential drop and is focused into a thin fan beam. Applying this model to the Crab Nebula, we find that the emerging synchrotron emission spectrum peaks above 100 MeV and is close to the spectral shape of a single electron. The flare inverse Compton emission is negligible and no detectable emission is expected at other wavelengths. This mechanism provides a plausible explanation for the gamma-ray flares in the Crab Nebula and could be at work in other astrophysical objects such as relativistic jets in active galactic nuclei.

The recent discovery of day-long gamma-ray flares in the Crab Nebula, presumed to be synchrotron emission by PeV (10^{15} eV) electrons in milligauss magnetic fields, presents a strong challenge to particle acceleration models. The observed photon energies exceed the upper limit (~100 MeV) obtained by balancing the acceleration rate and synchrotron radiation losses under standard conditions where the electric field is smaller than the magnetic field. We argue that a linear electric accelerator, operating at magnetic reconnection sites, is able to circumvent this difficulty. Sufficiently energetic electrons have gyroradii so large that their motion is insensitive to small-scale turbulent structures in the reconnection layer and is controlled only by large-scale fields. We show that such particles are guided into the reconnection layer by the reversing magnetic field as they are accelerated by the reconnection electric field. As these electrons become confined within the current sheet, they experience a decreasing perpendicular magnetic field that may drop below the accelerating electric field. This enables them to reach higher energies before suffering radiation losses and hence to emit synchrotron radiation in excess of the 100 MeV limit, providing a natural resolution to the Crab gamma-ray flare paradox.

The microquasar Cygnus X-3 was detected at high energies by the gamma-ray space telescopes AGILE and Fermi. The 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. In this study, we examine quantitatively the effect of pair production in Cygnus X-3 and put constraints on the location of the gamma-ray source. 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, 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. 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 due to 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. The lack of 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.

Cygnus X-3 is an accreting high-mass X-ray binary composed of a Wolf-Rayet star and an unknown compact object, possibly a black hole. The gamma-ray space telescope Fermi found definitive evidence that high-energy emission is produced in this system. We propose a scenario to explain the GeV gamma-ray emission in Cygnus X-3. In this model, energetic electron-positron pairs are accelerated at a specific location in the relativistic jet, possibly related to a recollimation shock, and upscatter the stellar photons to high energies. The comparison with Fermi observations shows that the jet should be inclined close to the line of sight and pairs should not be located within the system. Energetically speaking, a massive compact object is favored. We report also on our investigations of the gamma-ray absorption of GeV photons with the radiation emitted by a standard accretion disk in Cygnus X-3. This study shows that the gamma-ray source should not lie too close to the compact object.

LS 5039 is a Galactic binary system emitting high and very-high energy gamma rays. The gamma-ray flux is modulated on the orbital period and the TeV lightcurve shaped by photon-photon annihilation. The observed very-high energy modulation can be reproduced with a simple leptonic model but fails to explain the flux detected by HESS at superior conjunction, where gamma rays are fully absorbed. The contribution from an electron-positron pair cascade could be strong and prevail over the primary flux at superior conjunction. The created pairs can be isotropized by the magnetic field, resulting in a three-dimensional cascade. The aim of this article is to investigate the gamma ray radiation from this pair cascade in LS 5039. This additional component could account for HESS observations at superior conjunction in the system. A semi-analytical and a Monte Carlo method for computing three-dimensional cascade radiation are presented and applied in the context of binaries. Three-dimensional cascade radiation contributes significantly at every orbital phase in the TeV lightcurve, and dominates close to superior conjunction. The amplitude of the gamma-ray modulation is correctly reproduced for an inclination of the orbit of about 40 degrees. Primary pairs should be injected close to the compact object location, otherwise the shape of the modulation is not explained. In addition, synchrotron emission from the cascade in X-rays constrains the ambient magnetic field to below 10 G. The radiation from a three-dimensional pair cascade can account for the TeV flux detected by HESS at superior conjunction in LS 5039, but the very-high energy spectrum at low fluxes remains difficult to explain in this model.

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. The X-ray and the very high energy (>100 GeV, VHE) gamma-ray emissions 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. 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 (LS I +61 303, PSR B1259-63). 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. Doppler boosting might provide an explanation for the puzzling phasing of the VHE peak in this system.

High energy gamma-rays have been detected from Cygnus 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 hr 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. Modeling 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.

Gamma-ray binaries are orbital modulated gamma-ray sources in the Galaxy detected both at GeV and TeV energies. The high-energy radiation may come from the interaction of energetic electrons injected by a young pulsar and photons from the massive companion star. We present a model for the production of high-energy radiation where emission, absorption and pair cascading are considered. New observations of LS I +61 303 by the Fermi Space Telescope revealed an exponential cut-off in the spectrum at a few GeV, inconsistent with gamma-gamma absorption. Electrons radiating at GeV and TeV have probably two different origins. We investigate whether the emission from the unshocked pulsar wind explains the GeV component in LS I +61 303.

In gamma-ray binaries such as LS 5039 a large number of electron-positron pairs are created by the annihilation of primary very high energy (VHE) gamma-rays with photons from the massive star. The radiation from these particles contributes to the total high energy gamma-ray flux and can initiate a cascade, decreasing the effective gamma-ray opacity in the system. The aim of this paper is to model the cascade emission and investigate if it can account for the VHE gamma-ray flux detected by HESS from LS 5039 at superior conjunction, where the primary gamma-rays are expected to be fully absorbed. A one-dimensional cascade develops along the line-of-sight if the deflections of pairs induced by the surrounding magnetic field can be neglected. A semi-analytical approach can then be adopted, including the effects of the anisotropic seed radiation field from the companion star. Cascade equations are numerically solved, yielding the density of pairs and photons. In LS 5039, the cascade contribution to the total flux is large and anti-correlated with the orbital modulation of the primary VHE gamma-rays. The cascade emission dominates close to superior conjunction but is too strong to be compatible with HESS measurements. Positron annihilation does not produce detectable 511 keV emission. This study provides an upper limit to cascade emission in gamma-ray binaries at orbital phases where absorption is strong. The pairs are likely to be deflected or isotropized by the ambient magnetic field, which will reduce the resulting emission seen by the observer. Cascade emission remains a viable explanation for the detected gamma-rays at superior conjunction in LS 5039.

Gamma-ray binaries have been uncovered as a new class of Galactic objects in the very high energy sky (> 100 GeV). The three systems known today have hard X-ray spectra (photon index ~ 1.5), extended radio emission and a high luminosity in gamma-rays. Recent monitoring campaigns of LSI +61 303 in X-rays have confirmed variability in these systems and revealed a spectral hardening with increasing flux. In a generic one-zone leptonic model, the cooling of relativistic electrons accounts for the main spectral and temporal features observed at high energy. Persistent hard X-ray emission is expected to extend well beyond 10 keV. We explain how Simbol-X will constrain the existing models in connection with Fermi Space Telescope measurements. Because of its unprecedented sensitivity in hard X-rays, Simbol-X will also play a role in the discovery of new gamma-ray binaries, giving new insights into the evolution of compact binaries.

LS 5039 and LSI +61$\degr$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) pulsar wind. 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$\degr$303 by computing the unshocked wind emission and comparing it to observations. 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 gamma-gamma absorption and the shock geometry. 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 difficult to avoid in LS 5039 and LSI +61$\degr$303. If the particles in the pulsar are monoenergetic then the observations constrain their energy to roughly 10-100 GeV. For more complex particle distributions, the free pulsar wind emission will be difficult to distinguish from the shocked pulsar wind emission.

Gamma-ray binaries, composed of a massive star and compact object, have been established as a new class of sources of very high energy (VHE) photons. The gamma-rays are produced by inverse Compton scattering of the stellar light by VHE electrons accelerated in the vicinity of the compact object. The VHE emission from LS 5039 displays an orbital modulation. The inverse Compton spectrum depends on the angle between the incoming and outgoing photon in the electron rest frame. Since the angle at which an observer sees the star and electrons changes with the orbit, a phase dependence of the spectrum is expected. The phase-dependent spectrum of LS 5039 is calculated, assuming a continuous injection of electrons. The shape of the electron distribution depends on the injected power-law and on the magnetic field intensity. Anisotropic scattering produces hard emission at inferior conjunction, when attenuation due to pair production of the VHE gamma-rays on star light is minimum. The computed lightcurve and spectra provide good fits to the HESS and EGRET observations, except at phases of maximum attenuation where pair cascade emission may be significant for HESS. Detailed predictions are made for a modulation in the GLAST energy range. The magnetic field intensity at periastron is 0.8+-0.2 G. Anisotropic inverse Compton scattering plays a major role in LS 5039. The derived magnetic field intensity, injection energy and slope suggest a rotation-powered pulsar wind nebula. Gamma-ray binaries are promising sources to study the environment of pulsars on small scales.