Thomas M Michelitsch

Alejandro P Riascos

Continuous time random walk and diffusion with generalized fractional Poisson process

Keywords: Renewal process, fractional Poisson process and distribution, fractional Kolmogorov-Feller equation, continuous time random walk, generalized fractional diffusion

A non-Markovian counting process, the 'generalized fractional Poisson process' (GFPP) introduced by Cahoy and Polito in 2013 is analyzed. The GFPP contains two index parameters 0 < β ≤ 1, α > 0 and a time scale parameter. Generalizations to Laskin's fractional Poisson distribution and to the fractional Kolmogorov-Feller equation are derived. We develop a continuous time random walk subordinated to a GFPP in the infinite integer lattice Z d . For this stochastic motion, we deduce a 'generalized fractional diffusion equation'. In a well-scaled diffusion limit this motion is governed by the same type of fractional diffusion equation as with the fractional Poisson process exhibiting subdiffusive t β -power law for the mean-square displacement. In the special cases α = 1 with 0 < β < 1 the equations of the Laskin fractional Poisson process and for α = 1 with β = 1 the classical equations of the standard Poisson process are recovered. The remarkably rich dynamics introduced by the GFPP opens a wide field of applications in anomalous transport and in the dynamics of complex systems.

Introduction

In many complex systems, one has empirically observed asymptotic power laws instead of the classical exponential patterns. One finds such characteristic behavior for instance in anomalous transport and diffusion [START_REF] Zaslavsky | Chaos, fractional kinetics, and anomalous transport[END_REF][START_REF] Shlesinger | Origins and applications of the Montroll-Weiss continuous time random walk[END_REF][START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Gorenflo | From Power Laws to Fractional Diffusion: the Direct Way[END_REF][START_REF] Metzler | The Random Walk's Guide to Anomalous Diffusion : A Fractional Dynamics Approach[END_REF][START_REF] Metzler | The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics[END_REF][START_REF] Barkai | From continuous time random walks to the fractional Fokker-Planck equation[END_REF]. In the description of these phenomena the celebrated Montroll-Weiss continuous time random walk (CTRW) model [START_REF] Shlesinger | Origins and applications of the Montroll-Weiss continuous time random walk[END_REF][START_REF] Montroll | Random walks on lattices II[END_REF][START_REF] Scher | Stochastic Transport in a Disordered Solid. I. Theory[END_REF][START_REF] Kutner | The continuous time random walk, still trendy: fifty-year history, state of art and outlook[END_REF] has become a major approach where the empirical power laws lead to descriptions that involve fractional calculus.

Meanwhile, the CTRW of Montroll and Weiss has been applied to a vast number of problems in various contexts especially to model anomalous diffusion and transport phenomena [START_REF] Shlesinger | Origins and applications of the Montroll-Weiss continuous time random walk[END_REF][START_REF] Kutner | The continuous time random walk, still trendy: fifty-year history, state of art and outlook[END_REF][START_REF] Klafter | Beyond Brownian Motion[END_REF] including stochastic processes in protein folding [START_REF] Bryngelson | Intermediates and barrier crossing in a random energy model (with applications to protein folding)[END_REF], the dynamics of chemical reactions [START_REF] Sung | Fractional dynamics approach to diffusion-assisted reactions in disordered media[END_REF], transport phenomena in low dimensional chaotic systems [START_REF] Solomon | Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow[END_REF], 'Aging Continuous Time Random Walks' [START_REF] Monthus | Models of traps and glass phenomenology[END_REF][START_REF] Barkai | Aging Continuous Time Random Walks[END_REF], CTRW with correlated memory kernels have been introduced [START_REF] Chechkin | Continuous-time random walk with correlated waiting times[END_REF], and many further applications of the CTRW exist. In the classical picture, CTRWs are random walks subordinated to a Poisson renewal process with exponentially decaying waiting time densities and the Markovian property [START_REF] Metzler | The Random Walk's Guide to Anomalous Diffusion : A Fractional Dynamics Approach[END_REF][START_REF] Metzler | The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics[END_REF][START_REF] Feller | An introduction to probability theory and its applications[END_REF]. However, one has found that many complex systems are rather governed by non-Markovian fat-tailed waiting time densities not compatible with the classical Poisson process picture [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Metzler | The Random Walk's Guide to Anomalous Diffusion : A Fractional Dynamics Approach[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF][START_REF] Mainardi | A fractional generalization of the Poisson processes[END_REF].

Based on a renewal process to our present knowledge first introduced by Repin and Saichev [START_REF] Repin | Fractional Poisson law[END_REF], Laskin developed a non-Markovian fractional generalization of the Markovian Poisson process. He called this process 'fractional Poisson process' [START_REF] Laskin | Fractional Poisson process[END_REF] and demonstrated various important applications [START_REF] Laskin | Some applications of the fractional Poisson probability distribution[END_REF]. The fractional Poisson process is of utmost importance as it produces power-law patterns as observed in anomalous subdiffusion due to fat-tailed waiting time densities. Meanwhile, a large literature exists that analyze fractional Poisson processes and related random motions [START_REF] Gorenflo | From Power Laws to Fractional Diffusion: the Direct Way[END_REF][START_REF] Metzler | The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics[END_REF][START_REF] Beghin | Fractional Poisson processes and related random motions[END_REF][START_REF] Gorenflo | The asymptotic universality of the Mittag-Leffler waiting time law in continuous time random walks[END_REF][START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF][START_REF] Meerschaert | The Fractional Poisson Process and the Inverse Stable Subordinator[END_REF].

In the present paper, our aim is the development of a stochastic motion based on a generalization of the Laskin's fractional Poisson process. This process was introduced by Cahoy and Polito in 2013 [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF] and several features were analyzed in that paper. In the present paper, our goal is to complement that analysis and to investigate related stochastic motions. We call this renewal process 'generalized fractional Poisson process' (GFPP). Both the GFPP like the fractional Poisson process are non-Markovian and introduce long-memory effects. We will come back to this issue subsequently.

The present paper is organized as follows. In the first part (Section 2) we briefly outline general concepts of renewal theory and the CTRW approach by Montroll and Weiss [START_REF] Montroll | Random walks on lattices II[END_REF][START_REF] Scher | Stochastic Transport in a Disordered Solid. I. Theory[END_REF]. In Section 3 we briefly discuss essential features of Laskin's fractional Poisson counting process and of the fractional Poisson distribution. In Section 4 we generalize the Laskin process and define the 'generalized fractional Poisson process' (GFPP). We derive its waiting time density and related distributions. In Section 5 the probability distribution of n arrivals in a GFPP is deduced. We refer this distribution to as 'generalized fractional Poisson distribution' (GFPD).

We derive in Section 6 the expected number of arrivals in a GFPP and analyze its asymptotic power law features. The GFPP contains two index parameters 0 < β ≤ 1 and α > 0 and one additional parameter (ξ -1 β ) that defines a characteristic time scale of the process. The Laskin fractional Poisson process is a special case for α = 1 and 0 < β < 1 of the GFPP, as well as the standard Poisson process (recovered for α = 1 and β = 1). Further the case β = 1 and α > 0 produces an Erlang type process.

As an application, in Sections 7-9, we consider a random walk subordinated to a GFPP on undirected networks and multidimensional infinite lattices in the framework of Montroll-Weiss CTRW. We derive 'generalized fractional Kolmogorov-Feller equations' and for the resulting diffusion process a 'generalized fractional diffusion equation'. We show explicitly that for α = 1 (0 < β < 1) these equations reproduce their fractional counterparts [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF] of the Laskin fractional Poisson process and for α = 1, β = 1 the classical equations of standard Poisson. Finally we analyze the 'well-scaled' diffusion limit and obtain the same type of fractional diffusion equation as in walks subordinated to Laskin's fractional Poisson process.

As a major analytical tool in this paper, we utilize Laplace transforms where we deal with causal non-negative probability distributions (non-negative probability measures) in the distributional sense of Gelfand and Shilov generalized functions [START_REF] Gelfand | Generalized Functions[END_REF]. The Laplace transform f (s) of a causal function f (t) = Θ(t) f (t) is defined by, e.g. [START_REF] Gelfand | Generalized Functions[END_REF] 1

L{f (t)} =: f (s) = ∞ 0 - e -st Θ(t) f (t)dt, s = σ + iω, σ > σ 0 , (1)
where Θ(t) indicates the Heaviside step function defined by Θ(t) = 1 for t ≥ 0 and Θ(t) = 0 for t < 0, especially Θ(0) = 1. However, we emphasize that Θ(0 -) = lim →0+ Θ(-) = 0 ( > 0) is vanishing as the lower integration limit since 0 -= -( → 0+) is infinitesimally negative. The essential point is that the Laplace transform captures the entire nonzero contributions of a causal distribution f (t). We utilize throughout this paper definition [START_REF] Zaslavsky | Chaos, fractional kinetics, and anomalous transport[END_REF] for Laplace transforms since it especially convenient when dealing with causal probability measures and causal operators. The Laplace inversion is performed by

f (t) = f (t)Θ(t) = L -1 { f (s)} = 1 2πi σ+i∞ σ-i∞ e st f (s)ds, (s) = σ > σ 0 , = e σt 2π ∞ -∞ e iωt f (σ + iω)dω (2)
with suitably chosen (s) = σ. We only mention here the important property for causal functions

L d m dt m (Θ(t) f (t)) = s m f (s) - m-1 k=0 s m-1-k d k dt k (f (t)Θ(t)) t=0 - =0 L -1 {s m f (s)} = d m dt m Θ(t) f (t) , m ∈ N (3)
i.e. Laplace inversion of s m f (s) captures in the derivatives the Heaviside step function2 . This allows covering causal distributions that exhibit discontinuities especially when they are 'switched on' at t = 0. We emphasize property (3) as it allows us later to determine the appropriate fractional derivative kernels (See especially Appendix A). In our demonstration when no derivatives are involved we often skip the Heaviside Θ(t) step function in expressions of causal distributions. However, for clarity, we include the notation "Θ(t)" in all cases such as (3) where time derivatives are involved. We mention further with the above definition of Laplace transform [START_REF] Zaslavsky | Chaos, fractional kinetics, and anomalous transport[END_REF], that the Dirac's δ-function is entirely captured. This is expressed by L{δ(t)} = 1. Further properties are discussed in Appendix A. We avoided in the present paper too extensive mathematical derivations which can be found in a complementary paper [START_REF] Michelitsch | Generalized fractional Poisson process and related stochastic dynamics[END_REF].

Contrary to the version of the present paper to appear in Physica A, in the present version we have analyzed in Section 9 a 'well-scaled' diffusion limit which yields the time fractional diffusion equation (87.B).

Renewal process and CTRW

Let us briefly outline some basic features of renewal processes and waiting time distributions (See e.g. [START_REF] Gorenflo | From Power Laws to Fractional Diffusion: the Direct Way[END_REF]) as a basis of the CTRW approach of Montroll and Weiss [START_REF] Montroll | Random walks on lattices II[END_REF]. We assume that a walker makes jumps at random times 0 ≤ t 1 , t 2 , . . . t n , . . . ∞. The non-negative random times when jumps occur we also refer to as 'arrival times' t k . We start the observation of the walk at time t = 0 where the walker is sitting on his initial position until making the first jump at t = t 1 (arrival of the first jump event) and so forth. The waiting times ∆t k = t kt k-1 ≥ 0 between two successive (jump-) events are assumed to be drawn for each jump from the same probability density distribution (PDF) χ(t) which we refer to as 'waiting time density' or 'jump density'. The waiting times ∆t k then are independent and identically distributed (IID) random variables. χ(τ )dτ indicates the probability that the first jump event (first arrival) happens exactly at time τ . The probability that the walker at least has jumped once within the time interval [0, t] is given by the cumulative probability distribution

prob(∆t ≤ t) = Ψ(t) = t 0 χ(τ )dτ, t ≥ 0, lim t→∞ Ψ(t) = 1 -0 (4)
with initial condition Ψ(t = 0) = 0. The probability that the walker within [0, t] still is waiting on its departure site is given by

prob(∆t ≥ t) = Φ (0 (t) = 1 -Ψ(t) = ∞ t χ(τ )dτ. ( 5
)
The cumulative probability distribution (5) that the walker during [0, t] still is waiting on its initial site also is called 'survival probability' and indicates the probability that the waiting time ∆t ≥ t. We notice that the jump densities3 such as χ(t) have dimension [time] -1 whereas the cumulative distributions (4), ( 5) are dimensionless probabilities. In addition, the waiting time density χ(τ ) is normalized (see ( 4)) which is reflected by the property χ(s)| s=0 = 1 of its Laplace transform. It follows that the survival probability (5) tends to zero as time approaches infinity lim t→∞ Φ (0 (t) → 0. The PDF χ (n) (t) that the walker makes its nth jump exactly at time t then is4

χ (n) (t) = ∞ 0 . . . ∞ 0 χ(τ 1 ) . . . χ(τ n )δ t - n j=1 τ j dτ 1 . . . dτ n (6)
having the Laplace transform

χ(n) (s) = ∞ 0 dτ 1 . . . ∞ 0 dτ n χ(τ 1 ) . . . χ(τ n ) ∞ 0 - dte -st δ t - n j=1 τ j , = ∞ 0 . . . ∞ 0 χ(τ 1 ) . . . χ(τ n )dτ 1 . . . dτ n e -s n j=1 τ j = ∞ 0 χ(τ )e -sτ dτ n = ( χ(s)) n (7)
which by putting s = 0 yields the normalization condition ∞ 0 χ(n) (t)dt = 1 as a consequence of the normalization of the one-jump density χ(s)| s=0 = 1.

The probability that a walker within the time interval [0, t] has made n = 0, 1, 2, . . . steps then is related with (6) by the convolution

Φ (n) (t) = t 0 (1 -Ψ(t -τ )) χ (n) (τ )dτ = t 0 Φ (0) (t -τ )χ (n) (τ )dτ (8)
where we account for that the walker makes its nth jump at a time τ < t with probability χ (n) (τ )dτ and waits during tτ with probability Φ (0) (t) = 1 -Ψ(tτ ) where over all combinations from τ = 0 to t is integrated. The probabilities [START_REF] Montroll | Random walks on lattices II[END_REF] are dimensionless distributions whereas the n-jump densities [START_REF] Barkai | From continuous time random walks to the fractional Fokker-Planck equation[END_REF] have dimension [time] -1 . We mention the normalization condition shown later

∞ n=0 Φ (n) (t) = 1. (9)
We observe then the recursive relations

χ (n) (t) = t 0 χ (n-1) (τ )χ(t -τ )dτ, χ (0) (t) = δ(t) Φ (n) (t) = t 0 Φ (n-1) (τ )χ(t -τ )dτ, Φ (0) (t) = 1 -Ψ(t). ( 10
)
We then get for the Laplace transform of ( 8)

Φ(n) (s) = 1 -χ(s) s ( χ(s)) n = 1 s χ(n) (s) -χ(n+1) (s) , n = 0, 1, 2, . . . ∞, (11)
recovering Φ0 (s) = 1-χ(s) s when we account for χ 0 (t) = δ(t) with χ(0) (s) = L{δ(t)} = 1. Introducing a generating function in the Laplace domain as5

G(s, v) = ∞ n=0 Φ(n) (s) v n = 1 -χ(s) s 1 1 -v χ(s) , |v| ≤ 1 (12)
where G(s, 1) = s -1 proving then normalization condition [START_REF] Scher | Stochastic Transport in a Disordered Solid. I. Theory[END_REF]. A quantity of great interest in the following analysis is the (dimensionless) expected number of arrivals n(t) (expected number of jumps) that happen within the time interval [0, t]. This quantity is obtained by

n(t) = ∞ n=0 nΦ (n) (t) = d dv G(t, v) v=1 = L -1 χ(s) s(1 -χ(s)) . ( 13
)
For further discussions of renewal processes we refer to the references [START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF][START_REF] Cox | Renewal Theory[END_REF]. With these general remarks on renewal processes, let us briefly discuss in the following section Laskin's fractional Poisson process in order to generalize this process in Section 4.

Fractional Poisson distribution

The fractional Poisson process was introduced by Laskin and several aspects of this process were thoroughly analyzed in several papers [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF][START_REF] Mainardi | A fractional generalization of the Poisson processes[END_REF][START_REF] Laskin | Some applications of the fractional Poisson probability distribution[END_REF]. The Laskin fractional Poisson process is defined by a renewal process with a Mittag-Leffler waiting time density6

χ β (t) = ∞ n=1 (-1) n-1 ξ n t nβ-1 (nβ -1)! = ξt β-1 ∞ n=0 (-1) n ξ n t nβ Γ(βn + β) = ξt β-1 E β,β (-ξt β ), t > 0 = d dt [1 -E β (-ξt β )], ξ > 0 0 < β ≤ 1, ξ > 0 (14)
where we introduced the Mittag-Leffler function of order β [START_REF] Gorenflo | From Power Laws to Fractional Diffusion: the Direct Way[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF][START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | Mittag-Leffler Functions, Related Topics and Applications[END_REF]]

E β (z) = ∞ n=0 z n Γ(βn + 1) (15)
and the generalized Mittag-Leffler function

E β,γ (z) = ∞ n=0 z n Γ(βn + γ) (16)
where E β,1 (z) = E β (z). Since the argument ξt β of the Mittag-Leffler function is dimensionless, ξ is a characteristic parameter of physical dimension [time] -β . It is straight-forward to show that the Mittag-Leffler jump density [START_REF] Solomon | Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow[END_REF] has the Laplace transform

χβ (s) = ∞ n=1 (-1) n-1 ξ n ∞ 0 e -st t nβ-1 (nβ -1)! dt = ∞ n=1 (-1) n-1 ξ n s -nβ = ξ ξ + s β , s = σ + iω, σ > ξ 1 β (17)
converging for σ > ξ

β (t) = 1 -Ψ β (t) = L -1 s β-1 ξ + s β = E β (-ξt β ), 0 < β ≤ 1 ( 18
)
where for β = 1 the standard Poisson process with survival probability E 1 (-ξt) = e -ξt is recovered. The fractional Poisson distribution is defined as the probability Φ

(n) β (t) of exactly n arrivals within time interval [0, t] for a renewal process with jump density [START_REF] Solomon | Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow[END_REF] and is obtained as [START_REF] Laskin | Fractional Poisson process[END_REF] Φ (n)

β (t) = (ξt β ) n n! d n dτ n E β (τ ) τ =-ξt β = (ξt β ) n n! ∞ m=0 (m + n)! m! (-ξt β ) m Γ(β(m + n) + 1) , 0 < β ≤ 1, n = 0, 1, 2, . . . (19)
For β = 1 the fractional Poisson distribution and α = 1, β = 1 recovers the standard Poisson distribution

Φ (n) 1 (t) = Θ(t) (ξt) n n! e -ξt , n = 0, 1, 2, . . . ( 20
)
4 Generalization of the fractional Poisson process

We now analyze a counting process that generalizes the fractional Poisson process. The waiting time density of this generalized process has the Laplace transform

χβ,α (s) = ξ α (s β + ξ) α , 0 < β ≤ 1, α > 0, ξ > 0 (21)
and was introduced by Cahoy and Polito [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF]. The goal of the remaining parts is to complement their analysis of the process and to develop a CTRW that is governed by this process to derive generalizations of the fractional Kolmogorov-and diffusion equations. We see that [START_REF] Repin | Fractional Poisson law[END_REF] generalizes [START_REF] Chechkin | Continuous-time random walk with correlated waiting times[END_REF] of the Laskin fractional Poisson process where any non-negative power α ∈ R + is admissible in [START_REF] Repin | Fractional Poisson law[END_REF]. We call the renewal process with waiting time density defined by ( 21) the generalized fractional Poisson process (GFPP). The characteristic dimensional constant ξ in (21) has physical dimension [time] -β and defines a characteristic time scale. Per construction χβ,α (s)| s=0 = ∞ 0 χ β,α (t)dt = 1 reflecting normalization of the waiting time density. The GFPP recovers for α = 1 , 0 < β < 1 the fractional Poisson process, for α > 0, β = 1 an Erlang type process (See [START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF]) allowing also fractional powers α, and for α = 1, β = 1 the standard Poisson process.

We observe in ( 21) that the jump density has for 0 < β < 1 and for all α > 0 diverging mean (diverging expected time for first arrival), namely

-d ds χβ,α (s)| s=0 = ∞ 0 tχ β,α (t)
→ +∞, and a finite mean first arrival time in the Erlang regime αξ -1 for α > 0 with β = 1.

Beforehand we observe two limiting cases of 'extreme' behaviors: For infinitesimally small α we have lim α→0+ χβ,α (s) = 1 and hence lim α→0+ χ β,α (t) = L -1 {1} = δ(t), i.e. the smaller α the more the jump density becomes localized around t = 0+ and taking for α → 0+ the shape of a δ(t)-peak. In this limiting case the first arrival happens 'immediately' when the observation starts at t = 0+ with an extremely fast jump dynamics (Eqs. ( 29)). On the other hand in the limit α → ∞ the Laplace transform of the jump density (1 + ξ -1 s β ) -α is equal to one at s = 0 and nonzero extremely localized around s ≈ 0 dropping to zero immediately at s = 0. As a consequence in this limit the jump density becomes extremely delocalized with quasi-equally distributed waiting times χ β,α (t) ∼ → 0+ for all t > 0 thus extremely long waiting times may occur (shown subsequently by Eq. ( 35)). In order to determine the time domain representation of [START_REF] Repin | Fractional Poisson law[END_REF] it is convenient to introduce the Pochhammer symbol which is defined as [START_REF] Mathai | Some properties of Mittag-Leffler functions and matrix variant analogues: A statistical perspective[END_REF] (c

) m = Γ(c + m) Γ(c) =    1, m = 0 c(c + 1) . . . (c + m -1), m = 1, 2, . . . (22)
Despite Γ(c) is singular at c = 0 the Pochhammer symbol can be defined also for c = 0 by the limit (0) m = lim c→0+ (c) m = δ m0 which is also fulfilled by the right hand side of [START_REF] Laskin | Some applications of the fractional Poisson probability distribution[END_REF]. Then (c) m is defined for all c ∈ C. The jump density defined by ( 21) can be evaluated by accounting for

χβ,α (s) = ξ α s -βα (1 + ξs -β ) -α = ξ α ∞ m=0 (-1) m (α) m m! ξ m s -β(m+α) , σ = (s) > ξ 1 β . ( 23
)
The time domain representation of the jump density of the GFPP is then obtained by using Laplace inversion (106), namely (See also [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF])

χ β,α (t) = L -1 ξ α (s β + ξ) α = e σt 2π ∞ -∞ e iωt ξ α (ξ + (σ + iω) β ) α dω = ξ α ∞ m=0 (-1) m (α) m m! ξ m L -1 {s -β(m+α) }, t > 0, σ = {s} > ξ 1 β , 0 < β ≤ 1, α > 0 = ξ α t βα-1 ∞ m=0 (α) m m! (-ξt β ) m Γ(βm + αβ) = ξ α t βα-1 E α β,αβ (-ξt β ). (24)
In expression [START_REF] Gorenflo | The asymptotic universality of the Mittag-Leffler waiting time law in continuous time random walks[END_REF] appears a generalization of the Mittag-Leffler function which was introduced by Prabhakar [START_REF] Prabhakar | A singular integral equation with a generalized Mittag-Leffler function in the kernel[END_REF] defined by

E c a,b (z) = ∞ m=0 (c) m m! z m Γ(am + b) , {a} > 0, {b} > 0, c, z ∈ C ( 25
)
where with (0

) m = δ m0 we have E 0 a,b (z) = 1.
The generalized Mittag-Leffler function [START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF] was analyzed by several authors [START_REF] Mathai | Some properties of Mittag-Leffler functions and matrix variant analogues: A statistical perspective[END_REF][START_REF] Shukla | On a generalization of Mittag-Leffler function and its properties[END_REF][START_REF] Haubold | Mittag-Leffler functions and their applications[END_REF][START_REF] Gara | The Prabhakar or three parameter Mittag-Leffler function: Theory and application[END_REF]. This series converges absolutely in the entire complex z-plane. For t small the jump density behaves as

lim t→0+ χ β,α (t) = lim t→0 ξ α Γ(αβ) t αβ-1 →            0, αβ > 1, ξ 1 β , αβ = 1, ∞ αβ < 1. (26)
The power-law asymptotic behavior ∼ t αβ-1 occurring for t small is different from the fractional Poisson process which has ∼ t β-1 and is reproduced for α = 1 (See the order n = 0 in the Mittag-Leffler density ( 14)). The GFPP jump density allows exponents αβ > 1 where the jump density approaches zero for t → 0+. This behavior does not exist in the case α = 1 of the fractional Poisson process.

We see in the Figures 1 and2 that lim t→0+ χ β,α (t) ∼ t αβ-1 → +∞ occurs in the regime αβ < 1. The jump density in this regime is extremely high at short observation times, where this tendency is amplified for α → 0+ with δ(t)-peak localized jump densities. Further we observe in Figure 1 the 'phase transition' at αβ = 1 separating the behaviors of Eq. ( 26) for t → 0+. On the other hand for α → ∞ the jump density exhibits an extremely broad distribution of waiting times. For α = 1 the density (24) recovers the Mittag-Leffler type jump density of the fractional Poisson process [START_REF] Solomon | Observation of anomalous diffusion and Lévy flights in a two-dimensional rotating flow[END_REF], and for β = 1, α > 0 is an Erlang type density where non-integer α > 0 are allowed (See Eq. ( 33)). For α = 1, β = 1 the jump density χ 1,1 (t)) = ξe -ξt of the standard Poisson process is reproduced. For α = 1 the expression of the jump density [START_REF] Gorenflo | The asymptotic universality of the Mittag-Leffler waiting time law in continuous time random walks[END_REF] recovers the jump density ( 14) of the fractional Poisson process. Integrating (24) yields the cumulated distribution (4) of first arrivals

Ψ β,α (t) = t 0 χ β,α (τ )dτ = ξ α t αβ ∞ m=0 (α) m m! (-ξt β ) m Γ(βm + αβ + 1) = ξ α t αβ E α β,(αβ+1) (-ξt β ). ( 27
)
The probability distribution ( 27) is plotted for different values of 0 < β < 1, and α = 2 in Figure 1 and for different values α > 0 in Figure 2 for β = 0.5. For small ξ 1 β t, this distribution is dominated by order n = 0, namely 24) and ( 27). We represent with dashed lines the results for the particular case αβ = 1.

Ψ β,α (t) ≈ ξ α Γ(αβ + 1) t αβ , ξ 1 β t → 0. ( 28
)
We see in this relation that in the admissible range 0 < β ≤ 1, α > 0 the initial condition Ψ β,α (0) = 0 is always fulfilled. We hence observe the following distributional properties

lim α→0+ Ψ β,α (t) = lim α→0+ ξ α t αβ Γ(αβ + 1) = Θ(t), lim α→0+ χ β,α (t) = d dt Θ(t) = δ(t) (29)
and similar relations are obtained in the limit αβ → 0+. This behavior is consistent with the observation that in the Laplace domain Ψβ,α (s) = s -1 χβ,α (s) for α → 0 behaves as s -1 . The t αβ -power-law behavior can be observed in the Figures 1 and2. The smaller αβ → 0 the more (28) approaches a Heaviside step function shape where for small observation times an extremely high jump frequency occurs so that Ψ β,α (t) rapidly approaches maximum value one. The survival probability then is given by

Φ (0) β,α (t) = ∞ t χ β,α (τ )dτ = 1 -Ψ β,α (t) = 1 -ξ α t αβ E α β,(αβ+1) (-ξt β ) (30) with lim t→∞ Φ (0) β,α (t) = 0. For α = 1 (30) becomes Φ (0) β,1 (t) = 1 -ξt β E 1 β,(β+1) (-ξt β ) = 1 + t β ∞ m=0 (-1) m+1 ξ m+1 (-ξt β ) m Γ(βm + β + 1) = ∞ m=0 (-ξt β ) m Γ(βm + 1) = E β (-ξt β ) (31)
reproducing the Mittag-Leffler survival probability [START_REF] Feller | An introduction to probability theory and its applications[END_REF] of the fractional Poisson process Mittag-Leffler jump density

χ β,1 (t) = χ β (t) = ξt β-1 E β,β (-ξt β ) = d dt (1 -E β (-ξt β )) (32)
coinciding with the known results [START_REF] Laskin | Fractional Poisson process[END_REF]. On the other hand let us consider β = 1 and α > 0 which gives

χ 1,α (t) = ξ α t α-1 E α 1,α (-ξt) = ξ α t α-1 Γ(α) ∞ m=0 (-ξt) m m! = ξ α t α-1 Γ(α) e -ξt (33)
i.e. for β = 1 the jump density is light-tailed with exponentially evanescent behavior for t → ∞, but for t → 0 we get power-law behavior χ 1,α (t) ∼ t α-1 . The density [START_REF] Mathai | Some properties of Mittag-Leffler functions and matrix variant analogues: A statistical perspective[END_REF] is the so-called Erlang density and the renewal process generated with this density is referred to as Erlang process [START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF]. For α = 1 the Erlang density [START_REF] Mathai | Some properties of Mittag-Leffler functions and matrix variant analogues: A statistical perspective[END_REF] again recovers the density χ 1,1 (t) = ξe -ξt of the standard Poisson process.

Let us now consider the behavior for large times t → ∞. Expanding [START_REF] Repin | Fractional Poisson law[END_REF] for |s| small yields7

χβ,α (s

) = 1 + s β ξ -α = ∞ m=0 (α) m m! (-1) m ξ (-m) s mβ = 1 - α ξ s β + . . . (34)
which yields for large observation times for 0

< β < 1, α > 0 fat-tailed behavior 8 χ β,α (t) ≈ - α ξΓ(-β) t -β-1 , 0 < β < 1, α > 0, t → ∞ ( 35
)
where the exponent does not depend on α.

Generalized fractional Poisson distribution

In this section we derive the generalized counterpart Φ

(n)

β,α (t) to Laskin's fractional Poisson distribution [START_REF] Laskin | Fractional Poisson process[END_REF], i.e. the probabilities for n arrivals 9 within [0, t] in a GFPP. For the evaluation we utilize the general relation [START_REF] Klafter | Beyond Brownian Motion[END_REF] together with [START_REF] Repin | Fractional Poisson law[END_REF], namely

Φ (n) β,α (t) = L -1 1 s χn β,α (s) -χn+1 β,α (s) = L -1 1 s χβ,nα (s) -χβ,(n+1)α (s) , n = 0, 1, 2, . . . (36
) where χn β,α (s) = ξ nα (ξ+s β ) nα . We then obtain by using Eq. ( 24)

Ψ β,nα (t) = L -1 1 s χβ,nα (s) = ξ nα ∞ m=0 (-1) m (nα) m m! ξ m L -1 {s -βm-nαβ-1) }, σ = {s} > ξ 1 β = ξ nα ∞ m=0 (nα) m m! (-1) m ξ m t mβ+nαβ Γ(mβ + nαβ + 1) = ξ nα t nαβ E nα β,(nαβ+1) (-ξt β ), n = 0, 1, 2, . . . (37)
With this result, from [START_REF] Haubold | Mittag-Leffler functions and their applications[END_REF] we obtain the probability of exactly n arrivals within [0, t] as

