TY - JOUR A1 - Cherstvy, Andrey G. A1 - Nagel, Oliver A1 - Beta, Carsten A1 - Metzler, Ralf T1 - Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - What is the underlying diffusion process governing the spreading dynamics and search strategies employed by amoeboid cells? Based on the statistical analysis of experimental single-cell tracking data of the two-dimensional motion of the Dictyostelium discoideum amoeboid cells, we quantify their diffusive behaviour based on a number of standard and complementary statistical indicators. We compute the ensemble- and time-averaged mean-squared displacements (MSDs) of the diffusing amoebae cells and observe a pronounced spread of short-time diffusion coefficients and anomalous MSD-scaling exponents for individual cells. The distribution functions of the cell displacements, the long-tailed distribution of instantaneous speeds, and the velocity autocorrelations are also computed. In particular, we observe a systematic superdiffusive short-time behaviour for the ensemble- and time-averaged MSDs of the amoeboid cells. Also, a clear anti-correlation of scaling exponents and generalised diffusivity values for different cells is detected. Most significantly, we demonstrate that the distribution function of the cell displacements has a strongly non-Gaussian shape andusing a rescaled spatio-temporal variablethe cell-displacement data collapse onto a universal master curve. The current analysis of single-cell motions can be implemented for quantifying diffusive behaviours in other living-matter systems, in particular, when effects of active transport, non-Gaussian displacements, and heterogeneity of the population are involved in the dynamics. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp04254c SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 35 SP - 23034 EP - 23054 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Sposini, Vittoria A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - First passage statistics for diffusing diffusivity JF - Journal of physics : A, Mathematical and theoretical N2 - A rapidly increasing number of systems is identified in which the stochastic motion of tracer particles follows the Brownian law < r(2)(t)> similar or equal to Dt yet the distribution of particle displacements is strongly non-Gaussian. A central approach to describe this effect is the diffusing diffusivity (DD) model in which the diffusion coefficient itself is a stochastic quantity, mimicking heterogeneities of the environment encountered by the tracer particle on its path. We here quantify in terms of analytical and numerical approaches the first passage behaviour of the DD model. We observe significant modifications compared to Brownian-Gaussian diffusion, in particular that the DD model may have a faster first passage dynamics. Moreover we find a universal crossover point of the survival probability independent of the initial condition. KW - diffusion KW - superstatistics KW - first passage Y1 - 2018 U6 - https://doi.org/10.1088/1751-8121/aaf6ff SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 4 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Dybiec, Bartlomiej A1 - Capala, Karol A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Conservative random walks in confining potentials JF - Journal of physics : A, Mathematical and theoretical N2 - Levy walks are continuous time random walks with spatio-temporal coupling of jump lengths and waiting times, often used to model superdiffusive spreading processes such as animals searching for food, tracer motion in weakly chaotic systems, or even the dynamics in quantum systems such as cold atoms. In the simplest version Levy walks move with a finite speed. Here, we present an extension of the Levy walk scenario for the case when external force fields influence the motion. The resulting motion is a combination of the response to the deterministic force acting on the particle, changing its velocity according to the principle of total energy conservation, and random velocity reversals governed by the distribution of waiting times. For the fact that the motion stays conservative, that is, on a constant energy surface, our scenario is fundamentally different from thermal motion in the same external potentials. In particular, we present results for the velocity and position distributions for single well potentials of different steepness. The observed dynamics with its continuous velocity changes enriches the theory of Levy walk processes and will be of use in a variety of systems, for which the particles are externally confined. KW - Levy walk KW - conservative random walks KW - Levy flight Y1 - 2018 U6 - https://doi.org/10.1088/1751-8121/aaefc2 SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 1 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Hou, Ru A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Akimoto, Takuma T1 - Biased continuous-time random walks for ordinary and equilibrium cases BT - facilitation of diffusion, ergodicity breaking and ageing JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent alpha are considered for the WTD psi(t) similar to 1/t(1+alpha). We treat continuous-time random walks (CTRWs) with 0 < alpha < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < alpha < 2 and alpha > 2. We demonstrate that in the presence of a drift CTRWs with alpha < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < alpha < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with alpha > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < alpha < 2 and alpha > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with alpha > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent alpha. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp01863d SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 32 SP - 20827 EP - 20848 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averages and their statistical variation for the Ornstein-Uhlenbeck process BT - Role of initial particle distributions and relaxation to stationarity JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and outof-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022134 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Aydiner, Ekrem A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Wealth distribution, Pareto law, and stretched exponential decay of money BT - Computer simulations analysis of agent-based models JF - Physica : europhysics journal ; A, Statistical mechanics and its applications N2 - We study by Monte Carlo simulations a kinetic exchange trading model for both fixed and distributed saving propensities of the agents and rationalize the person and wealth distributions. We show that the newly introduced wealth distribution – that may be more amenable in certain situations – features a different power-law exponent, particularly for distributed saving propensities of the agents. For open agent-based systems, we analyze the person and wealth distributions and find that the presence of trap agents alters their amplitude, leaving however the scaling exponents nearly unaffected. For an open system, we show that the total wealth – for different trap agent densities and saving propensities of the agents – decreases in time according to the classical Kohlrausch–Williams–Watts stretched exponential law. Interestingly, this decay does not depend on the trap agent density, but rather on saving propensities. The system relaxation for fixed and distributed saving schemes are found to be different. KW - Econophysics KW - Wealth and income distribution KW - Pareto law KW - Scaling exponents Y1 - 2017 U6 - https://doi.org/10.1016/j.physa.2017.08.017 SN - 0378-4371 SN - 1873-2119 VL - 490 SP - 278 EP - 288 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Krapf, Diego A1 - Marinari, Enzo A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Xu, Xinran A1 - Squarcini, Alessio T1 - Power spectral density of a single Brownian trajectory BT - what one can and cannot learn from it JF - New journal of physics : the open-access journal for physics N2 - The power spectral density (PSD) of any time-dependent stochastic processX (t) is ameaningful feature of its spectral content. In its text-book definition, the PSD is the Fourier transform of the covariance function of X-t over an infinitely large observation timeT, that is, it is defined as an ensemble-averaged property taken in the limitT -> infinity. Alegitimate question is what information on the PSD can be reliably obtained from single-trajectory experiments, if one goes beyond the standard definition and analyzes the PSD of a single trajectory recorded for a finite observation timeT. In quest for this answer, for a d-dimensional Brownian motion (BM) we calculate the probability density function of a single-trajectory PSD for arbitrary frequency f, finite observation time T and arbitrary number k of projections of the trajectory on different axes. We show analytically that the scaling exponent for the frequency-dependence of the PSD specific to an ensemble of BM trajectories can be already obtained from a single trajectory, while the numerical amplitude in the relation between the ensemble-averaged and single-trajectory PSDs is afluctuating property which varies from realization to realization. The distribution of this amplitude is calculated exactly and is discussed in detail. Our results are confirmed by numerical simulations and single-particle tracking experiments, with remarkably good agreement. In addition we consider a truncated Wiener representation of BM, and the case of a discrete-time lattice random walk. We highlight some differences in the behavior of a single-trajectory PSD for BM and for the two latter situations. The framework developed herein will allow for meaningful physical analysis of experimental stochastic trajectories. KW - power spectral density KW - single-trajectory analysis KW - probability density function KW - exact results Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aaa67c SN - 1367-2630 VL - 20 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Estrada, Ernesto A1 - Delvenne, Jean-Charles A1 - Hatano, Naomichi A1 - Mateos, Jose L. A1 - Metzler, Ralf A1 - Riascos, Alejandro P. A1 - Schaub, Michael T. T1 - Random multi-hopper model BT - super-fast random walks on graphs JF - Journal of Complex Networks N2 - We develop a mathematical model considering a random walker with long-range hops on arbitrary graphs. The random multi-hopper can jump to any node of the graph from an initial position, with a probability that decays as a function of the shortest-path distance between the two nodes in the graph. We consider here two decaying functions in the form of Laplace and Mellin transforms of the shortest-path distances. We prove that when the parameters of these transforms approach zero asymptotically, the hitting time in the multi-hopper approaches the minimum possible value for a normal random walker. We show by computational experiments that the multi-hopper explores a graph with clusters or skewed degree distributions more efficiently than a normal random walker. We provide computational evidences of the advantages of the random multi-hopper model with respect to the normal random walk by studying deterministic, random and real-world networks. Y1 - 2018 U6 - https://doi.org/10.1093/comnet/cnx043 SN - 2051-1310 SN - 2051-1329 VL - 6 IS - 3 SP - 382 EP - 403 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Akimoto, Takuma A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity, rejuvenation, enhancement, and slow relaxation of diffusion in biased continuous-time random walks JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022105 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Mardoukhi, Yousof A1 - Jeon, Jae-Hyung A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Fluctuations of random walks in critical random environments JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Percolation networks have been widely used in the description of porous media but are now found to be relevant to understand the motion of particles in cellular membranes or the nucleus of biological cells. Random walks on the infinite cluster at criticality of a percolation network are asymptotically ergodic. On any finite size cluster of the network stationarity is reached at finite times, depending on the cluster's size. Despite of this we here demonstrate by combination of analytical calculations and simulations that at criticality the disorder and cluster size average of the ensemble of clusters leads to a non-vanishing variance of the time averaged mean squared displacement, regardless of the measurement time. Fluctuations of this relevant experimental quantity due to the disorder average of such ensembles are thus persistent and non-negligible. The relevance of our results for single particle tracking analysis in complex and biological systems is discussed. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp03212b SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 31 SP - 20427 EP - 20438 PB - Royal Society of Chemistry CY - Cambridge ER -