An analytical solution to the the three-dimensional telegraph equation is presented. This equation has recently received some attention but so far the treatment has been one-dimensional. By using the structural similarity to the Klein-Gordon equation, the telegraph equation can be solved in closed form. Illustrative examples are used to discuss the qualitative differences to the diffusion solution. The comparison with a numerical test-particle simulation reveals that some features of an intensity profile can be better explained using the telegraph approach.

Aims. The random walk of energetic charged particles in turbulent magnetic fields is investigated. Special focus is placed on transport across the mean magnetic field, which had been found to be subdiffusive on many occasions. Therefore, a characterization using the concept of ergodicity is attempted by noting the connection to the time evolution of the mean-square displacement. Methods. Based on the test-particle approach, a numerical Monte-Carlo simulation code is used to integrate the equation of motion for particles that are scattered by magnetic turbulence. The turbulent fields are generated by superposing plane waves with a Kolmogorov-type power spectrum. The individual particle trajectories are then used to calculate a variety of statistical quantities. Results. The simulation results clearly demonstrate how the heterogeneity of the particle ensemble causes the system to be weakly non-ergodic. In addition, it is shown how the step length distribution varies with the particle energy. In conclusion, cross-field transport is non-Gaussian but still almost diffusive.

Aims. Numerical test-particle simulations are a reliable and frequently used tool to test analytical transport theories and to predict mean-free paths. The comparison between solutions of the diffusion equation and the particle flux is used to critically judge the applicability of diffusion to the stochastic transport of energetic particles in magnetized turbulence. Methods. A Monte-Carlo simulation code is extended to allow for the generation of intensity profiles as well as anisotropy-time profiles. Due to the relatively low number density of computational particles, a kernel function has to be used to describe the spatial extent of each particle. Results. The obtained intensity profiles are interpreted as solutions of the diffusion equation by inserting the diffusion coefficients that have been directly determined from the mean-square displacements. The comparison shows that the time dependence of the diffusion coefficients needs to be considered, in particular the initial ballistic phase and the often sub-diffusive perpendicular coefficient. Conclusions. It is argued that the perpendicular component of the distribution function is essential if agreement between the diffusion solution and the simulated flux is to be obtained. In addition, time-dependent diffusion can provide a better description than the classic diffusion equation only after the initial ballistic phase.

Context. The diffusive shock acceleration mechanism has been widely accepted as the acceleration mechanism for galactic cosmic rays. While self-consistent hybrid simulations have shown how power-law spectra are produced, detailed information on the interplay of diffusive particle motion and the turbulent electromagnetic fields responsible for repeated shock crossings are still elusive. Aims. The framework of test-particle theory is applied to investigate the effect of diffusive shock acceleration by inspecting the obtained cosmic-ray energy spectra. The resulting energy spectra can be obtained this way from the particle motion and, depending on the prescribed turbulence model, the influence of stochastic acceleration through plasma waves can be studied. Methods. A numerical Monte-Carlo simulation code is extended to include collisionless shock waves. This allows one to trace the trajectories of test particle while they are being accelerated. In addition, the diffusion coefficients can be obtained directly from the particle motion, which allows for a detailed understanding of the acceleration process. Results. The classic result of an energy spectrum with $E^{-2}$ is only reproduced for parallel shocks, while, for all other cases, the energy spectral index is reduced depending on the shock obliqueness. Qualitatively, this can be explained in terms of the diffusion coefficients in the directions that are parallel and perpendicular to the shock front.

The problem of cosmic-ray scattering in the turbulent electromagnetic fields of the interstellar medium and the solar wind is of great importance due to the variety of applications of the resulting diffusion coefficients. Examples are diffusive shock acceleration, cosmic-ray observations, and, in the solar system, the propagation of coronal mass ejections. In recent years, it was found that the simple diffusive motion that had been assumed for decades is often in disagreement both with numerical and observational results. Here, an overview is given of the interaction processes of cosmic rays and turbulent electromagnetic fields. First, the formation of turbulent fields due to plasma instabilities is treated, where especially the non-linear behavior of the resulting unstable wave modes is discussed. Second, the analytical and the numerical side of high-energy particle propagation will be reviewed by presenting non-linear analytical theories and Monte-Carlo simulations. For the example of the solar wind, the impact of anisotropic and dynamical turbulence models will be discussed. In addition, it will be shown how further complications can be treated that arise from the large-scale magnetic field geometry and turbulent electric fields. The transport properties of energetic particles can thus be calculated for current turbulence models so that they withstand a comparison with measurements taken in the solar wind.