Φ (n) β,α (t) = ξ nα t nαβ E nα β,(nαβ+1) (-ξt β ) -ξ α t αβ E (n+1)α β,((n+1)αβ+1) (-ξt β ) , = Ψ β,nα (t) -Ψ β,(n+1)α (t), 0 < β ≤ 1, α > 0, n = 0, 1, 2, . . . (38)
where E c a,b (z) denotes the generalized Mittag-Leffler function defined in Eq. ( 25). This expression was also obtained in Ref. [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF] (Eq. (2.8) in that paper). We call the distribution Φ (n) β,α (t) in Eq. ( 38) the generalized fractional Poisson distribution (GFPD). In the second line of Eq. [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] we have taken into account that cumulated probabilities Ψ β,να (t) = t 0 χ β,να (t)dt indicate the probability for at least ν arrivals within [0, t]. The Ψ β,να (t) are determined by Eq. ( 27) by replacing there α → να. The case n = 0 is also covered in relation [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] and yields Φ (0) β,α (t) = 1 -Ψ β,α (t) of Eq. [START_REF] Michelitsch | Generalized fractional Poisson process and related stochastic dynamics[END_REF]. We can then write distribution [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] in the following compact form

Φ (n) β,α (t) = (ξt β ) nα 1 Γ(nαβ + 1) + ∞ m=1 A m,n α,β (t β ξ) (-ξt β ) m (39)
with

