The role of fractional time-derivative operators on anomalous diffusion

TL;DR

This paper explores how new fractional time-derivative operators with non-singular kernels can model complex anomalous diffusion processes, revealing diverse diffusive behaviors, non-Gaussian distributions, and connections to stochastic resetting.

Contribution

It introduces and analyzes the use of Caputo-Fabrizio and Atangana-Baleanu fractional operators in diffusion equations, expanding modeling capabilities for anomalous diffusion.

Findings

Derived exact probability distributions for various diffusion regimes

Identified crossovers between different diffusive behaviors

Linked new fractional operators to stochastic resetting processes

Abstract

The generalized diffusion equations with fractional order derivatives have shown be quite efficient to describe the diffusion in complex systems, with the advantage of producing exact expressions for the underlying diffusive properties. Recently, researchers have proposed different fractional-time operators (namely: the Caputo-Fabrizio and Atangana-Baleanu) which, differently from the well-known Riemann-Liouville operator, are defined by non-singular memory kernels. Here we proposed to use these new operators to generalize the usual diffusion equation. By analyzing the corresponding fractional diffusion equations within the continuous time random walk framework, we obtained waiting time distributions characterized by exponential, stretched exponential, and power-law functions, as well as a crossover between two behaviors. For the mean square displacement, we found crossovers between usual and confined diffusion, and between usual and sub-diffusion. We obtained the exact expressions for the probability distributions, where non-Gaussian and stationary distributions emerged. This former feature is remarkable because the fractional diffusion equation is solved without external forces and subjected to the free diffusion boundary conditions. We have further shown that these new fractional diffusion equations are related to diffusive processes with stochastic resetting, and to fractional diffusion equations with derivatives of distributed order. Thus, our results show that these new operators are a simple and efficient way for incorporating different structural aspects into the system, opening new possibilities for modeling and investigating anomalous diffusive processes.

Tables

Table 1. Table 1: Summary of the changes caused by the different fractional operators on the diffusion equation ( 1 ).

Fractional Operator	Waiting time distribution	Mean square displacement	Probability distribution
Riemann- Liouville	power-law	power-law and scale-invariant	non-Gaussian
Caputo- Fabrizio	exponential	crossover from usual to confined diffusion	crossover from Gaussian to non-Gaussian with steady state
Atangana- Baleanu	crossover from stretched exp. to power-law	crossover from usual to sub-diffusion	crossover from Gaussian to non-Gaussian

Equations

$$\frac{\partial }{\partial t}\rho (x,t)={}_0{\mathcal{D}}_t^{1-\alpha }\left(\mathcal{L}\{\rho (x,t)\}\right),$$

$$\mathcal{L}\{\rho (x,t)\}=D\frac{{\partial }^2}{\partial x^2}\rho (x,t)-\frac{\partial }{\partial x}[F(x,t)\rho (x,t)],$$

$${}_0{\mathcal{D}}_t^{1-\alpha }\rho (x,t)=\frac{1}{\Gamma \left(\alpha \right)}\frac{d}{dt}{\int }_0^t\frac{\rho (x,t^{\prime })}{(t-t^{\prime }{)}^{1-\alpha }}d{t}^{\prime },$$

$$\frac{\partial }{\partial t}\rho (x,t)=D{\mathcal{F}}_t^{\alpha }\left(\frac{{\partial }^2}{\partial x^2}\rho (x,t)\right),$$

$${\mathcal{F}}_t^{\alpha }\rho (x,t)=\frac{\partial }{\partial t}{\int }_0^t\rho (x,t^{\prime })\mathcal{K}(t-t^{\prime })d{t}^{\prime },$$

$$\mathcal{K}(t)=\frac{{\left(t/\tau \right)}^{\alpha -1}}{\Gamma \left(\alpha \right)},$$

$$\mathcal{K}(t)$$

$$\mathcal{K}(t)=b{E}_{\alpha }\left[-\frac{\alpha }{1-\alpha }{\left(\frac{t}{\tau }\right)}^{\alpha }\right],$$

$$\rho (k,s)=\frac{\phi (k)}{s+sD\mathcal{K}(s)k^2},$$

$$\omega (s)=\frac{\mathcal{K}(s)/{\tau }_c}{1+\mathcal{K}(s)/{\tau }_c}$$

