Turbulent diffusion of streaming cosmic rays in compressible, partially ionised plasma

Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfv\'en speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this paper, we use a large ensemble of magnetohydrodynamical turbulence simulations to quantify how the basic parameters describing ISM turbulence -- the sonic Mach number, $\mathcal{M}$ (plasma compressibility), Alfv\'en Mach number, $\mathcal{M}_{A0}$ (strength of the large-scale field with respect to the turbulence), and ionisation fraction by mass, $\chi$ -- affect the transport of streaming CRs. We show that the large-scale transport of CRs whose small-scale motion consists of streaming along field lines is well described as a combination of streaming along the mean field and superdiffusion both along (parallel to) and across (perpendicular to) it; $\mathcal{M}_{A0}$ drives the level of anisotropy between parallel and perpendicular diffusion and $\chi$ modulates the magnitude of the diffusion coefficients, while in our choice of units, $\mathcal{M}$ is unimportant except in the sub-Alfv\'enic ($\mathcal{M}_{A0} \lesssim 0.5$) regime. Our finding that superdiffusion is ubiquitous potentially explains the apparent discrepancy between CR diffusion coefficients inferred from measurements close to individual sources compared to those measured on larger, Galactic scales. Finally, we present empirical fits for the diffusion coefficients as a function of plasma parameters that may be used as sub-grid recipes for global interstellar medium, galaxy or cosmological simulations.