Diffusing diffusivity: a new derivation and comparison with simulations

Many experiments are now available where it has been shown that the probability distribution function (pdf) for the position of a Brownian particle diffusing in a heterogeneous medium is not Gaussian. However, in spite of this non-Gaussianity, the mean square displacement (MSD) still remains Fickian, i.e., . One possible explanation of this non-Gaussian yet Brownian behavior is that the diffusivity of the particle itself is ``diffusing''. Chubynsky and Slater (Phys. Rev. Lett. 113 098302 2014) proposed a model of ``diffusing diffusivity'' which they were able to solve analytically at small time scales, but simulations were performed for intermediate to large time scales. We present here a class of diffusing diffusivity models and show that the problem of calculating pdf for the position of diffusing particle is equivalent to calculating the survival probability of a particle undergoing Brownian motion in the presence of a sink. We give exact analytical results for all time scales and show that the pdf is non-Gaussian at short times which crosses over to a Gaussian at long times. The MSD is also shown to vary linearly with time at all times. We find that our results reproduce the numerical results of Chubynsky and Slater quite well