TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - A molecular relay race: sequential first-passage events to the terminal reaction centre in a cascade of diffusion controlled processes JF - New Journal of Physics (NJP) N2 - We consider a sequential cascade of molecular first-reaction events towards a terminal reaction centre in which each reaction step is controlled by diffusive motion of the particles. The model studied here represents a typical reaction setting encountered in diverse molecular biology systems, in which, e.g. a signal transduction proceeds via a series of consecutive 'messengers': the first messenger has to find its respective immobile target site triggering a launch of the second messenger, the second messenger seeks its own target site and provokes a launch of the third messenger and so on, resembling a relay race in human competitions. For such a molecular relay race taking place in infinite one-, two- and three-dimensional systems, we find exact expressions for the probability density function of the time instant of the terminal reaction event, conditioned on preceding successful reaction events on an ordered array of target sites. The obtained expressions pertain to the most general conditions: number of intermediate stages and the corresponding diffusion coefficients, the sizes of the target sites, the distances between them, as well as their reactivities are arbitrary. KW - diffusion KW - reaction cascade KW - first passage time Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/ac1e42 SN - 1367-2630 VL - 23 PB - IOP - Institute of Physics Publishing CY - Bristol ER - TY - JOUR A1 - Guggenberger, Tobias A1 - Chechkin, Aleksei A1 - Metzler, Ralf T1 - Absence of stationary states and non-Boltzmann distributions of fractional Brownian motion in shallow external potentials JF - New journal of physics : the open-access journal for physics N2 - We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials. KW - diffusion KW - Boltzmann distribution KW - fractional Brownian motion Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac7b3c SN - 1367-2630 VL - 24 IS - 7 PB - Dt. Physikalische Ges. CY - [Bad Honnef] ER - TY - JOUR A1 - Padash, Amin A1 - Sandev, Trifce A1 - Kantz, Holger A1 - Metzler, Ralf A1 - Chechkin, Aleksei T1 - Asymmetric Levy flights are more efficient in random search JF - Fractal and fractional N2 - We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta. We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times. We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient. When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity. KW - asymmetric Levy flights KW - first-arrival density KW - search efficiency Y1 - 2022 U6 - https://doi.org/10.3390/fractalfract6050260 SN - 2504-3110 VL - 6 IS - 5 PB - MDPI CY - Basel ER - TY - JOUR A1 - Seckler, Henrik A1 - Metzler, Ralf T1 - Bayesian deep learning for error estimation in the analysis of anomalous diffusion JF - Nature Communnications N2 - Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusionmodel and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a wellcalibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output. KW - random-walk KW - models Y1 - 2022 U6 - https://doi.org/10.1038/s41467-022-34305-6 SN - 2041-1723 VL - 13 PB - Nature Publishing Group UK CY - London ER - TY - JOUR A1 - Seckler, Henrik A1 - Metzler, Ralf T1 - Bayesian deep learning for error estimation in the analysis of anomalous diffusion JF - Nature Communications N2 - Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusion model and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a well-calibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output.
Diffusive motions in complex environments such as living biological cells or soft matter systems can be analyzed with single-particle-tracking approaches, where accuracy of output may vary. The authors involve a machine-learning technique for decoding anomalous-diffusion data and provide an uncertainty estimate together with predicted output. Y1 - 2022 U6 - https://doi.org/10.1038/s41467-022-34305-6 SN - 2041-1723 VL - 13 IS - 1 PB - Nature portfolio CY - Berlin ER - TY - JOUR A1 - Granado, Felipe Le Vot A1 - Abad, Enrique A1 - Metzler, Ralf A1 - Yuste, Santos B. T1 - Continuous time random walk in a velocity field BT - role of domain growth, Galilei-invariant advection-diffusion, and kinetics of particle mixing JF - New Journal of Physics N2 - We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed. KW - diffusion KW - expanding medium KW - continuous time random walk Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/ab9ae2 SN - 1367-2630 VL - 22 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Distribution of first-reaction times with target regions on boundaries of shell-like domains JF - New Journal of Physics (NJP) N2 - We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes. KW - diffusion KW - first-passage time KW - first-reaction time KW - shell-like geometries KW - approximate methods Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/ac4282 SN - 1367-2630 VL - 2021 SP - 1 EP - 23 PB - IOP Publishing CY - London ET - 23 ER - TY - JOUR A1 - Sarabadani, Jalal A1 - Metzler, Ralf A1 - Ala-Nissila, Tapio T1 - Driven polymer translocation into a channel: Isoflux tension propagation theory and Langevin dynamics simulations JF - Physical Review Research N2 - Isoflux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f βc for the translocation time, where β = −1 is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics and the force exponent β = −1 regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant and the resulting entropic force causes the force exponent to drop below unity. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.033003 SN - 2643-1564 VL - 4 SP - 033003-1 EP - 033003-14 PB - American Physical Society CY - College Park, Maryland, USA ET - 3 ER - TY - JOUR A1 - Mejia-Monasterio, Carlos A1 - Metzler, Ralf A1 - Vollmer, Jürgen T1 - Editorial: anomalous transport BT - applications, mathematical perspectives, and big data JF - Frontiers in Physics KW - anomalous (or non-Fickian) diffusion KW - anomalous heat conduction KW - stochastic dynamics KW - molecular overcrowding KW - dynamical systems Y1 - 2020 U6 - https://doi.org/10.3389/fphy.2020.622417 SN - 2296-424X VL - 8 PB - Frontiers Media CY - Lausanne ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Sposini, Vittoria A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Seno, Flavio T1 - Exact distributions of the maximum and range of random diffusivity processes JF - New Journal of Physics N2 - We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process. KW - random diffusivity KW - extremal values KW - maximum and range KW - diffusion KW - Brownian motion Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/abd313 SN - 1367-2630 VL - 23 PB - Dt. Physikalische Ges. CY - Bad Honnef ER -