We explore the influence of turbulence on the transport of energetic particles by using test-particle simulations. We compute parallel and perpendicular diffusion coefficients for two-component turbulence, isotropic turbulence, a model based on Goldreich-Sridhar scaling, noisy reduced magnetohydrodynamic turbulence, and a noisy slab model. We show that diffusion coefficients have a similar rigidity dependence regardless which turbulence model is used and, thus, we conclude that the influence of turbulence on particle transport is not as strong as originally thought. Only fundamental quantities such as particle rigidity and the Kubo number are relevant. In the current paper we also confirm the unified nonlinear transport theory for noisy slab turbulence. To double-check the validity and accuracy of our numerical results, we use a second test-particle code. We show that both codes provide very similar results confirming the validity of our conclusions.

Understanding turbulent transport of charged particles in magnetized plasmas often requires a model for the description of random variations in the particle's pitch angle. The Fokker-Planck coefficient of pitch-angle scattering, which is used to describe scattering parallel to the mean magnetic field, is therefore of central importance. Whereas quasi-linear theory assumes a homogeneous mean magnetic field, such a condition is often not fulfilled, especially for high-energy particles. Here, a new derivation of the quasi-linear approach is given that is based on the unperturbed orbit found for an adiabatically focused mean magnetic field. The results show that, depending on the ratio of the focusing length and the particle's Larmor radius, the Fokker-Planck coefficient is significantly modified but agrees with the classical expression in the limit of a homogeneous mean magnetic field.

Aims. The process of pitch-angle isotropization is important for many applications ranging from diffusive shock acceleration to large-scale cosmic-ray transport. Here, the basic analytical description is revisited on the basis of recent simulation results. Methods. Both an analytical and a numerical investigation were undertaken of the Fokker-Planck equation for pitch-angle scattering. Additional test-particle simulations obtained with the help of a Monte-Carlo code were used to verify the conclusions. Results. It is shown that the usual definition of the pitch-angle Fokker-Planck coefficient via the mean-square displacement is flawed. The reason can be traced back to the assumption of homogeneity in time which does not hold for pitch-angle scattering. Conclusions. Calculating the mean free path via the Fokker-Planck coefficient has often proven to give an accurate description. For numerical purposes, accordingly, it is the definition that has to be exchanged in favor of the pitch-angle correlation function.

Context. Spacecraft observations have motivated the need for a refined description of the phase-space distribution function. Of particular importance is the pitch-angle diffusion coefficient that occurs in the Fokker-Planck transport equation. Aims. Simulations and analytical test-particle theories are compared to verify the diffusion description of particle transport, which does not allow for non-Markovian behavior. Methods. A Monte-Carlo simulation code was used to trace the trajectories of test particles moving in turbulent magnetic fields. From the ensemble average, the pitch-angle Fokker-Planck coefficient is obtained via the mean square displacement. Results. It is shown that, while excellent agreement with analytical theories can be obtained for slab turbulence, considerable deviations are found for isotropic turbulence. In addition, all Fokker-Planck coefficients tend to zero for high time values.

The maximum magnetic field strength generated by Weibel-type plasma instabilities is estimated for typical conditions in the interstellar medium. The relevant kinetic dispersion relations are evaluated by conducting a parameter study both for Maxwellian and for suprathermal particle distributions showing that micro Gauss magnetic fields can be generated. It is shown that, depending on the streaming velocity and the plasma temperatures, either the longitudinal or a transverse instability will be dominant. In the presence of an ambient magnetic field, the filamentation instability is typically suppressed while the two-stream and the classic Weibel instability are retained.

Context: Typical space plasmas contain spatially and temporally variable turbulent electromagnetic fields. Understanding the transport of energetic particles and the acceleration mechanisms for charged particles is an important goal of today's astroparticle physics. Aims: To understand the acceleration mechanisms at the particle source, subsequent effects have to be known. Therefore, the modification of a particle energy distribution, due to stochastic acceleration, needs to be investigated. Methods: The diffusion in momentum space was investigated by using both a Monte-Carlo simulation code and by analytically solving the momentum-diffusion equation. For simplicity, the turbulence was assumed to consist of one-dimensional Alfven waves. Results: Using both methods, it is shown that, on average, all particles with velocities comparable to the Alfven speeds are accelerated. This influences the energy distribution by significantly increasing the energy spectral index. Conclusions: Because of electromagnetic turbulence, a particle energy spectrum measured at Earth can drastically deviate from its initial spectrum. However, for particles with velocities significantly above the Alfven speed, the effect becomes negligible.

