Ergodicity of perpendicular cosmic ray transport
Zentrum für Astronomie und Astrophysik, Technische Universität
Berlin,
Hardenbergstraße 36,
10623
Berlin,
Germany
email:
robert.c.tautz@gmail.com
Received: 26 February 2016
Accepted: 30 April 2016
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 meansquare displacement.
Methods. Based on the testparticle 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 Kolmogorovtype 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 nonergodic. In addition, it is shown how the step length distribution varies with the particle energy. In conclusion, crossfield transport is nonGaussian but still almost diffusive.
Key words: diffusion / turbulence / magnetic fields / cosmic rays / methods: numerical
© ESO, 2016
1. Introduction
Describing the stochastic motion of energetic charged particles due to the interactions with turbulent (electro)magnetic fields has been a longstanding problem in the fields of highenergy astrophysics and (laboratory) plasma physics (see, e.g., Schlickeiser 2002; and Shalchi 2009 for an introduction). A prominent example is that of cosmic rays being scattered in the fluctuating magnetic fields of the Milky Way and the heliosphere, depending on their kinetic energy. On smaller scales, particles interacting with planetary magnetospheres and the prediction of space weather are important, not least for the safety of electronic devices (e.g., Scherer et al. 2005; Bothmer & Daglis 2006).
By neglecting binary Coulomb collisions, a FokkerPlanck and subsequently a diffusionconvection equation can be derived from the Vlasov equation. The diffusion coefficients are then solely determined by the properties of the magnetic fields. Considerable effort has been put into the analytical derivation of the diffusion coefficients, and numerous approaches – both analytical and numerical – have been invoked (e.g., Shalchi 2009; Tautz & Dosch 2013). An important parameter for the classification of transport processes (Bakunin 2015) is the Kubo number (Kubo 1963; Zimbardo et al. 2000; Gustavsson & Mehlig 2011) which, for fluid turbulence, is defined as the relative turbulence strength multiplied by the ratio of the correlation lengths along and across the mean magnetic field R = (δB/B_{0})(ℓ_{∥}/ℓ_{⊥}).
In recent years, increasing evidence has been found that anomalous transport may play an important role in astrophysics (e.g., Zimbardo et al. 2006; Perri & Zimbardo 2008; Klages et al. 2008; Spatschek 2008; Tautz & Shalchi 2010). Accordingly, a theoretical description of the underlying physics and a detailed characterization of the data – obtained numerically and observationally – are required. A general proof that the deflections induced by turbulent magnetic fields indeed cause a diffusive behavior is still elusive. This problem is of interest both for theoretical investigations and for the analysis of data for instance taken in situ by spacecraft. A nondiffusive behavior results in timedependent (“running”) diffusion coefficients, in which case the solution of the diffusion equation becomes considerably more involved and may even be replaced by a fractional differential equation (e.g., Metzler & Klafter 2004; Tautz et al. 2016, and references therein), In addition, nonergodic behavior emphasizes the individuality of the particles, thus requiring more care when drawing conclusions based on a limited data set.
Here, the transport of energetic charged particles – in particular, cosmic rays – will be investigated based on the trajectories of individual particles (see also Metzler & Jeon 2012). Special focus is placed on crossfield transport, the importance of which had often been overlooked. Recently, the role of subdiffusive perpendicular diffusion (e.g., Tautz & Shalchi 2010, but cf. Qin et al. 2002; Xu & Yan 2013) was emphasized in the extraction of diffusion coefficients from intensity profiles (Tautz et al. 2016). Generally, the necessity of anisotropic diffusion is increasingly recognized in the community (e.g., Kissmann 2014; Girichidis et al. 2016).
This article is organized as follows. In Sect. 2, the meansquare displacement and its importance in the context of random walks and turbulent transport is introduced. The numerical simulations that are used to calculate the transport of charged energetic particles in turbulent magnetic fields are briefly described in Sect. 3. The results presented in Sect. 4 comprise several quantitative criteria and distributions that are based on the meansquare displacement. In Sect. 5, results are shown for the ensemble heterogeneity, which is used to characterize the typical behavior of random walks. Section 6 provides a short conclusion.
2. Turbulent transport processes
The mathematical framework developed to treat the random walk (see, e.g., Chandrasekhar 1943, for an introduction) of an ensemble of identical particles is based on the similarity to the diffusive mixing of gases and liquids (Fick 1855). For further applications such as the analytical or numerical treatment of shock acceleration, diffusivity is then often simply assumed.
Owing to the vanishing mean of the particle’s displacement, the second moment is required, which leads to the meansquare displacement (MSD). Based on the concept of ergodicity (Gustavsson & Mehlig 2011), which states that ensemble averages and time averages should be equal, two versions are widely used. The timeaveraged MSD is defined as (1)for each particle individually. It should be noted that the “running” time coordinate during the process is denoted as t in order to distinguish it from the total simulation time or time span of the measurement, T = max(t). In the simulations presented in Sect. 4, the measurement time T is always chosen such that parallel transport – with respect to the mean magnetic field – has already become diffusive. This is important since strictly speaking diffusion applies only in the infinitely long time limit.
For ergodic processes, the timeaveraged MSD for individual particles is equal to the ensembleaveraged MSD as (2)which are both proportional to a power law in t. If the two functions are not equal, the process is generally characterized as weakly nonergodic to distinguish it from strongly nonergodic processes, in which case the phase space of the particles is separated.
As is shown in Sect. 5, the timeaveraged MSD typically has a large variance for (weakly) nonergodic processes, which reflects the heterogeneity of the individual particles. Therefore, another variable can be introduced, which is the ensembleaveraged, timeaveraged MSD defined through (3)To describe the transport of energetic charged particles and to evaluate the time evolution of the probability distribution function, the diffusion coefficients are calculated via κ = ⟨ Δx(t)^{2} ⟩ /(2t) and are inserted in a transport equation such as the Parker (1965) or Roelof (1969) equation. In such cases a diffusive behavior – i.e., α = 1 in Eq. (2)– is implicitly assumed. There are, however, indications that anomalous diffusion may play an important role in many scenarios, including the propagation of solar energetic particles toward Earth (e.g., Ablaßmayer et al. 2016, and references therein).
Much effort has been put into characterizing the process in terms of a waiting time distribution (e.g., Meroz & Sokolov 2015, and references therein). Originally, a continuously distributed waiting time was used to describe the discrete motion charge carriers in amorphous semiconductors. However, the concept is difficult to apply for the turbulent transport of charged particles, which – at least in the absence of electric fields – move at a constant speed. Instead, the distribution of step lengths may be better suited as shown in Sect. 4.2.
3. Monte Carlo simulation
Measurements in the solar winds (see, e.g., Bruno & Carbone 2005, for an overview) have revealed that the Fourier spectrum of the turbulent fields is roughly in agreement with Kolmogorov’s prediction (Kolmogorov 1941). The turbulent magnetic fields therefore can be best modeled in Fourier space using a kappatype power spectrum (Shalchi & Weinhorst 2009) as (4)where q = 3 is the energy range spectral index (cf. Giacalone & Jokipii 1999) for isotropic turbulence. The turbulence bendover scale, ℓ_{0} ≈ 0.03 au, reflects the transition from the energy range G(k) ∝ k^{q} to the Kolmogorovtype inertial range, where G(k) ∝ k^{− s} with s = 5/3.
In the following, the numerical Monte Carlo code Padian will be used to evaluate the properties of perpendicular transport for the turbulence model described above. A general description of the code and the underlying numerical techniques can be found elsewhere (Tautz 2010, see also Michałek & Ostrowski 1996; Giacalone & Jokipii 1999; Laitinen et al. 2013). Specifically, the turbulent magnetic fields are generated via a superposition of N plane waves as (5)where the wavenumbers k_{n} are distributed logarithmically in the interval and where β is a random phase angle. The polarization vector is chosen so that with , respectively, which ensures the solenoidality condition. For isotropic turbulence, the primed coordinates are determined for each n by a rotation matrix with random angles.
From the integration of the equation of motion, the parallel diffusion coefficient, κ_{∥}, and the parallel meanfree path, λ_{∥}, can be calculated by averaging over an ensemble of particles and by determining the MSD in the direction parallel to the background magnetic field as κ_{∥} = (v/ 3)λ_{∥} = ⟨ (Δz)^{2} ⟩ /(2t). For the perpendicular diffusion coefficient, a similar procedure is performed, but the results are also averaged over the x and y directions. It is important to note that realistic and robust results can be obtained only if a number of turbulence realizations (i.e., sets of random numbers) is considered, over which the MSD needs to be averaged. This is discussed again in Sect. 5.
The dynamics of charged particles that interact with magnetic fields are determined through their rigidity pc/q, i.e., the momentum per charge. Therefore, the Padian code uses a normalized rigidity variable R = γv/Ωℓ_{0}, where γ = (1 + v^{2}/c^{2})^{− 1/2} is the relativistic Lorentz factor and Ω = qB/mc is the gyrofrequency with B the strength of the mean magnetic field. In addition, it should be noted that throughout this paper all lengths are normalized to the turbulence bendover scale, i.e., x = x_{phys}/ℓ_{0} and all times are normalized to the gyrofrequency, i.e., t = Ωt_{phys}. Here, the index “phys” denotes physical lengths and times given in centimeters and seconds. Consequently, diffusion coefficients are given as .
Fig. 1 Walk dimension as obtained from the timeaveraged, ensemble averaged MSD defined in Eq. (3). The black solid lines show the simulation results, while the red dotted lines show the power law given in Eq. (6)with the best fit values as summarized in Table 1. 

