No Arabic abstract
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.
We use particle-in-magnetohydrodynamics-cells to model particle acceleration and magnetic field amplification in a high Mach, parallel shock in three dimensions and compare the result to 2-D models. This allows us to determine whether 2-D simulations can be relied upon to yield accurate results in terms of particle acceleration, magnetic field amplification and the growth rate of instabilities. Our simulations show that the behaviour of the gas and the evolution of the instabilities are qualitatively similar for both the 2-D and 3-D models, with only minor quantitative differences that relate primarily to the growth speed of the instabilities. The main difference between 2-D and 3-D models can be found in the spectral energy distributions (SEDs) of the non-thermal particles. The 2-D simulations prove to be more efficient, accelerating a larger fraction of the particles and achieving higher velocities. We conclude that, while 2-D models are sufficient to investigate the instabilities in the gas, their results have to be treated with some caution when predicting the expected SED of a given shock.
We study diffusive shock acceleration (DSA) of electrons in non-relativistic quasi-perpendicular shocks using self-consistent one-dimensional particle-in-cell (PIC) simulations. By exploring the parameter space of sonic and Alfv{e}nic Mach numbers we find that high Mach number quasi-perpendicular shocks can efficiently accelerate electrons to power-law downstream spectra with slopes consistent with DSA prediction. Electrons are reflected by magnetic mirroring at the shock and drive non-resonant waves in the upstream. Reflected electrons are trapped between the shock front and upstream waves and undergo multiple cycles of shock drift acceleration before the injection into DSA. Strong current-driven waves also temporarily change the shock obliquity and cause mild proton pre-acceleration even in quasi-perpendicular shocks, which otherwise do not accelerate protons. These results can be used to understand nonthermal emission in supernova remnants and intracluster medium in galaxy clusters.
A new particle acceleration process in a developing Alfv{e}n turbulence in the course of successive parametric instabilities of a relativistic pair plasma is investigated by utilyzing one-dimensional electromagnetic full particle code. Coherent wave-particle interactions result in efficient particle acceleration leading to a power-law like energy distribution function. In the simulation high energy particles having large relativistic masses are preferentially accelerated as the turbulence spectrum evolves in time. Main acceleration mechanism is simultaneous relativistic resonance between a particle and two different waves. An analytical expression of maximum attainable energy in such wave-particle interactions is derived.
We develop a model of particle acceleration in explosive reconnection events in relativistic magnetically-dominated plasmas and apply it to explain gamma-ray flares from the Crab Nebula. The model relies on development of current-driven instabilities on macroscopic scales (not related to plasma skin depths). Using analytical and numerical methods (fluid and particle-in-cell simulations), we study a number of model problems in relativistic magnetically-dominated plasma: (i) we extend Syrovatskys classical model of explosive X-point collapse to magnetically-dominated plasmas; (ii) we consider instability of two-dimensional force-free system of magnetic flux tubes; (iii) we consider merger of two zero total poloidal current magnetic flux tubes. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization. The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced reconnection, further accelerating the particles, but proceeds at a slower reconnection rate. Crab flares result from the initial explosive stages of magnetic island mergers of magnetic flux tubes produced in the bulk of nebula at intermediate polar regions. The post-termination shock plasma flow in the wind sectors with mild magnetization naturally generates large-scale highly magnetized structures. Internal kink-like instabilities lead to the formation of macroscopic current-carrying magnetic flux tubes that merge explosively.
We develop a model of gamma-ray flares of the Crab Nebula resulting from the magnetic reconnection events in highly-magnetized relativistic plasma. We first discuss physical parameters of the Crab nebula and review the theory of pulsar winds and termination shocks. We also review the principle points of particle acceleration in explosive reconnection events (Lyutikov et al. 2017a,b). It is required that particles producing flares are accelerated in highly magnetized regions of the nebula. Flares originate from the poleward regions at the base of Crabs polar outflow, where both the magnetization and the magnetic field strength are sufficiently high. The post-termination shock flow develops macroscopic (not related to the plasma properties on the skin-depth scale) kink-type instabilities. The resulting large-scales magnetic stresses drive explosive reconnection events on the light-crossing time of the reconnection region. Flares are produced at the initial stage of the current sheet development, during the X-point collapse. The model has all the ingredients needed for Crab flares: natural formation of highly magnetized regions, explosive dynamics on light travel time, development of high electric fields on macroscopic scales and acceleration of particles to energies well exceeding the average magnetic energy per particle.