Normalized to: Rafighi, I.
[1]
oai:arXiv.org:1803.02990 [pdf] - 1670642
The electrostatic instability for realistic pair distributions in
blazar/EBL cascades
Submitted: 2018-03-08
This work revisits the electrostatic instability for blazar-induced pair
beams propagating through IGM with the methods of linear analysis and PIC
simulations. We study the impact of the realistic distribution function of
pairs resulting from interaction of high-energy gamma-rays with the
extragalactic background light. We present analytical and numerical
calculations of the linear growth rate of the instability for arbitrary
orientation of wave vectors. Our results explicitly demonstrate that the finite
angular spread of the beam dramatically affects the growth rate of the waves,
leading to fastest growth for wave vectors quasi-parallel to the beam direction
and a growth rate at oblique directions that is only by a factor of 2-4 smaller
compared to the maximum. To study the non-linear beam relaxation, we performed
PIC simulations that take into account a realistic wide-energy distribution of
beam particles. The parameters of the simulated beam-plasma system provide an
adequate physical picture that can be extrapolated to realistic blazar-induced
pairs. In our simulations the beam looses only 1\% percent of its energy, and
we analytically estimate that the beam would lose its total energy over about
$100$ simulation times. Analytical scaling is then used to extrapolate to the
parameters of realistic blazar-induced pair beams. We find that they can
dissipate their energy slightly faster by the electrostatic instability than
through inverse-Compton scattering. The uncertainties arising from, e.g.,
details of the primary gamma-ray spectrum are too large to make firm statements
for individual blazars, and an analysis based on their specific properties is
required.
[3]
oai:arXiv.org:1708.07797 [pdf] - 1593597
Plasma effects on relativistic pair beams from TeV blazars: PIC
simulations and analytical predictions
Submitted: 2017-08-25
Pair beams produced by very high-energy radiation from TeV blazars emit gamma
rays in the GeV band by inverse-Compton scattering of soft photons. The
observed GeV-band signal is smaller than that expected from the full
electromagnetic cascade. This means that the pair beams must be affected by
other physical processes reducing their energy flux. One possible loss
mechanism involves beam-plasma instabilities that we consider in the present
work. For realistic parameters the pair beams can not be simulated by modern
computers. Instead, we use a simple analytical model to find a range of the
beam parameters that (i) provides a physical picture similar to that of
realistic pair beams and (ii) at the same time can be handled by available
computational resources. Afterwards, we performed corresponding 2D PIC
simulations. We confirm that the beams experience only small changes in the
relevant parameter regime, and other processes such as deflection in magnetic
field must be at play.
[5]
oai:arXiv.org:1602.05064 [pdf] - 1378944
Non-relativistic perpendicular shocks modeling young supernova remnants:
nonstationary dynamics and particle acceleration at forward and reverse
shocks
Submitted: 2016-02-16
For parameters that are applicable to the conditions at young supernova
remnants, we present results of 2D3V particle-in-cell simulations of a
non-relativistic plasma shock with a large-scale perpendicular magnetic field
inclined at 45-deg angle to the simulation plane to approximate 3D physics. We
developed an improved clean setup that uses the collision of two plasma slabs
with different density and velocity, leading to the development of two
distinctive shocks and a contact discontinuity. The shock formation is mediated
by Weibel-type filamentation instabilities that generate magnetic turbulence.
Cyclic reformation is observed in both shocks with similar period, for which we
note global variations on account of shock rippling and local variations
arising from turbulent current filaments. The shock rippling occurs on spatial
and temporal scales given by gyro-motions of shock-reflected ions. The drift
motion of electrons and ions is not a gradient drift, but commensurates with E
x B drift. We observe a stable suprathermal tail in the ion spectra, but no
electron acceleration because the amplitude of Buneman modes in the shock foot
is insufficient for trapping relativistic electrons. We see no evidence of
turbulent reconnection. A comparison with other 2D simulation results suggests
that the plasma beta and the ion-to-electron mass ratio are not decisive for
efficient electron acceleration, but pre-acceleration efficacy might be reduced
with respect to the 2D results once three-dimensional effects are fully
accounted for. Other microphysical factors may also be at play to limit the
amplitude of Buneman waves or prevent return of electrons to the foot region.