Open with DEXTER 
The major importance of such simulations is found in the direct calculation of the transport parameters, as opposed to simulations that employ an often simplified model for the diffusion coefficients. In addition, the diffusion coefficients show a nontrivial time dependence, which means that at least three separate phases can be distinguished (see, e.g., Tautz & Shalchi 2011; Prosekin et al. 2015; Ablaßmayer et al. 2016).
4. Results for the MSD
In the following two sections, the results from the Monte Carlo simulation described in Sect. 3 is presented. Unless stated otherwise, three representative values are assumed for the particle rigidity, which are R = 10^{2}, 10^{1}, and 1.
4.1. AlexanderOrbach relation
For the random walk on clusters of fractals, Alexander & Orbach (1982) showed that there is a relation connecting three dimensions relevant for the random walk of particles in an arbitrary domain. A confirmation of this conjecture has proven to be challenging (e.g., Leyvraz & Stanley 1983; Nakayama & Yakubo 2003). For this reason, here its validity will be tested for an entirely different system than originally anticipated.
First, the relevant dimensions include the fractal dimension d_{f}, which in the present case is simply the number of spatial dimensions, i.e., d_{f} = 2 for perpendicular transport. Second, the walk dimension is obtained from the MSD via (6)so that d_{w} = 2 would indicate a diffusive walk. The fit to the simulation results of the perpendicular MSD is shown in Fig. 1 with the fit parameters listed in Table 1. To a good accuracy, the time evolution of the MSD can be described as a power law, not least after the initial transient effects (the socalled ballistic phase in which particles move freely until they are scattered for the first time). We also note that the timeaveraged, ensembleaveraged MSD is identical to the usual definition of the MSD. For low rigidities, we note that the particles have moved at most a few bendover scales, which are proportional to the turbulence correlation length (see, e.g., Shalchi 2009). Formally, however, a diffusive description can be applied only in the limit of large times when the particle has traveled several bendover scales so that its motion is no longer correlated with the initial conditions. Even though such a simulation is computationally very demanding, several tests were performed that did not reveal any later variations in the slope of the MSD for larger T, thus suggesting that a shorter simulation time is sufficient.
Third, the spectral dimension d_{s} is related to the probability density at the origin, i.e., the source, of perpendicular transport, which here is the zaxis. The relevant quantity to be determined from the simulated particle trajectories is thus (7)which can be obtained by recording the number of particles (density) in the form of an intensity profile^{1} (Tautz et al. 2016).
Fig. 2 Spectral dimension as obtained from the decreasing particle flux at the zaxis, which illustrates the turbulent transport in the perpendicular directions. For the three rigidities, the fit parameter values for the power law given in Eq. (7)are summarized in Table 1. 

