# Non-Linear Langevin and Fractional FokkerâPlanck Equations for Anomalous Diffusion by Lévy Stable Processes

###### Abstract

The numerical solutions to a non-linear Fractional Fokker–Planck (FFP) equation are studied estimating the generalized diffusion coefficients. The aim is to model anomalous diffusion using an FFP description with fractional velocity derivatives and Langevin dynamics where Lévy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space. Distribution functions are found using numerical means for varying degrees of fractionality of the stable Lévy distribution as solutions to the FFP equation. The statistical properties of the distribution functions are assessed by a generalized normalized expectation measure and entropy and modified transport coefficient. The transport coefficient significantly increases with decreasing fractality which is corroborated by analysis of experimental data.

###### pacs:

05.40Fb; 02.50Ey; 05.40-a## I Introduction

In magnetically confined (MC) plasma devices transport driven by turbulent fluctuations often severely limit the confinement time and thus impede the development of fusion as an alternative for electricity production. It is pertinent to understand and mitigate the effects of the turbulently driven transport where simplified models often are employed in order to elucidate the main features of the plasma turbulence. In magnetised plasmas, it is commonly accepted that turbulence is the primary cause of anomalous (i.e., elevated compared to collisional) transport horton (); krommes2002 (). It has been recognized that the nature of the anomalous transport processes is dominated by a significant ballistic or non-local component where a diffusive description is improper. The turbulence in MC tokamak plasmas is anisotropic in the parallel and perpendicular length scales to the magnetic field and taps free energy from the pressure gradient that can drive fluctuations in electrostatic potential and density horton (); krommes2002 (). The super-diffusive properties are often ubiquitously found in plasmas, such as the thermal and particle fluxes in the gradient region or in the Scrape-Off Layer (SOL) where the transport is dominated by the coherent structures (blobs) carreras1996 (); carreras1999 (); milligen2005 (); sanchez2008 (); negrete2005 (); sanchez2006 (); hahm () and inherently possess a non-local character Zweben (); Naulin (); kaye (); Lopez (); Gentle (); Mantica (); van-Milligen (). Moreover, there is a large quantity of experimental evidence that density and potential fluctuations measured by Langmuir probes at different fusion devices support the idea that these fluctuations are distributed according to Lévy statistics. This was illustrated for example in carreras1999 (), where probability density functions (PDFs) of the turbulence induced fluxes at the edge of the W7-AS stellarator were shown to exhibit power law characteristics in contrast to exponential decay expected from Gaussian statistics. Furthermore, the experimental evidence of the wave-number spectrum characterized by power laws over a wide range of wave-numbers can be directly linked to the values of Lévy index of the PDFs of the underlying turbulent processes. One widely used simplified model of a plasma is the Hasegawa–Wakatani model which was recently analysed by statistical methods in Reference anderson2017 (). It was concluded that even simplified models may have components of fractionality stemming from the non-linear interactions and the generation of large scale modes such as zonal and shear flows. The Hasegawa–Wakatani model allows for the electrons to dynamically and self-consistently determine the relationship between the density and the electrostatic potential through the turbulence. Moreover, fractal features in transport have been observed in many experiments in many different fields of research. In particular it has been found that there is strong evidence of non-local heat transport in JET plasmas moradi2018 (). In this paper, fractal features is synonymous to a system where power law statistics is found. Here it is important to keep in mind that, although a simplified fractional transport model is used, it indicates that there is a lack of physics in the current transport models based on the mean field theory, namely the super-diffusive character of heat transport. Finding a proper kinetic description of dynamical systems with chaotic behaviour is one of the main problems in classical physics schlesinger1993 (); sokolov2002 (); klafter2005 (); metzler2000 (); metzler2004 (); mandelbrot1982 (); anderson2010 (); kim2009 (); moradipop2011 (); moradipop2012 (); anderson2014 (); moradipop2016 (). Over the past two decades it has become obvious that behaviour much more complex than standard diffusion can occur in dynamical Hamiltonian chaotic systems. In principle, the orbits in dynamical systems are always theoretically predictable since they arise as solutions to simple system of equations such as Newton’s equations; however, these orbits are sensitive to initial conditions and thus very small changes in initial conditions may yield widely different outcomes. From the macroscopic point of view, the rapid mixing of orbits has been used as a motivation for assumptions of randomness of the motion and the random walk models schlesinger1993 (). In characterizing the diffusion processes in plasmas, the starting point is often Brownian motion where the mean value vanishes, whereas the second moment or variance grows linearly in time according to . However, taking into account the experimental data found in plasma experiments, it is evident that many phenomena exhibit anomalous diffusion where variance grows non-linearly in time such that . The reason an anomalous diffusion approach is needed is due to the restrictive assumptions of locality in space and time, and lack of long-range correlations that is the basis of Brownian motion. There are two limits of interest for where the first is super-diffusion with and the second is sub-diffusion with . A super-diffusive description is most often appropriate for fusion plasmas. Lévy statistics describing fractal processes (Lévy index where ) lie at the heart of complex processes such as anomalous diffusion. Lévy statistics can be generated by random processes that are scale-invariant. This means that a trajectory will possess many scales, but no single scale will be characteristic and dominate the process. Geometrically, this implies the fractal property that a trajectory, viewed at different resolutions, will look self-similar. Such strange kinetics schlesinger1993 (); mandelbrot1982 () may be generated by accelerated or sticky motions along the trajectory of the random walk krommes2002 (). Super-diffusivity may also occur as a result of variation in the step length of the motion, which breaks the assumption that a unique step length may, e.g., give rise to long-range correlations in the dynamics generated by the presence of anomalously large particle displacements connecting otherwise physically disjoint domains.