$$\omega (t)$$

$$\omega (t)=\frac{1}{{\tau }_c}{e}^{-t/{\tau }_c}.$$

$$\omega (t)=\frac{1}{{\tau }_c}{\left(\frac{t}{\tau }\right)}^{\alpha -1}{E}_{\alpha ,\alpha }\left[-\left(\frac{\tau }{{\tau }_c}\right){\left(\frac{t}{\tau }\right)}^{\alpha }\right],$$

$$\omega (t)=\frac{\xi b\gamma }{\pi {\tau }_c}\sin \left(\pi \gamma \right){\int }_0^{\infty }d\eta \frac{{\eta }^{\gamma }{e}^{-\eta b(1-\alpha )t/{\tau }_c}}{{\left(1-\eta \right)}^2+2\xi \left(1-\eta \right){\eta }^{\gamma }\cos \left(\pi \gamma \right)+{\eta }^{2\gamma }},$$

$$\omega (t)\sim E_{\alpha }\left[-\frac{\alpha }{1-\alpha }{\left(\frac{t}{\tau }\right)}^{\alpha }\right]\sim e^{-\frac{\alpha }{\Gamma (\alpha )(1-\alpha )}{\left(\frac{t}{\tau }\right)}^{\alpha }},$$

$$\mathcal{K}(s)={\tau }^{1-\alpha }{s}^{-\alpha }$$

$$\rho (x,t)={\int }_{-\infty }^{\infty }d{x}^{\prime }\mathcal{G}(x-x^{\prime },t)\phi (x^{\prime }),$$

$$\mathcal{G}(x,t)=\frac{1}{2|x|}{\text{mbox}H}_{1,1}^{1,0}\left[\frac{|x|}{\sqrt{D\tau {\left(t/\tau \right)}^{\alpha }}}|{}_{\left(1,1\right)}^{\left(1,\frac{\alpha }{2}\right)}\right].$$

$${\left(\Delta x\right)}^2=⟨{\left(x-⟨x⟩\right)}^2⟩=\frac{2D\tau }{\Gamma (1+\alpha )}{\left(\frac{t}{\tau }\right)}^{\alpha },$$

$$\mathcal{K}(s)$$

$$\rho (k,s)=\frac{\left(s+\frac{\alpha }{(1-\alpha )\tau }\right)\phi (k)}{s\left(s+\frac{\alpha }{(1-\alpha )\tau }+Db{k}^2\right)}.$$

$$\rho (x,t)$$

$$\mathcal{G}(x,t)=\frac{e^{-\frac{\alpha t}{(1-\alpha )\tau }}}{\sqrt{4\pi Dbt}}{e}^{-\frac{x^2}{4Dbt}},$$

$$\rho (x,t\to \infty )\sim {\int }_{-\infty }^{\infty }d{x}^{\prime }\phi (x^{\prime })e^{-\frac{\alpha |x-x^{\prime }|}{(1-\alpha )Db}},$$

$${\left(\Delta x\right)}^2=\frac{2Db\tau (1-\alpha )}{\alpha }\left(1-e^{-\frac{\alpha t}{(1-\alpha )\tau }}\right),$$

$$\frac{\partial }{\partial t}\rho (x,t)=Db\frac{{\partial }^2}{\partial x^2}\rho (x,t)-\frac{Db\alpha }{(1-\alpha )\tau }{\int }_0^{t}d{t}^{\prime }{e}^{-\frac{\alpha (t-t^{\prime })}{(1-\alpha )\tau }}\frac{{\partial }^2}{\partial x^2}\rho (x,t^{\prime }).$$

$$Db{\int }_0^{t}d{t}^{\prime }{e}^{-\frac{\alpha (t-t^{\prime })}{(1-\alpha )\tau }}\frac{{\partial }^2}{\partial x^2}\rho (x,t^{\prime })=\rho (x,t)-\phi (x),$$

$$\frac{\partial }{\partial t}\rho (x,t)=Db\frac{{\partial }^2}{\partial x^2}\rho (x,t)-\frac{\alpha }{(1-\alpha )\tau }\left[\rho (x,t)-\phi (x)\right].$$

$$\mathcal{J}(x,t)$$

