TY - JOUR A1 - Abdolvahab, Rouhollah Haji A1 - Metzler, Ralf A1 - Ejtehadi, Mohammad Reza T1 - First passage time distribution of chaperone driven polymer translocation through a nanopore homopolymer and heteropolymer cases JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - Combining the advection-diffusion equation approach with Monte Carlo simulations we study chaperone driven polymer translocation of a stiff polymer through a nanopore. We demonstrate that the probability density function of first passage times across the pore depends solely on the Peclet number, a dimensionless parameter comparing drift strength and diffusivity. Moreover it is shown that the characteristic exponent in the power-law dependence of the translocation time on the chain length, a function of the chaperone-polymer binding energy, the chaperone concentration, and the chain length, is also effectively determined by the Peclet number. We investigate the effect of the chaperone size on the translocation process. In particular, for large chaperone size, the translocation progress and the mean waiting time as function of the reaction coordinate exhibit pronounced sawtooth-shapes. The effects of a heterogeneous polymer sequence on the translocation dynamics is studied in terms of the translocation velocity, the probability distribution for the translocation progress, and the monomer waiting times. (C) 2011 American Institute of Physics. Y1 - 2011 U6 - https://doi.org/10.1063/1.3669427 SN - 0021-9606 VL - 135 IS - 24 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Adamcik, Jozef A1 - Jeon, Jae-Hyung A1 - Karczewski, Konrad J. A1 - Metzler, Ralf A1 - Dietler, Giovanni T1 - Quantifying supercoiling-induced denaturation bubbles in DNA JF - Soft matter N2 - In both eukaryotic and prokaryotic DNA sequences of 30-100 base-pairs rich in AT base-pairs have been identified at which the double helix preferentially unwinds. Such DNA unwinding elements are commonly associated with origins for DNA replication and transcription, and with chromosomal matrix attachment regions. Here we present a quantitative study of local DNA unwinding based on extensive single DNA plasmid imaging. We demonstrate that long-lived single-stranded denaturation bubbles exist in negatively supercoiled DNA, at the expense of partial twist release. Remarkably, we observe a linear relation between the degree of supercoiling and the bubble size, in excellent agreement with statistical modelling. Furthermore, we obtain the full distribution of bubble sizes and the opening probabilities at varying salt and temperature conditions. The results presented herein underline the important role of denaturation bubbles in negatively supercoiled DNA for biological processes such as transcription and replication initiation in vivo. Y1 - 2012 U6 - https://doi.org/10.1039/c2sm26089a SN - 1744-683X VL - 8 IS - 33 SP - 8651 EP - 8658 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Akimoto, Takuma A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity, rejuvenation, enhancement, and slow relaxation of diffusion in biased continuous-time random walks JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022105 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Awad, Emad A1 - Metzler, Ralf T1 - Crossover dynamics from superdiffusion to subdiffusion BT - models and solutions JF - Fractional calculus and applied analysis : an international journal for theory and applications N2 - 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. KW - Cattaneo equation KW - telegrapher's equation KW - crossover dynamics KW - fractional dynamic equations KW - anomalous diffusion KW - superdiffusion and KW - subdiffusion KW - Fox H-functions Y1 - 2020 U6 - https://doi.org/10.1515/fca-2020-0003 SN - 1311-0454 SN - 1314-2224 VL - 23 IS - 1 SP - 55 EP - 102 PB - De Gruyter CY - Berlin ; Boston 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 - TY - JOUR A1 - Aydiner, Ekrem A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Wealth distribution, Pareto law, and stretched exponential decay of money BT - Computer simulations analysis of agent-based models JF - Physica : europhysics journal ; A, Statistical mechanics and its applications N2 - We study by Monte Carlo simulations a kinetic exchange trading model for both fixed and distributed saving propensities of the agents and rationalize the person and wealth distributions. We show that the newly introduced wealth distribution – that may be more amenable in certain situations – features a different power-law exponent, particularly for distributed saving propensities of the agents. For open agent-based systems, we analyze the person and wealth distributions and find that the presence of trap agents alters their amplitude, leaving however the scaling exponents nearly unaffected. For an open system, we show that the total wealth – for different trap agent densities and saving propensities of the agents – decreases in time according to the classical Kohlrausch–Williams–Watts stretched exponential law. Interestingly, this decay does not depend on the trap agent density, but rather on saving propensities. The system relaxation for fixed and distributed saving schemes are found to be different. KW - Econophysics KW - Wealth and income distribution KW - Pareto law KW - Scaling exponents Y1 - 2017 U6 - https://doi.org/10.1016/j.physa.2017.08.017 SN - 0378-4371 SN - 1873-2119 VL - 490 SP - 278 EP - 288 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Aydiner, Ekrem A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Money distribution in agent-based models with position-exchange dynamics BT - the Pareto paradigm revisited JF - The European physical journal : B, Condensed matter and complex systems N2 - Wealth and income distributions are known to feature country-specific Pareto exponents for their long power-law tails. To propose a rationale for this, we introduce an agent-based dynamic model and use Monte Carlo simulations to unveil the wealth distributions in closed and open economical systems. The standard money-exchange scenario is supplemented with the position-exchange agent dynamics that vitally affects the Pareto law. Specifically, in closed systems with position-exchange dynamics the power law changes to an exponential shape, while for open systems with traps the Pareto law remains valid. KW - Statistical and Nonlinear Physics Y1 - 2019 U6 - https://doi.org/10.1140/epjb/e2019-90674-0 SN - 1434-6028 SN - 1434-6036 VL - 92 IS - 5 PB - Springer CY - New York ER - TY - JOUR A1 - Barkai, Eli A1 - Garini, Yuval A1 - Metzler, Ralf T1 - Strange Kinetics of single molecules in living cells JF - Physics today Y1 - 2012 SN - 0031-9228 VL - 65 IS - 8 SP - 29 EP - 35 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Bauer, Maximilian A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Diffusion of finite-size particles in two-dimensional channels with random wall configurations JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - 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 Y1 - 2014 U6 - https://doi.org/10.1039/c3cp55160a SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 13 SP - 6118 EP - 6128 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Bauer, Maximilian A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - Diffusion of finite-size particles in two-dimensional channels with random wall configurations JF - Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies N2 - 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]. KW - anomalous diffusion KW - fractional dynamics KW - transport KW - nonergodicity KW - coefficient Y1 - 2014 U6 - https://doi.org/10.1039/C3CP55160A SN - 1463-9084 SN - 1463-9076 VL - 16 IS - 13 SP - 6118 EP - 6128 PB - RSC Publications CY - Cambridge ER - TY - JOUR A1 - Bauer, Maximilian A1 - Metzler, Ralf T1 - In vivo facilitated diffusion model JF - PLoS one N2 - Under dilute in vitro conditions transcription factors rapidly locate their target sequence on DNA by using the facilitated diffusion mechanism. However, whether this strategy of alternating between three-dimensional bulk diffusion and one-dimensional sliding along the DNA contour is still beneficial in the crowded interior of cells is highly disputed. Here we use a simple model for the bacterial genome inside the cell and present a semi-analytical model for the in vivo target search of transcription factors within the facilitated diffusion framework. Without having to resort to extensive simulations we determine the mean search time of a lac repressor in a living E. coli cell by including parameters deduced from experimental measurements. The results agree very well with experimental findings, and thus the facilitated diffusion picture emerges as a quantitative approach to gene regulation in living bacteria cells. Furthermore we see that the search time is not very sensitive to the parameters characterizing the DNA configuration and that the cell seems to operate very close to optimal conditions for target localization. Local searches as implied by the colocalization mechanism are only found to mildly accelerate the mean search time within our model. Y1 - 2013 U6 - https://doi.org/10.1371/journal.pone.0053956 SN - 1932-6203 VL - 8 IS - 1 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Bauer, Maximilian A1 - Metzler, Ralf T1 - Generalized facilitated diffusion model for DNA-binding proteins with search and recognition states JF - Biophysical journal N2 - Transcription factors (TFs) such as the lac repressor find their target sequence on DNA at remarkably high rates. In the established Berg-von Hippel model for this search process, the TF alternates between three-dimensional diffusion in the bulk solution and one-dimensional sliding along the DNA chain. To overcome the so-called speed-stability paradox, in similar models the TF was considered as being present in two conformations (search state and recognition state) between which it switches stochastically. Combining both the facilitated diffusion model and alternating states, we obtain a generalized model. We explicitly treat bulk excursions for rodlike chains arranged in parallel and consider a simplified model for coiled DNA. Compared to previously considered facilitated diffusion models, corresponding to limiting cases of our generalized model, we surprisingly find a reduced target search rate. Moreover, at optimal conditions there is no longer an equipartition between the time spent by the protein on and off the DNA chain. Y1 - 2012 U6 - https://doi.org/10.1016/j.bpj.2012.04.008 SN - 0006-3495 VL - 102 IS - 10 SP - 2321 EP - 2330 PB - Cell Press CY - Cambridge ER - TY - JOUR A1 - Bauer, Maximilian A1 - Rasmussen, Emil S. A1 - Lomholt, Michael A. A1 - Metzler, Ralf T1 - Real sequence effects on the search dynamics of transcription factors on DNA JF - Scientific reports N2 - Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs. Y1 - 2015 U6 - https://doi.org/10.1038/srep10072 SN - 2045-2322 VL - 5 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Blavatska, Viktoria A1 - Metzler, Ralf T1 - Conformational properties of complex polymers: rosette versus star-like structures JF - Journal of physics : A, Mathematical and theoretical N2 - Multiple loop formation in polymer macromolecules is an important feature of the chromatin organization and DNA compactification in the nuclei. We analyse the size and shape characteristics of complex polymer structures, containing in general f(1) loops (petals) and f(2) linear chains (branches). Within the frames of continuous model of Gaussian macromolecule, we apply the path integration method and obtain the estimates for gyration radius R-g and asphericity (A) over cap of typical conformation as functions of parameters f(1), f(2). In particular, our results qualitatively reveal the extent of anisotropy of star-like topologies as compared to the rosette structures of the same total molecular weight. KW - polymers KW - path integration KW - conformational properties Y1 - 2015 U6 - https://doi.org/10.1088/1751-8113/48/13/135001 SN - 1751-8113 SN - 1751-8121 VL - 48 IS - 13 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bodrova, Anna A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Quantifying non-ergodic dynamics of force-free granular gases JF - Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies N2 - Brownianmotion is ergodic in the Boltzmann–Khinchin sense that long time averages of physical observables such as the mean squared displacement provide the same information as the corresponding ensemble average, even at out-of-equilibrium conditions. This property is the fundamental prerequisite for single particle tracking and its analysis in simple liquids. We study analytically and by event-driven molecular dynamics simulations the dynamics of force-free cooling granular gases and reveal a violation of ergodicity in this Boltzmann– Khinchin sense as well as distinct ageing of the system. Such granular gases comprise materials such as dilute gases of stones, sand, various types of powders, or large molecules, and their mixtures are ubiquitous in Nature and technology, in particular in Space. We treat—depending on the physical-chemical properties of the inter-particle interaction upon their pair collisions—both a constant and a velocity-dependent (viscoelastic) restitution coefficient e. Moreover we compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behaviour of the ensemble mean squared displacement (MSD) and the velocity correlations in the limit of weak dissipation. Qualitatively, the reported non-ergodic behaviour is generic for granular gases with any realistic dependence of e on the impact velocity of particles. Y1 - 2015 U6 - https://doi.org/10.1039/C5CP02824H SN - 1463-9084 IS - 17 SP - 21791 EP - 21798 ER - TY - JOUR A1 - Bodrova, Anna A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Quantifying non-ergodic dynamics of force-free granular gases JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Brownian motion is ergodic in the Boltzmann-Khinchin sense that long time averages of physical observables such as the mean squared displacement provide the same information as the corresponding ensemble average, even at out-of-equilibrium conditions. This property is the fundamental prerequisite for single particle tracking and its analysis in simple liquids. We study analytically and by event-driven molecular dynamics simulations the dynamics of force-free cooling granular gases and reveal a violation of ergodicity in this Boltzmann-Khinchin sense as well as distinct ageing of the system. Such granular gases comprise materials such as dilute gases of stones, sand, various types of powders, or large molecules, and their mixtures are ubiquitous in Nature and technology, in particular in Space. We treat-depending on the physical-chemical properties of the inter-particle interaction upon their pair collisions-both a constant and a velocity-dependent (viscoelastic) restitution coefficient epsilon. Moreover we compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behaviour of the ensemble mean squared displacement (MSD) and the velocity correlations in the limit of weak dissipation. Qualitatively, the reported non-ergodic behaviour is generic for granular gases with any realistic dependence of epsilon on the impact velocity of particles. Y1 - 2015 U6 - https://doi.org/10.1039/c5cp02824h SN - 1463-9076 SN - 1463-9084 VL - 17 IS - 34 SP - 21791 EP - 21798 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf 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 D(t) similar or equal to 1/t. 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 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion BT - (non-)existence of the overdamped limit in anomalous diffusion JF - Scientific reports N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. Y1 - 2016 U6 - https://doi.org/10.1038/srep30520 SN - 2045-2322 VL - 6 PB - Nature Publishing Group CY - London ER - TY - JOUR A1 - Bodrova, Anna S. A1 - Chechkin, Aleksei V. A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion JF - Scientific reports N2 - It is quite generally assumed that the overdamped Langevin equation provides a quantitative description of the dynamics of a classical Brownian particle in the long time limit. We establish and investigate a paradigm anomalous diffusion process governed by an underdamped Langevin equation with an explicit time dependence of the system temperature and thus the diffusion and damping coefficients. We show that for this underdamped scaled Brownian motion (UDSBM) the overdamped limit fails to describe the long time behaviour of the system and may practically even not exist at all for a certain range of the parameter values. Thus persistent inertial effects play a non-negligible role even at significantly long times. From this study a general questions on the applicability of the overdamped limit to describe the long time motion of an anomalously diffusing particle arises, with profound consequences for the relevance of overdamped anomalous diffusion models. We elucidate our results in view of analytical and simulations results for the anomalous diffusion of particles in free cooling granular gases. Y1 - 2016 U6 - https://doi.org/10.1038/srep30520 SN - 2045-2322 VL - 6 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Caetano, Daniel L. Z. A1 - Carvalho, Sidney Jurado de A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of multiple polyelectrolytes onto a nanosphere BT - splitting the adsorption-desorption transition boundary JF - Interface : journal of the Royal Society N2 - Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges. KW - nanoparticles KW - polyelectrolytes KW - electrostatics KW - critical adsorption KW - phase-transition boundary Y1 - 2020 U6 - https://doi.org/10.1098/rsif.2020.0199 SN - 1742-5689 SN - 1742-5662 VL - 17 IS - 167 PB - Royal Society CY - London ER - TY - JOUR A1 - Caetano, Daniel L. Z. A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of periodic and random polyampholytes onto charged surfaces JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - How different are the properties of critical adsorption of polyampholytes and polyelectrolytes onto charged surfaces? How important are the details of polyampholyte charge distribution on the onset of critical adsorption transition? What are the scaling relations governing the dependence of critical surface charge density on salt concentration in the surrounding solution? Here, we employ Metropolis Monte Carlo simulations and uncover the scaling relations for critical adsorption for quenched periodic and random charge distributions along the polyampholyte chains. We also evaluate and discuss the dependence of the adsorbed layer width on solution salinity and details of the charge distribution. We contrast our findings to the known results for polyelectrolyte adsorption onto oppositely charged surfaces, in particular, their dependence on electrolyte concentration. Y1 - 2017 U6 - https://doi.org/10.1039/c7cp04040g SN - 1463-9076 SN - 1463-9084 VL - 19 SP - 23397 EP - 23413 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Capała, Karol A1 - Padash, Amin A1 - Chechkin, Aleksei V. A1 - Shokri, Babak A1 - Metzler, Ralf A1 - Dybiec, Bartłomiej T1 - Levy noise-driven escape from arctangent potential wells JF - Chaos : an interdisciplinary journal of nonlinear science N2 - The escape from a potential well is an archetypal problem in the study of stochastic dynamical systems, representing real-world situations from chemical reactions to leaving an established home range in movement ecology. Concurrently, Levy noise is a well-established approach to model systems characterized by statistical outliers and diverging higher order moments, ranging from gene expression control to the movement patterns of animals and humans. Here, we study the problem of Levy noise-driven escape from an almost rectangular, arctangent potential well restricted by two absorbing boundaries, mostly under the action of the Cauchy noise. We unveil analogies of the observed transient dynamics to the general properties of stationary states of Levy processes in single-well potentials. The first-escape dynamics is shown to exhibit exponential tails. We examine the dependence of the escape on the shape parameters, steepness, and height of the arctangent potential. Finally, we explore in detail the behavior of the probability densities of the first-escape time and the last-hitting point. Y1 - 2020 U6 - https://doi.org/10.1063/5.0021795 SN - 1054-1500 SN - 1089-7682 VL - 30 IS - 12 PB - American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Kantz, Holger A1 - Metzler, Ralf T1 - Ageing effects in ultraslow continuous time random walks JF - The European physical journal : B, Condensed matter and complex systems N2 - In ageing systems physical observables explicitly depend on the time span elapsing between the original initiation of the system and the actual start of the recording of the particle motion. We here study the signatures of ageing in the framework of ultraslow continuous time random walk processes with super-heavy tailed waiting time densities. We derive the density for the forward or recurrent waiting time of the motion as function of the ageing time, generalise the Montroll-Weiss equation for this process, and analyse the ageing behaviour of the ensemble and time averaged mean squared displacements. Y1 - 2017 U6 - https://doi.org/10.1140/epjb/e2017-80270-9 SN - 1434-6028 SN - 1434-6036 VL - 90 PB - Springer CY - New York ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Seno, Flavio A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Brownian yet Non-Gaussian Diffusion: From Superstatistics to Subordination of Diffusing Diffusivities JF - Physical review : X, Expanding access N2 - A growing number of biological, soft, and active matter systems are observed to exhibit normal diffusive dynamics with a linear growth of the mean-squared displacement, yet with a non-Gaussian distribution of increments. Based on the Chubinsky-Slater idea of a diffusing diffusivity, we here establish and analyze a minimal model framework of diffusion processes with fluctuating diffusivity. In particular, we demonstrate the equivalence of the diffusing diffusivity process with a superstatistical approach with a distribution of diffusivities, at times shorter than the diffusivity correlation time. At longer times, a crossover to a Gaussian distribution with an effective diffusivity emerges. Specifically, we establish a subordination picture of Brownian but non-Gaussian diffusion processes, which can be used for a wide class of diffusivity fluctuation statistics. Our results are shown to be in excellent agreement with simulations and numerical evaluations. Y1 - 2017 U6 - https://doi.org/10.1103/PhysRevX.7.021002 SN - 2160-3308 VL - 7 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Zaid, I. M. A1 - Lomholt, M. A. A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Bulk-mediated surface diffusion on a cylinder in the fast exchange limit JF - Mathematical modelling of natural phenomena N2 - In various biological systems and small scale technological applications particles transiently bind to a cylindrical surface. Upon unbinding the particles diffuse in the vicinal bulk before rebinding to the surface. Such bulk-mediated excursions give rise to an effective surface translation, for which we here derive and discuss the dynamic equations, including additional surface diffusion. We discuss the time evolution of the number of surface-bound particles, the effective surface mean squared displacement, and the surface propagator. In particular, we observe sub- and superdiffusive regimes. A plateau of the surface mean-squared displacement reflects a stalling of the surface diffusion at longer times. Finally, the corresponding first passage problem for the cylindrical geometry is analysed. KW - Bulk-mediated diffusion KW - anomalous diffusion KW - Levy flights KW - stochastic processes Y1 - 2013 U6 - https://doi.org/10.1051/mmnp/20138208 SN - 0973-5348 VL - 8 IS - 2 SP - 114 EP - 126 PB - EDP Sciences CY - Les Ulis ER - TY - JOUR A1 - Chechkin, Aleksei V. A1 - Zaid, Irwin M. A1 - Lomholt, Michael A. A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Bulk-mediated diffusion on a planar surface full solution JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random-walk approach, we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations, we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover, we study the long-time dynamics for the surface motion. Y1 - 2012 U6 - https://doi.org/10.1103/PhysRevE.86.041101 SN - 1539-3755 VL - 86 IS - 4 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity JF - Soft matter N2 - We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells. KW - anomalous diffusion KW - intracellular-transport KW - adenoassociated virus KW - infection pathway KW - escherichia-coli KW - endosomal escape KW - living cells KW - trafficking KW - cytoplasm KW - models Y1 - 2014 U6 - https://doi.org/10.1039/c3sm52846d SN - 2046-2069 VL - 2014 IS - 10 SP - 1591 EP - 1601 PB - Royal Society of Chemistry ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes JF - New journal of physics : the open-access journal for physics N2 - We demonstrate the non-ergodicity of a simple Markovian stochastic process with space-dependent diffusion coefficient D(x). For power-law forms D(x) similar or equal to vertical bar x vertical bar(alpha), this process yields anomalous diffusion of the form < x(2)(t)> similar or equal to t(2/(2-alpha)). Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time-averaged mean-squared displacement <(delta(2)(Delta))over bar> remains linear in the lag time Delta and thus differs from the corresponding ensemble average < x(2)(t)>. We analyse the non-ergodic behaviour of this process in terms of the time-averaged mean- squared displacement (delta(2)) over bar and its random features, i.e. the statistical distribution of (delta(2)) over bar and the ergodicity breaking parameters. The heterogeneous diffusion model represents an alternative approach to non- ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media. Y1 - 2013 U6 - https://doi.org/10.1088/1367-2630/15/8/083039 SN - 1367-2630 VL - 15 IS - 15 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Ageing and confinement in non-ergodic heterogeneous diffusion processes JF - Journal of physics : A, Mathematical and theoretical N2 - We study the effects of ageing-the time delay between initiation of the physical process at t = 0 and start of observation at some time t(a) > 0-and spatial confinement on the properties of heterogeneous diffusion processes (HDPs) with deterministic power-law space-dependent diffusivities, D(x) = D-0 vertical bar x vertical bar(alpha). From analysis of the ensemble and time averaged mean squared displacements and the ergodicity breaking parameter quantifying the inherent degree of irreproducibility of individual realizations of the HDP we obtain striking similarities to ageing subdiffusive continuous time random walks with scale-free waiting time distributions. We also explore how both processes can be distinguished. For confined HDPs we study the long-time saturation of the ensemble and time averaged particle displacements as well as the magnitude of the inherent scatter of time averaged displacements and contrast the outcomes to the results known for other anomalous diffusion processes under confinement. KW - stochastic processes KW - anomalous diffusion KW - ageing KW - weak ergodicity breaking Y1 - 2014 U6 - https://doi.org/10.1088/1751-8113/47/48/485002 SN - 1751-8113 SN - 1751-8121 VL - 47 IS - 48 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity JF - Soft matter N2 - We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells. Y1 - 2014 U6 - https://doi.org/10.1039/c3sm52846d SN - 1744-683X SN - 1744-6848 VL - 10 IS - 10 SP - 1591 EP - 1601 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes JF - Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies N2 - We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion. Y1 - 2016 U6 - https://doi.org/10.1039/C6CP03101C SN - 1463-9084 SN - 1463-9076 VL - 18 SP - 23840 EP - 23852 PB - RSC Publ. CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability. Y1 - 2013 U6 - https://doi.org/10.1039/c3cp53056f SN - 1463-9076 SN - 1463-9084 VL - 15 IS - 46 SP - 20220 EP - 20235 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes N2 - We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability. Y1 - 2008 UR - http://pubs.rsc.org/en/content/articlehtml/2013/cp/c3cp53056f U6 - https://doi.org/10.1039/C3CP53056F ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity breaking, ageing, and confinement in generalized diffusion processes with position and time dependent diffusivity JF - Journal of statistical mechanics: theory and experiment N2 - We study generalized anomalous diffusion processes whose diffusion coefficient D(x, t) similar to D-0x(alpha)t(beta) depends on both the position x of the test particle and the process time t. This process thus combines the features of scaled Brownian motion and heterogeneous diffusion parent processes. We compute the ensemble and time averaged mean squared displacements of this generalized diffusion process. The scaling exponent of the ensemble averaged mean squared displacement is shown to be the product of the critical exponents of the parent processes, and describes both subdiffusive and superdiffusive systems. We quantify the amplitude fluctuations of the time averaged mean squared displacement as function of the length of the time series and the lag time. In particular, we observe a weak ergodicity breaking of this generalized diffusion process: even in the long time limit the ensemble and time averaged mean squared displacements are strictly disparate. When we start to observe this process some time after its initiation we observe distinct features of ageing. We derive a universal ageing factor for the time averaged mean squared displacement containing all information on the ageing time and the measurement time. External confinement is shown to alter the magnitudes and statistics of the ensemble and time averaged mean squared displacements. KW - diffusion Y1 - 2015 U6 - https://doi.org/10.1088/1742-5468/2015/05/P05010 SN - 1742-5468 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Ergodicity breaking and particle spreading in noisy heterogeneous diffusion processes JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We study noisy heterogeneous diffusion processes with a position dependent diffusivity of the form D(x) similar to D-0 vertical bar x vertical bar (alpha 0) in the presence of annealed and quenched disorder of the environment, corresponding to an effective variation of the exponent a in time and space. In the case of annealed disorder, for which effectively alpha(0) = alpha(0)(t), we show how the long time scaling of the ensemble mean squared displacement (MSD) and the amplitude variation of individual realizations of the time averaged MSD are affected by the disorder strength. For the case of quenched disorder, the long time behavior becomes effectively Brownian after a number of jumps between the domains of a stratified medium. In the latter situation, the averages are taken over both an ensemble of particles and different realizations of the disorder. As physical observables, we analyze in detail the ensemble and time averaged MSDs, the ergodicity breaking parameter, and higher order moments of the time averages. (C) 2015 AIP Publishing LLC. Y1 - 2015 U6 - https://doi.org/10.1063/1.4917077 SN - 0021-9606 SN - 1089-7690 VL - 142 IS - 14 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Nonergodicity, fluctuations, and criticality in heterogeneous diffusion processes JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We study the stochastic behavior of heterogeneous diffusion processes with the power-law dependence D(x) similar to vertical bar x vertical bar(alpha) of the generalized diffusion coefficient encompassing sub- and superdiffusive anomalous diffusion. Based on statistical measures such as the amplitude scatter of the time-averaged mean-squared displacement of individual realizations, the ergodicity breaking and non-Gaussianity parameters, as well as the probability density function P(x, t), we analyze the weakly nonergodic character of the heterogeneous diffusion process and, particularly, the degree of irreproducibility of individual realizations. As we show, the fluctuations between individual realizations increase with growing modulus vertical bar alpha vertical bar of the scaling exponent. The fluctuations appear to diverge when the critical value alpha = 2 is approached, while for even larger alpha the fluctuations decrease, again. At criticality, the power-law behavior of the mean-squared displacement changes to an exponentially fast growth, and the fluctuations of the time-averaged mean-squared displacement do not converge for increasing number of realizations. From a systematic comparison we observe some striking similarities of the heterogeneous diffusion process with the familiar subdiffusive continuous time random walk process with power-law waiting time distribution and diverging characteristic waiting time. Y1 - 2014 U6 - https://doi.org/10.1103/PhysRevE.90.012134 SN - 1539-3755 SN - 1550-2376 VL - 90 IS - 1 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We investigate the ensemble and time averaged mean squared displacements for particle diffusion in a simple model for disordered media by assuming that the local diffusivity is both fluctuating in time and has a deterministic average growth or decay in time. In this study we compare computer simulations of the stochastic Langevin equation for this random diffusion process with analytical results. We explore the regimes of normal Brownian motion as well as anomalous diffusion in the sub- and superdiffusive regimes. We also consider effects of the inertial term on the particle motion. The investigation of the resulting diffusion is performed for unconfined and confined motion. Y1 - 2016 U6 - https://doi.org/10.1039/c6cp03101c SN - 1463-9076 SN - 1463-9084 VL - 18 SP - 23840 EP - 23852 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Nagel, Oliver A1 - Beta, Carsten A1 - Metzler, Ralf T1 - Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - What is the underlying diffusion process governing the spreading dynamics and search strategies employed by amoeboid cells? Based on the statistical analysis of experimental single-cell tracking data of the two-dimensional motion of the Dictyostelium discoideum amoeboid cells, we quantify their diffusive behaviour based on a number of standard and complementary statistical indicators. We compute the ensemble- and time-averaged mean-squared displacements (MSDs) of the diffusing amoebae cells and observe a pronounced spread of short-time diffusion coefficients and anomalous MSD-scaling exponents for individual cells. The distribution functions of the cell displacements, the long-tailed distribution of instantaneous speeds, and the velocity autocorrelations are also computed. In particular, we observe a systematic superdiffusive short-time behaviour for the ensemble- and time-averaged MSDs of the amoeboid cells. Also, a clear anti-correlation of scaling exponents and generalised diffusivity values for different cells is detected. Most significantly, we demonstrate that the distribution function of the cell displacements has a strongly non-Gaussian shape andusing a rescaled spatio-temporal variablethe cell-displacement data collapse onto a universal master curve. The current analysis of single-cell motions can be implemented for quantifying diffusive behaviours in other living-matter systems, in particular, when effects of active transport, non-Gaussian displacements, and heterogeneity of the population are involved in the dynamics. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp04254c SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 35 SP - 23034 EP - 23054 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Safdari, Hadiseh A1 - Metzler, Ralf T1 - Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity BT - striking differences for massive versus massless particles JF - Journal of physics. D, Applied physics N2 - 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. KW - anomalous diffusion KW - scaled Brownian motion KW - stochastic processes KW - nonstationary diffusivity KW - water diffusion in the brain KW - nonergodicity Y1 - 2021 U6 - https://doi.org/10.1088/1361-6463/abdff0 SN - 0022-3727 SN - 1361-6463 VL - 54 IS - 19 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Mardoukhi, Yousof A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Time averages and their statistical variation for the Ornstein-Uhlenbeck process BT - Role of initial particle distributions and relaxation to stationarity JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How ergodic is diffusion under harmonic confinements? How strongly do ensemble- and time-averaged displacements differ for a thermally-agitated particle performing confined motion for different initial conditions? We here study these questions for the generic Ornstein-Uhlenbeck (OU) process and derive the analytical expressions for the second and fourth moment. These quantifiers are particularly relevant for the increasing number of single-particle tracking experiments using optical traps. For a fixed starting position, we discuss the definitions underlying the ensemble averages. We also quantify effects of equilibrium and nonequilibrium initial particle distributions onto the relaxation properties and emerging nonequivalence of the ensemble- and time-averaged displacements (even in the limit of long trajectories). We derive analytical expressions for the ergodicity breaking parameter quantifying the amplitude scatter of individual time-averaged trajectories, both for equilibrium and outof-equilibrium initial particle positions, in the entire range of lag times. Our analytical predictions are in excellent agreement with results of computer simulations of the Langevin equation in a parabolic potential. We also examine the validity of the Einstein relation for the ensemble- and time-averaged moments of the OU-particle. Some physical systems, in which the relaxation and nonergodic features we unveiled may be observable, are discussed. Y1 - 2018 U6 - https://doi.org/10.1103/PhysRevE.98.022134 SN - 2470-0045 SN - 2470-0053 VL - 98 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Thapa, Samudrajit A1 - Wagner, Caroline E. A1 - Metzler, Ralf T1 - Non-Gaussian, non-ergodic, and non-Fickian diffusion of tracers in mucin hydrogels JF - Soft matter N2 - Native mucus is polymer-based soft-matter material of paramount biological importance. How non-Gaussian and non-ergodic is the diffusive spreading of pathogens in mucus? We study the passive, thermally driven motion of micron-sized tracers in hydrogels of mucins, the main polymeric component of mucus. We report the results of the Bayesian analysis for ranking several diffusion models for a set of tracer trajectories [C. E. Wagner et al., Biomacromolecules, 2017, 18, 3654]. The models with "diffusing diffusivity', fractional and standard Brownian motion are used. The likelihood functions and evidences of each model are computed, ranking the significance of each model for individual traces. We find that viscoelastic anomalous diffusion is often most probable, followed by Brownian motion, while the model with a diffusing diffusion coefficient is only realised rarely. Our analysis also clarifies the distribution of time-averaged displacements, correlations of scaling exponents and diffusion coefficients, and the degree of non-Gaussianity of displacements at varying pH levels. Weak ergodicity breaking is also quantified. We conclude that-consistent with the original study-diffusion of tracers in the mucin gels is most non-Gaussian and non-ergodic at low pH that corresponds to the most heterogeneous networks. Using the Bayesian approach with the nested-sampling algorithm, together with the quantitative analysis of multiple statistical measures, we report new insights into possible physical mechanisms of diffusion in mucin gels. Y1 - 2019 U6 - https://doi.org/10.1039/c8sm02096e SN - 1744-683X SN - 1744-6848 VL - 15 IS - 12 SP - 2526 EP - 2551 PB - Royal Society of Chemistry CY - Cambridge 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 - 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 : the open-access journal for 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 - 135 EP - 147 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Vinod, Deepak A1 - Aghion, Erez A1 - Sokolov, Igor M. A1 - Metzler, Ralf T1 - Scaled geometric Brownian motion features sub- or superexponential ensemble-averaged, but linear time-averaged mean-squared displacements JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - Various mathematical Black-Scholes-Merton-like models of option pricing employ the paradigmatic stochastic process of geometric Brownian motion (GBM). The innate property of such models and of real stock-market prices is the roughly exponential growth of prices with time [on average, in crisis-free times]. We here explore the ensemble- and time averages of a multiplicative-noise stochastic process with power-law-like time-dependent volatility, sigma(t) similar to t(alpha), named scaled GBM (SGBM). For SGBM, the mean-squared displacement (MSD) computed for an ensemble of statistically equivalent trajectories can grow faster than exponentially in time, while the time-averaged MSD (TAMSD)-based on a sliding-window averaging along a single trajectory-is always linear at short lag times Delta. The proportionality factor between these the two averages of the time series is Delta/T at short lag times, where T is the trajectory length, similarly to GBM. This discrepancy of the scaling relations and pronounced nonequivalence of the MSD and TAMSD at Delta/T << 1 is a manifestation of weak ergodicity breaking for standard GBM and for SGBM with s (t)-modulation, the main focus of our analysis. The analytical predictions for the MSD and mean TAMSD for SGBM are in quantitative agreement with the results of stochastic computer simulations. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.103.062127 SN - 2470-0045 SN - 2470-0053 VL - 103 IS - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Cherstvy, Andrey G. A1 - Wang, Wei A1 - Metzler, Ralf A1 - Sokolov, Igor M. T1 - Inertia triggers nonergodicity of fractional Brownian motion JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM. Y1 - 2021 U6 - https://doi.org/10.1103/PhysRevE.104.024115 SN - 2470-0045 SN - 2470-0053 VL - 104 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces BT - the nonlinear Poisson–Boltzmann approach JF - New journal of physics : the open-access journal for physics N2 - We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems. KW - polyelectrolyte adsorption KW - electrostatic interactions KW - critical phenomena KW - Debye screening Y1 - 2016 U6 - https://doi.org/10.1088/1367-2630/18/8/083037 SN - 1367-2630 VL - 18 PB - IOP Publ. CY - London ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto charged Janus nanospheres JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Based on extensive Monte Carlo simulations and analytical considerations we study the electrostatically driven adsorption of flexible polyelectrolyte chains onto charged Janus nanospheres. These net-neutral colloids are composed of two equally but oppositely charged hemispheres. The critical binding conditions for polyelectrolyte chains are analysed as function of the radius of the Janus particle and its surface charge density, as well as the salt concentration in the ambient solution. Specifically for the adsorption of finite-length polyelectrolyte chains onto Janus nanoparticles, we demonstrate that the critical adsorption conditions drastically differ when the size of the Janus particle or the screening length of the electrolyte are varied. We compare the scaling laws obtained for the adsorption-desorption threshold to the known results for uniformly charged spherical particles, observing significant disparities. We also contrast the changes to the polyelectrolyte chain conformations close to the surface of the Janus nanoparticles as compared to those for simple spherical particles. Finally, we discuss experimentally relevant physicochemical systems for which our simulations results may become important. In particular, we observe similar trends with polyelectrolyte complexation with oppositely but heterogeneously charged proteins. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp02207f SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 29 SP - 15539 EP - 15550 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces: the nonlinear Poisson-Boltzmann approach JF - NEW JOURNAL OF PHYSICS N2 - We study the adsorption-desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces-are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye-Huckel approximation is often not feasible and the nonlinear Poisson-Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson-Boltzmann equation is smaller than the Debye-Huckel result, such that the required critical surface charge density for polyelectrolyte adsorption sigma(c) increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical-chemical and biophysical systems. KW - polyelectrolyte adsorption KW - electrostatic interactions KW - critical phenomena KW - Debye screening Y1 - 2016 U6 - https://doi.org/10.1088/1367-2630/18/8/083037 SN - 1367-2630 VL - 18 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Inverted critical adsorption of polyelectrolytes in confinement JF - Soft matter N2 - What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer–surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density sc—defining the adsorption–desorption transition—are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of sc for the concave interfaces versus the Debye screening length 1/k and the extent of confinement a for these three interfaces for small ka values. For large ka the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/k. We also rationalize how sc(k) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption—the effect often hard to tackle theoretically—putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids. Y1 - 2015 U6 - https://doi.org/10.1039/C5SM00635J SN - 1744-6848 SN - 1744-683X IS - 11 SP - 4430 EP - 4443 PB - Royal Society of Chemistry CY - London ER - TY - JOUR A1 - de Carvalho, Sidney J. A1 - Metzler, Ralf A1 - Cherstvy, Andrey G. T1 - Inverted critical adsorption of polyelectrolytes in confinement JF - Soft matter N2 - What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer-surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density sigma(c)-defining the adsorption-desorption transition-are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of sigma(c) for the concave interfaces versus the Debye screening length 1/kappa and the extent of confinement a for these three interfaces for small kappa a values. For large kappa a the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/kappa. We also rationalize how sigma(c)(kappa) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids. Y1 - 2015 U6 - https://doi.org/10.1039/c5sm00635j SN - 1744-683X SN - 1744-6848 VL - 11 IS - 22 SP - 4430 EP - 4443 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Doerries, Timo J. A1 - Chechkin, Aleksei A1 - Schumer, Rina A1 - Metzler, Ralf T1 - Rate equations, spatial moments, and concentration profiles for mobile-immobile models with power-law and mixed waiting time distributions JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We present a framework for systems in which diffusion-advection transport of a tracer substance in a mobile zone is interrupted by trapping in an immobile zone. Our model unifies different model approaches based on distributed-order diffusion equations, exciton diffusion rate models, and random-walk models for multirate mobile-immobile mass transport. We study various forms for the trapping time dynamics and their effects on the tracer mass in the mobile zone. Moreover, we find the associated breakthrough curves, the tracer density at a fixed point in space as a function of time, and the mobile and immobile concentration profiles and the respective moments of the transport. Specifically, we derive explicit forms for the anomalous transport dynamics and an asymptotic power-law decay of the mobile mass for a Mittag-Leffler trapping time distribution. In our analysis we point out that even for exponential trapping time densities, transient anomalous transport is observed. Our results have direct applications in geophysical contexts, but also in biological, soft matter, and solid state systems. Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevE.105.014105 SN - 2470-0045 SN - 2470-0053 VL - 105 IS - 1 PB - The American Institute of Physics CY - Woodbury, NY ER - TY - JOUR A1 - Doerries, Timo J. A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Apparent anomalous diffusion and non-Gaussian distributions in a simple mobile-immobile transport model with Poissonian switching JF - Interface : journal of the Royal Society N2 - We analyse mobile-immobile transport of particles that switch between the mobile and immobile phases with finite rates. Despite this seemingly simple assumption of Poissonian switching, we unveil a rich transport dynamics including significant transient anomalous diffusion and non-Gaussian displacement distributions. Our discussion is based on experimental parameters for tau proteins in neuronal cells, but the results obtained here are expected to be of relevance for a broad class of processes in complex systems. Specifically, we obtain that, when the mean binding time is significantly longer than the mean mobile time, transient anomalous diffusion is observed at short and intermediate time scales, with a strong dependence on the fraction of initially mobile and immobile particles. We unveil a Laplace distribution of particle displacements at relevant intermediate time scales. For any initial fraction of mobile particles, the respective mean squared displacement (MSD) displays a plateau. Moreover, we demonstrate a short-time cubic time dependence of the MSD for immobile tracers when initially all particles are immobile. KW - diffusion KW - mobile-immobile model KW - tau proteins Y1 - 2022 U6 - https://doi.org/10.1098/rsif.2022.0233 SN - 1742-5689 SN - 1742-5662 VL - 19 IS - 192 PB - Royal Society CY - London ER - TY - JOUR A1 - Dybiec, Bartlomiej A1 - Capala, Karol A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Conservative random walks in confining potentials JF - Journal of physics : A, Mathematical and theoretical N2 - Levy walks are continuous time random walks with spatio-temporal coupling of jump lengths and waiting times, often used to model superdiffusive spreading processes such as animals searching for food, tracer motion in weakly chaotic systems, or even the dynamics in quantum systems such as cold atoms. In the simplest version Levy walks move with a finite speed. Here, we present an extension of the Levy walk scenario for the case when external force fields influence the motion. The resulting motion is a combination of the response to the deterministic force acting on the particle, changing its velocity according to the principle of total energy conservation, and random velocity reversals governed by the distribution of waiting times. For the fact that the motion stays conservative, that is, on a constant energy surface, our scenario is fundamentally different from thermal motion in the same external potentials. In particular, we present results for the velocity and position distributions for single well potentials of different steepness. The observed dynamics with its continuous velocity changes enriches the theory of Levy walk processes and will be of use in a variety of systems, for which the particles are externally confined. KW - Levy walk KW - conservative random walks KW - Levy flight Y1 - 2018 U6 - https://doi.org/10.1088/1751-8121/aaefc2 SN - 1751-8113 SN - 1751-8121 VL - 52 IS - 1 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf T1 - Anomalous statistics of random relaxations in random environments JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - We comprehensively analyze the emergence of anomalous statistics in the context of the random relaxation ( RARE) model [Eliazar and Metzler, J. Chem. Phys. 137, 234106 ( 2012)], a recently introduced versatile model of random relaxations in random environments. The RARE model considers excitations scattered randomly across a metric space around a reaction center. The excitations react randomly with the center, the reaction rates depending on the excitations' distances from this center. Relaxation occurs upon the first reaction between an excitation and the center. Addressing both the relaxation time and the relaxation range, we explore when these random variables display anomalous statistics, namely, heavy tails at zero and at infinity that manifest, respectively, exceptionally high occurrence probabilities of very small and very large outliers. A cohesive set of closed-form analytic results is established, determining precisely when such anomalous statistics emerge. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevE.87.022141 SN - 1539-3755 SN - 1550-2376 VL - 87 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf T1 - The RARE model a generalized approach to random relaxation processes in disordered systems JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - This paper introduces and analyses a general statistical model, termed the RAndom RElaxations (RARE) model, of random relaxation processes in disordered systems. The model considers excitations that are randomly scattered around a reaction center in a general embedding space. The model's input quantities are the spatial scattering statistics of the excitations around the reaction center, and the chemical reaction rates between the excitations and the reaction center as a function of their mutual distance. The framework of the RARE model is versatile and a detailed stochastic analysis of the random relaxation processes is established. Analytic results regarding the duration and the range of the random relaxation processes, as well as the model's thermodynamic limit, are obtained in closed form. In particular, the case of power-law inputs, which turn out to yield stretched exponential relaxation patterns and asymptotically Paretian relaxation ranges, is addressed in detail. KW - chemical relaxation KW - Pareto analysis KW - reaction kinetics theory KW - reaction rate constants KW - stochastic processes Y1 - 2012 U6 - https://doi.org/10.1063/1.4770266 SN - 0021-9606 SN - 1089-7690 VL - 137 IS - 23 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf A1 - Reuveni, Shlomi T1 - Poisson-process limit laws yield Gumbel max-min and min-max JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - “A chain is only as strong as its weakest link” says the proverb. But what about a collection of statistically identical chains: How long till all chains fail? The answer to this question is given by the max-min of a matrix whose (i,j)entry is the failure time of link j of chain i: take the minimum of each row, and then the maximum of the rows' minima. The corresponding min-max is obtained by taking the maximum of each column, and then the minimum of the columns' maxima. The min-max applies to the storage of critical data. Indeed, consider multiple backup copies of a set of critical data items, and consider the (i,j) matrix entry to be the time at which item j on copy i is lost; then, the min-max is the time at which the first critical data item is lost. In this paper we address random matrices whose entries are independent and identically distributed random variables. We establish Poisson-process limit laws for the row's minima and for the columns' maxima. Then, we further establish Gumbel limit laws for the max-min and for the min-max. The limit laws hold whenever the entries' distribution has a density, and yield highly applicable approximation tools and design tools for the max-min and min-max of large random matrices. A brief of the results presented herein is given in: Gumbel central limit theorem for max-min and min-max Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.022129 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Eliazar, Iddo A1 - Metzler, Ralf A1 - Reuveni, Shlomi T1 - Gumbel central limit theorem for max-min and min-max JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - The max-min and min-max of matrices arise prevalently in science and engineering. However, in many real-world situations the computation of the max-min and min-max is challenging as matrices are large and full information about their entries is lacking. Here we take a statistical-physics approach and establish limit laws—akin to the central limit theorem—for the max-min and min-max of large random matrices. The limit laws intertwine random-matrix theory and extreme-value theory, couple the matrix dimensions geometrically, and assert that Gumbel statistics emerge irrespective of the matrix entries' distribution. Due to their generality and universality, as well as their practicality, these results are expected to have a host of applications in the physical sciences and beyond. Y1 - 2019 U6 - https://doi.org/10.1103/PhysRevE.100.020104 SN - 2470-0045 SN - 2470-0053 VL - 100 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Emanuel, Marc D. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Gompper, Gerhard T1 - Buckling transitions and soft-phase invasion of two-component icosahedral shells JF - Physical review / publ. by The American Physical Society. E, Statistical, nonlinear, and soft matter physics N2 - What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevE.102.062104 SN - 2470-0045 SN - 2470-0053 SN - 2470-0061 SN - 1538-4519 VL - 102 IS - 6 PB - Woodbury CY - New York ER - TY - JOUR A1 - Estrada, Ernesto A1 - Delvenne, Jean-Charles A1 - Hatano, Naomichi A1 - Mateos, Jose L. A1 - Metzler, Ralf A1 - Riascos, Alejandro P. A1 - Schaub, Michael T. T1 - Random multi-hopper model BT - super-fast random walks on graphs JF - Journal of Complex Networks N2 - We develop a mathematical model considering a random walker with long-range hops on arbitrary graphs. The random multi-hopper can jump to any node of the graph from an initial position, with a probability that decays as a function of the shortest-path distance between the two nodes in the graph. We consider here two decaying functions in the form of Laplace and Mellin transforms of the shortest-path distances. We prove that when the parameters of these transforms approach zero asymptotically, the hitting time in the multi-hopper approaches the minimum possible value for a normal random walker. We show by computational experiments that the multi-hopper explores a graph with clusters or skewed degree distributions more efficiently than a normal random walker. We provide computational evidences of the advantages of the random multi-hopper model with respect to the normal random walk by studying deterministic, random and real-world networks. Y1 - 2018 U6 - https://doi.org/10.1093/comnet/cnx043 SN - 2051-1310 SN - 2051-1329 VL - 6 IS - 3 SP - 382 EP - 403 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Fernandez, Amanda Diez A1 - Charchar, Patrick A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Finnis, Michael W. T1 - The diffusion of doxorubicin drug molecules in silica nanoslits is non-Gaussian, intermittent and anticorrelated JF - Physical chemistry, chemical physics N2 - In this study we investigate, using all-atom molecular-dynamics computer simulations, the in-plane diffusion of a doxorubicin drug molecule in a thin film of water confined between two silica surfaces. We find that the molecule diffuses along the channel in the manner of a Gaussian diffusion process, but with parameters that vary according to its varying transversal position. Our analysis identifies that four Gaussians, each describing particle motion in a given transversal region, are needed to adequately describe the data. Each of these processes by itself evolves with time at a rate slower than that associated with classical Brownian motion due to a predominance of anticorrelated displacements. Long adsorption events lead to ageing, a property observed when the diffusion is intermittently hindered for periods of time with an average duration which is theoretically infinite. This study presents a simple system in which many interesting features of anomalous diffusion can be explored. It exposes the complexity of diffusion in nanoconfinement and highlights the need to develop new understanding. Y1 - 2020 U6 - https://doi.org/10.1039/d0cp03849k SN - 1463-9076 SN - 1463-9084 VL - 22 IS - 48 SP - 27955 EP - 27965 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Grebenkov, Denis S. A1 - Metzler, Ralf T1 - Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments JF - NEW JOURNAL OF PHYSICS N2 - A topic of intense current investigation pursues the question of how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment made up of quenched distributions of monodisperse crowders of varying sizes in finite circular two-dimensional domains. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain or at the outer boundary, exhibiting distinctly different features of the observed anomalous diffusion for heterogeneous distributions of crowders. Specifically we reveal an asymmetric spreading of tracers even at moderate crowding. In addition to the mean squared displacement (MSD) and local diffusion exponent we investigate the magnitude and the amplitude scatter of the time averaged MSD of individual tracer trajectories, the non-Gaussianity parameter, and the van Hove correlation function. We also quantify how the average tracer diffusivity varies with the position in the domain with a heterogeneous radial distribution of crowders and examine the behaviour of the survival probability and the dynamics of the tracer survival probability. Inter alia, the systems we investigate are related to the passive transport of lipid molecules and proteins in two-dimensional crowded membranes or the motion in colloidal solutions or emulsions in effectively two-dimensional geometries, as well as inside supercrowded, surface adhered cells. KW - anomalous diffusion KW - crowded fluids KW - stochastic processes Y1 - 2016 U6 - https://doi.org/10.1088/1367-2630/18/1/013027 SN - 1367-2630 VL - 18 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf ED - Metzler, Ralf T1 - Non-universal tracer diffusion in crowded media of non-inert obstacles JF - Physical Chemistry Chemical Physics N2 - We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer–obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer–obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer–crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids. KW - fluorescence correlation spectroscopy KW - single-particle tracking KW - anomalous diffusion KW - living cells KW - physiological consequences KW - langevin equation KW - infection pathway KW - excluded volume KW - brownian-motion KW - random-walks Y1 - 2014 SN - 1463-9076 VL - 3 IS - 17 SP - 1847 EP - 1858 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Deformation propagation in responsive polymer network films JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - We study the elastic deformations in a cross-linked polymer network film triggered by the binding of submicron particles with a sticky surface, mimicking the interactions of viral pathogens with thin films of stimulus-responsive polymeric materials such as hydrogels. From extensive Langevin Dynamics simulations we quantify how far the network deformations propagate depending on the elasticity parameters of the network and the adhesion strength of the particles. We examine the dynamics of the collective area shrinkage of the network and obtain some simple relations for the associated characteristic decay lengths. A detailed analysis elucidates how the elastic energy of the network is distributed between stretching and compression modes in response to the particle binding. We also examine the force-distance curves of the repulsion or attraction interactions for a pair of sticky particles in the polymer network film as a function of the particle-particle separation. The results of this computational study provide new insight into collective phenomena in soft polymer network films and may, in particular, be applied to applications for visual detection of pathogens such as viruses via a macroscopic response of thin films of cross-linked hydrogels. (C) 2014 AIP Publishing LLC. Y1 - 2014 U6 - https://doi.org/10.1063/1.4893056 SN - 0021-9606 SN - 1089-7690 VL - 141 IS - 7 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf T1 - Non-universal tracer diffusion in crowded media of non-inert obstacles JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer-obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer-obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer-crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids. Y1 - 2015 U6 - https://doi.org/10.1039/c4cp03599b SN - 1463-9076 SN - 1463-9084 VL - 17 IS - 3 SP - 1847 EP - 1858 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Petrov, Eugene P. A1 - Metzler, Ralf T1 - Interactions of rod-like particles on responsive elastic sheets JF - Soft matter N2 - What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive–repulsive rod–rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed. Y1 - 2016 U6 - https://doi.org/10.1039/C6SM01522K SN - 1744-6848 SN - 1744-683X PB - RSC CY - London ER - TY - JOUR A1 - Ghosh, Surya K. A1 - Cherstvy, Andrey G. A1 - Petrov, Eugene P. A1 - Metzler, Ralf T1 - Interactions of rod-like particles on responsive elastic sheets JF - Soft matter N2 - What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive-repulsive rod-rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed. Y1 - 2016 U6 - https://doi.org/10.1039/c6sm01522k SN - 1744-683X SN - 1744-6848 VL - 12 SP - 7908 EP - 7919 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Godec, Aljaz A1 - Bauer, Maximilian A1 - Metzler, Ralf T1 - Collective dynamics effect transient subdiffusion of inert tracers in flexible gel networks JF - New journal of physics : the open-access journal for physics N2 - Based on extensive Brownian dynamics simulations we study the thermal motion of a tracer bead in a cross-linked, flexible gel in the limit when the tracer particle size is comparable to or even larger than the equilibrium mesh size of the gel. The analysis of long individual trajectories of the tracer demonstrates the existence of pronounced transient anomalous diffusion. From the time averaged mean squared displacement and the time averaged van Hove correlation functions we elucidate the many-body origin of the non-Brownian tracer bead dynamics. Our results shed new light onto the ongoing debate over the physical origin of steric tracer interactions with structured environments. KW - anomalous diffusion KW - gel network KW - van Hove correlation Y1 - 2014 U6 - https://doi.org/10.1088/1367-2630/16/9/092002 SN - 1367-2630 VL - 16 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Godec, Aljaz A1 - Chechkin, Aleksei V. A1 - Barkai, Eli A1 - Kantz, Holger A1 - Metzler, Ralf T1 - Localisation and universal fluctuations in ultraslow diffusion processes JF - Journal of physics : A, Mathematical and theoretical N2 - We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) < x(2)(t)> similar or equal to log(gamma)t. Comparison of annealed (renewal) continuous time random walks (CTRWs) with logarithmic waiting time distribution psi(tau) similar or equal to 1/(tau log(1+gamma)tau) and Sinai diffusion in quenched random landscapes reveals striking similarities, despite the great differences in their physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time-averaged and ensemble-averaged MSDs. Remarkably, for the CTRW we observe that the fluctuations of time averages become universal, with an exponential suppression of mobile trajectories. We discuss the fundamental connection between the Golosov localization effect and non-ergodicity in the sense of the disparity between ensemble-averaged MSD and time-averaged MSD. KW - Sinai diffusion KW - anomalous diffusion KW - quenched energy landscape Y1 - 2014 U6 - https://doi.org/10.1088/1751-8113/47/49/492002 SN - 1751-8113 SN - 1751-8121 VL - 47 IS - 49 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Finite-Time effects and ultraweak ergodicity breaking in superdiffusive dynamics JF - Physical review letters N2 - We study the ergodic properties of superdiffusive, spatiotemporally coupled Levy walk processes. For trajectories of finite duration, we reveal a distinct scatter of the scaling exponents of the time averaged mean squared displacement (delta x(2)) over bar around the ensemble value 3 - alpha (1 < alpha < 2) ranging from ballistic motion to subdiffusion, in strong contrast to the behavior of subdiffusive processes. In addition we find a significant dependence of the average of (delta x(2)) over bar over an ensemble of trajectories as a function of the finite measurement time. This so-called finite-time amplitude depression and the scatter of the scaling exponent is vital in the quantitative evaluation of superdiffusive processes. Comparing the long time average of the second moment with the ensemble mean squared displacement, these only differ by a constant factor, an ultraweak ergodicity breaking. Y1 - 2013 U6 - https://doi.org/10.1103/PhysRevLett.110.020603 SN - 0031-9007 VL - 110 IS - 2 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Optimization and universality of Brownian search in a basic model of quenched heterogeneous media JF - Physical review : E, Statistical, nonlinear and soft matter physics N2 - The kinetics of a variety of transport-controlled processes can be reduced to the problem of determining the mean time needed to arrive at a given location for the first time, the so-called mean first-passage time ( MFPT) problem. The occurrence of occasional large jumps or intermittent patterns combining various types of motion are known to outperform the standard random walk with respect to the MFPT, by reducing oversampling of space. Here we show that a regular but spatially heterogeneous random walk can significantly and universally enhance the search in any spatial dimension. In a generic minimal model we consider a spherically symmetric system comprising two concentric regions with piecewise constant diffusivity. The MFPT is analyzed under the constraint of conserved average dynamics, that is, the spatially averaged diffusivity is kept constant. Our analytical calculations and extensive numerical simulations demonstrate the existence of an optimal heterogeneity minimizing the MFPT to the target. We prove that the MFPT for a random walk is completely dominated by what we term direct trajectories towards the target and reveal a remarkable universality of the spatially heterogeneous search with respect to target size and system dimensionality. In contrast to intermittent strategies, which are most profitable in low spatial dimensions, the spatially inhomogeneous search performs best in higher dimensions. Discussing our results alongside recent experiments on single-particle tracking in living cells, we argue that the observed spatial heterogeneity may be beneficial for cellular signaling processes. Y1 - 2015 U6 - https://doi.org/10.1103/PhysRevE.91.052134 SN - 1539-3755 SN - 1550-2376 VL - 91 IS - 5 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time JF - Scientific reports N2 - The first passage is a generic concept for quantifying when a random quantity such as the position of a diffusing molecule or the value of a stock crosses a preset threshold (target) for the first time. The last decade saw an enlightening series of new results focusing mostly on the so-called mean and global first passage time (MFPT and GFPT, respectively) of such processes. Here we push the understanding of first passage processes one step further. For a simple heterogeneous system we derive rigorously the complete distribution of first passage times (FPTs). Our results demonstrate that the typical FPT significantly differs from the MFPT, which corresponds to the long time behaviour of the FPT distribution. Conversely, the short time behaviour is shown to correspond to trajectories connecting directly from the initial value to the target. Remarkably, we reveal a previously overlooked third characteristic time scale of the first passage dynamics mirroring brief excursion away from the target. Y1 - 2016 U6 - https://doi.org/10.1038/srep20349 SN - 2045-2322 VL - 6 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Universal Proximity Effect in Target Search Kinetics in the Few-Encounter Limit JF - Physical review : X, Expanding access N2 - When does a diffusing particle reach its target for the first time? This first-passage time (FPT) problem is central to the kinetics of molecular reactions in chemistry and molecular biology. Here, we explain the behavior of smooth FPT densities, for which all moments are finite, and demonstrate universal yet generally non-Poissonian long-time asymptotics for a broad variety of transport processes. While Poisson-like asymptotics arise generically in the presence of an effective repulsion in the immediate vicinity of the target, a time-scale separation between direct and reflected indirect trajectories gives rise to a universal proximity effect: Direct paths, heading more or less straight from the point of release to the target, become typical and focused, with a narrow spread of the corresponding first-passage times. Conversely, statistically dominant indirect paths exploring the entire system tend to be massively dissimilar. The initial distance to the target particularly impacts gene regulatory or competitive stochastic processes, for which few binding events often determine the regulatory outcome. The proximity effect is independent of details of the transport, highlighting the robust character of the FPT features uncovered here. Y1 - 2016 U6 - https://doi.org/10.1103/PhysRevX.6.041037 SN - 2160-3308 VL - 6 PB - American Physical Society CY - College Park ER - TY - JOUR A1 - Godec, Aljaz A1 - Metzler, Ralf T1 - Active transport improves the precision of linear long distance molecular signalling JF - Journal of physics : A, Mathematical and theoretical N2 - Molecular signalling in living cells occurs at low copy numbers and is thereby inherently limited by the noise imposed by thermal diffusion. The precision at which biochemical receptors can count signalling molecules is intimately related to the noise correlation time. In addition to passive thermal diffusion, messenger RNA and vesicle-engulfed signalling molecules can transiently bind to molecular motors and are actively transported across biological cells. Active transport is most beneficial when trafficking occurs over large distances, for instance up to the order of 1 metre in neurons. Here we explain how intermittent active transport allows for faster equilibration upon a change in concentration triggered by biochemical stimuli. Moreover, we show how intermittent active excursions induce qualitative changes in the noise in effectively one-dimensional systems such as dendrites. Thereby they allow for significantly improved signalling precision in the sense of a smaller relative deviation in the concentration read-out by the receptor. On the basis of linear response theory we derive the exact mean field precision limit for counting actively transported molecules. We explain how intermittent active excursions disrupt the recurrence in the molecular motion, thereby facilitating improved signalling accuracy. Our results provide a deeper understanding of how recurrence affects molecular signalling precision in biological cells and novel medical-diagnostic devices. KW - noise in biochemical signalling KW - Brownian motion KW - active transport KW - linear response theory KW - fluctuation-dissipation theorem KW - generalised Langevin equation KW - recurrence Y1 - 2016 U6 - https://doi.org/10.1088/1751-8113/49/36/364001 SN - 1751-8113 SN - 1751-8121 VL - 49 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Godec, Aljaž A1 - Metzler, Ralf T1 - First passage time statistics for two-channel diffusion JF - Journal of physics : A, Mathematical and theoretical N2 - We present rigorous results for the mean first passage time and first passage time statistics for two-channel Markov additive diffusion in a 3-dimensional spherical domain. Inspired by biophysical examples we assume that the particle can only recognise the target in one of the modes, which is shown to effect a non-trivial first passage behaviour. We also address the scenario of intermittent immobilisation. In both cases we prove that despite the perfectly non-recurrent motion of two-channel Markov additive diffusion in 3 dimensions the first passage statistics at long times do not display Poisson-like behaviour if none of the phases has a vanishing diffusion coefficient. This stands in stark contrast to the standard (one-channel) Markov diffusion counterpart. We also discuss the relevance of our results in the context of cellular signalling. KW - first passage time KW - Markov additive processes KW - Fokker-Planck equation KW - random search processes KW - coupled initial boundary value problem KW - cellular signalling KW - asymptotic analysis Y1 - 2017 U6 - https://doi.org/10.1088/1751-8121/aa5204 SN - 1751-8113 SN - 1751-8121 VL - 50 IS - 8 PB - IOP Publ. Ltd. CY - Bristol 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 - Goychuk, Igor A1 - Kharchenko, Vasyl O. A1 - Metzler, Ralf T1 - How Molecular Motors Work in the Crowded Environment of Living Cells: Coexistence and Efficiency of Normal and Anomalous Transport JF - PLoS one N2 - Recent experiments reveal both passive subdiffusion of various nanoparticles and anomalous active transport of such particles by molecular motors in the molecularly crowded environment of living biological cells. Passive and active microrheology reveals that the origin of this anomalous dynamics is due to the viscoelasticity of the intracellular fluid. How do molecular motors perform in such a highly viscous, dissipative environment? Can we explain the observed co-existence of the anomalous transport of relatively large particles of 100 to 500 nm in size by kinesin motors with the normal transport of smaller particles by the same molecular motors? What is the efficiency of molecular motors in the anomalous transport regime? Here we answer these seemingly conflicting questions and consistently explain experimental findings in a generalization of the well-known continuous diffusion model for molecular motors with two conformational states in which viscoelastic effects are included. Y1 - 2014 U6 - https://doi.org/10.1371/journal.pone.0091700 SN - 1932-6203 VL - 9 IS - 3 PB - PLoS CY - San Fransisco ER - TY - JOUR A1 - Goychuk, Igor 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 : a journal of European Chemical Societies 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. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp01234h SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 31 SP - 16524 EP - 16535 PB - Royal Society of Chemistry CY - Cambridge 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 - 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 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains JF - New journal of physics : the open-access journal for physics N2 - We study the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry-characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. A similar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA. We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters. We analyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations. KW - first passage time KW - cylindrical geometry KW - aspect ratio KW - protein search Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa8ed9 SN - 1367-2630 VL - 19 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - From single-particle stochastic kinetics to macroscopic reaction rates BT - fastest first-passage time of N random walkers JF - New Journal of Physics N2 - We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells. KW - diffusion KW - first-passage KW - fastest first-passage time of N walkers Y1 - 2020 U6 - https://doi.org/10.1088/1367-2630/abb1de 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 - Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains JF - New journal of physics N2 - Westudy the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry—characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. Asimilar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA.We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters.Weanalyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations. KW - first passage time KW - cylindrical geometry KW - aspect ratio KW - protein search Y1 - 2017 U6 - https://doi.org/10.1088/1367-2630/aa8ed9 SN - 1367-2630 VL - 19 SP - 1 EP - 11 PB - IOP CY - London ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Strong defocusing of molecular reaction times results from an interplay of geometry and reaction control JF - Communications Chemistry N2 - Textbook concepts of diffusion-versus kinetic-control are well-defined for reaction-kinetics involving macroscopic concentrations of diffusive reactants that are adequately described by rate-constants—the inverse of the mean-first-passage-time to the reaction-event. In contradiction, an open important question is whether the mean-first-passage-time alone is a sufficient measure for biochemical reactions that involve nanomolar reactant concentrations. Here, using a simple yet generic, exactly solvable model we study the effect of diffusion and chemical reaction-limitations on the full reaction-time distribution. We show that it has a complex structure with four distinct regimes delineated by three characteristic time scales spanning a window of several decades. Consequently, the reaction-times are defocused: no unique time-scale characterises the reaction-process, diffusion- and kinetic-control can no longer be disentangled, and it is imperative to know the full reaction-time distribution. We introduce the concepts of geometry- and reaction-control, and also quantify each regime by calculating the corresponding reaction depth. Y1 - 2018 U6 - https://doi.org/10.1038/s42004-018-0096-x SN - 2399-3669 VL - 1 PB - Macmillan Publishers Limited CY - London ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Full distribution of first exit times in the narrow escape problem JF - New Journal of Physics N2 - In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems. KW - narrow escape problem KW - first-passage time distribution KW - mean versus most probable reaction times KW - mixed boundary conditions Y1 - 2019 U6 - https://doi.org/10.1088/1367-2630/ab5de4 SN - 1367-2630 VL - 21 PB - Dt. Physikalische Ges. CY - Bad Honnef ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Towards a full quantitative description of single-molecule reaction kinetics in biological cells JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - The first-passage time (FPT), i.e., the moment when a stochastic process reaches a given threshold value for the first time, is a fundamental mathematical concept with immediate applications. In particular, it quantifies the statistics of instances when biomolecules in a biological cell reach their specific binding sites and trigger cellular regulation. Typically, the first-passage properties are given in terms of mean first-passage times. However, modern experiments now monitor single-molecular binding-processes in living cells and thus provide access to the full statistics of the underlying first-passage events, in particular, inherent cell-to-cell fluctuations. We here present a robust explicit approach for obtaining the distribution of FPTs to a small partially reactive target in cylindrical-annulus domains, which represent typical bacterial and neuronal cell shapes. We investigate various asymptotic behaviours of this FPT distribution and show that it is typically very broad in many biological situations, thus, the mean FPT can differ from the most probable FPT by orders of magnitude. The most probable FPT is shown to strongly depend only on the starting position within the geometry and to be almost independent of the target size and reactivity. These findings demonstrate the dramatic relevance of knowing the full distribution of FPTs and thus open new perspectives for a more reliable description of many intracellular processes initiated by the arrival of one or few biomolecules to a small, spatially localised region inside the cell. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp02043d SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 24 SP - 16393 EP - 16401 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb T1 - Search efficiency in the Adam-Delbruck reduction-of-dimensionality scenario versus direct diffusive search JF - New journal of physics : the open-access journal for physics N2 - The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search. KW - first-passage times KW - Adam-Delbruck scenario KW - dimensional reduction KW - bulk KW - and surface diffusion Y1 - 2022 U6 - https://doi.org/10.1088/1367-2630/ac8824 SN - 1367-2630 VL - 24 IS - 8 PB - IOP Publ. Ltd. CY - Bristol ER - TY - JOUR A1 - Grebenkov, Denis S. A1 - Metzler, Ralf A1 - Oshanin, Gleb A1 - Dagdug, Leonardo A1 - Berezhkovskii, Alexander M. A1 - Skvortsov, Alexei T. T1 - Trapping of diffusing particles by periodic absorbing rings on a cylindrical tube JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr Y1 - 2019 U6 - https://doi.org/10.1063/1.5098390 SN - 0021-9606 SN - 1089-7690 VL - 150 IS - 20 PB - American Institute of Physics CY - Melville 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 - 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 - Guggenberger, Tobias A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions JF - Journal of physics : A, Mathematical and theoretical N2 - 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. KW - anomalous diffusion KW - Boltzmann distribution KW - non-Gaussian distribution Y1 - 2021 U6 - https://doi.org/10.1088/1751-8121/ac019b SN - 1751-8113 SN - 1751-8121 VL - 54 IS - 29 PB - IOP Publ. Ltd. CY - Bristol 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 - 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 - Herrmann, Carl J. J. A1 - Metzler, Ralf A1 - Engbert, 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 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Hou, Ru A1 - Cherstvy, Andrey G. A1 - Metzler, Ralf A1 - Akimoto, Takuma T1 - Biased continuous-time random walks for ordinary and equilibrium cases BT - facilitation of diffusion, ergodicity breaking and ageing JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent alpha are considered for the WTD psi(t) similar to 1/t(1+alpha). We treat continuous-time random walks (CTRWs) with 0 < alpha < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < alpha < 2 and alpha > 2. We demonstrate that in the presence of a drift CTRWs with alpha < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < alpha < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with alpha > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < alpha < 2 and alpha > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with alpha > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent alpha. Y1 - 2018 U6 - https://doi.org/10.1039/c8cp01863d SN - 1463-9076 SN - 1463-9084 VL - 20 IS - 32 SP - 20827 EP - 20848 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Javanainen, Matti A1 - Hammaren, Henrik A1 - Monticelli, Luca A1 - Jeon, Jae-Hyung A1 - Miettinen, Markus S. A1 - Martinez-Seara, Hector A1 - Metzler, Ralf A1 - Vattulainen, Ilpo T1 - Anomalous and normal diffusion of proteins and lipids in crowded lipid membranes JF - Faraday discussions N2 - Lateral diffusion plays a crucial role in numerous processes that take place in cell membranes, yet it is quite poorly understood in native membranes characterized by, e.g., domain formation and large concentration of proteins. In this article, we use atomistic and coarse-grained simulations to consider how packing of membranes and crowding with proteins affect the lateral dynamics of lipids and membrane proteins. We find that both packing and protein crowding have a profound effect on lateral diffusion, slowing it down. Anomalous diffusion is observed to be an inherent property in both protein-free and protein-rich membranes, and the time scales of anomalous diffusion and the exponent associated with anomalous diffusion are found to strongly depend on packing and crowding. Crowding with proteins also has a striking effect on the decay rate of dynamical correlations associated with lateral single-particle motion, as the transition from anomalous to normal diffusion is found to take place at macroscopic time scales: while in protein-poor conditions normal diffusion is typically observed in hundreds of nanoseconds, in protein-rich conditions the onset of normal diffusion is tens of microseconds, and in the most crowded systems as large as milliseconds. The computational challenge which results from these time scales is not easy to deal with, not even in coarse-grained simulations. We also briefly discuss the physical limits of protein motion. Our results suggest that protein concentration is anything but constant in the plane of cell membranes. Instead, it is strongly dependent on proteins' preference for aggregation. Y1 - 2013 U6 - https://doi.org/10.1039/c2fd20085f SN - 1359-6640 VL - 161 IS - 1 SP - 397 EP - 417 PB - Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Javanainen, Matti A1 - Martinez-Seara, Hector A1 - Metzler, Ralf A1 - Vattulainen, Ilpo T1 - Diffusion of Integral Membrane Proteins in Protein-Rich Membranes JF - The journal of physical chemistry letters N2 - The lateral diffusion of embedded proteins along lipid membranes in protein-poor conditions has been successfully described in terms of the Saffman-Delbruck (SD) model, which predicts that the protein diffusion coefficient D is weakly dependent on its radius R as D proportional to ln(1/R). However, instead of being protein-poor, native cell membranes are extremely crowded with proteins. On the basis of extensive molecular simulations, we here demonstrate that protein crowding of the membrane at physiological levels leads to deviations from the SD relation and to the emergence of a stronger Stokes-like dependence D proportional to 1/R. We propose that this 1/R law mainly arises due to geometrical factors: smaller proteins are able to avoid confinement effects much better than their larger counterparts. The results highlight that the lateral dynamics in the crowded setting found in native membranes is radically different from protein-poor conditions and plays a significant role in formation of functional multiprotein complexes. Y1 - 2017 U6 - https://doi.org/10.1021/acs.jpclett.7b01758 SN - 1948-7185 VL - 8 SP - 4308 EP - 4313 PB - American Chemical Society CY - Washington ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Barkai, Eli A1 - Metzler, Ralf T1 - Noisy continuous time random walks JF - The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr N2 - Experimental studies of the diffusion of biomolecules within biological cells are routinely confronted with multiple sources of stochasticity, whose identification renders the detailed data analysis of single molecule trajectories quite intricate. Here, we consider subdiffusive continuous time random walks that represent a seminal model for the anomalous diffusion of tracer particles in complex environments. This motion is characterized by multiple trapping events with infinite mean sojourn time. In real physical situations, however, instead of the full immobilization predicted by the continuous time random walk model, the motion of the tracer particle shows additional jiggling, for instance, due to thermal agitation of the environment. We here present and analyze in detail an extension of the continuous time random walk model. Superimposing the multiple trapping behavior with additive Gaussian noise of variable strength, we demonstrate that the resulting process exhibits a rich variety of apparent dynamic regimes. In particular, such noisy continuous time random walks may appear ergodic, while the bare continuous time random walk exhibits weak ergodicity breaking. Detailed knowledge of this behavior will be useful for the truthful physical analysis of experimentally observed subdiffusion. Y1 - 2013 U6 - https://doi.org/10.1063/1.4816635 SN - 0021-9606 SN - 1089-7690 VL - 139 IS - 12 PB - American Institute of Physics CY - Melville ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion JF - Physical chemistry, chemical physics : PCCP N2 - Anomalous diffusion is frequently described by scaled Brownian motion (SBM){,} a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is ?x2(t)? [similar{,} equals] 2K(t)t with K(t) [similar{,} equals] t[small alpha]-1 for 0 < [small alpha] < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion{,} for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely{,} we demonstrate that under confinement{,} the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments{,} in particular{,} under confinement inside cellular compartments or when optical tweezers tracking methods are used. KW - single-particle tracking KW - living cells KW - random-walks KW - subdiffusion KW - dynamics KW - nonergodicity KW - coefficients KW - transport KW - membrane KW - behavior Y1 - 2014 U6 - https://doi.org/10.1039/C4CP02019G VL - 30 IS - 16 SP - 15811 EP - 15817 PB - The Royal Society of Chemistry CY - Cambridge ER - TY - JOUR A1 - Jeon, Jae-Hyung A1 - Chechkin, Aleksei V. A1 - Metzler, Ralf T1 - Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion JF - Physical chemistry, chemical physics : a journal of European Chemical Societies N2 - Anomalous diffusion is frequently described by scaled Brownian motion (SBM), a Gaussian process with a power-law time dependent diffusion coefficient. Its mean squared displacement is < x(2)(t) similar or equal to 2K(t)t with K(t) similar or equal to t(alpha-1) for 0 < alpha < 2. SBM may provide a seemingly adequate description in the case of unbounded diffusion, for which its probability density function coincides with that of fractional Brownian motion. Here we show that free SBM is weakly non-ergodic but does not exhibit a significant amplitude scatter of the time averaged mean squared displacement. More severely, we demonstrate that under confinement, the dynamics encoded by SBM is fundamentally different from both fractional Brownian motion and continuous time random walks. SBM is highly non-stationary and cannot provide a physical description for particles in a thermalised stationary system. Our findings have direct impact on the modelling of single particle tracking experiments, in particular, under confinement inside cellular compartments or when optical tweezers tracking methods are used. Y1 - 2014 U6 - https://doi.org/10.1039/c4cp02019g SN - 1463-9076 SN - 1463-9084 VL - 16 IS - 30 SP - 15811 EP - 15817 PB - Royal Society of Chemistry CY - Cambridge ER -