Time-fractional Caputo derivative versus other integrodifferential operators in generalized Fokker-Planck and generalized Langevin equations

Fractional diffusion and Fokker-Planck equations are widely used tools to describe anomalous diffusion in a large variety of complex systems. The equivalent formulations in terms of Caputo or Riemann-Liouville fractional derivatives can be derived as continuum limits of continuous-time random walks and are associated with the Mittag-Leffler relaxation of Fourier modes, interpolating between a short-time stretched exponential and a long-time inverse power-law scaling. More recently, a number of other integrodifferential operators have been proposed, including the Caputo-Fabrizio and Atangana-Baleanu forms. Moreover, the conformable derivative has been introduced. We study here the dynamics of the associated generalized Fokker-Planck equations from the perspective of the moments, the time-averaged mean-squared displacements, and the autocovariance functions. We also study generalized Langevin equations based on these generalized operators. The differences between the Fokker-Planck and Langevin equations with different integrodifferential operators are discussed and compared with the dynamic behavior of established models of scaled Brownian motion and fractional Brownian motion. We demonstrate that the integrodifferential operators with exponential and Mittag-Leffler kernels are not suitable to be introduced to Fokker-Planck and Langevin equations for the physically relevant diffusion scenarios discussed in our paper. The conformable and Caputo Langevin equations are unveiled to share similar properties with scaled and fractional Brownian motion, respectively.