$$\mathcal{K}(s)=\frac{b{s}^{\alpha -1}}{\left(s^{\alpha }+\frac{\alpha }{(1-\alpha ){\tau }^{\alpha }}\right)},$$

Full text

I Introduction

The random walk concept is one of the broadest and versatile paradigms to deal with statistical fluctuations. The term “random walk” was coined in 1905 by Pearson Pearson , but the fundamental relationship between this concept and the usual diffusion equation was reported earlier in the seminal works of Rayleigh Rayleigh1 ; Rayleigh2 ; Rayleigh3 in sound theory, Bachelier Bachelier in economics, Einstein Einstein and Smoluchowski Smoluchowski in the Brownian motion theory. Due to this intrinsic relation, the usual random walk is characterized by Gaussian, Markovian, and ergodic properties, which lead to a linear time dependence of the mean square displacement, $(\Delta x)^{2}\sim t$. The versatility of this concept relies on the possibility of generalizations and extensions to describe systems with one or more characteristics of anomalous diffusion: non-Gaussian distributions; long-range memory effects (non-Markovian); non-ergodicity; divergent mean square displacement (Lévy walks); and nonlinear mean square displacement, $(\Delta x)^{2}\sim t^{\alpha}$ (sub-diffusion: $\alpha<1$, superdiffusion: $\alpha>1$, confined or saturated diffusion: $\alpha=0$).

In the context of generalizations, the first landmark is the work proposed in 1965 by Montroll and Weiss Montroll , in which they introduce the continuous time random walk concept (see Ref. Kutner for a general overview). This framework is characterized by a joint distribution of jump length and waiting time $\psi(x,t)$, where $\lambda(x)=\int_{0}^{\infty}dt\psi(x,t)$ is the jump length distribution and $\omega(t)=\int_{-\infty}^{\infty}dx\psi(x,t)$ is the waiting time distribution. Subsequently, a connection with a generalized master equation is proposed to discuss memory effects in the continuous time random walk MontrollScher ; Kenkre ; Shlesinger ; the waiting time distribution $\omega(t)$ is strictly related to the memory kernel of the generalized master equation. Moreover, natural extensions of both random walk and continuous time random walk were proposed to study transport properties in systems with structural complexity, such as disordered, random, and fractal environments Havlin2000 .

The second landmark, “a modern era of the continuous time random walk” according to Kutner and Masoliver Kutner , is the development of the intrinsic relationship between this formalism and the fractional diffusion equations. Among the seminal works, we have found a simple mention in the work of Klafter et al. Klafter1987 on the possibility of having a fractional diffusion equation to describe anomalous transport. However, were Hilfer et al. Hilfer1995 that, ninety years after Einstein’s work Einstein , established a rigorous and precise connection between the continuous time random walk and the fractional master equation as well as with the fractional diffusion equation (a special case of the former). This result was later extended by Compte Compte1996 in the long-time limit, where it is shown that any decoupled continuous time random walk having no characteristic scale of time or space (power-law memories) corresponds to a time- or space-fractional diffusion equation, respectively with Riemann-Liouville time derivative or Riesz space derivate. The subsequent success and development of the fractional approach are well documented in two review articles by Meztler and Klafter 15 ; Restaurant , and in several articles by Barkai Barkai1 ; Barkai2 ; Barkai3 ; Barkai4 .

Since memory effects underlying a continuous time random walk are implicitly considered by the differential operators, the versatility of the fractional formalism is mainly related to two remarkable features. First, this formalism handles very well the physical requirements of a system by dealing boundary conditions and external forces in a simple manner. Second, it takes advantage of traditional tools from mathematical physics and statistics for obtaining exact expressions to describe complex systems with anomalous behaviors. For instance, in the following fractional differential equation

[TABLE]

where

[TABLE]

the nonusual relaxation can be associated with a continuous time random walk where the waiting time distribution is a power-law. This process is also strictly related to the Riemann-Liouville fractional operator BarkaiMetzlerKlafter

[TABLE]

