Information Geometric Investigation of Solutions to the Fractional Fokker–Planck Equation

A novel method for measuring distances between statistical states as represented by probability distribution functions (PDF) has been proposed, namely the information length. The information length enables the computation of the total number of statistically different states that a system evolves through in time. Anomalous transport can presumably be modeled fractional velocity derivatives and Langevin dynamics in a Fractional Fokker–Planck (FFP) approach. The numerical solutions or PDFs are found for varying degree of fractionality (α) of the stable Lévy distribution as solutions to the FFP equation. Specifically, the information length of time-dependent PDFs for a given fractional index α is computed.


Introduction
Anomalous transport processes is ubiquitus in many different fields where a diffusive description is improper. As an example, anomalous transport may be a result from turbulent processes, e.g., super-diffusive transport properties are often found in plasmas. In the gradient region or in the scrape-off layer (SOL), the thermal and particle fluxes can be dominated by coherent structures such as (blobs) [1][2][3][4][5][6][7][8] which inherently possess a non-local character [9][10][11][12][13][14][15][16]. The fluctuations under such plasma conditions are often distributed according to Lévy statistics in contrast to the Gaussian charactistics as was displayed in [17]. the turbulence induced fluxes at the edge of the W7-AS stellarator where shown to have probability density functions (PDFs) with power law characteristics however most often it is expected that Gaussian statistics dominate, which induces exponentially decaying tails of the distributions.
Brownian motion is often a starting point in characterizing diffusion processes. In Brownian motion the mean value vanishes, on the other hand the variance or the mean square displacement grows linearly in time according to δx 2 = 2Dt. Interestingly, there are many physical processes that deviates from a linear growth in time. These processes are often called anomalous, in mathematical terms this is whenever δx 2 = 2Dt α see for example, [18][19][20]. Here, there are in principle, two limits of interest with super-and sub-diffusive properties. the distinguishing property of a super-diffusive process is to have an α > 1 whereas for a sub-diffusive process an α < 1. Lévy statistics describe super-diffusive fractal processes in terms of the Lévy index α, where α ∈ R and 0 < α ≤ 2, lie at the heart of complex processes such as anomalous diffusion. Note that Brownian motion is the special case where α = 2.0. There are several paths to generate a process exhibiting stable Lévy statistics and can be generated by random processes that are scale-invariant. A scale invariant process is obtained where no scale is dominant (e.g., there is no typical mean free path) meaning that the process will possess many scales. Geometrically, this implies that the path of a test particle, viewed at different resolutions, will look self-similar.
One of the main goals of statistical mechanics is understanding the time evolution of in and out of equilibrium systems. One possible characterization of the time dependent process is to endow a metric for the thermodynamic length [21][22][23][24][25][26] and the information length [27][28][29][30][31][32] which is a generalization of the thermodynamic length to non-equilibrium systems. This yields a mathematical framework to compare two PDFs by a distance. One of the possible choices is to use the PDFs, e.g., of a turbulent process, to construct the Fisher information metric [27,28]. Using this methodology gives a novel way to measure distance in statistical space. Here it is is important to remember that the information length methodology is path dependent in statistical space yielding an unique possibility of classification of turbulent processes.
The present paper is limited to the study of super-diffusive properties and a more detailed study of sub-diffusive dynamics is beyond this scope. Moreover, sub-diffusive dynamics have been studied previously in different context where transport is often inhibited by sticky motion, see for example [18,19]. One important case of sub-diffusion is the dynamics of holes in amorphous semiconductors, where the waiting time distribution with extended tails has been identified [33]. A generalized Langevin equation with fractional Gaussian noise has been used to describe the sub-diffusive processes within a single protein molecule by [34]. In general, 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 information length concept pertaining to time-dependent solutions of the fractional Fokker-Planck (FFP) equation resulting from a Langevin description driven by Lévy stochastic force. the present work is based on previous efforts reported in Anderson et al. [35][36][37] and may provide new insights on the recent developments in the understanding of the anomalous transport and turbulent processes. One particular example is the turbulently driven transport of charged particles in plasmas such as the super-diffusive heat transport found in JET plasmas [38] and numerically generated data of the Hasegawa-Wakatani model where the fractionality is induced by the generation of large-scale events such as zonal flows and streamer-like structures [37]. the information length may give an opportunity to identify large-scale events by a sudden change in the statistical state described by the PDFs, yielding a sudden change in information, i.e., a sharp peak in the dynamic time. Furthermore, it is known from experiments that edge localized mode severely limits the confinement of the plasma in fusion experiments and the prediction and mitigation of these modes is a high priority, thus the information length concept may give indications and even predictive capabilities of these events.
The model used for the Fokker-Planck description is presented in Section 2, and the numerical results are shown and discussed in Section 3. the final section presents a discussion and conclusions.