A m,n α,β (t β ξ) = (nα) m m!Γ(nαβ + βm + 1) + (ξt β ) α-1 ((n + 1)α) m-1 (m -1)!Γ(nαβ + βm + (α -1)β + 1) (40)
for m = 1, 2, . . . and n = 0, 1, 2, . . .. The GFPD is a dimensionless probability distribution depending only on the dimensionless time tξ 1 β where ξ -1 β has dimension of time and defines a characteristic time scale in the GFPP. It follows that for small (dimensionless) times the GFPD behaves as

Φ (n) β,α (t) ≈ (ξt β ) nα Γ(nαβ + 1)
, tξ

1 β → 0, n = 0, 1, 2, . . . (41)
which is the inverse Laplace transform of ξ nα s -nαβ-1 , i.e. of the order m = 0 in the expansion [START_REF] Gara | The Prabhakar or three parameter Mittag-Leffler function: Theory and application[END_REF]. It follows that the general relation Φ

(n) β,α (t) t=0 = δ n0 (42)
is fulfilled, i.e. at t = 0 per construction the walker is on his departure site. Further of interest is the asymptotic behavior for large (dimensionless) times tξ 1 β . To this end, let us expand the Laplace transform for small s → 0 in [START_REF] Haubold | Mittag-Leffler functions and their applications[END_REF] up to the lowest non-vanishing order in s β ξ to arrive at

