2016ApJ...830...19L

Cosmic-ray anisotropy has been observed in a wide energy range and at different angular scales by a variety of experiments over the past decade. However, no comprehensive or satisfactory explanation has been put forth to date. The arrival distribution of cosmic rays at Earth is the convolution of the distribution of their sources and of the effects of geometry and properties of the magnetic field through which particles propagate. It is generally believed that the anisotropy topology at the largest angular scale is adiabatically shaped by diffusion in the structured interstellar magnetic field. On the contrary, the medium- and small-scale angular structure could be an effect of nondiffusive propagation of cosmic rays in perturbed magnetic fields. In particular, a possible explanation for the observed small-scale anisotropy observed at the TeV energy scale may be the effect of particle propagation in turbulent magnetized plasmas. We perform numerical integration of test particle trajectories in low-β compressible magnetohydrodynamic turbulence to study how the cosmic rays' arrival direction distribution is perturbed when they stream along the local turbulent magnetic field. We utilize Liouville's theorem for obtaining the anisotropy at Earth and provide the theoretical framework for the application of the theorem in the specific case of cosmic-ray arrival distribution. In this work, we discuss the effects on the anisotropy arising from propagation in this inhomogeneous and turbulent interstellar magnetic field.