The Fokker-Planck and Langevin Equations
Phenomenological models pertinent for anomalous plasma transport including a fractional derivative have recently been studied [4,39,40] and in particular the properties of fractional derivatives in velocity space [35,36,38]. Fractional kinetics is a powerful tool in modeling anomalous transport processes exhibiting Lévy statistics. Fractional kinetics provides a mathematical framework and description of the non-Gaussian self-similar nature of particle displacement PDFs, as well as the anomalous scaling of moments of the distribution function. Additionally, the non-local character of the transport processeses are inherently captured by the integro-differential nature of the fractional derivative operators. As an example it can be mentioned that, the local Fourier-Fick's law is in fractional diffusion replaced by an integral operator in which the flux at a given point in space depends globally on the whole parameter space.
A Langevin equation may be used to describe the Brownian motion of a colloidal particle of the form: where v is the speed of the particle, −νv is the friction, and A(t) is the white stochastic force such that A(t)A(t ) = 2Dδ(t − t ). In the case of Brownian motion it is assumed that A(t) is a Gaussian stochastic force, this ultimately leads to a Maxwellian velocity distribution. the standard Fokker-Planck (FP) equation for the evolution of the distribution function can be found to be: Here P is the distribution function, v is the velocity, F is an external force, e.g., the electromagnetic force, m is the mass, ν is the friction, and D is the diffusion coefficient. the model may be generalized by assuming that A(t) is a stochastic noise with the properties of a Lévy-stable process. However, the FP equation is also modified in order to accommodate for the power law tails of the distribution function of the form P(v) ∝ v −α−1 for a Lévy stable with fractional index α [35,36]. This results in a fractional FP equation: where 0 < α ≤ 2 and |v| < ∞. Here, the time-dependent solution is readily found in Fourier space where the fractional Riesz operator in 1+1D can be transformed to: where the Fourier transformed distribution function can be determined to be: The Fourier transform is used to define the fractional Riesz derivative −∞D [41] for more information. It will be seen later that the main effect of the time derivative is the introduction of a relaxation time which depends on the friction and the fractionality α. Here a smaller fractional index α yield a longer relaxation time. Note that estimating the parameters D, ν, and α have to be done for the situation at hand often by experimental means such as the global scaling properties found for heat flux in the JET Tokamak [38].