where $0<\alpha<1$ is the fractional order exponent (or the anomalous exponent), a quantity that can be interpreted as an index of memory in empirical systems Du2013 . The fractional operator is also responsible for introducing a nonlinear time dependence in the mean square displacement of the system 15 . Thus, a large class of complex phenomena can be effectively described by extending the standard differential operator to a non-integer order zahran ; zahran2 ; zahran3 ; villa ; burioni ; maringa ; iomin5 ; iomin6 ; maringa2 ; ChemPhys ; indeed, as pointed out by West WestBook , the fractional calculus provides a suitable framework to deal with complex systems.

Recently, researchers have made and promoted remarkable progress toward improving experimental techniques for investigating diffusive processes, mainly illustrated by the developments in the single-particle tracking technique Saxton1997 ; Wirtz2009 ; Gal2013 ; Hoze2017 . Such improvements yield novel insights into transport properties of biological systems Franosch2013 ; Manzo2015 ; Shen2017 and nanomaterials Zagato2014 ; Karger2012 ; Karger2016 , where the high-resolution of the experiments has found different diffusive behaviors depending on the time scale. In this context, an important question is whether other forms of fractional differential operators (replacing the Riemann-Liouville one) such as those recently-proposed with non-singular kernels nonsingular1 ; singularJordan ; GA ; nonsingular2 ; AB ; nonsingular3 are suitable to describe the aforementioned situations. To answer this question, we investigate an one-dimensional diffusive process described by the fractional diffusion equation

[TABLE]

where $D$ is the generalized diffusion coefficient. This equation is also subjected to the free diffusion boundary conditions $\rho(\pm\infty,t)=0$ and to the initial condition $\rho(x,0)=\varphi(x)$.

The fractional operator in Eq. (3) is defined as

[TABLE]

in order to consider situations with singular and non-singular kernels in a unified way. It is worth noting that ${\cal{K}}(t)=\delta(t)$ recovers the usual diffusion equation. Here we consider three different forms for the kernel ${\cal{K}}(t)$. The first one is

[TABLE]

which corresponds to the well-known Riemann-Liouville fractional operator Livro1 for $0<\alpha<1$. The second one is

[TABLE]

which corresponds to the fractional operator of Caputo-Fabrizio nonsingular1 . In the context of the diffusion equation, the use of this operator is related with a diffusion equation combined with first-order kinetics Sano .

Finally, the third one is

[TABLE]

where $E_{\alpha}(\dots)$ is the Mittag-Leffler function Livro1 . This kernel corresponds to the fractional operator of Atangana-Baleanu AB . Further possibilities for the kernel ${\cal{K}}(t)$ are discussed by Gómez-Aguilar et al. GA . We observe that the Riemann-Liouville operator have a singularity at the origin ($t=0$), while the recently-proposed Caputo-Fabrizio and Atangana-Baleanu are non-singular operators nonsingular1 ; singularJordan ; GA ; nonsingular2 ; AB ; nonsingular3 . In the previous definitions, the parameter $b$ is a normalization constant, $\alpha$ is the fractional order exponent, and $\tau$ is a characteristic time that controls the shape of the kernels.

Our main goal here is to verify how these different fractional operators modify the fractional diffusion equation (1) and what are the effect of these choices on the underlying diffusive properties of a system modeled by this equation. The rest of this manuscript is organized as follows. In Section II, we investigate general solutions and processes related with Eq. (3) when considering different choices (singular and non-singular) for the kernel ${\cal{K}}(t)$. In Section III, we present a summary of our results and some concluding remarks.

II Diffusion and Fractional Operators

We start by noting that the solution of the fractional diffusion equation (3) in the Fourier-Laplace space is

[TABLE]

where $\rho(k,s)$ is the Fourier-Laplace transformation of the probability distribution $\rho(x,t)$. This result can be related to different situations depending on the choice of the kernel ${\cal{K}}(s)$.

Within the continuous time random walk formalism and by following the works of Meztler and Klafter 15 , we can show that the waiting time $\omega(t)$ and the jump $\lambda(x)$ probability distribution associated with Eq. (3) are (in the Laplace and Fourier spaces)

[TABLE]

and $\lambda(k)=1-k^{2}D\tau_{c}\,$, where $\tau_{c}$ is a characteristic waiting time of the underlying continuous time random walk. We observe that the jump probability distribution is characterized by a Gaussian asymptotic behavior [$\lambda(x)\sim e^{-x^{2}/(2D\tau_{c})^{2}}$] and thus has a finite characteristic jump length, regardless of the choice for the kernel ${\cal{K}}(t)$. On the other hand, the inverse Laplace transform of the waiting time distribution is given by

