An Empirical Method to Characterize Displacement Distribution Functions for Anomalous and Transient Diffusion

We propose a practical empirical fitting function to characterize the non-Gaussian displacement distribution functions (DispD) often observed for heterogeneous diffusion problems. We first test this fitting function with the problem of a colloidal particle diffusing between two walls using Langevin Dynamics (LD) simulations of a raspberry particle coupled to a lattice Boltzmann (LB) fluid. We also test the function with a simple model of anomalous diffusion on a square lattice with obstacles. In both cases, the fitting parameters provide more physical information than just the Kurtosis (which is often the method used to quantify the degree of anomaly of the dynamics), including a length scale that marks where the tails of the DispD begin. In all cases, the fitting parameters smoothly converge towards the Gaussian values when the systems become less anomalous