We note that, although sub-diffusive processes are beyond the scope of the present work, its properties have been studied in many different contexts where transport is often inhibited by sticky motion. Among sub-diffusive phenomena are holes in amorphous semiconductors, where a waiting time distribution with long tails has been introduced montroll1973 (). The sub-diffusive processes within a single protein molecule have been described by generalized Langevin equation with fractional Gaussian noise kou2004 (). Turbulence and related anomalous diffusion phenomena have been observed in a wide variety of complex systems such as high energy plasmas, semiconductors, glassy materials, nanopores, biological cells, and epidemic proliferation.

The objective of the present paper is to explore the non-linear character of the fractional Fokker–Planck (FFP) equation resulting from a Langevin description driven by Lévy stochastic force with a non-linear interaction in the velocity. The present work is based on previous efforts reported in Reference anderson2014 () and may provide new insights on the recent developments in the modelling of the anomalous transport of charged particles in magnetised plasmas, such as the non-local heat transport found in JET plasmas.

## Ii The Fokker–Planck and Langevin Equations

Fractional kinetics is a powerful framework in describing anomalous transport processes exhibiting Lévy statistics. It is able to reproduce key aspects of anomalous transport including the non-Gaussian self-similar nature of the PDFs of particle displacement, and the anomalous scaling of the moments. It has been shown that the chaotic dynamics can be described by using the FFP equation with coordinate fractional derivatives as a possible tool for the description of anomalous diffusion zaslavsky (). Previous papers on plasma transport have used models including a fractional derivatives on phenomenological premises sanchez2008 (); del-Castillo-Negrete2004 (); del-Castillo-Negrete2010 (). Additionally, the integro-differential nature of the fractional derivative operators allows the description of spatiotemporal nonlocal transport processes. In particular, in fractional diffusion, the local Fourier–Fick’s law is replaced by an integral operator in which the flux at a given point in space depends globally on the spatial distribution of the transported scalar and on the time history of the transport process. Using fractional generalizations of the Liouville equation, kinetic descriptions have been developed zaslavsky2002 (); tarasov2005 (); tarasov2006 (). The currently applied model is based on the Langevin equation with a Lévy-stable noise term, where the applied noise exhibits a power law tail levy (); seshadri (). The generalized Central Limit Theorem for Lévy-stable processes is a particular weak-convergence theorem in probability theory. It expresses the fact that a sum of many independent and identically distributed (i.i.d.) random elements, or alternatively, random elements with specific types of dependence, will tend to be distributed according to one of a small set of attractor distributions. There are here two cases of special interest: the first is when the variance of the i.i.d. variables is finite and the attractor distribution is then a normal distribution, and the second is where the sum of a number of i.i.d. random elements with power law tail distributions decreasing as where (therefore having infinite variance) will tend to a Lévy-stable distribution with a fractality index of as the number of elements in the set increases.

The motion of a colloidal particle can be described by the Langevin equation in the case of Brownian motion and it will take the form

