Particle acceleration and wave excitation in quasi-parallel high-Mach-number collisionless shocks: Particle-in-cell simulation

We herein investigate shock formation and particle acceleration processes for both protons and electrons in a quasi-parallel high-Mach-number collisionless shock through a long-term, large-scale particle-in-cell simulation. We show that both protons and electrons are accelerated in the shock and that these accelerated particles generate large-amplitude Alfv{e}nic waves in the upstream region of the shock. After the upstream waves have grown sufficiently, the local structure of the collisionless shock becomes substantially similar to that of a quasi-perpendicular shock due to the large transverse magnetic field of the waves. A fraction of protons are accelerated in the shock with a power-law-like energy distribution. The rate of proton injection to the acceleration process is approximately constant, and in the injection process, the phase-trapping mechanism for the protons by the upstream waves can play an important role. The dominant acceleration process is a Fermi-like process through repeated shock crossings of the protons. This process is a `fast process in the sense that the time required for most of the accelerated protons to complete one cycle of the acceleration process is much shorter than the diffusion time. A fraction of the electrons is also accelerated by the same mechanism, and have a power-law-like energy distribution. However, the injection does not enter a steady state during the simulation, which may be related to the intermittent activity of the upstream waves. Upstream of the shock, a fraction of the electrons is pre-accelerated before reaching the shock, which may contribute to steady electron injection at a later time.

Energy Loss of High-Energy Particles in Particle-in-Cell Simulation

When a charged particle moves through a plasma at a speed much higher than the thermal velocity of the plasma, it is subjected to the force of the electrostatic field induced in the plasma by itself and loses its energy. This process is well-known as the stopping power of a plasma. In this paper we show that the same process works in particle-in-cell (PIC) simulations as well and the energy loss rate of fast particles due to this process is mainly determined by the number of plasma electrons contained in the electron skin depth volume. However, since there are generally very few particles in that volume in PIC simulations compared with real plasmas, the energy loss effect can be exaggerated significantly and can affect the results. Therefore, especially for the simulations that investigate the particle acceleration processes, the number of particles used in the simulations should be chosen large enough to avoid this artificial energy loss.

Nonrelativistic collisionless shocks in weakly magnetized electron--ion plasmas: two-dimensional particle-in-cell simulation of perpendicular shock

A two-dimensional particle-in-cell simulation is performed to investigate weakly magnetized perpendicular shocks with a magnetization parameter of 6 x 10^-5, which is equivalent to a high Alfven Mach number M_A of ~130. It is shown that current filaments form in the foot region of the shock due to the ion-beam--Weibel instability (or the ion filamentation instability) and that they generate a strong magnetic field there. In the downstream region, these current filaments also generate a tangled magnetic field that is typically 15 times stronger than the upstream magnetic field. The thermal energies of electrons and ions in the downstream region are not in equipartition and their temperature ratio is T_e / T_i ~ 0.3 - 0.4. Efficient electron acceleration was not observed in our simulation, although a fraction of the ions are accelerated slightly on reflection at the shock. The simulation results agree very well with the Rankine-Hugoniot relations. It is also shown that electrons and ions are heated in the foot region by the Buneman instability (for electrons) and the ion-acoustic instability (for both electrons and ions). However, the growth rate of the Buneman instability is significantly reduced due to the relatively high temperature of the reflected ions. For the same reason, ion-ion streaming instability does not grow in the foot region.

Electrostatic and electromagnetic instabilities associated with electrostatic shocks: two-dimensional particle-in-cell simulation

A two-dimensional electromagnetic particle-in-cell simulation with the realistic ion-to-electron mass ratio of 1836 is carried out to investigate the electrostatic collisionless shocks in relatively high-speed (~3000 km s^-1) plasma flows and also the influence of both electrostatic and electromagnetic instabilities, which can develop around the shocks, on the shock dynamics. It is shown that the electrostatic ion-ion instability can develop in front of the shocks, where the plasma is under counter-streaming condition, with highly oblique wave vectors as was shown previously. The electrostatic potential generated by the electrostatic ion-ion instability propagating obliquely to the shock surface becomes comparable with the shock potential and finally the shock structure is destroyed. It is also shown that in front of the shock the beam-Weibel instability gradually grows as well, consequently suggesting that the magnetic field generated by the beam-Weibel instability becomes important in long-term evolution of the shock and the Weibel-mediated shock forms long after the electrostatic shock vanished. It is also observed that the secondary electrostatic shock forms in the reflected ions in front of the primary electrostatic shock.

Non-relativistic Collisionless Shocks in Unmagnetized Electron-Ion Plasmas

We show that the Weibel-mediated collisionless shocks are driven at non-relativistic propagation speed (0.1c < V < 0.45c) in unmagnetized electron-ion plasmas by performing two-dimensional particle-in-cell simulations. It is shown that the profiles of the number density and the mean velocity in the vicinity of the shock transition region, which are normalized by the respective upstream values, are almost independent of the upstream bulk velocity, i.e., the shock velocity. In particular, the width of the shock transition region is ~100 ion inertial length independent of the shock velocity. For these shocks the energy density of the magnetic field generated by the Weibel-type instability within the shock transition region reaches typically 1-2% of the upstream bulk kinetic energy density. This mechanism probably explains the robust formation of collisionless shocks, for example, driven by young supernova remnants, with no assumption of external magnetic field in the universe.

Relativistic Collisionless Shocks in Unmagnetized Electron-Positron Plasmas

It is shown that collisionless shock waves can be driven in unmagnetized electron-positron plasmas by performing a two-dimensional particle-in-cell simulation. At the shock transition region, strong magnetic fields are generated by a Weibel-like instability. The generated magnetic fields are strong enough to deflect the incoming particles from upstream of the shock at a large angle and provide an effective dissipation mechanism for the shock. The structure of the collisionless shock propagates at an almost constant speed. There is no linear wave corresponding to the shock wave and therefore this can be regarded as a kind of ``instability-driven shock wave. The generated magnetic fields rapidly decay in the downstream region. It is also observed that a fraction of the thermalized particles in the downstream region return upstream through the shock transition region. These particles interact with the upstream incoming particles and cause the generation of charge-separated current filaments in the upstream of the shock as well as the electrostatic beam instability. As a result, electric and magnetic fields are generated even upstream of the shock transition region. No efficient acceleration processes of particles were observed in our simulation.