@article{BohdanNiemiecKobzaretal.2017, author = {Bohdan, Artem and Niemiec, Jacek and Kobzar, Oleh and Pohl, Martin}, title = {Electron Pre-acceleration at Nonrelativistic High-Mach-number Perpendicular Shocks}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {847}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa872a}, pages = {17}, year = {2017}, abstract = {We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit the extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta beta p =5.10(-4) and 0.5 beta p =), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane, while keeping it at 90 degrees to the shock normal. Ion reflection off of the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all of the cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field, the fraction of suprathermal electrons is much higher than for other configurations, because only in this case are the Buneman modes fully captured by the 2D simulation grid. Shocks in plasma with moderate bp provide more efficient pre-acceleration. The relevance of our results to the physics of fully 3D systems is discussed.}, language = {en} } @article{BohdanNiemiecPohletal.2019, author = {Bohdan, Artem and Niemiec, Jacek and Pohl, Martin and Matsumoto, Yosuke and Amano, Takanobu and Hoshino, Masahiro}, title = {Kinetic Simulations of Nonrelativistic Perpendicular Shocks of Young Supernova Remnants}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {878}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/ab1b6d}, pages = {11}, year = {2019}, abstract = {Electron injection at high Mach number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell simulations and tracing individual particles, we in detail analyze the shock-surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.}, language = {en} } @article{VafinDekaPohletal.2019, author = {Vafin, Sergei and Deka, Pranab Jyoti and Pohl, Martin and Bohdan, Artem}, title = {Revisit of Nonlinear Landau Damping for Electrostatic Instability Driven by Blazar-induced Pair Beams}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {873}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/ab017b}, pages = {12}, year = {2019}, abstract = {We revisit the effect of nonlinear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in k-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that nonlinear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wave numbers where wave driving is weak or absent. The relaxation time also strongly depends on the intergalactic medium temperature, T-IGM, and for T-IGM << 10 eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small k, that nonlinearly drains the wave energy at the resonant k of the pair-beam instability. Any other dissipation process operating at small k would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other nonlinear processes, such as the modulation instability, could provide additional damping of the nonresonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.}, language = {en} }