Φ (n) β,α (t) ≈ α ξ L -1 {s β-1 } = α Γ(1 -β) tξ 1 β -β , tξ 1 β → ∞, 0 < β < 1, α > 0, n = 0, 1, . . . , ∞ (43)
where this inverse power law holds universally for all α > 0 for 0 < β < 1 and is independent of the number of arrivals n and contains the Laskin case α = 1. That is all 'states' n exhibit for tξ 1 β the same (universal) t -β inverse power-law decay. We conjecture that the asymptotic equal-distribution (43) in a wider sense can be interpreted as quasi-ergodic property of the GFPP. Now let us consider the important limit α = 1 for the GFPD in more details. Then the functions (40) become time independent coefficients, namely

A m,n 1,β = (n + m -1)! Γ(β(m + n) + 1) 1 m!(n -1)! + 1 (m -1)!n! = (n + m)! n!m! 1 Γ(β(m + n) + 1) (44)
and for the oder m = 0 we have

1 Γ(nβ+1) = (n+m)! n!m! 1 Γ(β(m+n)+1)
m=0 thus we get for the GFPP (39) for α = 1 the expression

Φ (n) β,1 (t) = (ξt β ) n n! ∞ m=0 (n + m)! m! (-ξt β ) m Γ(β(m + n) + 1) (45)
which we identify with Laskin's fractional Poisson distribution of Eq. ( 19) [START_REF] Laskin | Fractional Poisson process[END_REF][START_REF] Laskin | Some applications of the fractional Poisson probability distribution[END_REF].

Expected number of arrivals in a GFPP

Here we analyze the asymptotic behavior of the average number of arrivals n(t) within the interval of observation [0, t]. To this end we consider the generating function

Gβ,α (v, s) = ∞ n=0 v n Φ(n) β,α (s) = 1 -χβ,α (s) s(1 -v χβ,α (s)) = 1 s (1 + s β ξ ) α -1 (1 + s β ξ ) α -v , |v| ≤ 1 (46) with Φ(n) β,α (s) = 1 -χβ,α (s) s χβ,(nα) (s). ( 47
)
The expected number nβ,α (t) the walker makes in [0, t] is obtained from nβ,α (t

) = d dv G β,α (v, t)| v=1 = ∞ n=0 nΦ (n)
β,α (t) with the Laplace transform [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF] ñβ,α (s

) = d dv Gβ,α (v, s) v=1 = χβ,α (s) s(1 -χβ,α (s)) = 1 s[( s β ξ + 1) α -1] . ( 48
)
The Laplace transform behaves as ñβ,α (s) ≈ ξ α s -β-1 for |s| → 0. We obtain hence the asymptotic behavior for tξ

1 β → ∞ large nβ,α (t) ≈ ξ αΓ(β + 1) t β , 0 < β ≤ 1, α > 0, tξ 1 β → ∞. ( 49
)
This result includes the Erlang regime β = 1 where for all α > 0 the average number of arrivals n1,α (t) ≈ ξ α t increases linearly for large tξ Poisson process, e.g. [START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF]. On the other hand, we obtain the asymptotic behavior for tξ 1 β → 0 small when we expand ñβ,α (s) with for |s| → ∞ and obtain

nβ,α (t) = L -1 1 s[( s β ξ + 1) α -1] ≈ L -1 ξ α s -αβ-1 = (ξt β ) α Γ(αβ + 1)
, tξ

