TY - JOUR A1 - Ślęzak, Jakub A1 - Metzler, Ralf A1 - Magdziarz, Marcin T1 - Codifference can detect ergodicity breaking and non-Gaussianity JF - New Journal of Physics N2 - We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement. KW - diffusion KW - stochastic time series KW - anomalous diffusion Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab13f3 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft CY - Bad Honnef ER - TY - JOUR A1 - Goychuk, Igor A. A1 - Kharchenko, Vasyl O. A1 - Metzler, Ralf T1 - Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion JF - Physical Chemistry Chemical Physics N2 - The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient. KW - royal soc chemistry KW - thomas graham house KW - science park KW - milton rd KW - cambridge cb4 0wf KW - cambs KW - england Y1 - 2014 SN - 1463-9076 IS - 16 SP - 16524 EP - 16535 PB - the Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Bodrova, Anna S. T1 - Ultraslow scaled Brownian motion JF - New journal of physics : the open-access journal for physics N2 - We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form . For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations. KW - anomalous diffusion KW - stochastic processes KW - ageing Y1 - 2015 U6 - https://doi.org/10.1088/1367-2630/17/6/063038 SN - 1367-2630 VL - 17 IS - 063038 PB - Dt. Physikalische Ges., IOP CY - Bad Honnef, London ER - TY - JOUR A1 - Herrmann, Carl J. J. A1 - Metzler, Ralf T1 - A self-avoiding walk with neural delays as a model of fixational eye movements JF - Scientific reports N2 - Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades. Y1 - 2017 U6 - https://doi.org/10.1038/s41598-017-13489-8 SN - 2045-2322 VL - 7 SP - 1 EP - 17 PB - Springer Nature CY - London ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averaging, ageing and delay analysis of financial time series JF - New journal of physics N2 - We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black–Scholes–Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics. KW - time averaging KW - diffusion KW - geometric Brownian motion KW - financial time series Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa7199 SN - 1367-2630 VL - 19 SP - 1 EP - 11 PB - IOP CY - London ER - TY - JOUR A1 - Schwarzl, Maria A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - Quantifying non-ergodicity of anomalous diffusion with higher order moments JF - Scientific reports N2 - Anomalous diffusion is being discovered in a fast growing number of systems. The exact nature of this anomalous diffusion provides important information on the physical laws governing the studied system. One of the central properties analysed for finite particle motion time series is the intrinsic variability of the apparent diffusivity, typically quantified by the ergodicity breaking parameter EB. Here we demonstrate that frequently EB is insufficient to provide a meaningful measure for the observed variability of the data. Instead, important additional information is provided by the higher order moments entering by the skewness and kurtosis. We analyse these quantities for three popular anomalous diffusion models. In particular, we find that even for the Gaussian fractional Brownian motion a significant skewness in the results of physical measurements occurs and needs to be taken into account. Interestingly, the kurtosis and skewness may also provide sensitive estimates of the anomalous diffusion exponent underlying the data. We also derive a new result for the EB parameter of fractional Brownian motion valid for the whole range of the anomalous diffusion parameter. Our results are important for the analysis of anomalous diffusion but also provide new insights into the theory of anomalous stochastic processes. Y1 - 2017 U6 - https://doi.org/10.1038/s41598-017-03712-x VL - 7 PB - Macmillan Publishers Limited CY - London ER - TY - JOUR A1 - Guggenberger, Tobias A1 - Pagnini, Gianni A1 - Vojta, Thomas A1 - Metzler, Ralf T1 - Fractional Brownian motion in a finite interval BT - correlations effect depletion or accretion zones of particles near boundaries JF - New Journal of Physics N2 - Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers. KW - anomalous diffusion KW - fractional Brownian motion KW - reflecting boundary conditions Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab075f SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Single-trajectory spectral analysis of scaled Brownian motion JF - New Journal of Physics N2 - Astandard approach to study time-dependent stochastic processes is the power spectral density (PSD), an ensemble-averaged property defined as the Fourier transform of the autocorrelation function of the process in the asymptotic limit of long observation times, T → ∞. In many experimental situations one is able to garner only relatively few stochastic time series of finite T, such that practically neither an ensemble average nor the asymptotic limit T → ∞ can be achieved. To accommodate for a meaningful analysis of such finite-length data we here develop the framework of single-trajectory spectral analysis for one of the standard models of anomalous diffusion, scaled Brownian motion.Wedemonstrate that the frequency dependence of the single-trajectory PSD is exactly the same as for standard Brownian motion, which may lead one to the erroneous conclusion that the observed motion is normal-diffusive. However, a distinctive feature is shown to be provided by the explicit dependence on the measurement time T, and this ageing phenomenon can be used to deduce the anomalous diffusion exponent.Wealso compare our results to the single-trajectory PSD behaviour of another standard anomalous diffusion process, fractional Brownian motion, and work out the commonalities and differences. Our results represent an important step in establishing singletrajectory PSDs as an alternative (or complement) to analyses based on the time-averaged mean squared displacement. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - single trajectory analysis Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab2f52 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Ślęzak, Jakub A1 - Burnecki, Krzysztof A1 - Metzler, Ralf T1 - Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems JF - New Journal of Physics N2 - Many studies on biological and soft matter systems report the joint presence of a linear mean-squared displacement and a non-Gaussian probability density exhibiting, for instance, exponential or stretched-Gaussian tails. This phenomenon is ascribed to the heterogeneity of the medium and is captured by random parameter models such as ‘superstatistics’ or ‘diffusing diffusivity’. Independently, scientists working in the area of time series analysis and statistics have studied a class of discrete-time processes with similar properties, namely, random coefficient autoregressive models. In this work we try to reconcile these two approaches and thus provide a bridge between physical stochastic processes and autoregressive models.Westart from the basic Langevin equation of motion with time-varying damping or diffusion coefficients and establish the link to random coefficient autoregressive processes. By exploring that link we gain access to efficient statistical methods which can help to identify data exhibiting Brownian yet non-Gaussian diffusion. KW - diffusion KW - Langevin equation KW - Brownian yet non-Gaussian diffusion KW - diffusing diffusivity KW - superstatistics KW - autoregressive models KW - time series analysis KW - codifference Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab3366 SN - 1367-2630 VL - 21 PB - Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics CY - Bad Honnef und London ER - TY - JOUR A1 - Awad, Emad A1 - Metzler, Ralf T1 - Closed-form multi-dimensional solutions and asymptotic behaviours for subdiffusive processes with crossovers: II. Accelerating case JF - Journal of physics : A, Mathematical and theoretical N2 - Anomalous diffusion with a power-law time dependence vertical bar R vertical bar(2)(t) similar or equal to t(alpha i) of the mean squared displacement occurs quite ubiquitously in numerous complex systems. Often, this anomalous diffusion is characterised by crossovers between regimes with different anomalous diffusion exponents alpha(i). Here we consider the case when such a crossover occurs from a first regime with alpha(1) to a second regime with alpha(2) such that alpha(2) > alpha(1), i.e., accelerating anomalous diffusion. A widely used framework to describe such crossovers in a one-dimensional setting is the bi-fractional diffusion equation of the so-called modified type, involving two time-fractional derivatives defined in the Riemann-Liouville sense. We here generalise this bi-fractional diffusion equation to higher dimensions and derive its multidimensional propagator (Green's function) for the general case when also a space fractional derivative is present, taking into consideration long-ranged jumps (Levy flights). We derive the asymptotic behaviours for this propagator in both the short- and long-time as well the short- and long-distance regimes. Finally, we also calculate the mean squared displacement, skewness and kurtosis in all dimensions, demonstrating that in the general case the non-Gaussian shape of the probability density function changes. KW - multidimensional fractional diffusion equation KW - continuous time random KW - walks KW - crossover anomalous diffusion dynamics KW - non-Gaussian probability KW - density Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac5a90 SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 20 PB - IOP Publ. Ltd. CY - Bristol ER -