TY - JOUR A1 - Varghese, Alan J. A1 - Chechkin, Aleksei A1 - Metzler, Ralf A1 - Sujith, Raman I. T1 - Capturing multifractality of pressure fluctuations in thermoacoustic systems using fractional-order derivatives JF - Chaos : an interdisciplinary journal of nonlinear science N2 - The stable operation of a turbulent combustor is not completely silent; instead, there is a background of small amplitude aperiodic acoustic fluctuations known as combustion noise. Pressure fluctuations during this state of combustion noise are multifractal due to the presence of multiple temporal scales that contribute to its dynamics. However, existing models are unable to capture the multifractality in the pressure fluctuations. We conjecture an underlying fractional dynamics for the thermoacoustic system and obtain a fractional-order model for pressure fluctuations. The data from this model has remarkable visual similarity to the experimental data and also has a wide multifractal spectrum during the state of combustion noise. Quantitative similarity with the experimental data in terms of the Hurst exponent and the multifractal spectrum is observed during the state of combustion noise. This model is also able to produce pressure fluctuations that are qualitatively similar to the experimental data acquired during intermittency and thermoacoustic instability. Furthermore, we argue that the fractional dynamics vanish as we approach the state of thermoacoustic instability. Y1 - 2021 U6 - https://doi.org/10.1063/5.0032585 SN - 1054-1500 SN - 1089-7682 VL - 31 IS - 3 PB - American Institute of Physics, AIP CY - Melville ER - TY - JOUR A1 - Vilk, Ohad A1 - Aghion, Erez A1 - Avgar, Tal A1 - Beta, Carsten A1 - Nagel, Oliver A1 - Sabri, Adal A1 - Sarfati, Raphael A1 - Schwartz, Daniel K. A1 - Weiß, Matthias A1 - Krapf, Diego A1 - Nathan, Ran A1 - Metzler, Ralf A1 - Assaf, Michael T1 - Unravelling the origins of anomalous diffusion BT - from molecules to migrating storks JF - Physical review research / American Physical Society N2 - Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations (“Joseph effect”), fat-tailed probability density of increments (“Noah effect”), and nonstationarity (“Moses effect”). We show that such a decomposition of real-life data allows us to infer nontrivial behavioral predictions and to resolve open questions in the fields of single-particle tracking in living cells and movement ecology. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.033055 SN - 2643-1564 VL - 4 IS - 3 PB - American Physical Society CY - College Park, MD ER - TY - JOUR A1 - Vilk, Ohad A1 - Aghion, Erez A1 - Nathan, Ran A1 - Toledo, Sivan A1 - Metzler, Ralf A1 - Assaf, Michael T1 - Classification of anomalous diffusion in animal movement data using power spectral analysis JF - Journal of physics : A, Mathematical and theoretical N2 - The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find 1/f noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations. KW - diffusion KW - anomalous diffusion KW - power spectral analysis KW - ecological KW - movement data Y1 - 2022 U6 - https://doi.org/10.1088/1751-8121/ac7e8f SN - 1751-8113 SN - 1751-8121 VL - 55 IS - 33 PB - IOP Publishing CY - Bristol ER - TY - JOUR A1 - Vinod, Deepak A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Time-averaging and nonergodicity of reset geometric Brownian motion with drift JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How do near-bankruptcy events in the past affect the dynamics of stock-market prices in the future? Specifically, what are the long-time properties of a time-local exponential growth of stock-market prices under the influence of stochastically occurring economic crashes? Here, we derive the ensemble- and time-averaged properties of the respective "economic" or geometric Brownian motion (GBM) with a nonzero drift exposed to a Poissonian constant-rate price-restarting process of "resetting." We examine-based both on thorough analytical calculations and on findings from systematic stochastic computer simulations-the general situation of reset GBM with a nonzero [positive] drift and for all special cases emerging for varying parameters of drift, volatility, and reset rate in the model. We derive and summarize all short- and long-time dependencies for the mean-squared displacement (MSD), the variance, and the mean time-averaged MSD (TAMSD) of the process of Poisson-reset GBM under the conditions of both rare and frequent resetting. We consider three main regions of model parameters and categorize the crossovers between different functional behaviors of the statistical quantifiers of this process. The analytical relations are fully supported by the results of computer simulations. In particular, we obtain that Poisson-reset GBM is a nonergodic stochastic process, with generally MSD(Delta) not equal TAMSD(Delta) and Variance(Delta) not equal TAMSD(Delta) at short lag times Delta and for long trajectory lengths T. We investigate the behavior of the ergodicity-breaking parameter in each of the three regions of parameters and examine its dependence on the rate of reset at Delta/T << 1. Applications of these theoretical results to the analysis of prices of reset-containing options are pertinent. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.106.034137 SN - 2470-0045 SN - 2470-0053 VL - 106 IS - 3 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Vinod, Deepak A1 - Cherstvy, Andrey G. A1 - Wang, Wei A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Nonergodicity of reset geometric Brownian motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We derive. the ensemble-and time-averaged mean-squared displacements (MSD, TAMSD) for Poisson-reset geometric Brownian motion (GBM), in agreement with simulations. We find MSD and TAMSD saturation for frequent resetting, quantify the spread of TAMSDs via the ergodicity-breaking parameter and compute distributions of prices. General MSD-TAMSD nonequivalence proves reset GBM nonergodic. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.105.L012106 SN - 2470-0045 SN - 2470-0053 VL - 105 IS - 1 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Vojta, Thomas A1 - Skinner, Sarah A1 - Metzler, Ralf T1 - Probability density of the fractional Langevin equation with reflecting walls JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We investigate anomalous diffusion processes governed by the fractional Langevin equation and confined to a finite or semi-infinite interval by reflecting potential barriers. As the random and damping forces in the fractional Langevin equation fulfill the appropriate fluctuation-dissipation relation, the probability density on a finite interval converges for long times towards the expected uniform distribution prescribed by thermal equilibrium. In contrast, on a semi-infinite interval with a reflecting wall at the origin, the probability density shows pronounced deviations from the Gaussian behavior observed for normal diffusion. If the correlations of the random force are persistent (positive), particles accumulate at the reflecting wall while antipersistent (negative) correlations lead to a depletion of particles near the wall. We compare and contrast these results with the strong accumulation and depletion effects recently observed for nonthermal fractional Brownian motion with reflecting walls, and we discuss broader implications. Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.042142 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Thapa, Samudrajit A1 - Seno, Flavio A1 - Liu, Xianbin A1 - Metzler, Ralf T1 - Fractional Brownian motion with random diffusivity BT - emerging residual nonergodicity below the correlation time JF - Journal of physics : A, Mathematical and theoretical N2 - Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments. KW - stochastic processes KW - anomalous diffusion KW - fractional Brownian motion KW - diffusing diffusivity KW - weak ergodicity breaking Y1 - 2020 U6 - https://doi.org/10.1088/1751-8121/aba467 SN - 1751-8113 SN - 1751-8121 VL - 53 IS - 47 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Kantz, Holger A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Time averaging and emerging nonergodicity upon resetting of fractional Brownian motion and heterogeneous diffusion processes JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x) = D0|x|gamma and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicitybreaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB similar to(1/r )-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.104.024105 SN - 2470-0045 SN - 2470-0053 VL - 104 IS - 2 PB - American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Liu, Xianbin A1 - Metzler, Ralf T1 - Anomalous diffusion and nonergodicity for heterogeneous diffusion processes with fractional Gaussian noise JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Heterogeneous diffusion processes (HDPs) feature a space-dependent diffusivity of the form D(x) = D-0|x|(alpha). Such processes yield anomalous diffusion and weak ergodicity breaking, the asymptotic disparity between ensemble and time averaged observables, such as the mean-squared displacement. Fractional Brownian motion (FBM) with its long-range correlated yet Gaussian increments gives rise to anomalous and ergodic diffusion. Here, we study a combined model of HDPs and FBM to describe the particle dynamics in complex systems with position-dependent diffusivity driven by fractional Gaussian noise. This type of motion is, inter alia, relevant for tracer-particle diffusion in biological cells or heterogeneous complex fluids. We show that the long-time scaling behavior predicted theoretically and by simulations for the ensemble-and time-averaged mean-squared displacements couple the scaling exponents alpha of HDPs and the Hurst exponent H of FBM in a characteristic way. Our analysis of the simulated data in terms of the rescaled variable y similar to |x|(1/(2/(2-alpha)))/t(H) coupling particle position x and time t yields a simple, Gaussian probability density function (PDF), PHDP-FBM(y) = e(-y2)/root pi. Its universal shape agrees well with theoretical predictions for both uni- and bimodal PDF distributions. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.012146 SN - 2470-0045 SN - 2470-0053 SN - 1063-651X SN - 1539-3755 SN - 2470-0061 VL - 102 IS - 1 SP - 012146-1 EP - 012146-16 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Wang, Wei A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Restoring ergodicity of stochastically reset anomalous-diffusion processes JF - Physical Review Research N2 - How do different reset protocols affect ergodicity of a diffusion process in single-particle-tracking experiments? We here address the problem of resetting of an arbitrary stochastic anomalous-diffusion process (ADP) from the general mathematical points of view and assess ergodicity of such reset ADPs for an arbitrary resetting protocol. The process of stochastic resetting describes the events of the instantaneous restart of a particle’s motion via randomly distributed returns to a preset initial position (or a set of those). The waiting times of such resetting events obey the Poissonian, Gamma, or more generic distributions with specified conditions regarding the existence of moments. Within these general approaches, we derive general analytical results and support them by computer simulations for the behavior of the reset mean-squared displacement (MSD), the new reset increment-MSD (iMSD), and the mean reset time-averaged MSD (TAMSD). For parental nonreset ADPs with the MSD(t)∝ tμ we find a generic behavior and a switch of the short-time growth of the reset iMSD and mean reset TAMSDs from ∝ _μ for subdiffusive to ∝ _1 for superdiffusive reset ADPs. The critical condition for a reset ADP that recovers its ergodicity is found to be more general than that for the nonequilibrium stationary state, where obviously the iMSD and the mean TAMSD are equal. The consideration of the new statistical quantifier, the iMSD—as compared to the standard MSD—restores the ergodicity of an arbitrary reset ADP in all situations when the μth moment of the waiting-time distribution of resetting events is finite. Potential applications of these new resetting results are, inter alia, in the area of biophysical and soft-matter systems. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.013161 SN - 2643-1564 VL - 4 SP - 013161-1 EP - 013161-13 PB - American Physical Society CY - College Park, Maryland, United States ET - 1 ER - TY - JOUR A1 - Wang, Wei A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Anomalous diffusion, aging, and nonergodicity of scaled Brownian motion with fractional Gaussian noise: overview of related experimental observations and models JF - Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies N2 - How does a systematic time-dependence of the diffusion coefficient D(t) affect the ergodic and statistical characteristics of fractional Brownian motion (FBM)? Here, we answer this question via studying the characteristics of a set of standard statistical quantifiers relevant to single-particle-tracking (SPT) experiments. We examine, for instance, how the behavior of the ensemble- and time-averaged mean-squared displacements-denoted as the standard MSD < x(2)(Delta)> and TAMSD <<(delta(2)(Delta))over bar>> quantifiers-of FBM featuring < x(2) (Delta >> = <<(delta(2)(Delta >)over bar>> proportional to Delta(2H) (where H is the Hurst exponent and Delta is the [lag] time) changes in the presence of a power-law deterministically varying diffusivity D-proportional to(t) proportional to t(alpha-1) -germane to the process of scaled Brownian motion (SBM)-determining the strength of fractional Gaussian noise. The resulting compound "scaled-fractional" Brownian motion or FBM-SBM is found to be nonergodic, with < x(2)(Delta >> proportional to Delta(alpha+)(2H)(-1) and <(delta 2(Delta >) over bar > proportional to Delta(2H). We also detect a stalling behavior of the MSDs for very subdiffusive SBM and FBM, when alpha + 2H - 1 < 0. The distribution of particle displacements for FBM-SBM remains Gaussian, as that for the parent processes of FBM and SBM, in the entire region of scaling exponents (0 < alpha < 2 and 0 < H < 1). The FBM-SBM process is aging in a manner similar to SBM. The velocity autocorrelation function (ACF) of particle increments of FBM-SBM exhibits a dip when the parent FBM process is subdiffusive. Both for sub- and superdiffusive FBM contributions to the FBM-SBM process, the SBM exponent affects the long-time decay exponent of the ACF. Applications of the FBM-SBM-amalgamated process to the analysis of SPT data are discussed. A comparative tabulated overview of recent experimental (mainly SPT) and computational datasets amenable for interpretation in terms of FBM-, SBM-, and FBM-SBM-like models of diffusion culminates the presentation. The statistical aspects of the dynamics of a wide range of biological systems is compared in the table, from nanosized beads in living cells, to chromosomal loci, to water diffusion in the brain, and, finally, to patterns of animal movements. Y1 - 2022 U6 - https://doi.org/10.1039/d2cp01741e SN - 1463-9076 SN - 1463-9084 VL - 24 IS - 31 SP - 18482 EP - 18504 PB - RSC Publ. CY - Cambridge ER - TY - JOUR A1 - Wang, Wei A1 - Seno, Flavio A1 - Sokolov, Igor M. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Unexpected crossovers in correlated random-diffusivity processes JF - New Journal of Physics N2 - The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by 'viscoelastic' anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of 'viscoelastic yet non-Gaussian' diffusion. KW - diffusion KW - anomalous diffusion KW - non-Gaussianity KW - fractional Brownian motion Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/aba390 SN - 1367-2630 VL - 22 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Xu, Pengbo A1 - Metzler, Ralf A1 - Wang, Wanli T1 - Infinite density and relaxation for Levy walks in an external potential BT - Hermite polynomial approach JF - Physical review N2 - Levy walks are continuous-time random-walk processes with a spatiotemporal coupling of jump lengths and waiting times. We here apply the Hermite polynomial method to study the behavior of LWs with power-law walking time density for four different cases. First we show that the known result for the infinite density of an unconfined, unbiased LW is consistently recovered. We then derive the asymptotic behavior of the probability density function (PDF) for LWs in a constant force field, and we obtain the corresponding qth-order moments. In a harmonic external potential we derive the relaxation dynamic of the LW. For the case of a Poissonian walking time an exponential relaxation behavior is shown to emerge. Conversely, a power-law decay is obtained when the mean walking time diverges. Finally, we consider the case of an unconfined, unbiased LW with decaying speed v(r ) = v0/./r. When the mean walking time is finite, a universal Gaussian law for the position-PDF of the walker is obtained explicitly. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.105.044118 SN - 2470-0045 SN - 2470-0053 VL - 105 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Xu, Pengbo A1 - Zhou, Tian A1 - Metzler, Ralf A1 - Deng, Weihua T1 - Lévy walk dynamics in an external harmonic potential JF - Physical review : E, Statistical, nonlinear, and soft matter physics N2 - Levy walks (LWs) are spatiotemporally coupled random-walk processes describing superdiffusive heat conduction in solids, propagation of light in disordered optical materials, motion of molecular motors in living cells, or motion of animals, humans, robots, and viruses. We here investigate a key feature of LWs-their response to an external harmonic potential. In this generic setting for confined motion we demonstrate that LWs equilibrate exponentially and may assume a bimodal stationary distribution. We also show that the stationary distribution has a horizontal slope next to a reflecting boundary placed at the origin, in contrast to correlated superdiffusive processes. Our results generalize LWs to confining forces and settle some longstanding puzzles around LWs. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.101.062127 SN - 2470-0045 SN - 2470-0053 SN - 1550-2376 SN - 1063-651X SN - 1539-3755 VL - 101 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Xu, Pengbo A1 - Zhou, Tian A1 - Metzler, Ralf A1 - Deng, Weihua T1 - Stochastic harmonic trapping of a Lévy walk BT - transport and first-passage dynamics under soft resetting strategies JF - New journal of physics : the open-access journal for physics / Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics N2 - We introduce and study a Lévy walk (LW) model of particle spreading with a finite propagation speed combined with soft resets, stochastically occurring periods in which an harmonic external potential is switched on and forces the particle towards a specific position. Soft resets avoid instantaneous relocation of particles that in certain physical settings may be considered unphysical. Moreover, soft resets do not have a specific resetting point but lead the particle towards a resetting point by a restoring Hookean force. Depending on the exact choice for the LW waiting time density and the probability density of the periods when the harmonic potential is switched on, we demonstrate a rich emerging response behaviour including ballistic motion and superdiffusion. When the confinement periods of the soft-reset events are dominant, we observe a particle localisation with an associated non-equilibrium steady state. In this case the stationary particle probability density function turns out to acquire multimodal states. Our derivations are based on Markov chain ideas and LWs with multiple internal states, an approach that may be useful and flexible for the investigation of other generalised random walks with soft and hard resets. The spreading efficiency of soft-rest LWs is characterised by the first-passage time statistic. KW - diffusion KW - anomalous diffusion KW - stochastic resetting KW - Levy walks Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac5282 SN - 1367-2630 VL - 24 IS - 3 SP - 1 EP - 28 PB - Deutsche Physikalische Gesellschaft CY - Bad Honnef ER - TY - JOUR A1 - Xu, Yong A1 - Liu, Xuemei A1 - Li, Yongge A1 - Metzler, Ralf T1 - Heterogeneous diffusion processes and nonergodicity with Gaussian colored noise in layered diffusivity landscapes JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent a. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent alpha is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.062106 SN - 2470-0045 SN - 2470-0053 VL - 102 IS - 6 PB - American Physical Society CY - College Park 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 - Ślęzak, Jakub A1 - Metzler, Ralf A1 - Magdziarz, Marcin T1 - Superstatistical generalised Langevin equation BT - non-Gaussian viscoelastic anomalous diffusion JF - New Journal of Physics N2 - Recent advances in single particle tracking and supercomputing techniques demonstrate the emergence of normal or anomalous, viscoelastic diffusion in conjunction with non-Gaussian distributions in soft, biological, and active matter systems. We here formulate a stochastic model based on a generalised Langevin equation in which non-Gaussian shapes of the probability density function and normal or anomalous diffusion have a common origin, namely a random parametrisation of the stochastic force. We perform a detailed analysis demonstrating how various types of parameter distributions for the memory kernel result in exponential, power law, or power-log law tails of the memory functions. The studied system is also shown to exhibit a further unusual property: the velocity has a Gaussian one point probability density but non-Gaussian joint distributions. This behaviour is reflected in the relaxation from a Gaussian to a non-Gaussian distribution observed for the position variable. We show that our theoretical results are in excellent agreement with stochastic simulations. KW - anomalous diffusion KW - generalised langevin equation KW - superstatistics KW - non-Gaussian diffusion Y1 - 2018 U6 - https://doi.org/10.1088/1367-2630/aaa3d4 SN - 1367-2630 VL - 20 IS - 023026 SP - 1 EP - 25 PB - Deutsche Physikalische Gesellschaft / Institute of Physics CY - Bad Honnef und London ER - 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 -