1 β → 0. ( 50
)
It is important to note that both asymptotic expressions (49) and (50) coincide for α = 1 recovering the known expression of the fractional Poisson process.

Montroll-Weiss CTRW on undirected networks

As an application we consider a random walk on an undirected network subordinated to a GFPP and derive the generalization of the fractional Kolmogorov-Feller equation. Then in subsequent Section 8 we develop this process for the infinite d-dimensional integer lattice Z d and analyze the resulting 'generalized fractional diffusion'. To this end let us briefly outline the Montroll-Weiss CTRW on undirected networks. We assume an undirected connected network with N nodes which we denote with p = 1, . . . , N . For our convenience we employ Dirac's |bra ket| notation. In an undirected network the positive-semidefinite N × N Laplacian matrix L = (L pq ) characterizes the network topology and is defined by [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF][START_REF] Noh | Random walks on complex networks[END_REF][START_REF] Riascos | Fractional dynamics on networks: Emergence of anomalous diffusion and Lévy flights[END_REF]]

L pq = K p δ pq -A pq (51)
where A = (A pq ) denotes the adjacency matrix having elements one if a pair of nodes is connected and zero otherwise. Further we forbid that nodes are connected with themselves thus A pp = 0. In an undirected network adjacency and Laplacian matrices are symmetric. The diagonal elements L pp = K p of the Laplacian matrix are referred to as the degrees of the nodes p counting the number of neighbor nodes of a node p with K p = N q=1 A pq . The one-step transition matrix W = (W pq ) relating the network topology with the random walk then is defined by [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF][START_REF] Noh | Random walks on complex networks[END_REF]

W pq = 1 K p A pq = δ pq - 1 K p L pq ( 52
)
which generally is a non-symmetric matrix for networks with variable degrees K i = K j (i = j). The one-step transition matrix W pq defines the conditional probability that the walker which is on node p jumps in one step to node q where in one step only neighbor nodes with equal probability 1 Kp can be reached.

We consider now a random walker that performs a CTRW with IID random steps at random times 0 ≤ t 1 , t 2 , . . . t n , . . . ∞ on the network where the observation starts at t = 0. Each step of the walker from one to another node is associated with a jump event or arrival in a CTRW with identical transition probability (W pq ) for a step from node p to node q. We introduce the transition matrix P(t) = (P ij (t)) indicating the probability to find the walker at time t on node j under the condition that the walker was sitting at node i at time t = 0 when the observation starts.

The transition matrix fulfills the normalization condition N j=1 P ij (t) = 1 and stochasticity of the transition matrix implies that 0 ≤ P ij (t) ≤ 1. We restrict us on connected undirected networks and allow variable degrees. In such networks the transition matrix is non-symmetric P ij (t) = P ji (t) (See [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] for a detailed analysis). Assuming the initial condition P(t = 0), then the probability to find the walker on node j at time t is given by a series of the form [START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF][START_REF] Cox | Renewal Theory[END_REF]]

P(t) = P(0) ∞ n=0 Φ (n) (t)W n (53)
where the Φ (n) (t) indicate the probabilities of n arrivals, i.e. that the walker made n steps within interval [0, t] of relation [START_REF] Montroll | Random walks on lattices II[END_REF]. It is straight-forward to see in (53) that the normalization condition

N j=1 P ij (t) = ∞ n=0 Φ (n) (t) =
1 is fulfilled during the entire observation time t ≥ 0. The convergence of this series is proved by using that W has uniquely eigenvalues |λ m | ≤ 1 [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] together with | χ(s)| ≤ 1 for the Laplace transforms of the jump density. Let us assume that at t = 0 the walker is sitting on its departure node P ij (0) = (δ ij ), then the Laplace transform of (53) writes

P(s) = (1 -χ(s)) s {1 -χ(s)W} -1 . (54)
Let us take into account the canonic representation of the one-step transition matrix [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF] where P(s) has the eigenvalues