(1) |

where is the speed of the particle, is the friction, and is the white stochastic force such that . Moreover, by assuming that is a Gaussian stochastic force, a Maxwellian velocity distribution may be obtained and would lead to the standard Fokker–Planck (FP) equation describing the evolution of the distribution function:

(2) |

Here is the distribution function, is the velocity, is an external force, e.g., the electromagnetic force, is the mass, is the friction, and is the diffusion coefficient. The corresponding reduced FP equation , where the Lorentz force is neglected, to the Langevin equation is

(3) |

which yields to the stationary state Maxwellian velocity distribution for fogedby1994 (); fogedby19942 (). However, if is a stochastic noise with the properties of a Lévy-stable process, the FP equation has to be modified in order to accommodate for power law tails of the form for a Lévy stable with fractional index . The FFP equation becomes

(4) |

where and . The time-dependent solution is readily found in the Fourier space where the fractional Riesz operator in 1 + 1D can be transformed to

(5) |

where the Fourier transformed distribution function can be determined to be

(6) |

The fractional Riesz derivative is defined through its Fourier transform , see, e.g., metzler2000 () for more information. Here it should be noted that the time derivative only introduces a relaxation time dependent on the friction and the fractionality , where a smaller yield a longer relaxation time.

In Figure 1, the exponentially fast relaxation of the velocity PDFs with time is displayed. The PDFs of a Gaussian () and for a PDF with fractional index for times , and are computed numerically. We note that, at , the PDFs are close to the stationary state PDF, whereas the time evolution of the PDF depends on the fractional index such that the relaxation process is slower for a PDF with a lower fractional index. In general, the distributions found for the have more pronounced tails and sharper peaks, whereas, in the case, the system has a shorter relaxation time.

## Iii Results

The aim of the present paper is to look into the effects of a non-linear interaction in the Langevin equation, but it is here assumed that we can neglect the time dependence , i.e., the stationary state PDF (), and the FFP Equation can be written as

(7) |

Here , , and are constants. The equation is obtained by inclusion of the quartic potential, leading to the addition of a term of third order of the form . The main effect of retaining the temporal dynamics is to introduce a relaxation time. In the current model, square and quartic terms will be retained. The properties of the current non-linear terms are analogous to a potential well with square and quartic terms. Note that even terms in the potential provide proper stable equilibria, whereas odd terms yield an unstable equilibrium; thus, the third and fifth order terms are neglected. The Equation (7) is directly integrated by using a predictor—corrector method according to Adams-Bashforth-Moulton diethelm ().

To find an analytical solution of the original Langevin equation, the Fourier transform can be used;

(8) |

The found equation is a third order ordinary differential equation with variable coefficients. The general solutions to Equation (III) can only be determined by numerical means however a similar system was investigated in Reference kim2009 () suggesting a PDF proportional to , where is a constant. Furthermore, it is also possible to find an analytical solution for the tail of the PDF to leading order by using the WentzelâKramersâBrillouin (WKB) approximation for small values of . The WKB anzats is to assume a series solution to the Fourier transformed equation (III), of the form , here will be taken small and to be determined in terms of . It is then found that, the leading order tail contribution corroborates the findings in kim2009 () for . We note that the real space distribution function is convergent for and can only in general be obtained by numerical integration, and is here solved by using method described above. We note that there are three different interesting regimes: the first is where the diffusion is much larger than the quartic potential strength , the second is where the diffusion is comparable to the quartic potrential strength , and the third is where the diffusion is negligible to the quartic potential strength . In the third regime, the PDFs become may be expected to have similarities to the results found in kim2009 () for where for some constant . The values used in this study are chosen to illustrate these three regimes of interest. Note that the non-linear interaction, i.e., the term introduces three different regimes with richer dynamics which is in contrast to what was found in Reference anderson2014 (). In any linear fractal model based on the Lévy statistics the power law tails of the velocity PDF will be . Even more interestingly, in non-linear models the precise scaling of the PDF tails are still open.