Open with DEXTER 
The result is shown in Fig. 2. While for ergodic transport and from the diffusion equation a value of d_{s} = 2 is expected, it can be seen that this value is not exhibited for low rigidities. Instead, a reduced slope is observed, which is indicative of a subdiffusive process and reflects the fact that the particles’ residence time near the origin has increased.
The AlexanderOrbach relation connects these three dimensions and states that (8)By comparing Figs. 1 and 2, it can be seen that the simulation results confirm the validity of Eq. (8). Accordingly, two independent measurements confirm that α decreases with increasing rigidity (see Table 1).
Fig. 3 Distribution of steps, ψ(Δx_{⊥}), for a simulation with three rigidity values. The gray areas and the red lines illustrate the distribution of step lengths in the x and y directions, respectively. In addition, the mean (positive) step length is depicted by the blue dotted lines with the values given directly in the figure. 

Open with DEXTER 
The comparison of the walk and spectral dimensions is important to judge the general behavior: if – as is the case here – then the exploration of space is compact, which means that the particle will visit each point in space multiple times for an infinite duration of the random walk. In the direction perpendicular to the mean magnetic field, therefore, particles tend to remain confined instead of diffusively filling space.
4.2. Step lengths distribution
MSD as derived from the step length distribution.
In addition to the waiting times, an equally important distribution is that of the step lengths, hereafter denoted ψ(Δx_{⊥}). This distribution describes the probability for a certain displacement within the constant time step τ. It should be noted that this time step does not correspond to the adaptively refined time step used by the differential equation solver; instead, τ is the arbitrarily chosen spacing of the simulation output.
Fig. 4 Sample results for the timeaveraged MSD as a function of the time lag, t, which is the time variable during the simulation. The dashed black line shows the ensembleaveraged, timeaveraged MSD, which can be fitted to a power law ∝ t^{0.7235} for t ≫ t_{ball}, where t_{ball} is the ballistic time scale. The colored lines show the ensembleaveraged, timeaveraged MSD for individual turbulence realizations. 