W = |Φ 1 Φ1 | + N m=2 λ m |Φ m Φm | (
P (m, s) = (1 -χ(s)) s 1 (1 -λ m χ(s)) . ( 57
)
This expression is the celebrated Montroll-Weiss formula [START_REF] Montroll | Random walks on lattices II[END_REF]. Since λ 1 = 1 we see that P (1, s) = 1 s thus the stationary amplitude always is of the form P (1, t) = L -1 { 1 s } = Θ(t).

Generalized fractional Kolmogorov-Feller equation

With the general remarks on CTRWs of the previous section we now analyze a random walk on an undirected network which is subordinated to a GFPP. For this diffusional process we derive a generalization of the fractional Kolmogorov-Feller equation [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF]. To this end let us first consider the Laplace transforms of the GFPD Φ

(n) β,α (t) (See Eq. ( 36
)) Φ(n) β,α (s) = Φ(0) β,α (s) ξ nα (s β + ξ) nα , α > 0, 0 < β ≤ 1, n = 0, 1, 2, . . . (58)
where

χ(n) β,α (s) = χβ,(nα) (s) = ξ nα (s β + ξ) nα (59)
and especially for n = 0 we have

Φ(0) β,α (s) = 1 -χβ,α (s) s = (s β + ξ) α -ξ α s(s β + ξ) α . ( 60
)
10 reflecting row-normalization 11 In this demonstration we ignore the cases of bipartite graphs where a unique eigenvalue -1 occurs.

N q=1 W pq = 1.
These Laplace transforms fulfill

(s β + ξ) α Φ(n) β,α (s) = ξ α Φ(n-1) β,α (s), n = 1, 2, . . . , (61)
and for n = 0 we have

(s β + ξ) α Φ(0) β,α (s) = (s β + ξ) α -ξ α s . ( 62
)
In the time domain these equations write

0 D β,α t Φ (n) β,α (t) = ξ α Φ (n-1) β,α (t), n = 1, 2, . . . , (63)
and

0 D β,α t (ξ)Φ (0) β,α (t) = L -1 (s β + ξ) α -ξ α s = K (0) β,α (t) -ξ α Θ(t) ( 64
)
where 0 D β,α t Φ(t) =:

t 0 Φ(τ )D β,α (t -τ
)dτ has to be read as a convolution. Now the goal is to derive the generalized counterpart to the fractional Kolmogorov-Feller equation [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF]. A major role is played by the kernel

D β,α (t) = L -1 (s β + ξ) α . ( 65
)
This kernel has the explicit representation [START_REF] Michelitsch | Generalized fractional Poisson process and related stochastic dynamics[END_REF] (and see also Appendix A)12

D β,α (t) = L -1 (s β + ξ) α = d αβ dt αβ Θ(t)d β,α (t) =        d αβ dt αβ Θ(t)t αβ -βα-1 ∞ m=0 α! (α-m)!m! (ξt β ) m Γ(βm+ αβ -βα) , αβ / ∈ N d αβ dt αβ δ(t) + Θ(t) d dt ∞ m=0 α! (α-m)!m! (ξt β ) m Γ(mβ+1) , αβ ∈ N =    d αβ dt αβ Θ(t)t αβ -βα-1 E α,β,( αβ -αβ) (ξt β ) , αβ / ∈ N d αβ dt αβ δ(t) + Θ(t) d dt E α,β,1 (ξt β ) , αβ ∈ N. (66)
In this expression we introduced the ceiling function γ indicating the smallest integer greater or equal to γ. In (66) occurs a generalized Mittag-Leffler type function E c,a,b (z) which is related with the generalized Prabhakar Mittag-Leffler function [START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF] and given by

E c,a,b (z) = ∞ m=0 c! (c -m)!m! z m Γ(am + b) = ∞ m=0 (-c) m m! (-z) m Γ(am + b) = E -c a,b (-z), {a} > 0, {b} > 0.
(67) The second kernel K (0) β,α (t) of equation ( 64) is then explicitly obtained as [START_REF] Michelitsch | Generalized fractional Poisson process and related stochastic dynamics[END_REF]

K (0) β,α (t) = L -1 (s β + ξ) α s =            Θ(t)t -αβ E α,β,1-αβ (ξt β ), , 0 < αβ < 1 d αβ -1 dt αβ -1 Θ(t)t αβ -βα-1 E α,β,( αβ -αβ) (ξt β ) αβ > 1, αβ / ∈ N d αβ-1 dt αβ-1 δ(t) + Θ(t) d dt E α,β,1 (ξt β ) , αβ ≥ 1, αβ ∈ N. ( 68
)
The transition probability matrix (53) of the GFPP walk then writes

P β,α (t) = P(0)Φ (0) β,α (t) + P(0) ∞ n=1 Φ (n) β,α (t)W n . ( 69
)
Then we obtain with relations (63), (64) the convolution equation 13 0 D β,α t

• P β,α (t) = P(0) K (0) β,α (t) -ξ α Θ(t) + ξ α P β,α (t)W, t 0 D β,α (t -τ )P β,α (τ )dτ = P(0) K (0) β,α (t) -ξ α Θ(t) + ξ α P β,α (t)W. ( 70
)
We call equation ( 70) the generalized fractional Kolmogorov-Feller equation where

K (0)
β,α (t)ξ α Θ(t) can be conceived as memory kernel. Eq. ( 70) has for αβ / ∈ N the explicit representation

d αβ dt αβ t 0 (t -τ ) αβ -βα-1 E α,β,( αβ -αβ) (ξ(t -τ ) β )P β,α (τ )dτ, αβ / ∈ N = P(0) K (0) β,α (t) -ξ α Θ(t) + ξ α P β,α (t)W, 0 < β < 1, α > 0. ( 71
)
Let us next consider the cases α = 1 of Laskin's fractional Poisson process, and of α = 1, β = 1 representing the standard Poisson process.

We first consider (71) in the fractional Poisson process limit α = 1 and 0 < β < 1. Then we get for the kernel ( 68)

K (0) β,1 (t) = t -β E 1,β,1-β (ξt β ) = Θ(t) t -β Γ(1 -β) + ξΘ(t) (72)
with

E 1,β,1-β (ξt β ) = 1 Γ(1 -β) + ξt β Γ(1) = 1 Γ(1 -β) + ξt β . ( 73
)
Only the orders m = 0 and m = 1 in ξt β are nonzero. Then with αβ = β = 1 we get for (71)

d dt t 0 (t -τ ) -β E 1,β,1-β (ξ(t -τ ) β )P β,1 (τ )dτ = P(0) K (0) β,1 (t) -ξΘ(t) + ξP β,1 (t)W ( 74
)
which takes with (72) and (73) the form (t > 0)

d dt t 0 (t -τ ) -β 1 Γ(1 -β) + ξ(t -τ ) β P β,1 (τ )dτ = t -β Γ(1 -β) P(0) + ξP β,1 (t)W. ( 75
)
On the left hand side we identify 1

Γ(1-β) d dt t 0 (t-τ ) -β P β,1 (τ )dτ = 0 D β t P β,1 (
t) with the Riemann-Liouville fractional derivative of order β. The second term on the left hand side yields ξ d dt t 0 P β,1 (τ )dτ = ξP β,1 (t). Hence we can write (75) in the form

0 D β t P β,1 (t) = t -β Γ(1 -β) P(0) -ξP β,1 (t) + ξP β,1 (t)W, 0 < β < 1. ( 76
)
13 Read 0 D β,α t

• P β,α (t) = t 0 D β,α (t -τ )P β,α (τ )dτ =: d αβ dt αβ t 0 d β,α (t -τ )P β,α (τ )dτ .
We identify this equation with the fractional Kolmogorov-Feller equation [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF][START_REF] Laskin | Fractional Poisson process[END_REF]. In this way we have proved that the generalized fractional Kolmogorov-Feller equation (70) for α = 1 recovers the fractional Kolmogorov-Feller equation introduced in [START_REF] Saichev | Fractional kinetic equations: solutions and applications[END_REF] and see also the references [START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF][START_REF] Meerschaert | The Fractional Poisson Process and the Inverse Stable Subordinator[END_REF][START_REF] Hilfer | Fractional master equations and fractal time random walks[END_REF].

In the fractional Kolmogorov-Feller equation ( 76) occurs the 'memory term' t -β Γ(1-β) P(0) which reflects the slow power-law decay of the influence of the initial condition. Generally for the GFPP the long-time memory (non-markovianity) is a consequence that the jump density is fat-tailed ∼ t -1-β , see Eq. [START_REF] Shukla | On a generalization of Mittag-Leffler function and its properties[END_REF]. Let us demonstrate briefly that the same memory effect occurs in walks subordinated to GFPPs for any α > 0. To this end it is instructive to consider the generalized fractional Kolmogorov-Feller equation (70) in the Laplace domain

(s β + ξ) α Pβ,α (s) = (s β + ξ) α -ξ α s P(0) + ξ α W Pβ,α (s). (77)
For |s| → 0 small this equation takes the representation

s β Pβ,α (s) = s β-1 P(0) - ξ α Pβ,α (s) + ξ α W Pβ,α (s), 0 < β < 1. (78)
Transforming Eq. ( 78) into the time domain yields fractional Kolmogorov-Feller equation (76) as the the asymptotic limit for ξ 1 β t → ∞ of the generalized fractional fractional Kolmogorov-Feller equation (70) (with ξ → ξ = ξ α ). It follows that the GFPP generates asymptotically the same memory effect as the fractional Poisson process with the memory term t -β Γ(1-β) P(0) (0 < β < 1) which is independent of α. The memory effect becomes especially pronounced when β → 0+ with extremely slow decay. On the other hand in the limit β → 1 -0 the memory term takes δ(t)-representation thus the walk then becomes memoryless.

Finally we consider the limit α = 1 and β = 1 of standard Poisson. Then Eqs. ( 61) and (62) take in the time domain the form

d dt + ξ Φ (n) 1 (t) = ξΦ (n-1) 1 (t) (79)
and

d dt + ξ Φ (0) 1 (t) = δ(t) (80)
where these relations define the standard Poisson process and are easily confirmed to be fulfilled by the Poisson distribution [START_REF] Mainardi | A fractional generalization of the Poisson processes[END_REF]. Plugging the last two equations into (70) for β = 1 and α = 1 yields relation

d dt P 1,1 (t) = P(0)δ(t) -ξP 1,1 (t) + ξP 1,1 (t)W (81)
which is known as the Kolmogorov-Feller equation of the standard Poisson process [START_REF] Laskin | Fractional Poisson process[END_REF][START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF][START_REF] Gorenflo | On the Fractional Poisson Process and the Discretized Stable Subordinator[END_REF]. One can also recover the Kolmogorov Feller equation (81) from Eq. (76) in the limit β → 1 -0.

Generalized fractional diffusion in Z d

In this section we analyze the features of a random walk subordinated to a GFPP in the infinite d-dimensional integer lattice Z d . We refer the resulting diffusion process to as 'generalized fractional diffusion'. The eigenvalues of the one-step transition matrix (55) for the d-dimensional infinite lattice are given by the Fourier transforms [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF][START_REF] Michelitsch | On recurrence of random walks with long-range steps generated by fractional Laplacian matrices on regular networks and simple cubic lattices[END_REF]

λ( k) = 1 - 1 2d µ( k), µ( k) = 2d -2 d j=1 cos(k j ) λ( k) = 1 d d j cos(k j ), -π ≤ k j ≤ π. (82)
Here µ( k) denote the eigenvalues of the Laplacian matrix (51) of the lattice. In this walk the walker in one jump can reach only connected neighbor nodes. Such walk refers to as 'normal walk' [START_REF] Noh | Random walks on complex networks[END_REF]. In the present section we subordinate a normal walk on the Z d to a GFPP. In this lattice the eigenvectors become Bloch-waves |Φ m = | k . The Laplace transform (56) of the transition matrix is then obtained as

P(s, p -q) = (1 -χ(s)) s 1 (2π) d π -π dk 1 . . . π -π dk d e -i k•( p-q) (1 -λ( k) χ(s)) . (83)
The eigenvalues (82) for k → 0 (where k = | k|) then take the form

λ( k) ≈ 1 - 1 2d k 2 . ( 84
)
The Laplace transformed eigenvalues of the transition matrix then are obtained from canonic representation (83). For k small with (84) the Montroll-Weiss equation takes the form

P (k, s) = (1 -χ(s)) s 1 (1 -χ(s)λ(k)) = s -1 1 + χ(s) (1-χ(s)) µ( k) 2d ≈ s -1 1 + χ(s) (1-χ(s)) k 2 2d . (85)
We employ now the Montroll-Weiss equation (85) as point of departure to derive a diffusional equation governing a CTRW subordinated to a GFPP with the waiting time density described by Laplace transform [START_REF] Repin | Fractional Poisson law[END_REF]. After some elementary manipulations we can rewrite the Montroll-Weiss Eq. (85) in the Fourier-Laplace domain in the form

- ξ α 2d µ( k) Pβ,α (k, s) = (s β + ξ) α Pβ,α (k, s) - Φ(0) β,α (s) -ξ α Pβ,α (k, s) - h 2 2d k2 Pβ,α (h k, s) ≈ (1 + 1 ξ s β ) α -1 Pβ,α (h k, s) + 1 s 1 -(1 + 1 ξ s β ) α , h → 0 (86)
The first line is an exact equation which leads to the second line holding asymptotically for small k = kh. We introduced in the second line a new wave vector k = ( k1 , .., kd ). Then for any finite k we can choose h sufficiently small that µ( k) ≈ k 2 = ( kh) 2 → 0 thus the left-hand side of the second line tends to zero when h → 0. To maintain the second line of relation (86) 'small' requires that ξ(h) → ∞ for h → 0. Hence we can expand (1

+ 1 ξ s β ) α ≈ 1 + α ξ s β to arrive at - h 2 ξ 2dα k2 Pβ,α (h k, s) ≈ s β Pβ,α (h k, s) -s β-1 , (86.

B).

Assuming that h 2 ξ is kept constant when h → 0 leads to the scaling ξ(h) ∼ h -2 → ∞. Eq. (86.B) describes the diffusion limit and reflects the asymptotic GFPP behavior for large dimensionless times ξ 1 β (h)t → ∞ discussed previously (See Eq. ( 78)). The constant h 2 ξ 2dα has physical dimension cm 2 sec -β and can be interpreted as a generalized fractional diffusion constant (recovering for β = 1 the units cm 2 sec -1 of normal diffusion). We observe that the index α enters in Eq. (86.B) only as a scaling parameter.

Transforming the first line of Eq. ( 86) into the causal space-time domain yields (See also general relation (54))

- ξ α 2d L • P (β,α) (t) = 0 D β,α t • P β,α (t) -ξ α P β,α (t) + I ξ α Θ(t) -K (0) β,α (t) . ( 87
)
This relation is written in matrix form where L denotes the Laplacian and I the unit matrix in Z d indicating the initial condition that the walker at t = 0 is sitting in the origin q = 0. Eq. ( 87) describes generalized fractional diffusion in Z d . Its diffusion limit (h → 0 and ξ(h) ∼ h -2 ) is given by the space-time representation of relation (86.B) and yields 14

h 2 ξ 2dα ∆P (β,α) ( x, t) = 0 D β t • P (β,α) ( x, t) -δ (d) ( x) t -β Γ(1 -β) , (87.B)
where δ (d) ( x) denotes Dirac's δ-function in R d and x = h q ∈ R d are the rescaled quasi-continuous coordinates. 0 D β t • P β,1 (q, t) indicates the Riemann-Liouville fractional derivative of order β (See Appendix A). We refer the exact Eq. (87) to as generalized fractional diffusion equation having Eq. (87.B) as 'well-scaled' diffusion limit. Rewriting Eq. (87) in the Fourier-time domain yields

0 D β,α t • Pβ,α (k, t) = ξ α 1 - µ( k) 2d P(β,α) (k, t) + δ q, 0 K (0) β,α (t) -ξ α Θ(t) (88)
showing equivalence with the generalized fractional Kolmogorov-Feller equation ( 70) by accounting for the eigenvalues λ( k) = 1 -µ( k) 2d of W where we used initial condition (δ q, 0 ) = P(t = 0) that the walker sits on the origin at t = 0. The kernels D β,α (t) = L -1 {(s β + ξ) α } and K (0) β,α (t) = L -1 { (s β +ξ) α s } were given explicitly in Eqs. (66) and (68), respectively. Now let us consider the fractional Poisson limit α = 1 with 0 < β < 1. Then Eq. (87.B) takes the form

h 2 ξ 2d ∆P (β,1) ( x, t) = 0 D β t • P (β,1) ( x, t) -δ (d) ( x) t -β Γ(1 -β) 0 < β < 1 (89)
Eqs. ( 89) and (87.B) coincide with the fractional diffusion equation given by Metzler and Klafter [START_REF] Metzler | The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics[END_REF] (Eq. [START_REF] Kutner | The continuous time random walk, still trendy: fifty-year history, state of art and outlook[END_REF] in that paper). The fractional diffusion equation of the diffusion limit (87.B) represents also the asymptotic limit tξ (where ξ → ξ = ξ α ) describing the slow power-law decay of the contribution of the initial condition (See also Eq. ( 86) for |s| small).

Finally we consider the standard Poisson limit α = 1 and β = 1. Then Eq. (87) yields 15

(β,α) ( x, t) = lim h→0 1 (2π) d πh -1 2 -πh -1 2 d k1 • • • πh -1 2 -πh -1
2 d kd P (h k, s)e -i k• x having physical units cm -d . For more details about diffusive or long-wave limits, see e.g. [START_REF] Michelitsch | Fractional Dynamics on Networks and Lattices[END_REF][START_REF] Gorenflo | Mittag-Leffler Waiting Time, Power Laws, Rarefaction, Continuous Time Random Walk, Diffusion Limit[END_REF]. 15 The Fourier-Laplace domain representation of Eq. ( 90) is with (86) given by: -

h 2 ξ 2d k2 P1,1 ( kh, s) = s P1,1 (h k, s) -1. h 2 ξ 2d ∆P (1,1) ( x, t) = ∂ ∂t Θ(t)P (1,1) ( x, t) -δ d ( x)δ(t) (90)
where h 2 ξ 2d indicates the diffusion constant of standard (normal-) diffusion. It is important here that we take into account the Θ(t)-function 'under the derivative'. Then by using

∂ ∂t Θ(t)P (1,1) ( x, t) = δ(t)P (1,1) ( x, t) + Θ(t) ∂ ∂t P (1,1) ( x, t) = δ(t)δ (d) ( x) + Θ(t) ∂ ∂t P (1,1) ( x, t)
with (P 1,1 (q, 0) = δ q, 0 ), Eq. ( 90) recovers then Fick's second law of normal diffusion

h 2 ξ 2d ∆P (1,1) ( x, t) = ∂ ∂t P (1,1) ( x, t). (91)
Finally let us consider the mean squared displacement in generalized fractional diffusion. Then we obtain for the Laplace-Fourier transform of the mean squared displacement σ2 (s) from the general relation σ2

β,α (s) = -∆ k P (h k, s) k=0 = h 2 ñβ,α (s) = h 2 χβ,α (s) s(1 -χβ,α (s)) = 1 s[( s β ξ + 1) α -1] ≈ ξh 2 α s -β-1 , ξ(h) ∼ h -2 → ∞ (92)
where nβ,α (t) indicates the expected number of arrivals of Eq. ( 48). For the GFPP the asymptotic behavior of the expected arrival number nβ,α (t) was also obtained in relation (49). We hence obtain for the mean square displacement the behavior

σ 2 β,α (t) ≈ h 2 ξ αΓ(β + 1) t β , 0 < β ≤ 1 α > 0. (93)
We obtain for 0 < β < 1 and for all α > 0 including fractional Poisson α = 1 sublinear t β power-law behavior corresponding to fat-tailed jump density [START_REF] Shukla | On a generalization of Mittag-Leffler function and its properties[END_REF]. The sublinear power law σ 2 (t) ∼ t β (0 < β < 1) represents the universal limit for tξ 1 β → ∞ of fat-tailed jump PDFs. Such behavior is well known in the literature for anomalous subdiffusion (0 < β < 1) [START_REF] Metzler | The Random Walk's Guide to Anomalous Diffusion : A Fractional Dynamics Approach[END_REF].

In contrast for β = 1 for all α > 0 (Erlang regime which includes the standard Poisson α = 1, β = 1) we obtain normal diffusive behavior with linear increase σ 2 β,α (t) ≈ h 2 ξ α t of the mean squared displacement. This behavior can be attributed to normal diffusion [START_REF] Metzler | The Random Walk's Guide to Anomalous Diffusion : A Fractional Dynamics Approach[END_REF]. The occurrence of universal scaling laws for large tξ 1 β → ∞ reflects the asymptotic universality for time fractional dynamics shown by Gorenflo and Mainardi [START_REF] Gorenflo | The asymptotic universality of the Mittag-Leffler waiting time law in continuous time random walks[END_REF].

Conclusion

We have analyzed a generalization of the Laskin fractional Poisson process, the so called generalized fractional Poisson process introduced for the first time to our knowledge by Polito and Cahoy in 2013 [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF]. We developed the stochastic motions on undirected networks and d-dimensional infinite integer lattices for this renewal process. We derived the probability distribution for n arrivals within time interval [0, t] (Eq. ( 38)) which we call the 'generalized fractional Poisson distribution' (GFPD).

The GFPP is non-Markovian and introduces long-time memory effects. The GFPP contains two index parameters α > 0 and 0 < β ≤ 1 and a characteristic time scale ξ -1 β controlling the dynamic behavior. The GFPP recovers for α = 1 and 0 < β < 1 the Laskin fractional Poisson process, for β = 1, α > 0 the Erlang process, and for α = 1, β = 1 the standard Poisson process. We showed that for α = 1 and 0 < β < 1 the GFPD recovers Laskin's fractional Poisson distribution [START_REF] Laskin | Fractional Poisson process[END_REF] and for α = 1 with β = 1 the standard Poisson distribution. For small dimensionless observation times ξ 1 β t → 0 two distinct regimes emerge, αβ > 1 with 'slow' jump dynamics and αβ < 1 with 'fast' jump dynamics with a transition at αβ = 1 (See Eq. ( 26), and Figures 1 and2).

Based on the Montroll-Weiss CTRW approach we analyzed the GFPP on undirected networks and derived a 'generalized fractional Kolmogorov-Feller equation'. We analyzed for the resulting stochastic motion in d-dimensional infinite lattices the 'well-scaled' diffusion limit of the 'generalized fractional diffusion equation' representing a network variant of the 'generalized fractional Kolmogorov-Feller equation'.

We showed that the generalized Kolmogorov-Feller equation for α = 1 and 0 < β < 1 recovers the fractional Kolmogorov-Feller equation. For α = 1, β = 1 the classical Kolmogorov-Feller equations with Fick's second law as diffusion limit of the standard Poisson process are reproduced.

An essential feature of generalized fractional diffusion is the emergence of subdiffusive behavior σ 2 (t) ∼ t β in the diffusion limit governed by the same type of fractional diffusion equation (87.B) as occurring with the purely fractional Poisson process (α = 1 and 0 < β < 1). Generalized fractional diffusion turns into normal diffusion in the Erlang regime β = 1, α > 0 with σ 2 (t) ∼ t (Eq. ( 93) for β = 1).

We mention that further generalizations of the fractional Poisson process are of interest. For instance we can define a renewal processes which generalize the GFPP waiting time density ( 21) by [43] χγ,β,α (s) =

(ξ α 1 + ξ 2 ) γ {(s β + ξ 1 ) α + ξ 2 } γ 0 < β ≤ 1, α > 0, γ > 0, ξ 1 , ξ 2 > 0, (94)
where χγ,β,α (s)| s=0 = 1 (normalization) and recovers for α = 1 a GFPP with χ β,γ (t) waiting time density, then for α, γ = 1 with 0 < β < 1 the fractional Poisson process and finally for α, β, γ = 1 the standard Poisson process. The advantage of generalizations such as (94) is that it introduces three index parameters α, β, γ, and two time scale parameters, ξ 1,2 (ξ 1 having physical dimension sec -β , and ξ 2 is of dimension sec -βα ). Such generalizations offer a great flexibility to fit various real-world situations such as occurring in various problems of anomalous diffusion, e.g. [START_REF] Goychuk | Fractional diffusion modeling of ion channel gating[END_REF] where power law asymptotic features occur. The analysis of the process defined by (94) is beyond the scope of the present article and will be presented in a follow-up paper [43].

Our results suggest that the GFPP and further generalizations have huge potentials of applications in anomalous diffusion and transport phenomena including turbulence, non-Markovian dynamics on networks, and in the dynamics of complex systems.

Acknowledgment

We thank Federico Polito for his valuable comments and to have drawn our attention to Ref. [START_REF] Cahoy | Renewal processes based on generalized Mittag-Leffler waiting times[END_REF].

A Appendix: Laplace transforms of causal distributions and fractional operators

Here we discuss some properties of Laplace transforms that we use in the present paper. The Laplace transform of the mth time derivative of convolution of two causal functions h(t)Θ(t) and k(t)Θ(t) is given by

L d m dt m t 0 h(v)k(t -v)dv = ∞ -∞ ∞ -∞ dudvΘ(u)Θ(v)h(u)k(v) ∞ 0- e -st d m dt m δ(t -u -v)ds = s m ∞ -∞ ∞ -∞
dudvΘ(u)Θ(v)h(u)k(v)e -su e -sv = s m k(s) h(s), m = 0, 1, 2, . . . ∈ N 0 .

(95)

The last relation is crucial to define in our analysis 'good' fractional derivatives and integrals.

Let us consider here property (3) in more details, namely

L d m dt m (Θ(t)f (t)) = s m f (s) = L Θ(t) d m dt m f (t) + m m=1 m! (m -k)!k! ∞ 0 - e -st δ k-1 (t) d m-k dt m-k f (t)dt (96)
where is denoted δ s (t) = d s dt s δ(t). Then integrating by parts these terms and using d n dt n (e -st y(t)) = e -st ( d dts) n y(t) we arrive at

∞ 0 - e -st δ k-1 (t) d m-k dt m-k f (t)dt = ∞ 0 - δ(t)e -st (-1) k-1 d dt -s k-1 d m-k dt m-k f (t)dt = s - d dt k-1 d m-k dt m-k f (t) t=0 . (97)
We emphasize that the δ(t)-function is entirely captured thus in (97) occur only initial values at t = 0 (not at t = 0 -). Then by substituting (97) into (96) yields the well known relation

L Θ(t) d m dt m f (t) = s m f (s) -s - d dt -1 s m - d m dt m f (t) t=0 . ( 98
)
We often will meet in our analysis the inverse Laplace transform L -1 {s γ } ,(γ > 0) which we evaluate here in details. To this end it is useful to introduce the ceiling function γ indicating the smallest integer greater or equal to γ, i.e. for γ / ∈ N we have -1 < γγ < 0 and γ = γ if γ ∈ N. We can represent L -1 {s γ } as Fourier integral, namely We see that for integer γ ∈ N 0 integral (99) yields d γ dt γ δ(t) i.e. the convolution operator of the integer order derivative γ. Let us now consider the case γ > 0 and γ / ∈ N. Then we take into account that

L -1 {s γ } = 1 2πi
s γ = s γ s γ-γ =            s γ Γ( γ -γ) ∞ -∞ Θ(t)e -sτ τ γ -γ-1 dτ, γ / ∈ N s γ ∞ -∞ e -sτ δ(τ )dτ, γ ∈ N ( 100
)
where it is important to notice that we write s γ in the form (100) since the τ -integral converges as 0 < γγ < 1 and the Heaviside Θ(t)-function guarantees that the integral starts from τ = 0. Then we further use the simple relation s γ e -sτ = (-1) γ d γ dτ γ e -sτ . (101)

Then by substituting this relation into integral (100) yields

s γ =            (-1) γ Γ( γ -γ) ∞ -∞ Θ(t)τ γ -γ-1 d γ dτ γ e -sτ dτ, γ / ∈ N (-1) γ ∞ -∞ δ(τ ) d γ dτ γ e -sτ dτ, γ ∈ N (σ + iω) γ =            ∞ -∞
e -(σ+iω)τ d γ dτ γ Θ(τ )

τ γ -γ-1 Γ( γ -γ) dτ, γ / ∈ N ∞ -∞ e -(σ+iω)τ d γ dτ γ δ(τ )dτ, γ ∈ N (102)
where in the second line γ partial integrations have been performed and the boundary terms

d k dτ k (Θ(τ )τ γ -γ-1 . . .)
∞ -∞ are all vanishing. In this relation is important that the Heaviside Θ(τ ) is included within the differentiation. With representations in (102) the Laplace inversion (99) then yields

L -1 {s γ } =                      e σt ∞ -∞ ∞ -∞ e iω(t-τ )
2π dω =δ(t-τ )

e -στ d γ dτ γ Θ(τ )

τ γ -γ-1 Γ( γ -γ) dτ, γ / ∈ N e σt ∞ -∞ ∞ -∞
e iω(t-τ ) 2π dω =δ(t-τ ) e -στ d γ dτ γ δ(τ )dτ, γ ∈ N

L -1 {s γ } =          d γ dt γ Θ(t) t γ -γ-1 Γ( γ -γ) , γ / ∈ N d γ dt γ δ(t), γ ∈ N. (103)
For fractional γ / ∈ N we identify the kernel of the Riemann-Liouville fractional derivative of order γ, e.g. [START_REF] Oldham | The fractional calculus[END_REF][START_REF] Miller | An Introduction to the Fractional Calculus and Fractional Differential Equations[END_REF] and many others, and for γ ∈ N we get the distributional kernel producing integer derivatives of order γ. The Riemann-Liouville fractional derivative acts on causal functions f (t) = Θ(t)f (t) as

L -1 {s γ } • f (t)Θ(t) =: 0 D γ t f (t) = d γ dt γ ∞ -∞ Θ(t -τ ) (t -τ ) γ -γ-1 Γ( γ -γ) f (τ )Θ(t)dτ 0 D γ t f (t) = 1 Γ( γ -γ) d γ dτ γ t 0 (t -τ ) γ -γ-1 f (τ )dτ. ( 104
)
We emphasize that (103) requires causality and definition of the Laplace transform (1) that captures the entire non-zero contributions of the causal distribution with lower integration limit 0 -= lim →0+where all boundary terms at 0 -are vanishing. The same result (103) for the Riemann-Liouville fractional derivative kernel is obtained in the following short way (105)

In the same way as above one can derive the Laplace inversion of s -γ , γ > 0 which yields

L -1 {s -γ } = Θ(t) t γ-1 Γ(γ) , γ > 0 (106)
as a fractional generalization of integration and indeed can be identified with the Riemann-Liouville fractional integral kernel of order γ, e.g. [START_REF] Oldham | The fractional calculus[END_REF][START_REF] Miller | An Introduction to the Fractional Calculus and Fractional Differential Equations[END_REF].

1 β.

 1 For β = 1, (14) recovers the exponential jump density χ 1 (t) = ξe -ξt of the standard Poisson process. The survival probability in the fractional Poisson process is then of the form of a Mittag-Leffler function Φ (0)

Figure 1 :Figure 2 :

 12 Figure 1: (a) Jump density χ β,α (t) and (b) cumulated probability Ψ β,α (t) as a function of t. We explore the case α = 2.0 for different values 0 < β ≤ 1 codified in the colorbar. The results were obtained numerically using ξ = 1 and Eqs. (24) and (27). We depict with dashed lines the results for β = 0.5 (αβ = 1).

1 β

 1 and α = 1 recovers the well-known result of the standard

  55) where we have always a unique eigenvalue λ 1 = 1 10 |λ m | < 1 (m = 2, . . . , N ) [38] 11 . In Eq. (55) we have introduced the right-and left right eigenvectors |Φ j and Φj | of W, respectively. The first term that corresponds to λ 1 = 1 in |Φ 1 Φ1 | indicates the stationary distribution. We can then write the canonic representation of (54) in the form P(s) = N m=1 P (m, s)|Φ m Φm | (56)

