Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic (MHD) simulations of the ISM usually incorporate CR dynamics as an advection-diffusion equation for CR energy density, with anisotropic, magnetic field-aligned diffusion with the diffusive flux assumed to obey Fick's law. We compare test-particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time which appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.

Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight correlation to hold across all scales. We use test particle simulations, tracing the propagation of charged particles (protons) through a random magnetic field, to study the cosmic ray distribution at scales comparable to the correlation scale of the turbulent flow in the interstellar medium ($\simeq 100\,{\rm pc}$ in spiral galaxies). In these simulations, we find that there is no spatial correlation between the cosmic ray number density and the magnetic field energy density. In fact, their distributions are approximately statistically independent. We find that low-energy cosmic rays can become trapped between magnetic mirrors, whose location depends more on the structure of the field lines than on the field strength.

The propagation of cosmic rays in turbulent magnetic fields is a diffusive process driven by the scattering of the charged particles by random magnetic fluctuations. Such fields are usually highly intermittent, consisting of intense magnetic filaments and ribbons surrounded by weaker, unstructured fluctuations. Studies of cosmic ray propagation have largely overlooked intermittency, instead relying on Gaussian random magnetic fields. Using test particle simulations, we investigate cosmic ray diffusivity in intermittent, dynamo-generated magnetic fields. The results are compared with those obtained from non-intermittent magnetic fields having identical power spectra. The presence of magnetic intermittency significantly enhances cosmic ray diffusion over a wide range of particle energies. We demonstrate that the results can be interpreted in terms of a correlated random walk.

The propagation of charged particles, including cosmic rays, in a partially ordered magnetic field is characterized by a diffusion tensor whose components depend on the particle's Larmor radius $R_L$ and the degree of order in the magnetic field. Most studies of the particle diffusion presuppose a scale separation between the mean and random magnetic fields (e.g., there being a pronounced minimum in the magnetic power spectrum at intermediate scales). Scale separation is often a good approximation in laboratory plasmas, but not in most astrophysical environments such as the interstellar medium (ISM). Modern simulations of the ISM have numerical resolution of order 1 pc, so the Larmor radius of the cosmic rays that dominate in energy density is at least $10^{6}$ times smaller than the resolved scales. Large-scale simulations of cosmic ray propagation in the ISM thus rely on oversimplified forms of the diffusion tensor. We take the first steps towards a more realistic description of cosmic ray diffusion for such simulations, obtaining direct estimates of the diffusion tensor from test particle simulations in random magnetic fields (with the Larmor radius scale being fully resolved), for a range of particle energies corresponding to $10^{-2}\lesssim R_L/l_c \lesssim 10^{3}$, where $l_c$ is the magnetic correlation length. We obtain explicit expressions for the cosmic ray diffusion tensor for $R_L/l_c \ll 1$, that might be used in a sub-grid model of cosmic ray diffusion. The diffusion coefficients obtained are closely connected with existing transport theories that include the random walk of magnetic lines.

We present a model of the global magnetic field in the barred galaxy NGC 1365 based jointly on the large-scale velocity field of interstellar gas fitted to HI and CO observations of this galaxy and on mean-field dynamo theory. The aim of the paper is to present a detailed quantitative comparison of a galactic dynamo model with independent radio observations. We consider several gas dynamical and nonlinear dynamo models that include plausible variations of parameters that are poorly known. Models of cosmic ray distribution in the galaxy are introduced to produce synthetic radio polarization maps allowing direct comparison with those observed at 3.5cm and 6.2cm. We show that the dynamo model is robust in that the most important magnetic features are controlled by the relatively well established properties of the density distribution and gas velocity field. The optimal agreement between the synthetic polarization maps and observations is obtained when a uniform cosmic ray distribution is adopted. We find some indirect evidence for enhanced turbulence in the regions of strong velocity shear and within 1-2kpc of the galactic centre. We confirm that magnetic stresses can drive an inflow of gas into the inner 1kpc of the galaxy at a rate of a few solar masses per year. The dynamo models are successful to some extent in modelling the large scale regular magnetic field in this galaxy. Our results demonstrate that dynamo models and synthetic polarization maps can provide information about both the gas dynamical models and conditions in the interstellar medium. We demonstrate that the dynamical effects of magnetic fields cannot be everywhere ignored in galaxy modelling.

The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference magnetohydrodynamic code is discussed. It is shown that the standard implementation fails near singular X-points of the magnetic field, which are common if the field is random. A modification to the diffusion model for cosmic rays is described and the resulting telegraph equation (implemented by solving a dynamic equation for the diffusive flux of cosmic rays) is used; it is argued that this modification may better describe the physics of cosmic ray diffusion. The present model reproduces several processes important for the propagation and local confinement of cosmic rays, including spreading perpendicular to the local large-scale magnetic field, controlled by the random-to-total magnetic field ratio, and the balance between cosmic ray pressure and magnetic tension. Cosmic ray diffusion is discussed in the context of a random magnetic field produced by turbulent dynamo action. It is argued that energy equipartition between cosmic rays and other constituents of the interstellar medium do not necessarily imply that cosmic rays play a significant role in the balance of forces.

We present 3.5cm and 6.2cm radio continuum maps in total and polarized intensity of the barred galaxies NGC 1097 and NGC 1365. Both galaxies exhibit radio ridges roughly overlapping with the massive dust lanes in the bar region. The contrast in total intensity across the radio ridges is compatible with compression and shear of an isotropic random magnetic field. The contrast in polarized intensity is significantly smaller than that expected from compression and shearing of the regular magnetic field; this could be the result of decoupling of the regular field from the dense molecular clouds. The regular field in the ridge is probably strong enough to reduce significantly shear in the diffuse gas (to which it is coupled) and hence to reduce magnetic field amplification by shearing. This contributes to the misalignment of the observed field orientation with respect to the velocity vectors of the dense gas. Our observations, for the first time, indicate that magnetic forces can control the flow of the diffuse interstellar gas at kiloparsec scales. The total radio intensity reaches its maximum in the circumnuclear starburst regions, where the equipartition field strength is about 60\mu G, amongst the strongest fields detected in spiral galaxies so far. The regular field in the inner region has a spiral shape with large pitch angle, indicating the action of a dynamo. Magnetic stress leads to mass inflow towards the centre, sufficient to feed the active nucleus in NGC 1097. We detected diffuse X-ray emission, possibly forming a halo of hot gas around NGC 1097.