[TABLE]

yielding different situations that depends on ${\cal{K}}(t)$.

The choice ${\cal{K}}(t)=\delta(t)$ leads to usual diffusion and an exponential distribution for the waiting times

[TABLE]

For the fractional operator of Riemann-Liouville, we find

[TABLE]

where $E_{\alpha,\alpha}(\dots)$ is the generalized Mittag-Leffler function Livro1 whose asymptotic behavior is described by a power-law, $\omega(t)\sim 1/t^{1+\alpha}$ for $t\rightarrow\infty$. For the Atangana-Baleanu operator [Eq. (7)], the waiting time distribution is given by

[TABLE]

where $\gamma=1-\alpha$ and $\xi=(\alpha b)/(\tau_{c}\tau^{\alpha})$. This expression is very interesting because for small times we have a stretched exponential, that is,

[TABLE]

while for long times we have the same power-law behavior of the Riemann-Liouville operator. Thus, the Atangana-Baleanu operator yields a crossover between a stretched exponential and a power-law distribution.

In the case of the Caputo-Fabrizio operator, the connection with the continuous time random walk is more complex and not compatible with its standard interpretation. As we shall discuss later on, the diffusion equation associated with this operator is connected to a diffusive process with stochastic resetting Mendez2016 ; Shkilev , where the waiting time distribution is exponential.

Figure 1 depicts the behavior of the waiting time distribution $\omega(t)$ for the different kernels previously-discussed. For long times, we confirm that the operators of Riemann-Liouville and Atangana-Baleanu yield the same power-law decay for $\omega(t)$. We further note that the Atangana-Baleanu operator yields a non-divergent $\omega(t)$, an interesting feature that is not observed for the singular kernel of Riemann-Liouville.

We now focus on finding the formal solutions for the fractional diffusion equation (1) when considering the three different fractional operators. These solutions are obtained by performing the inverse of Fourier and Laplace transforms of the $\rho(k,s)$ expressed in Eq. (8), where the Laplace transform of the kernel ${\cal K}(s)$ appears. In the well-known case of the Riemann-Liouville operator 15 , we have

[TABLE]

and consequently

[TABLE]

where the Green function is

[TABLE]

Here ${\mbox{\large{H}}}(\dots)$ stands for the Fox H-function Fox . Having found the probability distribution, we can show that the mean square displacement is

[TABLE]

which corresponds to the typical case of anomalous diffusion, where $\alpha<1$ represents sub-diffusion and $\alpha\to 1$ recovers the usual diffusion. The time-dependent behavior of a typical probability distribution $\rho(x,t)$ (with $\alpha=1/2$) is shown in Fig. 2(a). We observe that the Riemann-Liouville operator leads to a tent-shaped distribution, whose tails are longer than the Gaussian distribution of the usual diffusion [Fig. 2(d)]. Figure 3 shows the corresponding behavior for mean square displacement of Eq. (18), which is a power-law function of the time $t$ with an exponent $\alpha$.

For the Caputo-Fabrizio operator, the Laplace transform of the kernel in Eq. (6) is

[TABLE]

which substituted into Eq. (8) yields

[TABLE]

By performing the inverse Fourier and Laplace transforms, we have

[TABLE]

where the Green function is

[TABLE]

A typical shape of this distribution is shown in Fig. 2(b). We observe that this distribution is very similar to a Gaussian for small times, and exhibits a tent-shape to long times. However, differently from the distribution obtained for the Riemann-Liouville operator [Eqs. (16) and (17)], the distribution obtained from Eqs. (21) and (22) displays a stationary behavior for $t\to\infty$, that is,

[TABLE]

a result that corresponds to confined diffusion. Figure 2(b) also shows this stationary solution (dashed line); in particular, we observe that the shape of $\rho(x,t)$ is practically constant for $t\gtrsim 5$ in that case. This behavior also appears in the mean square displacement

[TABLE]

