@article{CherstvySafdariMetzler2021,
  author    = {Cherstvy, Andrey G. and Safdari, Hadiseh and Metzler, Ralf},
  title     = {Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity},
  series = {Journal of physics. D, Applied physics},
  volume    = {54},
  journal   = {Journal of physics. D, Applied physics},
  number    = {19},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0022-3727},
  doi       = {10.1088/1361-6463/abdff0},
  pages     = {18},
  year      = {2021},
  abstract  = {We investigate a diffusion process with a time-dependent diffusion coefficient, both exponentially increasing and decreasing in time, D(t)=D-0(e +/- 2 alpha t). For this (hypothetical) nonstationary diffusion process we compute-both analytically and from extensive stochastic simulations-the behavior of the ensemble- and time-averaged mean-squared displacements (MSDs) of the particles, both in the over- and underdamped limits. Simple asymptotic relations derived for the short- and long-time behaviors are shown to be in excellent agreement with the results of simulations. The diffusive characteristics in the presence of ageing are also considered, with dramatic differences of the over- versus underdamped regime. Our results for D(t)=D-0(e +/- 2 alpha t) extend and generalize the class of diffusive systems obeying scaled Brownian motion featuring a power-law-like variation of the diffusivity with time, D(t) similar to t(alpha-1). We also examine the logarithmically increasing diffusivity, D(t)=D(0)log[t/tau(0)], as another fundamental functional dependence (in addition to the power-law and exponential) and as an example of diffusivity slowly varying in time. One of the main conclusions is that the behavior of the massive particles is predominantly ergodic, while weak ergodicity breaking is repeatedly found for the time-dependent diffusion of the massless particles at short times. The latter manifests itself in the nonequivalence of the (both nonaged and aged) MSD and the mean time-averaged MSD. The current findings are potentially applicable to a class of physical systems out of thermal equilibrium where a rapid increase or decrease of the particles' diffusivity is inherently realized. One biological system potentially featuring all three types of time-dependent diffusion (power-law-like, exponential, and logarithmic) is water diffusion in the brain tissues, as we thoroughly discuss in the end.},
  language  = {en}
}
@article{MutothyaXuLietal.2021,
  author    = {Mutothya, Nicholas Mwilu and Xu, Yong and Li, Yongge and Metzler, Ralf},
  title     = {Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {54},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {29},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/abfba6},
  pages     = {31},
  year      = {2021},
  abstract  = {We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments.},
  language  = {en}
}
@article{VilkAghionNathanetal.2022,
  author    = {Vilk, Ohad and Aghion, Erez and Nathan, Ran and Toledo, Sivan and Metzler, Ralf and Assaf, Michael},
  title     = {Classification of anomalous diffusion in animal movement data using power spectral analysis},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {55},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {33},
  publisher = {IOP Publishing},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/ac7e8f},
  pages     = {16},
  year      = {2022},
  abstract  = {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.},
  language  = {en}
}
@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{AwadMetzler2020,
  author    = {Awad, Emad and Metzler, Ralf},
  title     = {Crossover dynamics from superdiffusion to subdiffusion},
  series = {Fractional calculus and applied analysis : an international journal for theory and applications},
  volume    = {23},
  journal   = {Fractional calculus and applied analysis : an international journal for theory and applications},
  number    = {1},
  publisher = {De Gruyter},
  address   = {Berlin ; Boston},
  issn      = {1311-0454},
  doi       = {10.1515/fca-2020-0003},
  pages     = {55 -- 102},
  year      = {2020},
  abstract  = {The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function.},
  language  = {en}
}
@article{MolinaGarciaSandevSafdarietal.2018,
  author    = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf},
  title     = {Crossover from anomalous to normal diffusion},
  series = {New Journal of Physics},
  volume    = {20},
  journal   = {New Journal of Physics},
  publisher = {IOP Publishing Ltd},
  address   = {London und Bad Honnef},
  issn      = {1367-2630},
  doi       = {10.1088/1367-2630/aae4b2},
  pages     = {28},
  year      = {2018},
  abstract  = {Abstract The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive-diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive-diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion-diffusion and subdiffusion-subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling.},
  language  = {en}
}
@article{BauerGodecMetzler2014,
  author    = {Bauer, Maximilian and Godec, Aljaž and Metzler, Ralf},
  title     = {Diffusion of finite-size particles in two-dimensional channels with random wall configurations},
  series = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies},
  volume    = {16},
  journal   = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies},
  number    = {13},
  publisher = {RSC Publications},
  address   = {Cambridge},
  issn      = {1463-9084},
  doi       = {10.1039/C3CP55160A},
  pages     = {6118 -- 6128},
  year      = {2014},
  abstract  = {Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].},
  language  = {en}
}
@article{SliusarenkoVitaliSposinietal.2019,
  author    = {Sliusarenko, Oleksii Yu and Vitali, Silvia and Sposini, Vittoria and Paradisi, Paolo and Chechkin, Aleksei V. and Castellani, Gastone and Pagnini, Gianni},
  title     = {Finite-energy Levy-type motion through heterogeneous ensemble of Brownian particles},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {52},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {9},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/aafe90},
  pages     = {27},
  year      = {2019},
  abstract  = {Complex systems are known to display anomalous diffusion, whose signature is a space/time scaling x similar to t(delta) with delta not equal 1/2 in the probability density function (PDF). Anomalous diffusion can emerge jointly with both Gaussian, e.g. fractional Brownian motion, and power-law decaying distributions, e.g. Levy Flights or Levy Walks (LWs). Levy flights get anomalous scaling, but, being jumps of any size allowed even at short times, have infinite position variance, infinite energy and discontinuous paths. LWs, which are based on random trapping events, overcome these limitations: they resemble a Levy-type power-law distribution that is truncated in the large displacement range and have finite moments, finite energy and, even with discontinuous velocity, they are continuous. However, LWs do not take into account the role of strong heterogeneity in many complex systems, such as biological transport in the crowded cell environment. In this work we propose and discuss a model describing a heterogeneous ensemble of Brownian particles (HEBP). Velocity of each single particle obeys a standard underdamped Langevin equation for the velocity, with linear friction term and additive Gaussian noise. Each particle is characterized by its own relaxation time and velocity diffusivity. We show that, for proper distributions of relaxation time and velocity diffusivity, the HEBP resembles some LW statistical features, in particular power-law decaying PDF, long-range correlations and anomalous diffusion, at the same time keeping finite position moments and finite energy. The main differences between the HEBP model and two different LWs are investigated, finding that, even when both velocity and position PDFs are similar, they differ in four main aspects: (i) LWs are biscaling, while HEBP is monoscaling; (ii) a transition from anomalous (delta = 1/2) to normal (delta = 1/2) diffusion in the long-time regime is seen in the HEBP and not in LWs; (iii) the power-law index of the position PDF and the space/time diffusion scaling are independent in the HEBP, while they both depend on the scaling of the interevent time PDF in LWs; (iv) at variance with LWs, our HEBP model obeys a fluctuation-dissipation theorem.},
  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{GuggenbergerChechkinMetzler2021,
  author    = {Guggenberger, Tobias and Chechkin, Aleksei V. and Metzler, Ralf},
  title     = {Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions},
  series = {Journal of physics : A, Mathematical and theoretical},
  volume    = {54},
  journal   = {Journal of physics : A, Mathematical and theoretical},
  number    = {29},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {1751-8113},
  doi       = {10.1088/1751-8121/ac019b},
  pages     = {17},
  year      = {2021},
  abstract  = {We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.},
  language  = {en}
}