Open with DEXTER 
Together with the mean step length (positive) , the step distributions are shown in Fig. 3. The variations in the results for different particle rigidities are caused by the turbulence power spectrum, which introduces an absolute scale – the bendover scale, ℓ_{0} – that has to be related to the particle Larmor radius, R_{L}. For higher rigidities, it is significantly less likely for a particle to become trapped in smallscale structures and, at the same time, the particle motion becomes increasingly random and unpredictable. These findings are reflected in the fact that the distribution ψ is flat over a large range of step lengths, as illustrated in Fig. 3.
From the step length distribution, the diffusion constant can be inferred as (e.g., Klafter et al. 1987) (9)where ⟨ τ ⟩ would normally be the characteristic waiting time. For a continuously moving particle – particularly if its speed is constant – this quantity can be simply the time step at which the displacement is recorded. The results are summarized in Table 2 and confirm the validity of Eq. (9), provided that the maximum of the timedependent diffusion coefficient is taken into account.
5. Ensemble heterogeneity
Fig. 5 Distribution of the MSD randomness variable ξ as defined in Eq. (11). For the three rigidity values, the histograms of the numerical simulation are shown in logarithmic units. We note that there are significantly larger outliers that are omitted here for clarity. The red solid lines show power law fits ∝ ξ^{− β} with β = 2.226 (upper panel), β = 2.404 (middle panel), and β = 2.486 (lower panel). The black dashed line in the upper panel shows the distribution φ(ξ) as defined in Eq. (11). 

Open with DEXTER 
As stated in Sect. 2, ergodic behavior requires an ensemble of individual particles to move in a statistically independent manner in entirely random directions (Metzler et al. 2014). Reversing the argument, ergodicity breaking is caused by particles with a time average that deviates from that of the ensemble. In Fig. 4, the timeaveraged and ensembleaveraged MSDs are shown for charged particles being scattered by a stochastic magnetic field. Notably, the ensembleaveraged MSD agrees with the timeaveraged MSD as defined in Eq. (1)only if the latter is also ensembleaveraged. Therefore, Eq. (2)is not fulfilled for perpendicular cosmicray transport.
Fig. 6 Ergodicity breaking parameter as defined in Eq. (15)for the three rigidity values. The black solid and blue dashed lines (the latter only marginally visible in the bottom panel) show EB(t) for the x and y directions, respectively. For comparison purposes, the red dotdashed line shows EB(t) for the z direction. 