Results
Here the information geometry and the resulting information lengths will be explored. The methodology is to consider a time-dependent PDF (p(x, t)) for a stochastic variable x. The Fisher-Information metric g ij can be computed when the control parameters λ i that determines the PDF are known, yielding: The distribution function p(x, t) determines the probability of the system to be in state x at time t. Th metric tensor g ij , yield the information length [27][28][29][30] as: The information length or distance in Equation (7) is a measure of the statistical distance between consecutive PDFs. the statistical distance can then be a measure of the time evolution of the system. Note that the information length is path dependent and proportional to the time integral of the square root of the infinitesimal relative entropy [31].
In general, the information length can be computed by finding the dynamic time unit τ(t) and the total time in this unit even when we do not know control parameters that govern PDFs [27][28][29][30][31].
Here the dynamic time is the typical time-scale over which the PDF (p(x, t)) temporally changes on average at time t and then determine the total elapsed time in units of τ(t). The dynamic time τ(t) is directly linked to the second moment E (as can be inferred from combining Equations (7) and (6)) and can be computed as: τ in Equation (8) quantifies the correlation time over which the (dimensionless) information changes. The information length L(t) then follows [27], The numerical solutions by inverse Fourier transform to Equation (5) is found for positive times t > 0. It should the be noted that the PDFs are non-zero, which would lead to a singularity in Equation (9).
The information length is a dimensionless quantity representing the total different number of states between the initial and final times, t i and t f , respectively. It establishes a statistical space where distances can be measured between the initial and final PDFs. An example is the Gaussian process (α = 2.0) where statistically distinguishable states are determined by the standard deviation, which increases with the level of fluctuations. In this work, the numerical solutions to the linear FFP in Equation (3) is utilized as the continously changing PDF as time progresses. This determines a specific path in statistical space and thus ultimately leads to a linear increase in the information. However in general terms the information length between two PDFs is dependent on the path between the PDFs and can thus take an arbitrary value depending on the total number of different statistical states that a system passes through in time [27][28][29][30][31][32].
In comparison with this work Heseltine and Kim [32] compared the information length, relative entropy, and Jensen divergence, and showed that it was only the information length that captures the linear geometry of a linear Ornstein-Uhlenbeck process by a linear relation between the information length (in the long time limit) and the mean position of an initial Gaussian PDF [32]. the information length constitutes a geometric methodology to understand stochastic processes such as dynamic phase transition [29]. This study extends this to the general FFP as desribed by Anderson et al. in [35,36].
In Figure 1, the time evolution of the solutions to the FFP by the inverse Fourier transform of Equation (5) for α = 2.0 (left) and α = 1.5 (right) are displayed. the time evolution is also dependent of the diffusion D = 1.0 and viscosity ν = 1.0. A decreased diffusion results in peaked PDFs with decreased variance whereas increased diffusion has the opposite effect however the time dependency is not changed. By changing the viscosity both variance and the time evolution is changed. the information length in Equation (9) will be computed for different α and the effect of D and ν will be discussed.  In Figure 2, the dynamic time computed by Equation (8) and the information length Equation (9) are shown for D = 1.0 and ν = 1.0. A convergence to a constant dynamic time in the long time limit is evident and the information length thus increases linearly with time, which is in accordance with the results found in earlier papers [32,37]. Next a discussion of the effect of viscosity on the dynamic time and the information length. In Equation (5), it can be seen that a change in viscosity ν changes both the time evolution and the variance of the PDF. Figure 3, display the dynamic time (Equation (8), left) and information length (Equation (9), left) as a function of the viscosity ν = 1.0, 2.0, 5.0, 10.0 where the other parameters are fractional index α = 1.5 and the diffusion coefficient D = 1.0. Here the rapid change in the time evolution with increased viscosity is clearly visible as well as the change in the variance of the PDFs. the variance decreases with increasing ν in addition to a more rapid time-relaxation of the PDF is visible as a change in the instantaneous dynamic time, however the information length is still increasing linearly.

Summary and Conclusions
The present paper elucidates the dynamic time and information length for the solutions to the Fractional Fokker-Planck Equation. the solutions are based on previous work by Anderson et al. in [35,36]. the Fractional Fokker-Planck Equation (FFP) is a unique description of anomalous transport which is an ubiquitous phenomenon in fusion plasma dynamics and thus may give a deeper understanding in plasma transport far from equilibrium. the FFP is also mathematically interesting due to the inherent non-local nature of the solutions, thus it is a very broad and interesting topic crossing many different fields of research. the dynamic time, Equation (8) and information length, Equation (9), are novel tools in characterizing dynamical systems, in particular systems that exhibit anomalous transport. It is found that the dynamical time converge to a constant value in the long time limit thus the information length increases linearly with time. Here it is pertinent to point out that the model is in contrast to previous models due to the fractional velocity derivative and thus it is expected and also shown that the slope is rather different in comparison to Gaussian models. Furthermore, it is shown that the general linear increase of information holds in non-linear simulations, as is seen in [37], although fluctuations in the overall increase is due to a sudden generation of structures impacting the transport. However, it should be noted that the models used here are linear in time and space and thus affects the full effects of non-linearities that are not taken into account.
The current work opens up for new research on the statistical characterization of anomalous transport, in particular continuing the analysis of the experimental data obtained in the JET Tokamak, see [38], as was started with an estimation of a global scaling of heat flux where the fractional α∼0.8 was found. Further characterization of transport events and modes causing the anomalous transport in Tokamaks using the current results employing a fractional model would be of great interest. Note that the information length approach is path dependent in statistical space capturing the creation and interaction of modes and events, and thus may give a predictive capability of disastrous events that reduces confinement.