In Figures 2–4, the numerically found PDFs, by solving Equation (7) by the Adams-Bashforth-Moulton method diethelm (), are shown for the three different regimes: and , and , and and , respectively. The resolution in is except in the case of for and where the resolution is increased to , however increased resolution would not change the in any significant way for smaller —v—. As expected, in Figure 2 almost independently of where , the PDFs exhibit power law tails, although in the case of some exponential behaviour is observed at the tail of the distribution (for large modulus of velocity ). In Figures 3 and 4, we find more pronounced tails in particular in the low- case. In the low- case, the fractal term dominates the dynamics. We note that the PDFs displayed in Figure 3 exhibit a hybrid between fractal and Gaussian behaviour (), where in some cases the PDF is retains some Gaussian behaviour, which is particularly visible for small . In the regime the non-Gaussian effects of the PDFs are clearly visible in Figure 4. The PDFs are used to evaluate the dynamics of the system in terms of Tsallis’ statistical mechanics where -entropy and -energy determines the properties of the system for the three different regimes: and , and , and and in Figures 5 and 6, respectively. The , -entropy, and -energy values are determined by the following relations (see References tsallis1995lnp (); tsallis1996 (); tsallis1998 (); hamza (); barkai (); angulo (); balasis2011 (); pavlos2012 (); pavlos20122 (); tsallis1995 (); prato1999 (); anderson2014 ()):

(9) | |||||

(10) | |||||

(11) |

The entropy and -entropy are displayed in Figure 5. A maximum in the entropy, computed according to , is found, whereas almost constant -entropy, computed according to Equation (10), is found with increasing , where is determined according to Equation (9). An increasing trend with importance of non-local effects are also visible where the entropy and -entropy is increased in the cases with , which is in the regime where the fractal nature is more prominent. In Figure 6, a decreased energy and increased -energy with increasing for and small case is found.

The interpretation of this strange kinetics has to be based on the results from experimental data since there is no first principle method to compute the value of and thus is indeterminable. However, recent findings suggest that JET plasmas have a significant degree of super-diffusive transport with an , and it was found that this super-diffusive transport is slightly different for the ion and electron channels moradi2018 (). The analysis is based on a power balance where a large set of JET shots are used whereby a distribution in can be obtained with a mean value of approximately , suggesting that a convective model would be more appropriate with . The diffusion coefficient can be estimated by the velocity autocorrelation functions according to the Kubo formula, but such an estimate looks at the ratio of the generalized diffusion coefficient () and the Brownian diffusion coefficient , using the tempered -velocity correlators, computed by Equation (11).

We find that the ratio of the diffusion coefficients increases with smaller and significantly increases in the regime where fractality is pronounced, as shown in Figure 7. Interestingly enough, in the analysis presented in moradi2018 () it is evident that in the cases with increased transport a lower value of is obtained, indicating a strong non-diffusive component or equivalently, a significantly increased transport where processes following Lévy statistics dominate the transport. The qualitative increase in the generalized transport coefficient is thus qualitatively corroborated by what is seen experimentally using the power balance analysis.

## Iv Discussion and Conclusions

Understanding anomalous transport in MC plasmas is an outstanding issue in controlled fusion research. It is commonly accepted that, in these plasmas, turbulence is the primary cause of anomalous (i.e., elevated compared to collisional) transport. It has also been recognized that the nature of the anomalous transport processes is dominated by a significant ballistic or non-local component where a diffusive description is improper. A satisfactorily understanding of the non-local features as well as the non-Gaussian PDFs found in experimental measurements of particle and heat fluxes is still lacking anderson2010 (); kim2009 (), but there has been some recent progress in this direction. Fractional kinetics has been put forward for building a more physically relevant kinetic description for such dynamics. In these situations, kinetic descriptions, which arise as a consequence of averaging over the well-known Gaussian and Poissonian statistics (for diffusion in space and temporal measures, respectively), seem to fall short in describing the apparent randomness of dynamical chaotic systems schlesinger1993 (). This is due to the restrictive assumptions of locality in space and time, and the lack of long-range correlations that is the basis of these descriptions.

In magnetised plasma experiments, a predator–prey system exists with avalanches (strong driver of transport) and zonal flows (sheared flows that decorrelate turbulent eddies reducing transport). It has been suggested that an Fractional Fokker–Planck Equation on a comb-like potential background can be applied where meso-scale transport events (avalanching) occurs in between regions of strong zonal flow activity (see Milovanov and Rasmussen milovanov ()). This method is straightforward for applications in this setting; by assuming the used potential in between the zonal flow regions, it is suggested that the potential should be of degree 4 (or higher), as has been used here.

