Shock acceleration is considered one of the most important mechanisms for the acceleration of astrophysical energetic particles. In this work, we calculate the trajectories of a large number of test charged particles accurately in a parallel shock with magnetic turbulence. We investigate the time evolution of the accelerated-particle energy spectrum in the downstream of the shock in order to understand the acceleration mechanism of energetic particles. From simulation results we obtain power-law energy spectra with a bend-over energy, $E_0$, increasing with time. With the particle mean acceleration time and mean momentum change during each cycle of the shock crossing from diffusive shock acceleration model (following Drury), a time-dependent differential equation for the maximum energy, $E_{acc}$, of particles accelerated at the shock, can be approximately obtained. We assume the theoretical bend-over energy as $E_{acc}$. It is found that the bend-over energy from simulations agrees well with the theoretical bend-over energy using the non-linear diffusion theory, NLGCE-F, in contrast to that using the classic quasi-linear theory (QLT).

Using test particle simulations we study electron acceleration at collisionless shocks with a two-component model turbulent magnetic field with slab component including dissipation range. We investigate the importance of shock normal angle $\theta_{Bn}$, magnetic turbulence level $\left(b/B_0\right)^2$, and shock thickness on the acceleration efficiency of electrons. It is shown that at perpendicular shocks the electron acceleration efficiency is enhanced with the decreasing of $\left(b/B_0\right)^2$, and at $\left(b/B_0\right)^2=0.01$ the acceleration becomes significant due to strong drift electric field with long time particles staying near the shock front for shock drift acceleration (SDA). In addition, at parallel shocks the electron acceleration efficiency is increasing with the increasing of $\left(b/B_0\right)^2$, and at $\left(b/B_0\right)^2=10.0$ the acceleration is very strong due to sufficient pitch-angle scattering for first-order Fermi acceleration, as well as due to large local component of magnetic field perpendicular to shock normal angle for SDA. On the other hand, the high perpendicular shock acceleration with $\left(b/B_0\right)^2=0.01$ is stronger than the high parallel shock acceleration with ($\left(b/B_0\right)^2=10.0$), the reason might be the assumption that SDA is more efficient than first-order Fermi acceleration. Furthermore, for oblique shocks, the acceleration efficiency is small no matter the turbulence level is low or high. Moreover, for the effect of shock thickness on electron acceleration at perpendicular shocks, we show that there exists the bend-over thickness, $L_{\text{diff,b}}$. The acceleration efficiency does not change evidently if the shock thickness is much smaller than $L_{\text{diff,b}}$. However, if the shock thickness is much larger than $L_{\text{diff,b}}$, the acceleration efficiency starts to drop abruptly.

Using test particle simulations we study particle acceleration at highly perpendicular ($\theta_{Bn}\geq 75^\circ$) shocks under conditions of modeling magnetic turbulence. We adopt a backward-in-time method to solve the Newton-Lorentz equation using the observed shock parameters for quasi-perpendicular interplanetary shocks, and compare the simulation results with $ACE$/EPAM observations to obtain the injection energy and timescale of particle acceleration. With our modeling and observations we find that a large upstream speed is responsible for efficient particle acceleration. Our results also show that the quasi-perpendicular shocks are capable of accelerating thermal particles to high energies of the order of MeV for both kappa and Maxwellian upstream distributions, which may originate from the fact that in our model the local background magnetic field has a component parallel to the shock normal.