A modified method is presented to generate artificial magnetic turbulence that is used for test-particle simulations. Such turbulent fields are obtained from the superposition of a set of wave modes with random polarizations and random directions of propagation. First, it is shown that the new method simultaneously fulfils requirements of isotropy, equal mean amplitude and variance for all field components, and vanishing divergence. Second, the number of wave modes required for a stochastic particle behavior is investigated by using a Lyapunov approach. For the special case of slab turbulence, it is shown that already for 16 wave modes the particle behavior agrees with that shown for considerably larger numbers of wave modes.

Reproducing measurements of the scattering mean free paths for energetic particles propagating through the solar system has been a major problem in space physics. The pioneering work of Bieber et al. [Astrophys. J. 420, 294 (1994)] provided a theoretical explanation of such observations, which, however, was based on assumptions such as the questionable hypothesis that quasi-linear theory is correct for parallel diffusion. By employing a hybrid plasma-wave/magnetostatic turbulence model, a test-particle code is used to investigate the scattering of energetic particles. The results show excellent agreement with solar wind observations.

Using the conventional approach of superposing plane waves, it is not possible to create a strictly isotropic turbulent magnetic field structure that obeys all physical constraints, which are (i) equal mean of all magnetic field components; (ii) isotropy of the wave vectors; and (iii) vanishing divergence of the magnetic field. Such magnetic fields are widely implemented in test-particle Monte-Carlo simulations, which are used to obtain (i) scattering mean free paths of charged particles; (ii) field line diffusion coefficients. It is shown that, while the turbulent magnetic field strength plays an important role for the results, such does not seem to be the case for a non-zero magnetic field divergence and/or the isotropy of the wave vectors.

Discussion is given of non-linear soliton behavior including coupling between electrostatic and electromagnetic potentials for non-relativistic, weakly relativistic, and fully relativistic plasmas. For plasma distribution functions that are independent of the canonical momenta perpendicular to the soliton spatial structure direction there are, in fact, no soliton behaviors allowed because transverse currents are zero. Dependence on the transverse canonical momenta is necessary. When such is the case, it is shown that the presence or absence of a soliton is intimately connected to the functional form assumed for the particle distribution functions. Except for simple situations, the coupled non-linear equations for the electrostatic and electromagnetic potentials would seem to require numerical solution procedures. Examples are given to illustrate all of these points for non-relativistic, weakly relativistic, and fully relativistic plasmas.

Cosmic Ray transport in curved background magnetic fields is investigated using numerical Monte-Carlo simulation techniques. Special emphasis is laid on the Solar system, where the curvature of the magnetic field can be described in terms of the Parker spiral. Using such geometries, parallel and perpendicular diffusion coefficients have to be re-defined using the arc length of the field lines as the parallel displacement and the distance between field lines as the perpendicular displacement. Furthermore, the turbulent magnetic field is incorporated using a WKB approach for the field strength. Using a test-particle simulation, the diffusion coefficients are then calculated by averaging over a large number of particles starting at the same radial distance from the Sun and over a large number of turbulence realizations, thus enabling one to infer the effects due to the curvature of the magnetic fields and associated drift motions.

The particle distribution function that describes two interpenetrating plasma streams is re-investigated. It is shown how, based on the Maxwell-Boltzmann-J\"uttner distribution function that has been derived almost a century ago, a counterstreaming distribution function can be derived that uses velocity space. Such is necessary for various analytical calculations and numerical simulations that are reliant on velocity coordinates rather than momentum space. The application to the electrostatic two-stream instability illustrates the differences caused by the use of the relativistic distribution function.

In this letter, an initially unmagnetized pair plasma with asymmetric velocity distributions is investigated where any unstable Weibel mode must be isolated, with discrete values for the growth rates and the unstable wavenumbers. For both a non-relativistic distribution with thermal spread and a high-relativistic two-stream distribution it is shown that isolated modes are excited and that, as the asymmetry tends to zero, the growth rate remains finite, as long as the distribution function is not precisely symmetric.