Although there has previously been some criticism on the appropriateness of using the Tsallis method in describing processes with Lévy statistics, this is mainly concerning descriptions based on fractality in coordinate space not in velocity space. However, the aim of the present work was to shed light on the non-extensive properties of the velocity space statistics and characterization of the fractal processes by estimating the generalized diffusion coefficients of the FFP equation in terms of Tsallis statistics. Jespersen et al. jespersen () showed an example of the Langevin equation with a harmonic potential, and the Tsallis -statistics had limited usefulness. The reason for this is that, using the variational calculus of Equation (10) with the appropriate constraints, the relation between and is , which is different from Equation (9) and thus cannot reproduce the correct scaling. They then concluded that the Tsallis entropy was not the appropriate framework for Lévy flights in a harmonic potential described by the generalized FP equation. However, this limitation seems not to impede the usefulness of the application of Tsallis entropy on this Langevin equation where the correct scaling is obtained.

In summary, we have employed an FFP equation to find the PDFs and studied the -entropy and -energies in this system with a non-linear interaction in the FFP equation. We found a significantly elevated diffusion coefficient, which is qualitatively similar to what was expected in light of the analysis of the experimental data.

## References

- (1) Horton, W. Turbulent Transport in Magnetized Plasmas; World Scientific: Danvers, MA, USA, 2017; ISBN 978-981-3225-88-6.
- (2) Krommes, J.A. Fundamental statistical descriptions of plasma turbulence in magnetic fields. Phys. Rep. 2002, 360, 1–352.
- (3) Carreras, B.A.; Hidalgo, C.; Sanchez, E.; Pedrosa, M.A.; Balbin, R.; Garcia, C.I.; van Milligen, B.; Newman, D.E.; Lynch, V.E. Fluctuationâinduced flux at the plasma edge in toroidal devices. Phys. Plasmas 1996, 3, 2664–2672.
- (4) Carreras, B.A.; van Milligen, B.; Hidalgo, C.; Balbin, R.; Sanchez, E.; Cortes, I.G.; Pedrosa, M.A.; Bleuel, J.; Endler, M. Self-similarity properties of the probability distribution function of turbulence-induced particle fluxes at the plasma edge. Phys. Rev. Lett. 1999, 83, 3653–3656.
- (5) Van Milligen, B.P.; Sanchez, R.; Carreras, B.A.; Lynch, V.E.; LaBombard, B.; Pedrosa, M.A.; Hidalgo, C.; GonÃ§alves, B.; BalbÃn, R.; The W7-AS Team. Additional evidence for the universality of the probability distribution of turbulent fluctuations and fluxes in the scrape-off layer region of fusion plasmas. Phys. Plasmas 2005, 12, 52501–52507.
- (6) Sanchez, R.; Newman, D.E.; Leboeuf, J.N.; Decyk, V.K.; Carreras, B.A. Nature of Transport across Sheared Zonal Flows in Electrostatic Ion-Temperature-Gradient Gyrokinetic Plasma Turbulence. Phys. Rev. Lett. 2008, 101, 205002–205004.
- (7) del-Castillo-Negrete, D.; Carreras B.A.; Lynch, V.E. Front Dynamics in Reaction-Diffusion Systems with Levy Flights: A Fractional Diffusion Approach. Phys. Rev. Lett. 2005, 94, 18302–18304.
- (8) Sanchez, R.; Carreras, B.A.; Newman, D.E.; Lynch, V.E.; van Milligen, B.P. Renormalization of tracer turbulence leading to fractional differential equations. Phys. Rev. E 2006, 74, 16305–16311.
- (9) Hahm, T.S. Nonlinear gyrokinetic equations for tokamak microturbulence. Phys. Fluids 1988, 31, 2670–2673.
- (10) Zweben, S.J. Search for coherent structure within tokamak plasma turbulence. Phys. Fluids 2007, 28, 974–982.
- (11) Naulin V. Turbulent transport and the plasma edge. J. Nucl. Mater. 2007, 363-365, 24-31.
- (12) Kaye, S.M.; Barnes, C.W.; Bell, M.G.; DeBoo, J.C.; Greenwald, M.; Riedel, K.; Sigmar, D.; Uckan, N.; Waltz, N. Status of global energy confinement studies. Phys. Plasmas 1990, 2, 2926–2940.
- (13) Cardozo, N.J.L. Perturbative transport studies in fusion plasmas. Plasma Phys. Control. Fusion 1995, 37, 799–852.
- (14) Gentle, K.W.; Bravenec, R.V.; Cima, G.; Gasquet, H.; Hallock, G.A.; Phillips, P.E.; Ross, D.W.; Rowan, W.L.; Wootton, A.J. An experimental counterâexample to the local transport paradigm. Phys. Plasmas 1995, 2, 2292–2298.
- (15) Mantica, P.; Galli, P.; Gorini, G.; Hogeweij, G.M.D.; de Kloe, J.; Cardozo, N.J.L.; RTP Team. Nonlocal transient transport and thermal barriers in rijnhuizen tokamak project plasmas. Phys. Rev. Lett. 1999, 82, 5048–5051.
- (16) Van-Milligen, B.P.; de la Luna, E.; Tabars, F.L.; AscasÃbar, E.; Estrada, T.; CastejÃ³n, F.; Castellano, J.; CortÃ©s, I.G.; Herranz, J.; Hidalgo, C.; et al. Ballistic transport phenomena in TJ-II. Nucl. Fusion 2002, 42, 787–795.
- (17) Anderson, J.; Hnat, B. Statistical analysis of Hasegawa-Wakatani turbulence. Phys. Plasmas 2017, 24, 62301–62308.
- (18) Moradi, S.; Anderson, J.; Romanelli, M. Evidence of non-local heat transport model in JET plasmas. Presented at EU-US Transport Task Force Meeting, Seville, Spain, 11–14 September 2018.
- (19) Schlesinger, M.F.; Zaslavsky, G.M.; Klafter, J. Strange kinetics. Nature 1993, 363, 31–37.
- (20) Sokolov, I.M.; Klafter, J.; Blumen, A. Fractional kinetics. Phys. Today 2002, 55, 48–54.
- (21) Klafter, J.; Sokolov, I.M. Anomalous diffusion spreads its wings. Phys. World 2005 18, 29–32.
- (22) Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77.
- (23) Metzler, R.; Klafter, J. The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A 2004, 37, 161–208.
- (24) Mandelbrot, B.B. Fractals and Geometry of Nature; W. H. Freeman and Company: San Francisco, CA, USA, 1982; pp. 170–180.
- (25) Anderson, J.; Xanthopoulos, P. Signature of a universal statistical description for drift-wave plasma turbulence. Phys. Plasmas 2010, 17, 110702–110704.
- (26) Kim, E.; Liu, H.L.; Anderson, J. Probability distribution function for self-organization of shear flows. Phys. Plasmas 2009 16, 52301–52304.
- (27) Moradi, S.; Anderson J.; Weyssow, B. A theory of non-local linear drift wave transport. Phys. Plasmas 2011, 18, 062101–062106.
- (28) Moradi, S.; Anderson, J. Non-local gyrokinetic model of linear ion-temperature-gradient modes. Phys. Plasmas 2012, 19, 82301–82307.
- (29) Anderson, J.; Kim, E.; Moradi, S. A fractional Fokker–Planck model for anomalous diffusion. Phys. Plasmas 2014, 21, 122101–122108.
- (30) Moradi, S.; del Castillo, N.D.; Anderson, J. Charged particle dynamics in the presence of non-Gaussian LÃ©vy electrostatic fluctuations. Phys. Plasmas 2016, 23, 907041–907045.
- (31) Montroll, E.W.; Scher, H. Random walks on lattices. IV. Continuous-time walks and influence of absorbing boundaries. J. Stat. Phys. 1973, 9, 101–135.
- (32) Kou, S.C.; Sunney X. Generalized langevin equation with fractional Gaussian noise: Subdiffusion within a single protein molecule. Phys. Rev. Lett. 2004, 93, 1806031–1806034.
- (33) Combescure, M. Hamiltonian Chaos and Fractional Dynamics; Oxford University Press: Oxford, UK, 2005; pp. 0305–4470.
- (34) del Castillo, N.D.; Carreras, B.A.; Lynch, V.E. Fractional diffusion in plasma turbulence. Phys. Plasmas 2004, 11, 3854–3864.
- (35) Del Castillo, N.D. Non-diffusive, non-local transport in fluids and plasmas. Nonlinear Proc. Geophys. 2010, 17, 795–807.
- (36) Zaslavsky, G.M. Chaos, fractional kinetics, and anomalous transport. Phys. Rep. 2002, 371, 461–580.
- (37) Tarasov, V.E. Fractional Liouville and BBGKI equations. J. Phys. 2005, 7, 17–33.
- (38) Tarasov, V.E. Fractional statistical mechanics. Chaos 2006, 16, 331081–331087.
- (39) Lévy, P. ThÃ©orie de L’addition des Variables AlÃ©atoires; Gauthier-Villiers: Paris, France, 1937.
- (40) West, B.J.; Seshadri, V. Linear systems with LÃ©vy fluctuations. Physical A 1982, 113, 203–216.
- (41) Fogedby, H.C. Langevin equations for continuous time LÃ©vy flights. Phys. Rev. E 1994, 50, 1657–1660.
- (42) Fogedby, H.C. LÃ©vy Flights in Random Environments. Phys. Rev. Lett. 1994, 73, 2517–2520.
- (43) Diethelm, K.; Freed, A.D. The Fractional PECE Subroutine for the numerical solution of differential equations of fractional order. In Forschung und Wissenschaftliches Rechnen; Heinzel, S., Plesser, T., Eds.; Gessellschaft fur Wissenschaftliche Datenverarbeitung: Gottingen, Germany, 1999; pp. 57–71.
- (44) Tsallis, C.; de Souza, A.M.C.; Maynard, R. Derivation of LÃ©vy-type anomalous superdiffusion from generalized statistical mechanics. In LÃ©vy Flights and Related Topics in Physics; Springer: New York, NY, USA, 1995; Volume 450, pp. 269–289.
- (45) Tsallis, C.; Bukman, D.J. Anomalous diffusion in the presence of external forces: Exact time-dependent solutions and their thermostatistical basis. Phys. Rev. E 1996, 54, 2197–2200.
- (46) Tsallis, C.; Mendes, R.S.; Plastino, A.R. The role of constraints within generalized nonextensive statistics. Physical A 1998, 261, 534–554.
- (47) Hamza, A.B.; Krim, H. Jensen-Rényi divergence measure: Theoretical and Computational Perspectives. In Proceeding of the 2003 IEEE International Symposium on Information Theory, Yokohama, Japan, 29 June–4 July 2003.
- (48) Barkai, E. Stable equilibrium based on LÃ©vy statistics: Stochastic collision models approach. Phys. Rev. E. 2003, 68, 551041–551044.
- (49) Angulo, J.M.; Esquivel, F.J. Multifractal dimensional dependence assessment based on Tsallis mutual information. Entropy 2015, 17, 5382–5401.
- (50) Balasis, G.; Daglis, I.A.; Anastasiadis, A.; Papadimitriou, C.; Mandea, M.; Eftaxiasb, K. Universality in solar flare, magnetic storm and earthquake dynamics using Tsallis statistical mechanics. Physical A 2011, 390, 341–346.
- (51) Pavlos, G.P.; Karkatsanis, L.P.; Xenakis, M.N.; Sarafopoulos D.; Pavlos, E.G. Tsallis statistics and magnetospheric self-organization. Physical A 2012, 391, 3069–3080.
- (52) Pavlos, G.P.; Karkatsanis, L.P.; Xenakis, M.N. Tsallis non-extensive statistics, intermittent turbulence, SOC and chaos in the solar plasma, Part one: Sunspot dynamics. Physical A 2012, 391, 6287–6319.
- (53) Tsallis, C.; Lévy, S.V.F.; Souza, A.M.C.; Maynard, R. Statistical-mechanical foundation of the ubiquity of LÃ©vy distributions in nature. Phys. Rev. Lett. 1995, 75, 3589–3593.
- (54) Prato, D.; Tsallis, C. Nonextensive foundation of LÃ©vy distributions. Phys. Rev. E 1999, 60, 2398–2401.
- (55) Milovanov, A.V.; Rasmussen, J.J. LÃ©vy flights on a comb and the plasma staircase. Phys. Rev. E 2018, 98, 222081–222092.
- (56) Jespersen, S.; Metzler, R.; Fogedby, H.C. LÃ©vy flights in external force fields: Langevin and fractional Fokker–Planck equations and their solutions. Phys. Rev. E 1999, 59, 2736–2745.