In situ observations of energetic particles at the Earth's bow-shock attainable by the satellite missions have long created an opinion that electrons are most efficiently accelerated in a quasi-perpendicular shock geometry. However, shocks that deemed responsible for the production of cosmic ray electrons and their radiation from sources such as supernova remnants are much more powerful and larger than the Earth's bow-shock. Their remote observations suggest that electrons are accelerated very efficiently in the quasi-parallel shocks as well. In this paper we investigate the possibility that protons accelerated to high energies create sufficient wave turbulence necessary for the electron preheating and subsequent injection into the diffusive shock acceleration in a quasi-parallel shock geometry. An additional test-particle-electron population, meant as a low-density addition to the electron core-distribution the hybrid simulation operates on, is introduced. Its purpose is to investigate how these electrons are energized by the "hybrid" electromagnetic field. The reduced spatial dimensionality allowed us to dramatically increase the number of macro-ions per numerical cell and achieve the converged results for the velocity distributions of test-electrons. We discuss the electron preheating mechanisms which can make a significant part of thermal electrons accessible to the ion-driven waves observed in hybrid simulations. We find that the precursor wave field supplied by ions has considerable potential to preheat the electrons before they are shocked at the subshock. Our results indicate that a downstream thermal equilibration of the hot test-electrons and protons does not occur. Instead, the resulting electron-to-proton temperature ratio is a decreasing function of the shock Mach number, $M_\mathrm{A},$ tending to a saturation at high $M_\mathrm{A}.$

Cosmic ray (CR) spectra, both measured upon their arrival at the Earth's atmosphere and inferred from the emission in supernova remnants (SNRs), appear to be significantly steeper than the "standard" diffusive shock acceleration (DSA) theory predicts. Although the reconstruction of the primary spectra introduces an additional steepening due to propagation effects, there is a growing consensus in the CR community that these corrections fall short to explain the newest high-precision data. Using 2D hybrid simulations, we investigate a new mechanism that may steepen the spectrum during the acceleration in SNR shocks. Most of the DSA treatments are limited to homogeneous shock environments. To investigate whether inhomogeneity effects can produce the necessary extra steepening, we assume that the magnetic field changes its angle along the shock front. The rationale behind this approach is the strong dependence of the DSA efficiency upon the field angle, $\theta_\mathrm{Bn}$. Our results show that the variation of shock obliquity along its face results in a noticeable steepening of the DSA spectrum. Compared to simulations of quasi-parallel shocks, we observe an increase of the spectral index by $\Delta q=0.1-0.15$. Possible extrapolation of the limited simulation results to more realistic SNR conditions are briefly considered.

Precise measurements of galactic cosmic rays revealed a significant difference between the rigidity spectral indices of protons and helium ions. This finding is a notable contrast to the commonly accepted theoretical prediction that supernova remnant (SNR) shocks accelerate protons and helium ions with the same rigidity alike. Most of the earlier explanations for the "paradox" appealed to SNR environmental factors, such as inhomogeneous $p$/He mixes in the shock upstream medium, variable ionization states of He, or a multi-SNR origin of the observed spectra. The newest observations, however, are in tension with most of them. In this paper, we show by self-consistent hybrid simulations that such special conditions are not vital for the explanation of the cosmic ray rigidity spectra. In particular, our simulations prove that an SNR shock can modify the chemical composition of accelerated cosmic rays by preferentially extracting them from a homogeneous background plasma without additional, largely untestable assumptions. Our results confirm the earlier theoretical predictions of how the efficiency of injection depends on the shock Mach number $M.$ Its increase with the charge-to-mass ratio saturates at a level that grows with $M.$ We have convolved the time-dependent injection rates of protons and helium ions, obtained from the simulations, with a decreasing shock strength over the active life of SNRs. The integrated SNR rigidity spectrum for $p$/He ratio compares well with the AMS-02 and PAMELA data.

High precision spectrometry of galactic cosmic rays (CR) has revealed the lack of our understanding of how different CR elements are extracted from the supernova environments to be further accelerated in their shocks. Comparing the spectra of accelerated particles with different mass to charge ratios is a powerful tool for studying the physics of particle injection into the diffusive shock acceleration (DSA). Recent AMS-02 demonstration of the similarity of He/$p$, C/$p$, and O/$p$ rigidity spectra has provided new evidence that injection is a mass-to-charge dependent process. We performed hybrid simulations of collisionless shocks and analyzed a joint injection of $p$ and He$^{2+}$ in conjunction with upstream waves they generate. By implication, our results equally apply to C and O fully ionized ions, since they have similar mass to charge ratios. By convolving the time-dependent injection rates of $p$ and He, obtained from the simulations, with a decreasing shock strength over the active life of SNRs, we generate the integrated SNR spectra for $p$ and He. These spectra are consistent with the AMS-02 and Pamela data and earlier theoretical predictions. Our interpretation of the elemental anomaly is therefore intrinsic to collisionless shock mechanisms and does not require additional assumptions, such as the contributions from several different SNRs, their inhomogeneous environments or acceleration from grains.