Open with DEXTER 
It should be noted that the various curves in Fig. 4 do not show the timeaveraged MSD for individual particles. Instead, all curves represent averages over different subensembles, each for a given turbulence realization. This result therefore illustrates the strong effect of a particular turbulence realization, which has been pointed out by Mertsch & Funk (2015).
5.1. Heterogeneity distribution
To account for the ensemble heterogeneity in a more quantitative way, a dimensionless variable can be introduced as (10)which is based on the MSDs defined in Eqs. (1)and (3). Its probability distribution has been derived analytically (He et al. 2008) to be (11)where α is the slope of the MSD as given in Eq. (2).
The function g_{α}(x) is defined through a Laplace transformation as (12)for p> 0 and α ∈ (0,1). For rational values α = ℓ/k and , the solution is (Penson & Górska 2010) (13)where the Meijer G function^{2} (Beals & Szmigielski 2013) takes as parameters two special lists of elements defined as Σ(k,a) = a/k,(a + 1) /a,...,(a + k − 1) /k. Other representations of g_{α}(x) and simpler forms for special values of α have been given by Penson & Górska (2010).
The numerically obtained distribution is shown in Fig. 5 together with the histogram for the heterogeneity parameter, ξ, as obtained from the numerical simulations. Clearly, the analytically derived distribution φ(ξ) from Eq. (11)shows no agreement with the numerical result. In particular, it should be noted that for higher rigidities the slope of the MSD, α, increases, which causes the distribution φ(ξ) to be even more narrow. Instead, a simple power law provides a reasonable fit, but lacks theoretical justification.
5.2. Ergodicity breaking
In this section, the nonergodic particle behavior is quantified in terms of the socalled ergodicity breaking parameter. According to Cherstvy et al. (2013), a necessary condition for ergodicity can be formulated by stating that the ratio of time and ensembleaveraged MSD equals unity, i.e., (14)which, as emphasized earlier, is the case for cosmicray transport. The sufficient condition states that, in the limit of large times, the ergodicity breaking parameter (15)should vanish in the limit of large times.
Figure 6 shows that this is not the case. In fact, for perpendicular transport the parameter EB steadily increases instead of decreasing and is on the order of 10 for large simulation times. The corresponding result for the z direction, in contrast, typically remains one order of magnitude smaller. This result illustrates that for parallel cosmicray scattering, the ergodicity is far less broken than for perpendicular transport.
We note that for consistency reasons, some authors prefer to call EB a heterogeneity parameter (e.g., Thiel & Sokolov 2014). A connection of the ergodicity breaking parameter to the diffusivity power index α has been given by He et al. (2008) by requiring that . There is, however, no reason that the variance of ξ should be limited to one, and in fact it is not.
5.3. NonGaussianity
Fig. 7 Results for the nonGaussianity parameter α_{2} as defined in Eq. (16). In each plot, α_{2} is shown as a function of the simulation time, t, for the x and y direction (black solid and blue dashed lines, respectively). For comparison purposes, α_{2} is also shown for the z direction (red dotdashed line). 

Open with DEXTER 
As shown in Sect. 4.2, the distribution of step lengths does not follow a normal distribution. This deviation from Gaussianity can be characterized quantitatively (see Rahman 1964; Meroz & Sokolov 2015) through a nonGaussianity parameter (16)with σ = 3 for one dimension and σ = 2 for two dimensions. For a Gaussian distribution, α_{2} would be zero.
The result for the nonGaussianity parameter α_{2} is shown in Fig. 7, where for comparison the corresponding parameter for parallel transport is also shown. While the turbulent transport of charged particles in stochastic magnetic fields is essentially nonGaussian in all directions, it turns out that parallel transport tends to become regular in the limit of large times, as shown by the approximate limit α_{2} → 0 even for high rigidities.
We note that parts of the results show a deviation between the x and y directions, which to date can be explained only by the relative importance of outliers in the calculation of the fourth moment. The uncertainties on α_{2} are thus large, while the results nevertheless indicate that the perpendicular transport is nonGaussian.
Fig. 8 Timeaveraged MSD as a function of the measurement duration, T. We note that the simulation setup is such that for lower rigidities, the longer measurement time also entails a lower timeresolution. Accordingly, the data seems to begin at later times because the first time step cannot be placed at T = 0. 