1 β

 1 → ∞ of the generalized fractional diffusion equation (87) with the memory term -δ(d) ( x) t -β Γ(1-β)

2 j

 2 denotes the Laplace operator of the R d , x = h p ∈ R d indicate rescaled quasicontinuous spatial coordinates, and P (β,α) ( q, t) = h d P (β,α) ( x, t) with the transition probability kernel P

e

  iωt (σ + iω) γ dω, s = σ + iω, σ = (s) > 0. (99)

L - 1

 1 {s γ } = L -1 {s γ s γ-γ } = e σt σ + e iωt (σ + iω) γ-γ = e σt σ + d dt γ e -σt Θ(t) t γ -γ-1 ( γγ -1)! = d γ dt γ Θ(t) t γ -γ-1 Γ( γγ).

We denote as 0-(. . .) = lim →0+ -(. . .).

We notice the distributional representation (See Appendix A):

We abbreviate probability density functions or short 'densities' by 'PDF'.

We refer this density also to as n-jump density.

This series is converging for |v| ≤ 1 when s = 0 and for |v| < 1 when s = 0 as consequence of | χ(s)| ≤ 1 (equality for s = 0).

We employ throughout this paper as equivalent notations ζ! = Γ(ζ + 1) for the Γ-function.

The large time behavior emerges for ξ → ∞ (or equivalently for small |s| → 0) in[START_REF] Repin | Fractional Poisson law[END_REF].

Where -Γ(-β) = β -1 Γ(1β) > 0.

See Eq. (8).

Where we introduce d β,α (t) = L -1 s -αβ (s β + ξ) α .