Extreme value statistics and arcsine laws for heterogeneous diffusion processes
Prashant Singh
Abstract
Heterogeneous diffusion with a spatially changing diffusion coefficient arises in many experimental systems such as protein dynamics in the cell cytoplasm, mobility of cajal bodies, and confined hard-sphere fluids. Here, we showcase a simple model of heterogeneous diffusion where the diffusion coefficient $D(x)$ varies in a power-law way, i.e., $D(x)\ensuremath{\sim}{|x|}^{\ensuremath{-}\ensuremath{\alpha}}$ with the exponent $\ensuremath{\alpha}>\ensuremath{-}1$. This model is known to exhibit anomalous scaling of the mean-squared displacement (MSD) of the form $\ensuremath{\sim}{t}^{\frac{2}{2+\ensuremath{\alpha}}}$ and weak ergodicity breaking in the sense that ensemble averaged and time averaged MSDs do not converge. In this paper, we look at the extreme value statistics of this model and derive, for all $\ensuremath{\alpha}$, the exact probability distributions of the maximum spatial displacement $M(t)$ and arg-maximum ${t}_{m}(t)$ (i.e., the time at which this maximum is reached) till duration $t$. In the second part of our paper, we analyze the statistical properties of the residence time ${t}_{r}(t)$ and the last-passage time ${t}_{\ensuremath{\ell}}(t)$ and compute their distributions exactly for all values of $\ensuremath{\alpha}$. Our study unravels that the heterogeneous version $(\ensuremath{\alpha}\ensuremath{\ne}0)$ displays many rich and contrasting features compared to that of the standard Brownian motion (BM). For example, while for BM $(\ensuremath{\alpha}=0)$, the distributions of ${t}_{m}(t),\phantom{\rule{0.28em}{0ex}}{t}_{r}(t)$, and ${t}_{\ensuremath{\ell}}(t)$ are all identical (\'a la ``arcsine laws'' due to L\'evy), they turn out to be significantly different for nonzero $\ensuremath{\alpha}$. Another interesting property of ${t}_{r}(t)$ is the existence of a critical $\ensuremath{\alpha}$ (which we denote by ${\ensuremath{\alpha}}_{c}=\ensuremath{-}0.3182$) such that the distribution exhibits a local maximum at ${t}_{r}=t/2$ for $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$ whereas it has minima at ${t}_{r}=t/2$ for $\ensuremath{\alpha}\ensuremath{\ge}{\ensuremath{\alpha}}_{c}$. The underlying reasoning for this difference hints at the very contrasting natures of the process for $\ensuremath{\alpha}\ensuremath{\ge}{\ensuremath{\alpha}}_{c}$ and $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$ which we thoroughly examine in our paper. All our analytical results are backed by extensive numerical simulations.