Open with DEXTER 
5.4. Aging
If the timeaveraged MSD is taken to be a function of the measurement duration, T, then a variation could indicate the presence of an aging process (Metzler et al. 2014). Such a variation is typically observed for a broad distribution of waiting times τ, since the measurement duration can only include waiting times up to the measurement time, i.e., .
In Fig. 8, the timeaveraged MSD is shown as a function of T. For all rigidities, the results are compatible with a constant value. This result is to be expected because, as mentioned earlier, the concept of a waiting time is not applicable for continuously moving particles.
6. Summary and conclusion
The ergodicity of cosmicray transport in turbulent magnetic fields has been investigated by means of numerical Monte Carlo simulations. In particular, the transport across the mean magnetic field has been studied, which had previously been found to be subdiffusive on many occasions. In view of the controversy if, and under what circumstances, perpendicular transport behaves “normally” merits special attention.
In biology and chemistry, the characterization of particle transport processes in terms of waiting time distributions, Lévy flights, and aging systems has a certain tradition. In highenergy astrophysics, in contrast, a diffusive – and thus ergodic – behavior is often simply assumed. Here it has been shown that crossfield transport is, in general, weakly nonergodic but still almost diffusive, as confirmed through several mathematical criteria. The transport is nonGaussian and obeys the AlexanderOrbach relation, and the underlying system is nonaging as expected. The most prominent deviation from a simple ergodic process is the heterogeneity of the particle ensemble, which is characterized both by the diversity of the MSD for different turbulence realizations and by the heterogeneity parameter. An additional interesting finding is the step length distribution, which is flat (peaked) for particles with high (low) rigidities.
For future work it may prove interesting to relate the nonergodic behavior to the diffusion coefficients. The plethora of criteria that exist in the literature for the classification of random walks may be useful in numerical simulations and in observations. For instance, the validity of the AlexanderOrbach relation allows the intensity decay at a fixed point in time to be connected to the random walk through space. While beyond the scope of the present paper, it will be interesting to see if such relations also hold for more realistic turbulence models involving propagating plasma waves and transient structures.
Unlike in previous work (Tautz et al. 2016; Ablaßmayer et al. 2016), here no kernel function is used because (i) a sufficiently large number of particles will be registered at all times near the zaxis and (ii) potential inaccuracies due to the specific choice of a kernel will be avoided altogether.
Acknowledgments
I thank J. Pratt and the organizers of the XXXV Dynamics Days Europe conference in Exeter for inviting me to an inspiring meeting.
References
 Ablaßmayer, J., Tautz, R. C., & Dresing, N. 2016, Phys. Plasmas, 23, 012901 [NASA ADS] [CrossRef] [Google Scholar]
 Alexander, S., & Orbach, R. 1982, J. Phys., 43, L625 [Google Scholar]
 Bakunin, O. G. 2015, Uspekhi Fiziologicheskikh Nauk, 185, 271 [CrossRef] [Google Scholar]
 Beals, R., & Szmigielski, J. 2013, Notices of the AMS, 60, 866 [NASA ADS] [CrossRef] [Google Scholar]
 Bothmer, V., & Daglis, I. A. 2006, Space Weather: Physics and Effects (Berlin: Springer) [Google Scholar]
 Bruno, R., & Carbone, V. 2005, Liv. Rev. Sol. Phys., 2, 1 [Google Scholar]
 Chandrasekhar, S. 1943, Rev. Mod. Phys., 15, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Cherstvy, A. G., Chechkin, A. V., & Metzler, R. 2013, New J. Phys., 15, 083039 [NASA ADS] [CrossRef] [Google Scholar]
 Fick, A. 1855, Annalen der Physik, 170, 59 [NASA ADS] [CrossRef] [Google Scholar]
 Giacalone, J., & Jokipii, J. R. 1999, ApJ, 520, 204 [NASA ADS] [CrossRef] [Google Scholar]
 Girichidis, P., Naab, T., Walch, S., & Hanasz, M. 2016, MNRAS, submitted [arXiv:1406.4861] [Google Scholar]
 Gustavsson, K., & Mehlig, B. 2011, Europhys. Lett., 96, 60012 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 He, Y., Burov, S., Metzler, R., & Barkai, E. 2008, Phys. Rev. Lett., 101, 058101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Kissmann, R. 2014, Astropart. Phys., 55, 37 [NASA ADS] [CrossRef] [Google Scholar]
 Klafter, J., Blumen, A., & Shlesinger, M. F. 1987, Phys. Rev. A, 35, 3081 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Klages, R., Radons, G., & Sokolov, I. M. 2008, Anomalous Transport (Weinheim: WileyVCH) [Google Scholar]
 Kolmogorov, A. N. 1941, Proc. USSR Acad. Sci., 30, 299 [Google Scholar]
 Kubo, R. 1963, J. Math. Phys., 4, 174 [NASA ADS] [CrossRef] [Google Scholar]
 Laitinen, T., Dalla, S., & Marsh, M. S. 2013, ApJ, 773, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Leyvraz, F., & Stanley, H. E. 1983, Phys. Rev. Lett., 51, 2048 [NASA ADS] [CrossRef] [Google Scholar]
 Meroz, Y., & Sokolov, I. M. 2015, Phys. Rep., 573, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Mertsch, P., & Funk, S. 2015, Phys. Rev. Lett., 114, 021101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Metzler, R., & Jeon, J.H. 2012, Phys. Scr., 86, 058510 [CrossRef] [Google Scholar]
 Metzler, R., & Klafter, J. 2004, J. Phys. A: Math. Gen., 37, R161 [NASA ADS] [CrossRef] [Google Scholar]
 Metzler, R., Jeon, J.H., Cherstvya, A. G., & Barkaid, E. 2014, Phys. Chem. Chem. Phys., 16, 24128 [CrossRef] [PubMed] [Google Scholar]
 Michałek, G., & Ostrowski, M. 1996, Nonlin. Processes Geophys., 3, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Nakayama, T., & Yakubo, K. 2003, Fractal Concepts in Condensed Matter Physics (Berlin: Springer) [Google Scholar]
 Parker, E. N. 1965, Planet. Space Sci., 13, 9 [NASA ADS] [CrossRef] [Google Scholar]
 Penson, K. A., & Górska, K. 2010, Phys. Rev. Lett., 105, 210604 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Perri, S., & Zimbardo, G. 2008, J. Geophys. Res., 113, A03107 [NASA ADS] [Google Scholar]
 Prosekin, A. Y., Kelner, S. R., & Aharonian, F. A. 2015, Phys. Rev. D, 92, 083003 [NASA ADS] [CrossRef] [Google Scholar]
 Qin, G., Matthaeus, W. H., & Bieber, J. W. 2002, ApJ, 578, L117 [NASA ADS] [CrossRef] [Google Scholar]
 Rahman, A. 1964, Phys. Rev., 136, A405 [NASA ADS] [CrossRef] [Google Scholar]
 Roelof, E. C. 1969, in Lectures in high energy astrophysics, eds. H. B. Ögelmann, & J. R. Wayland Jr., NASA SP199, 111 [Google Scholar]
 Scherer, K., Fichtner, H., Heber, B., & Mall, U. 2005, Space Weather: The Physics Behind a Slogan (Berlin: Springer) [Google Scholar]
 Schlickeiser, R. 2002, Cosmic Ray Astrophysics (Berlin: Springer) [Google Scholar]
 Shalchi, A. 2009, Nonlinear Cosmic Ray Diffusion Theories (Berlin: Springer) [Google Scholar]
 Shalchi, A., & Weinhorst, B. 2009, Adv. Space Res., 43, 1429 [NASA ADS] [CrossRef] [Google Scholar]
 Spatschek, K. H. 2008, Plasma Phys. Contr. Fusion, 50, 124027 [NASA ADS] [CrossRef] [Google Scholar]
 Tautz, R. C. 2010, Comp. Phys. Commun., 181, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Tautz, R. C., & Dosch, A. 2013, Phys. Plasmas, 20, 022302 [NASA ADS] [CrossRef] [Google Scholar]
 Tautz, R. C., & Shalchi, A. 2010, J. Geophys. Res., 115, A03104 [NASA ADS] [CrossRef] [Google Scholar]
 Tautz, R. C., & Shalchi, A. 2011, ApJ, 735, 92 [NASA ADS] [CrossRef] [Google Scholar]
 Tautz, R. C., Bolte, J., & Shalchi, A. 2016, A&A, 586, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Thiel, F., & Sokolov, I. M. 2014, Phys. Rev. E, 89, 012136 [NASA ADS] [CrossRef] [Google Scholar]
 Xu, S., & Yan, H. 2013, ApJ, 779, 140 [NASA ADS] [CrossRef] [Google Scholar]
 Zimbardo, G., Pommois, P., & Veltri, P. 2000, Physica A, 280, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Zimbardo, G., Pommois, P., & Veltri, P. 2006, ApJ, 639, L91 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Walk dimension as obtained from the timeaveraged, ensemble averaged MSD defined in Eq. (3). The black solid lines show the simulation results, while the red dotted lines show the power law given in Eq. (6)with the best fit values as summarized in Table 1. 

