{'fraction_non_alphanumeric': 0.03793269937233663, 'fraction_numerical': 0.015765056573671724, 'mean_word_length': 4.438748715606804, 'pattern_counts': {'":': 0, '<': 19, '<?xml version=': 0, '>': 8, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 11, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Fractal Properties of Anomalous Diffusion in Intermittent Maps

An intermittent nonlinear map generating subdiffusion is investigated. Computer simulations show that the generalized diffusion coefficient of this map has a fractal, discontinuous dependence on control parameters. An amended continuous time random walk theory well approximates the coarse behavior of this quantity in terms of a continuous function. This theory also reproduces a full suppression of the strength of diffusion, which occurs at the dynamical phase transition from one type of diffusive behavior to another. Similarly, the probability density function of this map exhibits a nontrivial fine structure while its coarse functional form is governed by a time fractional diffusion equation. A more detailed understanding of the irregular structure of the generalized diffusion coefficient is provided by an anomalous Taylor-Green-Kubo formula establishing a relation to de Rham-type fractal functions.


I. INTRODUCTION
Since more than a century normal diffusion, characterized by a linear increase in time of the mean square displacement (MSD) of an ensemble of moving particles, provides a paradigmatic example of a stochastic process.
Here the MSD can be expressed by the ensemble average x 2 ∼ t β with an exponent β = 1, where x holds for the position of a particle on the real line at time t ≥ 0. However, exponents β = 1 are also possible yielding two important examples of anomalous diffusion: superdiffusion with β > 1 and subdiffusion with 0 < β < 1. More recently, the importance of anomalous diffusive regimes was realized not only in physics but also in chemistry, biology and economics. These regimes were found in theoretical models and experiments related, among others, to turbulence, amorphous semiconductors, porous media, surface diffusion, glasses, granular matter, reactiondiffusion processes, plasmas and biological cell motility, see Refs. [1,2,3,4,5,6,7,8,9] for reviews. Parallel to this development low dimensional deterministic maps attracted much attention as simple models of anomalous dynamics which can be understood analytically [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. A characteristic feature of the stochastic processes generated by these maps is provided by the probability density functions (PDFs) of the dynamical variables. In case of normal diffusion these PDFs exhibit Gaussian forms, whereas for subdiffusion they yield tails with stretched exponential decay and for superdiffusion Lévy power laws [17].
Several theoretical approaches have been worked out over the past few decades in order to explain anomalous diffusion. Perhaps the most famous one is continuous time random walk (CTRW) theory, pioneered by the work of Montroll, Weiss and Scher [27]. Their stochastic approach was later on adapted to sub-and superdiffu-sive deterministic maps [11,12,13,17,18,19,22,24]. Related maps were originally proposed by Pomeau and Manneville as simple models of intermittency [10]. A classification of their dynamics was provided particularly by Gaspard and Wang combining stochastic with dynamical systems theory [14]. Exploiting the deterministic properties of these maps, anomalous diffusion was further studied in the framework of the thermodynamic formalism [15,20], by means of periodic orbit theory [16,21,25] and by spectral decomposition techniques [23]. Currently aging phenomena [24,26], non-ergodic behavior [28] and infinite invariant measures [29] are in the focus of investigations demonstrating that these simple models provide ongoing inspiration for important new research.
However, all of the above studies focused on specific values of control parameters only for which these maps are to a large extent analytically tractable. That generally the dynamics is more intricate was shown by calculating the parameter-dependent diffusion coefficient of a piecewise linear one-dimensional map, which turned out to be fractal [30,31,32]. Similar behavior was detected in more complicated models like the climbing sine map [33], in experimentally accessible systems like particles bouncing on corrugated surfaces [34] and in models of Josephson junctions [35]. These findings can be explained by the existence of deterministic dynamical correlations that are topologically unstable under parameter variation [30,31,36]; for a review see Part 1 of Ref. [37].
In this paper we show that a fractal parameter dependence of physical quantities is not only typical for lowdimensional periodic deterministic dynamical systems exhibiting normal transport laws but also for anomalous dynamics. As an example we study parameter-dependent subdiffusion in a simple deterministic map, however, our main arguments hold for a broad class of anomalously diffusive systems. A specific feature of our analysis is that it employs a blending of techniques from stochastic theory and the theory of dynamical systems. That way we wish to contribute to the microscopic foundations of a general theory of anomalous deterministic transport, which does not yet seem to be fully developed.
Our paper is organized as follows: The model is introduced in Section II. As a rate of diffusivity we choose a generalized diffusion coefficient (GDC), which generalizes the diffusion constant known from normal diffusion. In Section III we present results from computer simulations showing that the GDC has a nonmonotonic, irregular dependence on control parameters of the map. We give a qualitative explanation of this phenomenon by arguing that the GDC is a self similar-like, fractal function. The GDC is furthermore conjectured to be everywhere discontinuous. In Section IV we briefly review CTRW theory. In Section V we show that an amended version of it correctly describes the coarse functional form of the GDC. Apart from fractal parameter dependencies, the map under consideration exhibits an interesting dynamical transition, which is studied in detail in Section VI. In Section VII we compute the PDF of this map from simulations and analyze it by means of a time fractional diffusion equation. Section VIII finally outlines a deterministic approach for analyzing the fractal GDC, which is based on a Taylor-Green-Kubo (TGK) formula for anomalous diffusion. This theory enables us to relate anomalous diffusion processes to fractal de Rham-type functions. Section IX summarizes our results. For an outline of this work we refer to Ref. [38].

II. THE MODEL
We focus on a subdiffusive map which, restricted onto the unit interval, was introduced by Pomeau and Manneville to describe intermittency [10], where the parameter z ≥ 1 holds for the degree of nonlinearity and a ≥ 1 is a second control parameter. In the following we will omit these two indices for convenience, M z,a ≡ M . The variable x n ∈ R denotes the position of a point particle at discrete time n ∈ N 0 . Translation symmetry, M (x + 1) = M (x) + 1, and reflection symmetry, , complete the definition of the map on the real line. The iterated dynamics of this model is illustrated in Fig. 1 for particular values of the two control parameters: Given some initial condition x 0 , the equations of motion Eq. (1) produce the next value x 1 , x 1 determines x 2 , and so on. A typical trajectory of the system, displayed in Fig. 1 in a cobweb plot, is therefore represented in form of jumps x 0 → x 1 → x 2 → . . . on the real line. For 1 ≤ z < 2 the diffusion process generated by this map is normal, whereas for z ≥ 2 it is anomalous [11,14,17]. This anomaly is due to the existence of marginally stable fixed points located at all integer values of x n around which a particle gets 'trapped' for a long time [15]. However, in regions with larger local derivative, where the particle can 'jump' to different unit intervals, the dynamics becomes more irregular, see Fig.  1. Thus, a typical trajectory of the spatially continued map Eq. (1) consists of long laminar pieces interrupted by short chaotic bursts. Such a behavior is the hallmark of what is called intermittency [10]. An example of a typical intermittent trajectory generated by this map is shown at the bottom of Fig. 1.
In the groundbreaking works of Thomae, Geisel [11], Zumofen and Klafter [17], main attention was paid to the time dependence of the MSD for specific fixed values of the parameter a. In contrast, our aim is to study how a suitably defined GDC behaves under variation of the two control parameters a and z. In Refs. [5,39] an anomalous diffusion coefficient was introduced by where Γ is the gamma function and β is some constant that will be specified in the following. For convenience we define the GDC by For numerical simulations of the spatially extended map Eq. (1) we have typically used an ensemble of 10 6 particles, where the initial conditions are uniformly distributed on the unit interval (0, 1). If not said otherwise, each trajectory was calculated up to 10 4 time steps. For different values of a we have confirmed that the numerically computed power law dependence of the MSD is in full agreement with the CTRW solution [17] We remark that at the transition point between normal and anomalous diffusion a logarithmically corrected dependence of the MSD is obtained, which again is in complete agreement with CTRW theory. Hence, for the rest of this paper β is defined by the theoretically correctly predicted values Eq. (4).

III. NUMERICAL RESULTS
We first fix the parameter a and study the dependence of K as a function of the degree of nonlinearity z. Numerical results are presented in Figs. 2 and 3. K as a function of z appears to be highly irregular exhibiting a lot of structure on fine scales. Particularly Fig. 2 shows how the fine structure as a function of z evolves under variation of a.
In Fig. 3 the calculations were performed both for the normal and for the anomalous diffusive regime. Eq. (5) suggests that at z = 2 the GDC exhibits a local minimum in form of a total suppression of the strength of anomalous diffusion, that is, K(2, a) = 0 for all values of a. Interestingly, this transition from normal to anomalous diffusion appears to be approached continuously by K as a function of z. The reason is that for z = 2 logarithmic terms are still present in the dynamics, but they contribute for transient times only. This peculiar behavior does not only diminish the GDC but also significantly slows down the convergence of the simulation results. The inset of Fig. 3 shows several representative values of K in the transition region, where different symbols correspond to different computation times and different numbers of trajectories for fixed a. Due to the slow convergence, even for the largest computation times the results for K are quantitatively still apart from the CTRW prediction, which holds in the limit of time to infinity. However, qualitatively there is a tendency towards zero at z = 2. This specific problem will be analyzed in full detail in Section V.
From normal diffusive maps it is known that iterations of the critical points of a map, which here are the points of discontinuity at x = 1/2+m, m ∈ Z, play a crucial role in order to understand the complicated parameter dependence of a diffusion process [30,31]. Accordingly, variations of the height of the map h := M ( 1 2 ) = 1 2 + a( 1 2 ) z strongly affect the parameter-dependent GDC, as we will discuss in detail below. Such variations are achieved both by changing z and a. In order to more clearly represent the impact of variations of z only on the GDC, we study K as a function of z for fixed height h. We emphasize that the topology of the orbit of the critical point is still affected by variations of z only, however, it is less sensitive to z than to varying h. Respective simulation results are shown in Fig. 4. As expected, K(z, h) for fixed h is considerably smoother compared to Figs. 2 and 3. Indeed, the local minimum in the transition from normal to anomalous diffusion is now less obscured by a non-trivial fine structure. However, note that very small irregularities forming a self similar-like pattern are still present, as is suggested by Fig. 4 and the inset, which depicts a magnification of a small characteristic region. The figure also shows the convergence of the data with respect to simulation time and number of trajectories from which the formation of the irregularities is clearly seen. Now we focus on the dependence of K on the parameter h with the parameter of nonlinearity z kept fixed. In Fig. 5 the behavior of the GDC as a function of h is shown for different values of z which, as expected, is more non-monotonous than that of K(z). The structure of K(h) for different z looks somewhat similar, however, we could not detect a simple scaling. For z = 1 the map Eq. (1) reduces to a piecewise linear map for which it was shown analytically and numerically that diffusion is normal and that the diffusion coefficient is a fractal function of the slope a [30,31,32,36,37,40], see the upper curve in Fig. 5. For z = 3 and over a larger range of a, K(a) as obtained from simulations is presented in Fig.  6 (a). Magnifications of the initial region, see Fig. 6 (b) and (c), show a self similar-like behavior of K on finer and finer scales indicating a fractal structure.
A qualitative explanation for the fractality of all these curves can be given in terms of turnstile dynamics, which has successfully been applied in order to understand the fractality of the normal diffusion coefficient in piecewise linear maps. Here we outline the basic idea of this approach only, for technical details we refer to Refs. [30,31].
The key ingredient of this method is the orbit generated by iterations of the critical point x c = 1/2 of the map M restricted onto the unit interval, x n = M n (x c ) mod 1. If this orbit is periodic it defines a Markov partition with certain partition parts representing coupling regions (turnstiles) where, in the associated periodically continued map, particles can jump from one unit interval to another one. Specific parameter values define Markov partitions with specific turnstile couplings which, in more physical terms, generate specific sequences of forwardand backward scattering of particles that start nearby the critical point. It turns out that these Markov parti- tions are topologically highly unstable under parameter variation, in the present case both for varying z and a.
In Fig. 6 (a) three parameter values are identified representing Markov partitions that correspond to two local minima and a local maximum of K(a) thus highlighting a specific fine structure in the whole curve. Higher-order Markov partitions can now be constructed (generated by higher-order iterations of the critical point) which are of the same type as the initial three identified in Fig.  6 (a), however, yielding new sets of parameter values. Such higher-order parameter triples are shown in Fig. 6 (b), (c). They identify the same type of fine structure on finer and finer scales. A similar analysis can be performed for other local peaks of K(a) qualitatively explaining the fractal structure of the GDC. More quantitatively, one may wish to calculate some fractal dimension for these curves. However, this is a hard task even if analytical formulas for a diffusion coefficient are available [40] and is left for further studies.
We conclude this section with an unexpected detail of the GDC. For the piecewise linear map at z = 1 the diffusion coefficient as a function of slope of the map is continuous [32]. For z = 3, however, numerical results suggest that the dependence of K on a is discontinuous. In the inset of Fig. 6 we have magnified a small region of K around a = 20 (analogous behavior is found around a = 12, 28, 36, . . .). Note that for a > 20 particles can jump for the first time from the first unit interval to the fourth one. Increasing the numerical precision by increasing the computation time, K(a) near a = 20 gradually approaches a step function, which suggests that a discontinuity develops. Similar discontinuities have already been reported for Lyapunov exponents as functions of a bias in piecewise linear maps [41]. More complicated ones are known to exist for diffusion coefficients of nonlinear maps exhibiting bifurcation scenarios [33] and have recently been highlighted in research on transport in polygonal billiard channels [42]. An explanation why for our model the GDC is discontinuous for z = 3 while it is continuous for z = 1 may be sought in the completely different character of the PDF of the map for z = 3, which develops non-integrable singularities at all marginally stable fixed points. Details of this peculiar behavior will be discussed later, see analysis after Eq. (52). Note that at suspected points of discontinuity of the GDC like a = 20 the critical point is getting mapped right onto these singular regions of the PDF at the very first iteration. Simultaneously, the orbit of the critical point exhibits a transition from forward to backward scattering under variation of a, which leads to a local maximum of the GDC in terms of turnstile dynamics. For z = 3 this change in the microscopic scattering process is drastically amplified by the singular behavior of the PDF probably leading to a discontinuity in the GDC. In contrast, for z = 1 the PDF is not singular but a non-differentiable step function [31], thus the GDC is locally maximal but still continuous. It would be valuable to have a mathematical proof of this conjectured discontinuity phenomenon.
The discontinuity at a = 20 is furthermore a boundary point of the specific fine structure identified in Fig. 6, which was previously discussed in terms of turnstile dynamics. As we have just emphasized, this parameter value is determined by a specific periodic orbit of the critical point of the map. However, as was outlined before, one can now construct infinitely many higher-order iterates of the same type of periodic orbits yielding new sets of parameter values. These parameter values are related to the same type of structure and, hence, identify the same type of discontinuity on finer and finer scales. Furthermore, such parameter values are typically densely distributed on the parameter axis [30,31]. This leads us to conjecture that K(a) shown in Fig. 6 is discontinuous on a dense set of parameter values (probably being of measure zero). Note that the argument carries over to K(z) in Fig. 4, where a careful study of the largest point of irregularity at z ≃ 3.4 reveals the same type of discontinuity.
So far we have focused on the fine structure of the GDC only. The next section proceeds with an understanding of its coarse functional form. As far as we know there is only one analytical result in the literature trying to predict the whole parameter dependence of the GDC of this map [11]. However, as we will show the respective calculation needs to be modified in order to match to computer simulation results. This necessitates to briefly review the whole approach by explaining our corrections.

IV. CONTINUOUS TIME RANDOM WALK THEORY FOR MAPS
The CTRW theory of Montroll, Weiss and Scher [27] has become a standard tool to model diffusion in intermittent maps like Eq. (1) [11,17]. This approach assumes that diffusion can be decomposed into two stochastic processes characterized by waiting times and jumps. Thus one has two sequences of independent identically distributed random variables, namely a sequence of positive random waiting times T 1 , T 2 , T 3 , . . . with PDF w(t), ∞ 0 w(t)dt = 1, and a sequence of random jumps ζ 1 , ζ 2 , ζ 3 , . . . with a PDF λ(x), ∞ −∞ λ(x)dx = 1. For example, if a particle starts at point x = 0 at time t 0 = 0 and makes a jump of length ζ n at time t n = T 1 + T 2 + ... + T n , its position is x = 0 for 0 ≤ t < T 1 = t 1 and x = ζ 1 + ζ 2 + ... + ζ n for t n ≤ t < t n+1 . The probability that at least one jump is performed within the time interval [0, t) is then t 0 dt ′ w(t ′ ) while the probability for no jump during this time interval reads . The master equation for the PDF P (x, t) to find a particle at position x and time t is then It has the following probabilistic meaning: The PDF to find a particle at position x at time t is equal to the PDF to find it at point x ′ at some previous time t ′ multiplied with the transition probability to get from (x ′ , t ′ ) to (x, t) integrated over all possible values of x ′ and t ′ . The second term accounts for the probability to remain at the initial position x = 0. It is easy to check that this equation yields a normalized PDF P (x, t). The most convenient representation of this equation is obtained in terms of the Fourier-Laplace transform of the PDF, where the hat stands for the Fourier transform and the tilde for the Laplace transform. Respectively, this function obeys the Fourier-Laplace transform of Eq. (6), which is called the Montroll-Weiss equation [27], The Laplace transform of the MSD can be readily obtained by differentiating the Fourier-Laplace transform of the PDF, In order to calculate the MSD within this theory, it thus suffices to know λ(x) and w(t) generating the stochastic process. For one-dimensional maps of the type of Eq. (1), using the symmetry of the map the waiting time distribution can be calculated from the approximation where we have introduced the continuous time t ≥ 0. Solving this equation for initial condition x 0 yields We now define that a particle makes a jump when it leaves the unit interval. Note that this definition is different from the one of Ref. [11], which used the interval [− 1 2 , 1 2 ]. That the unit interval is the right choice for calculating the parameter-dependent random walk diffusion coefficient was shown for z = 1 in Refs. [30,31,36,43,44]. From Eq. (11) one can then obtain the time t that a particle spends on the unit interval before making a jump according to the condition x t = 1. The waiting time thus becomes a function of the injection point x 0 , Accordingly, the waiting time PDF w(t) is related to the yet unknown PDF of injection points, Making the assumption that the PDF of injection points is uniform, P in ≃ 1, the waiting time PDF is calculated to for long enough times t, where normalization is used to obtain the prefactor. If jumps between neighboring cells only are taken into account, the jump PDF and its Fourier transform may be assumed as [17] λ(x) = δ(|x| − 1),λ(k) = cos(k) .
For the following calculations it is useful to define Note that for z ≥ 2, γ is identical with β defined in Eq. (4). For γ = 2 the Laplace transform of Eq. (14) reads where Here we are only interested in the long time behavior of the MSD, which corresponds to taking the limit s → 0 in Eqs. (18), (19) with γ being constant. For 0 < γ < 1, taking the lowest order of the expansion we get Γ(1 − γ, bs) ≃ Γ(1 − γ) and Eq. (16) then yields and its inverse Laplace transform is Similarly, for the normal diffusive case γ > 1 we get from Eq. (18) as the small s-asymptotics of the waiting time distributionw Using Eq. (16) and the inverse Laplace transform yields In the remaining case γ = 1, the expansion of the exponential integral is ei(−x) = ln(x) +γ, whereγ is Euler's constant. For small s Eq. (19) reads w(s) ≃ 1 + bs ln(bs) (25) and Eq. (16) gives˜ The inverse Laplace transform of the first term can be obtained by applying Karamata's Tauberian theorem [45], which relates the Laplace transform of a function given in the form f (s) ∼ A Γ(1+γ) , γ = 1 .
In summary, we obtain for the GDC with γ = 1 z−1 . For γ = 1 we get K = 0 because of the logarithmic term in Eq. (27). Now we compare this result with our numerical simulations. As an example, K as a function of a determined by Eq. (28) is shown in Fig. 6 (a), where z = 3. For z > 2, the difference to the result of Ref. [11] is by a factor of two. This yields approximately the right slope for small a, however, it does not reproduce the trivial value of K(4) = 0 when particles cannot escape from the unit interval anymore. Hence, further modifications are necessary.

V. MODIFIED CTRW THEORY
The easiest way to modify standard CTRW theory is by varying either the waiting time or the jump PDF. However, the waiting time distribution is straightforwardly determined by the model as a power law for all intermittent values of z. This suggests to change the jump PDF Eq. (15) only, which in a first attempt one may write as Here the second term reflects the fact that the particle can stay on the unit interval with probability (1 − p). By assuming that the density of particles is uniform on the  For normal diffusion it is well-known that K rw = p(a) provides a random walk approximation for the diffusion coefficient which is asymptotically correct in case of nearest neighbor jumps, p(a) ≪ 1 [30,31,36,43]. For farther than nearest neighbor jumps this approximation is straightforwardly generalized by weighting the integer jump distance squared with the respective probability, based on the escape region, to perform such a jump [44]. Assuming that one knows about the power law time dependence of the MSD Eq. (2) for this map, one can apply the same approximation to anomalous diffusion. In Fig. 7 the result is shown for K(z) at a = 5. As one can see, this simple argument indeed roughly reproduces the coarse functional form of K(z) and yields the exact value at z = 1.
Accordingly, the modification of standard CTRW theory Eq. (29) works for nearest neighbor jumps only. We thus further amend the jump pdf Eq. (29) by introducing a typical jump length l, with Fourier transformλ(k) = p cos(lk)+1−p. This jump length we define either by the actual mean displacement or, under the simplifying assumption of integer jumps, by the coarse-grained displacement where the square brackets hold for the largest integer less than the argument. Here the averages, denoted by curly brackets, are performed starting from a uniform initial ensemble of walkers with respect to the conditional probability density that particles leave the unit interval. The latter ensemble averages are then averaged over time.
Repeating the calculations from the previous section with Eq. (30), we obtain the MSD   Fig. 6 (b). It gives an asymptotically exact approximation of the GDC for small values of a, however, it also holds well for larger parameters. A quantitative comparison moreover indicates (within numerical accuracy) that this combination yields results which exactly correspond to our numerical findings for K at all integer values of h = M ( 1 2 ) > 0, which for z = 3 are equivalent to a = 4, 12, 20, .... This generalizes results obtained previously for the normal diffusion coefficient at z = 1 [31,43].
We now focus on K as a function of z. Fig. 7 shows the modified CTRW approximation Eqs. (32), (34) in comparison with simulation results. In the strongly anomalous regime of large z the approximation reproduces the GDC from simulations very well, however, for smaller values of z there are obvious deviations. Let us first explain the problem for z → 1, the region around z = 2 will be discussed in the following section.
Let us recall that according to Eqs. (2),(4), for z < 2 the map exhibits normal diffusion in the long-time limit. One can further split this regime into 3/2 ≤ z < 2, where the MSD shows transient anomalous diffusion which becomes normal for long times while for 1 ≤ z < 3/2 it represents purely normal dynamics [14]. However, in the latter case it is well-known that the waiting time distribution is exponential [31], whereγ is the escape rate, and not a power law such as the CTRW approximation Eq. (14). In fact, in the limit of z → 1 Eq. (34) leads to K = apl 2 i , whereas the correct random walk result for normal diffusion reads [31,43] Note that for z → 1 and by assuming a uniform density on the unit interval, in case of l 2 ≪ 2 this equation yields K = p, whereas for l 2 ≫ 2 we get K = l 2 2 thus recovering the two simple random walk results mentioned earlier. Eq. (36) is indeed straightforwardly obtained from CTRW theory if one repeats the calculations with the exponential waiting time distribution Eq. (35) andγ = 1, i.e., by assuming that on average a particle leaves the unit cell after one time-discrete jump. We thus conclude that the CTRW result Eq. (34) gives only a reliable approximation for the GDC if γ ≤ 2, whereas for γ > 2 the simple random walk result for normal diffusion Eq. (36) should be used. This is in line with Fig. 7.

VI. SUPPRESSION OF THE GDC AT A DYNAMICAL TRANSITION
Let us now understand the behavior of K(z) around z = 2, where the map exhibits a dynamical transition [15,47] from normal to anomalous diffusion, see Eqs.
(2),(4). Interestingly, this transition is clearly visible in the strength of the GDC, which to our knowledge has not been reported before. We now explain this phenomenon in detail.
As mentioned in Section III, around the transition point there are significant deviations between CTRW theory and the simulation results. These differences are displayed in Fig. 7. At z = 2 the CTRW approximation forms a non-differentiable little wedge with K(2) = 0 at the minimum, whereas the simulations yield K(2) > 0. We recall that by increasing the computation time there is some slow convergence of the simulation data towards the CTRW solution, see the inset of Fig. 3, but that we were not able to achieve quantitative agreement with the CTRW prediction of K(2) = 0.
This very slow convergence of the simulation results as well as the logarithmic term in the MSD at z = 2, see Eq. (33), suggest that the deviations between simulations and CTRW theory are due to logarithmic corrections in the MSD for parameter values around z = 2. Curiously, such terms are not present in Eq. (33) for z = 2 (γ = 1). By refining our CTRW analysis, we will now show that such logarithmic terms indeed exist. However, for z = 2 they hold for large but finite times only, whereas for z = 2 they persist in the limit of infinite time. As is obvious from Eqs. (33), (34), for z = 2 the surviving logarithmic term leads to a full suppression of the GDC, whereas close to the transition point finite-time logarithmic corrections yield a gradual suppression of the strength of diffusion.
We start again from the exact expression for the Laplace transform of the waiting time PDFw(s), Eq. (18). Combining this equation with the jump PDF Eq. (30) and Eqs. (8), (9) we obtain the Laplace transform of the MSD We immediately drop the second term, since its inverse Laplace transform is simply constant. Let us focus now on the limits s → 0 and γ → 1. In case of 0 < γ ≤ 1 the Γ function in Eq. (37) can be expanded to (38) From this equation it follows that as γ approaches 1 from below one may consider not only the first term as in Eq. (20) but also the second one, since both terms are divergent for γ → 1 − . All other terms are nonsingular and can be safely neglected. We further expand the second term of Eq. (38) to Note that here we have obtained a series of logarithmic terms which, as we shall see, provides the logarithmic corrections we are looking for. By imposing the condition that |(1 − γ) ln(bs)| < 1, we may keep only the first two terms of Eq. (39), Finally, we expand the first term of Eq. (38) to whereγ is again Euler's constant. Because we are interested in the small γ limit, we keep only the first term of this expansion. In summary, we obtain for the Laplace transform of the MSD Eq. (37) in the limits of |(1 − γ) ln(bs)| < 1 and γ → 1 − Inversion of the Laplace transform yields as the final result with t cr := be 1/(1−γ) . For t ≫ t cr and 0 < γ < 1 one can drop the second term in the Γ function expansion Eq. (38), and the asymptotic CTRW result Eq. (33) is recovered. Interestingly, t cr diverges when γ → 1 − , and at γ = 1 one thus arrives at the asymptotic t/ ln t dependence of Eq. (33). Analogous corrections are obtained for 1 < γ near the dynamical transition point, γ → 1 + . By applying the same arguments as above we get wheret cr := be 1/(γ−1) . In the long time limit we again recover the CTRW asymptotics Eq. (33) [48]. These analytical findings are in agreement with results from computer simulations. In the inset of Fig. 7, n β / x 2 is shown as a function of ln n close to the transition point at z = 2 for different values of z and a = 3.14. The logarithmic corrections are less obvious if z is sufficiently different from 2, where n β / x 2 quickly converges to 1/K. However, the logarithmic corrections are getting significant when z approaches the dynamical transition value. We have thus identified the precise dynamical origin of the suppression of the GDC at z = 2.
Dynamical transitions are quite ubiquitous in intermittent maps and have been widely discussed in the literature in terms of the time dependence of the MSD. However, it appears that so far no attention has been paid to a possible critical behavior of the associated GDC. Our results lead us to conjecture that suppression and enhancement of the GDC are typical for dynamical transitions in anomalous dynamics altogether [18].

VII. TIME-FRACTIONAL EQUATION FOR SUBDIFFUSION
The great success of the CTRW approach is related to the fact that it not only predicts the power γ but also the form of the coarse grained PDF P (x, t) of displacements [17]. Correspondingly the anomalous diffusion process generated by our model is not described by an ordinary diffusion equation but by a fractional generalization of it. Starting from the Montroll-Weiss equation and making use of the expressions for the jump and waiting time PDFs Eqs. (14), (15), we rewrite Eq. (8) in the long-time and -space asymptotic form with c = Γ(1 − γ) and b = γ/a. For the initial condition P (x, 0) = δ(x) of the PDF we haveP (k, 0) = 1. It is now helpful to recall the definition of the Caputo fractional derivative of a function G ∂ γ G ∂t γ := and its Laplace transform [49,50], Noticing that the left part of Eq. (45) precisely coincides with the Laplace transform of the Caputo derivative of the PDF and turning back to real space, we arrive at the time-fractional diffusion equation with D given by the modified CTRW theory Eqs. (3), (34). Time-fractional diffusion equations of such a form have already been extensively studied by mathematicians [51]. Note that in application to our model another version of such an equation was proposed in Ref. [24], which uses a Riemann-Liouville fractional derivative. It can be easily shown that both forms of time-fractional diffusion equations are equivalent under rather weak assumptions [52]. Yet two other forms of subdiffusive fractional equations (however, not with applications to maps) were proposed in Refs. [53] and [5]. Again, after some recasting these equations yield Eq. (48) [54]. The solution of Eq. (48) is expressed in terms of an M-function of Wright type [49], with ξ := |x|/ √ Dt γ/2 , where we use the representation of the M-function This solution gives exactly the same asymptotics that was obtained in Ref. [17] for small and large values of ξ, where a 1 , a 2 , b 1 are some constants that are given explicitly in Ref. [17].
Similar to the modified CTRW approximation of the parameter-dependent GDC we now find that the analytical PDF Eq. (49) describes only the coarse scale of the PDF obtained from simulations but not its fine structure. This mismatch can be understood related to the fact that our model consists of a periodic lattice defined by specific elementary cells. Any elementary cell develops a characteristic "microscopic" PDF, characterized by the map Eq. (1) restricted onto the unit interval, which exhibits singularities at the marginally stable fixed points. The inset of Fig. 8 schematically depicts such a PDF of an elementary cell after a large number of iterations. Precisely this functional form yields the building block for the PDF of the periodically continued map, see the main part of Fig. 8. If one eliminates the fine structure by averaging over whole unit intervals, one obtains a coarse-grained PDF that is in excellent agreement with the analytical solution of the fractional diffusion equation Eq. (48), see Fig. 8. Such an interplay between the "microscopic" PDF of a single scatterer and the "macroscopic" PDF of the spatially extended system was already reported for the normal diffusive model at z = 1 [31,37,44] and was recently also studied in billiards [55]. However, in case of anomalous diffusion the PDF exhibits a further remarkable property: Fig. 8 shows that the oscillatory structure of the whole PDF is bounded by two different functions. The upper curve is of M-function type and connects all local maxima of the microscopic PDF. These maxima are situated in regions of the map where the dynamics is regular, due to the marginally stable fixed points. The lower curve, on the other hand, is Gaussian and is determined by all local minima of the microscopic PDF. These minima are generated in regions of the map being far away from the marginally stable fixed points, where the dynamics is locally most strongly chaotic [56]. Fig. 8 thus nicely exemplifies the microscopic origin of anomalous diffusion in terms of intermittency.

VIII. TAYLOR-GREEN-KUBO APPROACH AND FRACTAL FUNCTIONS
We now turn back to the parameter dependence of the GDC. In Section V we have shown that an amended CTRW theory correctly describes the coarse dependence of the model's GDC, whereas the fractal fine structure is not captured by this approach. This reflects the fact that CTRW theory is a purely stochastic approach involving a randomness assumption between jumps, see the approximations involved as outlined in Section IV and V, whereas the origin of the fractality of the GDC lies in the existence of long-range deterministic dynamical correlations, see the analysis in Section III.
This motivates us to propose an alternative approach for analyzing the GDC, which is based on the Taylor-Green-Kubo (TGK) formula [57,58] for diffusion in maps. This theory has successfully been applied to the fractal diffusion coefficient of the normal diffusive map at z = 1, where it could partly be worked out analytically [31], to a nonlinear map with a complicated bifurcation scenario [33] and to normal diffusion in billiards [36]. An advantage of this analysis is that it systematically incorporates dynamical correlations. A disadvantage is that, in contrast to CTRW theory, for the map under consideration we could implement it only numerically.
The basic idea is to generalize the TGK formula for maps to anomalous diffusion. Following the usual derivation [36,57], the first step is to transform the MSD in Eq. (2) into sums over increments, or velocities, v k := For normal diffusion the second step requires to define the ensemble average expressed by the angular brackets in terms of the invariant PDF obtained for the map restricted onto the unit interval. In case of the map Eq. (1) this still works for z < 2, where the dynamics of the spatially extended map is normal diffusive. However, for z ≥ 2 it is well-known that the map of the elementary cell does not possess a non-trivial invariant PDF anymore. Mathematical analysis [29] shows that in this case there exist two physically relevant invariant measures, one which is concentrated on the marginally stable fixed points and one that lies in-between on the interval (0, 1). The first measure has still nice, so-called SRB properties and yields a PDF in form of a δ-distribution on the marginally stable fixed points. The second one, however, is not normalizable anymore and hence is called an infinite invariant measure. In other words, if one starts a computer simulation from an ensemble of points uniformly distributed on the unit interval, the underlying stochastic process is not stationary and the PDF of an elementary cell, computed by a histogram method, does not converge in time to a well-defined invariant PDF [14,18].
Consequently, in contrast to normal diffusion Eq. (52) cannot further be simplified by using time-translational invariance, and for the GDC Eq. (3) we have to stop at In numerical simulations we find that the first term alone is already proportional to n β . If the system is ergodic, as for z < 2 in our model, this term boils down to n < v 2 0 >, and by neglecting any higher-order terms in Eq. (53) we recover the random walk result for normal diffusion [31,36,43]. For z ≥ 2 the situation is again more complicated. Here only generalized ergodic theorems may hold [28,29], which is intimately related to the existence of infinite invariant measures as outlined above. Whether in this case a CTRW result such as Eq. (34) can be extracted from the first term is a non-trivial open question.
However, considering the first term only as an approximation of K for all z, the numerical results are depicted in Fig. 9. The comparison of this approximation with the previous simulation results shows that this term provides a first little step beyond a modified CTRW approach, since it reproduces the major irregularities of K as a function of a and even follows the suspected discontinuities discussed in Section III. However, our numerical precision is not sufficient to conclude whether it yields exact values for K at integer heights.
Eq. (53) thus provides a suitable starting point for a systematic evaluation of the fractal structure of the GDC. Since the series expansion in Eq. (53) is exact, working out further terms one should recover more and more structure in the GDC on fine scales [33,36]. Interestingly, the second term in Eq. (53) can be understood as a velocity autocorrelation function that, in contrast normal diffusion, depends on two time scales. The existence of a second time scale points to the phenomenon of aging in dynamical systems [24,26]. In our case this can be understood as a consequence of infinite invariant measures.
Here we do not further pursue these questions but focus instead onto a direct link between the GDC and fractal functions as provided by the TGK formula, see Refs. [31,33,36] for normal diffusion. As was shown in Ref. [36], Eq. (53) also holds if the velocities v k are replaced by the integer velocities j k := [x k+1 ] − [x k ]. It is now useful to start from K expressed again by Eq. (52), trivially rewritten as The two sums suggest to define the jump function which satisfies the recursion relation [31] J n (x) = j 0 (x) + J n−1 (M (x)) .
The jump function gives the integer value of the displacement of a particle after n time steps, which started at some initial position x ≡ x 0 . It thus contains essential information about the microscopic scattering process of a particle by showing how sensitively the displacement depends on initial conditions. Eqs. (54), (55) imply Fig. 10 displays the product of jump functions J 2 n (x), which governs the anomalous diffusion process, for a representative parameter value. Clearly, as time evolves the structure of this function is getting more complicated. The GDC, in turn, is determined by the cumulative function where we fix the integration constants by the condition that T 2 n (0) = 0 and by requiring that the whole function is continuous. Integration of Eq. (56) would yield a recursive functional equation for T 2 n (x), which is of the same type as the one derived in Ref. [33]. The solutions of such equations are generalized de Rham-, respectively generalized Takagi functions [31,33]. However, at present such generalized de Rham-equations cannot be solved analytically, hence we do not work out the details but instead compute T 2 n (x) directly from simulations, according to the definition Eqs. (55), (58). Results for the first four iterations of T 2 n (x) are shown in Fig. 11. With respect to the construction of T 2 n (x), it may not come too much as a surprise that in the limit of n → ∞ this function exhibits a fractal structure. Interestingly, the GDC is obtained by integrating over this structure. If one starts from a uniform ensemble of initial conditions the result reads K = lim n→∞ 1 n β T 2 n (1) .
This equation relates the GDC to a fractal function representing the microscopic scattering process of our model. The numerical result for the GDC Eq. (59) is in good agreement with the one employing the MSD. In the light of Fig. 11 and Eq. (59), it may not be too surprising anymore that under parameter variation a highly non-trivial GDC is obtained from this model.

IX. CONCLUSIONS
In this paper we have studied subdiffusion generated by a paradigmatic one-dimensional intermittent map. In contrast to previous works here we focused onto the parameter dependence of the GDC, which is an anomalously diffusive generalization of the diffusion coefficient known for normal diffusion. Computer simulations suggested that this GDC is a fractal function of control parameters. This finding was corroborated by a qualitative explanation of the fractal structure in terms of complicated sequences of forward-and backward scattering (turnstile dynamics), which are topologically unstable under parameter variation. Our analysis furthermore led us to conjecture that the GDC is a discontinuous function of control parameters.
In trying to understand the coarse functional form of the GDC, we applied standard CTRW theory to our model. By suitably amending previous calculations, we arrived at analytical approximations that enabled us to reproduce the whole coarse functional form of the GDC yielding asymptotically exact results in the limit of large and small parameter values. However, there are clear deviations between this theory and simulations around a dynamical transition from normal to anomalous diffusion. By refining our amended theory, we were able to explain these deviations in terms of logarithmic corrections leading to ultra-slow convergence of our simulation results and eventually yielding a full suppression of the GDC right at the transition point. These findings were confirmed by simulations of the MSD.
We then studied in detail the PDF of our model. We first derived a time-fractional subdiffusion equation from CTRW theory. The coarse-grained PDF of our model turned out to be in excellent agreement with the non-Gaussian solution of our fractional diffusion equation. On a fine scale, however, the simulation results yielded an oscillatory structure reflecting the microscopic details of the intermittent scattering process. This structure was generated by the invariant density of the single scatterers hence revealing an interesting interplay between microscopic scattering and macroscopic diffusion.
A more detailed understanding of the GDC was finally provided in terms of an anomalous TGK formula, which again is a generalization of the TGK formula for normal diffusion. The structure of this formula points at intimate relations to infinite invariant measures, aging, and generalized de Rham-type fractal functions. All these terms define very active topics of research. We thus hope that, along these lines, in future work more detailed relations between mathematical infinite ergodic theory and the physics of anomalous diffusive processes can be established.
In the present paper we have only studied a onedimensional subdiffusive map, however, we expect these findings to be typical for spatially extended, lowdimensional, anomalous deterministic dynamical systems altogether. Further studies may also focus on the relation between the anomalous TGK formula and CTRW theory, on a spectral analysis of the Frobenius-Perron operator of this model, and on analyzing a superdiffusive map by similar methods.