@article{ŚlęzakMetzlerMagdziarz2019, author = {Ślęzak, Jakub and Metzler, Ralf and Magdziarz, Marcin}, title = {Codifference can detect ergodicity breaking and non-Gaussianity}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab13f3}, pages = {25}, year = {2019}, abstract = {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.}, language = {en} } @article{GoychukKharchenkoMetzler2014, author = {Goychuk, Igor A. and Kharchenko, Vasyl O. and Metzler, Ralf}, title = {Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion}, series = {Physical Chemistry Chemical Physics}, journal = {Physical Chemistry Chemical Physics}, number = {16}, publisher = {the Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, pages = {16524 -- 16535}, year = {2014}, abstract = {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.}, language = {en} } @article{MetzlerCherstvyChechkinetal.2015, author = {Metzler, Ralf and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Bodrova, Anna S.}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, number = {063038}, publisher = {Dt. Physikalische Ges., IOP}, address = {Bad Honnef, London}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, year = {2015}, abstract = {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.}, language = {en} } @article{HerrmannMetzler2017, author = {Herrmann, Carl J. J. and Metzler, Ralf}, title = {A self-avoiding walk with neural delays as a model of fixational eye movements}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-13489-8}, pages = {1 -- 17}, year = {2017}, abstract = {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.}, language = {en} } @article{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, series = {New journal of physics}, volume = {19}, journal = {New journal of physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa7199}, pages = {1 -- 11}, year = {2017}, abstract = {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.}, language = {en} } @article{SchwarzlGodecMetzler2017, author = {Schwarzl, Maria and Godec, Aljaž and Metzler, Ralf}, title = {Quantifying non-ergodicity of anomalous diffusion with higher order moments}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Macmillan Publishers Limited}, address = {London}, doi = {10.1038/s41598-017-03712-x}, pages = {18}, year = {2017}, abstract = {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.}, language = {en} } @article{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab075f}, pages = {13}, year = {2019}, abstract = {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.}, language = {en} } @article{SposiniMetzlerOshanin2019, author = {Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb}, title = {Single-trajectory spectral analysis of scaled Brownian motion}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab2f52}, pages = {16}, year = {2019}, abstract = {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.}, language = {en} } @article{ŚlęzakBurneckiMetzler2019, author = {Ślęzak, Jakub and Burnecki, Krzysztof and Metzler, Ralf}, title = {Random coefficient autoregressive processes describe Brownian yet non-Gaussian diffusion in heterogeneous systems}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab3366}, pages = {18}, year = {2019}, abstract = {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.}, language = {en} } @article{AwadMetzler2022, author = {Awad, Emad and Metzler, Ralf}, title = {Closed-form multi-dimensional solutions and asymptotic behaviours for subdiffusive processes with crossovers: II. Accelerating case}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {20}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac5a90}, pages = {29}, year = {2022}, abstract = {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.}, language = {en} }