Open with DEXTER  
In the text 
Fig. 2 Spectral dimension as obtained from the decreasing particle flux at the zaxis, which illustrates the turbulent transport in the perpendicular directions. For the three rigidities, the fit parameter values for the power law given in Eq. (7)are summarized in Table 1. 

Open with DEXTER  
In the text 
Fig. 3 Distribution of steps, ψ(Δx_{⊥}), for a simulation with three rigidity values. The gray areas and the red lines illustrate the distribution of step lengths in the x and y directions, respectively. In addition, the mean (positive) step length is depicted by the blue dotted lines with the values given directly in the figure. 

Open with DEXTER  
In the text 
Fig. 4 Sample results for the timeaveraged MSD as a function of the time lag, t, which is the time variable during the simulation. The dashed black line shows the ensembleaveraged, timeaveraged MSD, which can be fitted to a power law ∝ t^{0.7235} for t ≫ t_{ball}, where t_{ball} is the ballistic time scale. The colored lines show the ensembleaveraged, timeaveraged MSD for individual turbulence realizations. 

Open with DEXTER  
In the text 
Fig. 5 Distribution of the MSD randomness variable ξ as defined in Eq. (11). For the three rigidity values, the histograms of the numerical simulation are shown in logarithmic units. We note that there are significantly larger outliers that are omitted here for clarity. The red solid lines show power law fits ∝ ξ^{− β} with β = 2.226 (upper panel), β = 2.404 (middle panel), and β = 2.486 (lower panel). The black dashed line in the upper panel shows the distribution φ(ξ) as defined in Eq. (11). 

Open with DEXTER  
In the text 
Fig. 6 Ergodicity breaking parameter as defined in Eq. (15)for the three rigidity values. The black solid and blue dashed lines (the latter only marginally visible in the bottom panel) show EB(t) for the x and y directions, respectively. For comparison purposes, the red dotdashed line shows EB(t) for the z direction. 

Open with DEXTER  
In the text 
Fig. 7 Results for the nonGaussianity parameter α_{2} as defined in Eq. (16). In each plot, α_{2} is shown as a function of the simulation time, t, for the x and y direction (black solid and blue dashed lines, respectively). For comparison purposes, α_{2} is also shown for the z direction (red dotdashed line). 

Open with DEXTER  
In the text 
Fig. 8 Timeaveraged MSD as a function of the measurement duration, T. We note that the simulation setup is such that for lower rigidities, the longer measurement time also entails a lower timeresolution. Accordingly, the data seems to begin at later times because the first time step cannot be placed at T = 0. 

Open with DEXTER  
In the text 