which behaves linearly in time for small times and saturates in ${2D\,b\,\tau(1-\alpha)}/{\alpha}$ for long times. Figure (3) illustrates this crossover, a common feature of systems where diffusion is confined or hindered Manzo2015 ; Mo2015 . In particular, the same crossover between usual and confined diffusion is observed in simulations of diffusion with immobile obstacles or obstacles moving according to an Ornstein-Uhlenbeck process Berry2014 ; Spakowitz2017 .

An intriguing feature of the diffusion equation with the Caputo-Fabrizio fractional operator is that it can be related to a diffusion with stochastic resetting Resetting . Indeed, we find out that the fractional diffusion equation (3) with the kernel of Eq. (6) can be rewritten as

[TABLE]

Also, by integrating both sides of the fractional diffusion equation (3), we obtain

[TABLE]

which after substituting into Eq. (25) yields

[TABLE]

Equation (27) with $\varphi(x)=\delta(x-x_{0})$ is the same obtained by Evans and Majumdar Resetting when studying a random walker whose position is redefined to the position $x_{0}$ with a rate $r=\frac{\alpha}{(1-\alpha)\tau}$. Thus, the fractional exponent $\alpha$ (as well as the parameter $\tau$) in the fractional diffusion equation of Caputo-Fabrizio can be related to a well-defined physical quantity (resetting rate).

Also, the mean square displacement of Eq. (24) is analogous to results obtained from a random walk description of a diffusive process with stochastic resetting, subjected to an exponential waiting time distribution Mendez2016 ; Shkilev . As discussed in these works, a suitable continuous time random walk formulation is established by considering a density of particles ${\cal{J}}(x,t)$ whose dynamics is governed by

[TABLE]

when particles start the random walk at the origin ($x=0$) with $\rho(x,t)=\int_{0}^{t}\Phi(t^{\prime}){\cal{J}}(x,t-t^{\prime})dt^{\prime}$ and $\Phi(t)=\int_{t}^{\infty}\omega(t^{\prime})dt^{\prime}$. In Eq. (28), $r$ is a resetting rate, $\psi(x,t)$ is joint distribution of jump length and waiting time, $\lambda(x)=\int_{0}^{\infty}dt\psi(x,t)$ is the jump length distribution, and $\omega(t)=\int_{-\infty}^{\infty}dx\psi(x,t)$ is waiting time distribution. By considering $\lambda(x)$ Gaussian and $\omega(t)$ exponentially distributed, we can show that this formalism leads to Eq. (27). It is worth remarking that by comparison with this framework, we can infer that the diffusion equation with the Caputo-Fabrizio operator leads to the same waiting time distribution of the usual diffusion, that is, an exponential.

Finally, for the Atangana-Baleanu operator, the Laplace transform of the kernel in Eq. (7) is

[TABLE]

which substituted into Eq. (8) yields

[TABLE]

the solution for the fractional diffusion equation (3) in the Fourier-Laplace space. By evaluating the inverse Fourier and Laplace transforms, we obtain

[TABLE]

where the Green function is

[TABLE]

Once again, ${\mbox{\large{H}}}(\dots)$ stands for the Fox H-function Fox . We can also show that for $|x|\rightarrow\infty$, Eq. (32) is approximated by

[TABLE]

where

[TABLE]

A typical behavior for the distribution $\rho(x,t)$ for this operator is shown in Fig. 2(c). Similarly to the Caputo-Fabrizio operator, the profile of $\rho(x,t)$ resembles a Gaussian for small times, while exhibits a tent-shape for long times. However, the distribution does not have a stationary solution for the Atangana-Baleanu operator. This crossover between two behaviors for $\rho(x,t)$ is also present in Eq. (34), and can be better quantified by analyzing the mean square displacement. For this operator, we have

[TABLE]

where $E_{\alpha,2}(\dots)$ is the generalized Mittag-Leffler function Livro1 . By considering the asymptotic limits of this function, we can show that $\left(\Delta x\right)^{2}\sim t$ for small times, and $\left(\Delta x\right)^{2}\sim t^{1-\alpha}$ for long times.

