We obtain exact results for fractional equations of Fokker-Planck type using evolution operator method. We employ exact forms of one-sided Levy stable distributions to generate a set of self-reproducing solutions. Explicit cases are reported and studied for various fractional order of derivatives, different initial conditions, and for differe
It is pointed out that, for the fractional Fokker-Planck equation for subdiffusion proposed by Metzler, Barkai, and Klafter [Phys. Rev. Lett. 82 (1999) 3563], there are four types of infinitely many exact solutions associated with the newly discovered exceptional orthogonal polynomials. They represent fractionally deform
We have derived a fractional Fokker-Planck equation for subdiffusion in a general space-and- time-dependent force field from power law waiting time continuous time random walks biased by Boltzmann weights. The governing equation is derived from a generalized master equation and is shown to be equivalent to a subordinated stochastic Langevin equation.
Anomalous dynamics characterized by non-Gaussian probability distributions (PDFs) and/or temporal long-range correlations can cause subtle modifications of conventional fluctuation relations. As prototypes we study three variants of a generic time-fractional Fokker-Planck equation with constant force. Type A generates superdiffusion, type B subdiffusion and type C both super- and subdiffusion depending on parameter variation. Furthermore type C obeys a fluctuation-dissipation relation whereas A and B do not. We calculate analytically the position PDFs for all three cases and explore numerically their strongly non-Gaussian shapes. While for type C we obtain the conventional transient work fluctuation relation, type A and type B both yield deviations by featuring a coefficient that depends on time and by a nonlinear dependence on the work. We discuss possible applications of these types of dynamics and fluctuation relations to experiments.
We study the connection between the parameters of the fractional Fokker-Planck equation, which is associated with the overdamped Langevin equation driven by noise with heavy-tailed increments, and the transition probability density of the noise generating process. Explicit expressions for these parameters are derived both for finite and infinite variance of the rescaled transition probability density.
Fokker-Planck equations are extensively employed in various scientific fields as they characterise the behaviour of stochastic systems at the level of probability density functions. Although broadly used, they allow for analytical treatment only in limited settings, and often is inevitable to resort to numerical solutions. Here, we develop a computational approach for simulating the time evolution of Fokker-Planck solutions in terms of a mean field limit of an interacting particle system. The interactions between particles are determined by the gradient of the logarithm of the particle density, approximated here by a novel statistical estimator. The performance of our method shows promising results, with more accurate and less fluctuating statistics compared to direct stochastic simulations of comparable particle number. Taken together, our framework allows for effortless and reliable particle-based simulations of Fokker-Planck equations in low and moderate dimensions. The proposed gradient-log-density estimator is also of independent interest, for example, in the context of optimal control.