This crossover between usual and sub-diffusion is present in several biological systems Skalski2000 ; Sawadac2001 ; Mieruszynsky2015 ; ribeiro2014investigating ; Ribeiro2016 and is also illustrated in Fig. 3 for $\alpha=1/2$. A similar situation appears in simulations of diffusion with obstacles moving according to a usual random walk Berry2014 ; Spakowitz2017 , where the same crossover between usual and sub-diffusion with $\alpha=1/2$ is observed. It is worth mentioning that crossovers between diffusive regimes can also be described by generalized Langevin equations Tateishi2012 and fractional (with the Riemann-Liouville operator) Kramers equations Barkai1 , among other approaches DifferentRegimes2 ; DifferentRegimes1 . In particular, the usual Langevin equation Bian2016 predicts a crossover between ballistic and usual diffusion, which has been experimentally observed only in 2011 HuangRaizen2011 . However, the diffusion equation in terms of these new operators lead to these crossovers without explicitly considering external forces, inertial effects, and reaction terms.

The fractional diffusion equation with the Atangana-Baleanu operator can be further related to fractional derivatives of distributed order as proposed by Caputo caputo1995mean ; caputo2001distributed and worked out in Refs. DifferentRegimes2 ; DifferentRegimes1 , that is,

[TABLE]

where $w(\nu)$ is the distribution of the fractional order exponent $\nu$ and

[TABLE]

is the fractional time derivative of Caputo. Indeed, by substituting the kernel of Eq. (7) into Eq. (3) and taking the Laplace transform, we have

[TABLE]

which can be rewritten as

[TABLE]

By calculating the inverse Laplace transform of the previous equation, we find

[TABLE]

which can also be written as

[TABLE]

We note that Eq. (41) is a special case of Eq. (36) with $w(\nu)=\delta(\nu-1)+\frac{\alpha}{(1-\alpha)\tau^{\alpha}}\delta(\nu+\alpha-1)$. Analogously to results reported here, the solutions of Eq. (41) are also characterized by two diffusive regimes DifferentRegimes2 ; DifferentRegimes1 .

III Discussion and Conclusions

We presented a detailed investigation of the changes in the fractional diffusion equation when the well-established Riemann-Liouville operator is replaced by the recently-proposed operators of Caputo-Fabrizio and Atangana-Baleanu. These changes are summarized in Table 1. Within the context of the continuous time random walk, we verified that these new fractional operators modify the behavior of the waiting time distribution. In the Caputo-Fabrizio case, we found that the waiting time distribution is described by an exponential distribution; while the Atangana-Baleanu operator yields a distribution that decays as a stretched exponential for small times and as a power-law (with the same exponent of the Riemann-Liouville operator) for long times.

We obtained the exact solutions of the fractional diffusion equation and the time dependence of the mean square displacement when considering these different fractional operators. Our results reveal that these new operators lead to non-Gaussian distributions and different diffusive regimes depending on the time scale. For the Caputo-Fabrizio operator, the probability distribution $\rho(x,t)$ displays a stationary state as well as saturated diffusion for long times. This is a remarkable feature because the fractional diffusion equation is solved without external forces and subjected to the free diffusion boundary conditions. For the Atangana-Baleanu operator, we found a crossover between two diffusive regimes: a usual for small times and a sub-diffusive for long times, a feature observed in several empirical systems.

By properly manipulating the fractional diffusion equations, we demonstrated that the results obtained with these new fractional operators could be connected with other diffusive models. The fractional diffusion equation with the Caputo-Fabrizio operator recovers a diffusive process with stochastic resetting, where the fractional order exponent is directly related to the resetting rate. Also, the equation with the Atangana-Baleanu operator can be associated with a fractional diffusion equation with derivatives of distributed order. Our results thus show that these new fractional operators are a simple and efficient way for incorporating different structural aspects into the system, opening new possibilities for modeling and investigating the interplay of different physical mechanisms of anomalous diffusion.

Data Availability

There no datasets involved in this work. All graphs were produced by using matplotlib, an open source plotting library for the Python programming language.

Competing Interests

We have no competing interests.

Authors’ Contributions

All the authors were involved in the design of the study, calculations, preparation of the manuscript, and gave final approval for publication.

Funding

AAT thanks the financial support of the CNPq under Grant No. 462067/2014-9. HVR thanks the financial support of the CNPq under Grant No. 440650/2014-3 and CAPES. EKL thanks the financial support of the CNPq under Grant No. 303642/2014-9.