@article{AbdolvahabMetzlerEjtehadi2011, author = {Abdolvahab, Rouhollah Haji and Metzler, Ralf and Ejtehadi, Mohammad Reza}, title = {First passage time distribution of chaperone driven polymer translocation through a nanopore homopolymer and heteropolymer cases}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {135}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {24}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.3669427}, pages = {8}, year = {2011}, abstract = {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.}, language = {en} } @article{AdamcikJeonKarczewskietal.2012, author = {Adamcik, Jozef and Jeon, Jae-Hyung and Karczewski, Konrad J. and Metzler, Ralf and Dietler, Giovanni}, title = {Quantifying supercoiling-induced denaturation bubbles in DNA}, series = {Soft matter}, volume = {8}, journal = {Soft matter}, number = {33}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c2sm26089a}, pages = {8651 -- 8658}, year = {2012}, abstract = {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.}, language = {en} } @article{AkimotoCherstvyMetzler2018, author = {Akimoto, Takuma and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ergodicity, rejuvenation, enhancement, and slow relaxation of diffusion in biased continuous-time random walks}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {98}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.98.022105}, pages = {6}, year = {2018}, abstract = {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.}, language = {en} } @article{AwadMetzler2020, author = {Awad, Emad and Metzler, Ralf}, title = {Crossover dynamics from superdiffusion to subdiffusion}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {23}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {1}, publisher = {De Gruyter}, address = {Berlin ; Boston}, issn = {1311-0454}, doi = {10.1515/fca-2020-0003}, pages = {55 -- 102}, year = {2020}, abstract = {The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function.}, language = {en} } @article{AwadMetzler2022, author = {Awad, Emad and Metzler, Ralf}, title = {Closed-form multi-dimensional solutions and asymptotic behaviours for subdiffusive processes with crossovers: II. Accelerating case}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {20}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac5a90}, pages = {29}, year = {2022}, abstract = {Anomalous diffusion with a power-law time dependence vertical bar R vertical bar(2)(t) similar or equal to t(alpha i) of the mean squared displacement occurs quite ubiquitously in numerous complex systems. Often, this anomalous diffusion is characterised by crossovers between regimes with different anomalous diffusion exponents alpha(i). Here we consider the case when such a crossover occurs from a first regime with alpha(1) to a second regime with alpha(2) such that alpha(2) > alpha(1), i.e., accelerating anomalous diffusion. A widely used framework to describe such crossovers in a one-dimensional setting is the bi-fractional diffusion equation of the so-called modified type, involving two time-fractional derivatives defined in the Riemann-Liouville sense. We here generalise this bi-fractional diffusion equation to higher dimensions and derive its multidimensional propagator (Green's function) for the general case when also a space fractional derivative is present, taking into consideration long-ranged jumps (Levy flights). We derive the asymptotic behaviours for this propagator in both the short- and long-time as well the short- and long-distance regimes. Finally, we also calculate the mean squared displacement, skewness and kurtosis in all dimensions, demonstrating that in the general case the non-Gaussian shape of the probability density function changes.}, language = {en} } @article{AydinerCherstvyMetzler2018, author = {Aydiner, Ekrem and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Wealth distribution, Pareto law, and stretched exponential decay of money}, series = {Physica : europhysics journal ; A, Statistical mechanics and its applications}, volume = {490}, journal = {Physica : europhysics journal ; A, Statistical mechanics and its applications}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0378-4371}, doi = {10.1016/j.physa.2017.08.017}, pages = {278 -- 288}, year = {2018}, abstract = {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.}, language = {en} } @article{AydinerCherstvyMetzler2019, author = {Aydiner, Ekrem and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Money distribution in agent-based models with position-exchange dynamics}, series = {The European physical journal : B, Condensed matter and complex systems}, volume = {92}, journal = {The European physical journal : B, Condensed matter and complex systems}, number = {5}, publisher = {Springer}, address = {New York}, issn = {1434-6028}, doi = {10.1140/epjb/e2019-90674-0}, pages = {4}, year = {2019}, abstract = {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.}, language = {en} } @article{BarkaiGariniMetzler2012, author = {Barkai, Eli and Garini, Yuval and Metzler, Ralf}, title = {Strange Kinetics of single molecules in living cells}, series = {Physics today}, volume = {65}, journal = {Physics today}, number = {8}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0031-9228}, pages = {29 -- 35}, year = {2012}, language = {en} } @article{BauerGodecMetzler2014, author = {Bauer, Maximilian and Godec, Aljaz and Metzler, Ralf}, title = {Diffusion of finite-size particles in two-dimensional channels with random wall configurations}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {16}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {13}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c3cp55160a}, pages = {6118 -- 6128}, year = {2014}, abstract = {Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda}, language = {en} } @article{BauerGodecMetzler2014, author = {Bauer, Maximilian and Godec, Aljaž and Metzler, Ralf}, title = {Diffusion of finite-size particles in two-dimensional channels with random wall configurations}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, volume = {16}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, number = {13}, publisher = {RSC Publications}, address = {Cambridge}, issn = {1463-9084}, doi = {10.1039/C3CP55160A}, pages = {6118 -- 6128}, year = {2014}, abstract = {Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].}, language = {en} } @article{BauerMetzler2013, author = {Bauer, Maximilian and Metzler, Ralf}, title = {In vivo facilitated diffusion model}, series = {PLoS one}, volume = {8}, journal = {PLoS one}, number = {1}, publisher = {PLoS}, address = {San Fransisco}, issn = {1932-6203}, doi = {10.1371/journal.pone.0053956}, pages = {8}, year = {2013}, abstract = {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.}, language = {en} } @article{BauerMetzler2012, author = {Bauer, Maximilian and Metzler, Ralf}, title = {Generalized facilitated diffusion model for DNA-binding proteins with search and recognition states}, series = {Biophysical journal}, volume = {102}, journal = {Biophysical journal}, number = {10}, publisher = {Cell Press}, address = {Cambridge}, issn = {0006-3495}, doi = {10.1016/j.bpj.2012.04.008}, pages = {2321 -- 2330}, year = {2012}, abstract = {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.}, language = {en} } @article{BauerRasmussenLomholtetal.2015, author = {Bauer, Maximilian and Rasmussen, Emil S. and Lomholt, Michael A. and Metzler, Ralf}, title = {Real sequence effects on the search dynamics of transcription factors on DNA}, series = {Scientific reports}, volume = {5}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep10072}, pages = {13}, year = {2015}, abstract = {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.}, language = {en} } @article{BlavatskaMetzler2015, author = {Blavatska, Viktoria and Metzler, Ralf}, title = {Conformational properties of complex polymers: rosette versus star-like structures}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {48}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {13}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/48/13/135001}, pages = {14}, year = {2015}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Quantifying non-ergodic dynamics of force-free granular gases}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, number = {17}, issn = {1463-9084}, doi = {10.1039/C5CP02824H}, pages = {21791 -- 21798}, year = {2015}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Quantifying non-ergodic dynamics of force-free granular gases}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {17}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {34}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c5cp02824h}, pages = {21791 -- 21798}, year = {2015}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2015, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ultraslow scaled Brownian motion}, series = {New journal of physics : the open-access journal for physics}, volume = {17}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/17/6/063038}, pages = {16}, year = {2015}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2016, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Safdari, Hadiseh and Sokolov, Igor M. and Metzler, Ralf}, title = {Underdamped scaled Brownian motion}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publishing Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep30520}, year = {2016}, abstract = {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.}, language = {en} } @article{BodrovaChechkinCherstvyetal.2016, author = {Bodrova, Anna S. and Chechkin, Aleksei V. and Cherstvy, Andrey G. and Safdari, Hadiseh and Sokolov, Igor M. and Metzler, Ralf}, title = {Underdamped scaled Brownian motion: (non-)existence of the overdamped limit in anomalous diffusion}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep30520}, pages = {16}, year = {2016}, abstract = {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.}, language = {en} } @article{CaetanoCarvalhoMetzleretal.2020, author = {Caetano, Daniel L. Z. and Carvalho, Sidney Jurado de and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of multiple polyelectrolytes onto a nanosphere}, series = {Interface : journal of the Royal Society}, volume = {17}, journal = {Interface : journal of the Royal Society}, number = {167}, publisher = {Royal Society}, address = {London}, issn = {1742-5689}, doi = {10.1098/rsif.2020.0199}, pages = {10}, year = {2020}, abstract = {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.}, language = {en} } @article{CaetanodeCarvalhoMetzleretal.2017, author = {Caetano, Daniel L. Z. and de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of periodic and random polyampholytes onto charged surfaces}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {19}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c7cp04040g}, pages = {23397 -- 23413}, year = {2017}, abstract = {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.}, language = {en} } @article{CapałaPadashChechkinetal.2020, author = {Capała, Karol and Padash, Amin and Chechkin, Aleksei V. and Shokri, Babak and Metzler, Ralf and Dybiec, Bartłomiej}, title = {Levy noise-driven escape from arctangent potential wells}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {30}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {12}, publisher = {American Institute of Physics}, address = {Woodbury, NY}, issn = {1054-1500}, doi = {10.1063/5.0021795}, pages = {15}, year = {2020}, abstract = {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.}, language = {en} } @article{ChechkinKantzMetzler2017, author = {Chechkin, Aleksei V. and Kantz, Holger and Metzler, Ralf}, title = {Ageing effects in ultraslow continuous time random walks}, series = {The European physical journal : B, Condensed matter and complex systems}, volume = {90}, journal = {The European physical journal : B, Condensed matter and complex systems}, publisher = {Springer}, address = {New York}, issn = {1434-6028}, doi = {10.1140/epjb/e2017-80270-9}, pages = {12}, year = {2017}, abstract = {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.}, language = {en} } @article{ChechkinSenoMetzleretal.2017, author = {Chechkin, Aleksei V. and Seno, Flavio and Metzler, Ralf and Sokolov, Igor M.}, title = {Brownian yet Non-Gaussian Diffusion: From Superstatistics to Subordination of Diffusing Diffusivities}, series = {Physical review : X, Expanding access}, volume = {7}, journal = {Physical review : X, Expanding access}, publisher = {American Physical Society}, address = {College Park}, issn = {2160-3308}, doi = {10.1103/PhysRevX.7.021002}, pages = {20}, year = {2017}, abstract = {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.}, language = {en} } @article{ChechkinZaidLomholtetal.2013, author = {Chechkin, Aleksei V. and Zaid, I. M. and Lomholt, M. A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated surface diffusion on a cylinder in the fast exchange limit}, series = {Mathematical modelling of natural phenomena}, volume = {8}, journal = {Mathematical modelling of natural phenomena}, number = {2}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0973-5348}, doi = {10.1051/mmnp/20138208}, pages = {114 -- 126}, year = {2013}, abstract = {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.}, language = {en} } @article{ChechkinZaidLomholtetal.2012, author = {Chechkin, Aleksei V. and Zaid, Irwin M. and Lomholt, Michael A. and Sokolov, Igor M. and Metzler, Ralf}, title = {Bulk-mediated diffusion on a planar surface full solution}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {86}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.86.041101}, pages = {11}, year = {2012}, abstract = {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.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, series = {Soft matter}, volume = {2014}, journal = {Soft matter}, number = {10}, publisher = {Royal Society of Chemistry}, issn = {2046-2069}, doi = {10.1039/c3sm52846d}, pages = {1591 -- 1601}, year = {2014}, abstract = {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.}, language = {en} } @article{CherstvyChechkinMetzler2013, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes}, series = {New journal of physics : the open-access journal for physics}, volume = {15}, journal = {New journal of physics : the open-access journal for physics}, number = {15}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/15/8/083039}, pages = {13}, year = {2013}, abstract = {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.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Ageing and confinement in non-ergodic heterogeneous diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {48}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/48/485002}, pages = {18}, year = {2014}, abstract = {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.}, language = {en} } @article{CherstvyChechkinMetzler2014, author = {Cherstvy, Andrey G. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity}, series = {Soft matter}, volume = {10}, journal = {Soft matter}, number = {10}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c3sm52846d}, pages = {1591 -- 1601}, year = {2014}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2016, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, volume = {18}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European chemical societies}, publisher = {RSC Publ.}, address = {Cambridge}, issn = {1463-9084}, doi = {10.1039/C6CP03101C}, pages = {23840 -- 23852}, year = {2016}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2013, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {15}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {46}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c3cp53056f}, pages = {20220 -- 20235}, year = {2013}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2008, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Population splitting, trapping, and non-ergodicity in heterogeneous diffusion processes}, doi = {10.1039/C3CP53056F}, year = {2008}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2015, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ergodicity breaking, ageing, and confinement in generalized diffusion processes with position and time dependent diffusivity}, series = {Journal of statistical mechanics: theory and experiment}, journal = {Journal of statistical mechanics: theory and experiment}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1742-5468}, doi = {10.1088/1742-5468/2015/05/P05010}, pages = {20}, year = {2015}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2015, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Ergodicity breaking and particle spreading in noisy heterogeneous diffusion processes}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {142}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {14}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4917077}, pages = {11}, year = {2015}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2014, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Nonergodicity, fluctuations, and criticality in heterogeneous diffusion processes}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {90}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {1}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.90.012134}, pages = {11}, year = {2014}, abstract = {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.}, language = {en} } @article{CherstvyMetzler2016, author = {Cherstvy, Andrey G. and Metzler, Ralf}, title = {Anomalous diffusion in time-fluctuating non-stationary diffusivity landscapes}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {18}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c6cp03101c}, pages = {23840 -- 23852}, year = {2016}, abstract = {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.}, language = {en} } @article{CherstvyNagelBetaetal.2018, author = {Cherstvy, Andrey G. and Nagel, Oliver and Beta, Carsten and Metzler, Ralf}, title = {Non-Gaussianity, population heterogeneity, and transient superdiffusion in the spreading dynamics of amoeboid cells}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {20}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {35}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp04254c}, pages = {23034 -- 23054}, year = {2018}, abstract = {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.}, language = {en} } @article{CherstvySafdariMetzler2021, author = {Cherstvy, Andrey G. and Safdari, Hadiseh and Metzler, Ralf}, title = {Anomalous diffusion, nonergodicity, and ageing for exponentially and logarithmically time-dependent diffusivity}, series = {Journal of physics. D, Applied physics}, volume = {54}, journal = {Journal of physics. D, Applied physics}, number = {19}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0022-3727}, doi = {10.1088/1361-6463/abdff0}, pages = {18}, year = {2021}, abstract = {We investigate a diffusion process with a time-dependent diffusion coefficient, both exponentially increasing and decreasing in time, D(t)=D-0(e +/- 2 alpha t). For this (hypothetical) nonstationary diffusion process we compute-both analytically and from extensive stochastic simulations-the behavior of the ensemble- and time-averaged mean-squared displacements (MSDs) of the particles, both in the over- and underdamped limits. Simple asymptotic relations derived for the short- and long-time behaviors are shown to be in excellent agreement with the results of simulations. The diffusive characteristics in the presence of ageing are also considered, with dramatic differences of the over- versus underdamped regime. Our results for D(t)=D-0(e +/- 2 alpha t) extend and generalize the class of diffusive systems obeying scaled Brownian motion featuring a power-law-like variation of the diffusivity with time, D(t) similar to t(alpha-1). We also examine the logarithmically increasing diffusivity, D(t)=D(0)log[t/tau(0)], as another fundamental functional dependence (in addition to the power-law and exponential) and as an example of diffusivity slowly varying in time. One of the main conclusions is that the behavior of the massive particles is predominantly ergodic, while weak ergodicity breaking is repeatedly found for the time-dependent diffusion of the massless particles at short times. The latter manifests itself in the nonequivalence of the (both nonaged and aged) MSD and the mean time-averaged MSD. The current findings are potentially applicable to a class of physical systems out of thermal equilibrium where a rapid increase or decrease of the particles' diffusivity is inherently realized. One biological system potentially featuring all three types of time-dependent diffusion (power-law-like, exponential, and logarithmic) is water diffusion in the brain tissues, as we thoroughly discuss in the end.}, language = {en} } @article{CherstvyThapaMardoukhietal.2018, author = {Cherstvy, Andrey G. and Thapa, Samudrajit and Mardoukhi, Yousof and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averages and their statistical variation for the Ornstein-Uhlenbeck process}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {98}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.98.022134}, pages = {15}, year = {2018}, abstract = {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.}, language = {en} } @article{CherstvyThapaWagneretal.2019, author = {Cherstvy, Andrey G. and Thapa, Samudrajit and Wagner, Caroline E. and Metzler, Ralf}, title = {Non-Gaussian, non-ergodic, and non-Fickian diffusion of tracers in mucin hydrogels}, series = {Soft matter}, volume = {15}, journal = {Soft matter}, number = {12}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c8sm02096e}, pages = {2526 -- 2551}, year = {2019}, abstract = {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.}, language = {en} } @article{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, series = {New journal of physics}, volume = {19}, journal = {New journal of physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa7199}, pages = {1 -- 11}, year = {2017}, abstract = {We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black-Scholes-Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.}, language = {en} } @article{CherstvyVinodAghionetal.2017, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Time averaging, ageing and delay analysis of financial time series}, series = {New journal of physics : the open-access journal for physics}, volume = {19}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa7199}, pages = {135 -- 147}, year = {2017}, abstract = {We introduce three strategies for the analysis of financial time series based on time averaged observables. These comprise the time averaged mean squared displacement (MSD) as well as the ageing and delay time methods for varying fractions of the financial time series. We explore these concepts via statistical analysis of historic time series for several Dow Jones Industrial indices for the period from the 1960s to 2015. Remarkably, we discover a simple universal law for the delay time averaged MSD. The observed features of the financial time series dynamics agree well with our analytical results for the time averaged measurables for geometric Brownian motion, underlying the famed Black-Scholes-Merton model. The concepts we promote here are shown to be useful for financial data analysis and enable one to unveil new universal features of stock market dynamics.}, language = {en} } @article{CherstvyVinodAghionetal.2021, author = {Cherstvy, Andrey G. and Vinod, Deepak and Aghion, Erez and Sokolov, Igor M. and Metzler, Ralf}, title = {Scaled geometric Brownian motion features sub- or superexponential ensemble-averaged, but linear time-averaged mean-squared displacements}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {103}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {6}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.103.062127}, pages = {11}, year = {2021}, abstract = {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.}, language = {en} } @article{CherstvyWangMetzleretal.2021, author = {Cherstvy, Andrey G. and Wang, Wei and Metzler, Ralf and Sokolov, Igor M.}, title = {Inertia triggers nonergodicity of fractional Brownian motion}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {104}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.104.024115}, pages = {12}, year = {2021}, abstract = {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.}, language = {en} } @article{deCarvalhoMetzlerCherstvy2016, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces}, series = {New journal of physics : the open-access journal for physics}, volume = {18}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ.}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/18/8/083037}, year = {2016}, abstract = {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{\"u}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{\"u}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.}, language = {en} } @article{deCarvalhoMetzlerCherstvy2014, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of polyelectrolytes onto charged Janus nanospheres}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {16}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {29}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c4cp02207f}, pages = {15539 -- 15550}, year = {2014}, abstract = {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.}, language = {en} } @article{deCarvalhoMetzlerCherstvy2016, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Critical adsorption of polyelectrolytes onto planar and convex highly charged surfaces: the nonlinear Poisson-Boltzmann approach}, series = {NEW JOURNAL OF PHYSICS}, volume = {18}, journal = {NEW JOURNAL OF PHYSICS}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/18/8/083037}, pages = {17}, year = {2016}, abstract = {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.}, language = {en} } @article{deCarvalhoMetzlerCherstvy2015, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Inverted critical adsorption of polyelectrolytes in confinement}, series = {Soft matter}, journal = {Soft matter}, number = {11}, publisher = {Royal Society of Chemistry}, address = {London}, issn = {1744-6848}, doi = {10.1039/C5SM00635J}, pages = {4430 -- 4443}, year = {2015}, abstract = {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.}, language = {en} } @article{deCarvalhoMetzlerCherstvy2015, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Inverted critical adsorption of polyelectrolytes in confinement}, series = {Soft matter}, volume = {11}, journal = {Soft matter}, number = {22}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c5sm00635j}, pages = {4430 -- 4443}, year = {2015}, abstract = {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.}, language = {en} } @article{DoerriesChechkinSchumeretal.2022, author = {Doerries, Timo J. and Chechkin, Aleksei and Schumer, Rina and Metzler, Ralf}, title = {Rate equations, spatial moments, and concentration profiles for mobile-immobile models with power-law and mixed waiting time distributions}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {105}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {1}, publisher = {The American Institute of Physics}, address = {Woodbury, NY}, issn = {2470-0045}, doi = {10.1103/PhysRevE.105.014105}, pages = {24}, year = {2022}, abstract = {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.}, language = {en} } @article{DoerriesChechkinMetzler2022, author = {Doerries, Timo J. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Apparent anomalous diffusion and non-Gaussian distributions in a simple mobile-immobile transport model with Poissonian switching}, series = {Interface : journal of the Royal Society}, volume = {19}, journal = {Interface : journal of the Royal Society}, number = {192}, publisher = {Royal Society}, address = {London}, issn = {1742-5689}, doi = {10.1098/rsif.2022.0233}, pages = {14}, year = {2022}, abstract = {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.}, language = {en} } @article{DybiecCapalaChechkinetal.2018, author = {Dybiec, Bartlomiej and Capala, Karol and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Conservative random walks in confining potentials}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaefc2}, pages = {25}, year = {2018}, abstract = {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.}, language = {en} } @article{EliazarMetzler2013, author = {Eliazar, Iddo and Metzler, Ralf}, title = {Anomalous statistics of random relaxations in random environments}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {87}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.87.022141}, pages = {12}, year = {2013}, abstract = {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.}, language = {en} } @article{EliazarMetzler2012, author = {Eliazar, Iddo and Metzler, Ralf}, title = {The RARE model a generalized approach to random relaxation processes in disordered systems}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {137}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {23}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4770266}, pages = {9}, year = {2012}, abstract = {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.}, language = {en} } @article{EliazarMetzlerReuveni2019, author = {Eliazar, Iddo and Metzler, Ralf and Reuveni, Shlomi}, title = {Poisson-process limit laws yield Gumbel max-min and min-max}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {100}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.100.022129}, pages = {12}, year = {2019}, abstract = {"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}, language = {en} } @article{EliazarMetzlerReuveni2019, author = {Eliazar, Iddo and Metzler, Ralf and Reuveni, Shlomi}, title = {Gumbel central limit theorem for max-min and min-max}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {100}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.100.020104}, pages = {6}, year = {2019}, abstract = {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.}, language = {en} } @article{EmanuelCherstvyMetzleretal.2020, author = {Emanuel, Marc D. and Cherstvy, Andrey G. and Metzler, Ralf and Gompper, Gerhard}, title = {Buckling transitions and soft-phase invasion of two-component icosahedral shells}, series = {Physical review / publ. by The American Physical Society. E, Statistical, nonlinear, and soft matter physics}, volume = {102}, journal = {Physical review / publ. by The American Physical Society. E, Statistical, nonlinear, and soft matter physics}, number = {6}, publisher = {Woodbury}, address = {New York}, issn = {2470-0045}, doi = {10.1103/PhysRevE.102.062104}, pages = {26}, year = {2020}, abstract = {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.}, language = {en} } @article{EstradaDelvenneHatanoetal.2018, author = {Estrada, Ernesto and Delvenne, Jean-Charles and Hatano, Naomichi and Mateos, Jose L. and Metzler, Ralf and Riascos, Alejandro P. and Schaub, Michael T.}, title = {Random multi-hopper model}, series = {Journal of Complex Networks}, volume = {6}, journal = {Journal of Complex Networks}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {2051-1310}, doi = {10.1093/comnet/cnx043}, pages = {382 -- 403}, year = {2018}, abstract = {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.}, language = {en} } @article{FernandezCharcharCherstvyetal.2020, author = {Fernandez, Amanda Diez and Charchar, Patrick and Cherstvy, Andrey G. and Metzler, Ralf and Finnis, Michael W.}, title = {The diffusion of doxorubicin drug molecules in silica nanoslits is non-Gaussian, intermittent and anticorrelated}, series = {Physical chemistry, chemical physics}, volume = {22}, journal = {Physical chemistry, chemical physics}, number = {48}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/d0cp03849k}, pages = {27955 -- 27965}, year = {2020}, abstract = {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.}, language = {en} } @article{GhoshCherstvyGrebenkovetal.2016, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Grebenkov, Denis S. and Metzler, Ralf}, title = {Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments}, series = {NEW JOURNAL OF PHYSICS}, volume = {18}, journal = {NEW JOURNAL OF PHYSICS}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/18/1/013027}, pages = {16}, year = {2016}, abstract = {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.}, language = {en} } @article{GhoshCherstvyMetzler2014, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Non-universal tracer diffusion in crowded media of non-inert obstacles}, series = {Physical Chemistry Chemical Physics}, volume = {3}, journal = {Physical Chemistry Chemical Physics}, number = {17}, editor = {Metzler, Ralf}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, pages = {1847 -- 1858}, year = {2014}, abstract = {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.}, language = {en} } @article{GhoshCherstvyMetzler2014, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Deformation propagation in responsive polymer network films}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {141}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {7}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4893056}, pages = {9}, year = {2014}, abstract = {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.}, language = {en} } @article{GhoshCherstvyMetzler2015, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Non-universal tracer diffusion in crowded media of non-inert obstacles}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {17}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {3}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c4cp03599b}, pages = {1847 -- 1858}, year = {2015}, abstract = {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.}, language = {en} } @article{GhoshCherstvyPetrovetal.2016, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Petrov, Eugene P. and Metzler, Ralf}, title = {Interactions of rod-like particles on responsive elastic sheets}, series = {Soft matter}, journal = {Soft matter}, publisher = {RSC}, address = {London}, issn = {1744-6848}, doi = {10.1039/C6SM01522K}, year = {2016}, abstract = {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.}, language = {en} } @article{GhoshCherstvyPetrovetal.2016, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Petrov, Eugene P. and Metzler, Ralf}, title = {Interactions of rod-like particles on responsive elastic sheets}, series = {Soft matter}, volume = {12}, journal = {Soft matter}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c6sm01522k}, pages = {7908 -- 7919}, year = {2016}, abstract = {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.}, language = {en} } @article{GodecBauerMetzler2014, author = {Godec, Aljaz and Bauer, Maximilian and Metzler, Ralf}, title = {Collective dynamics effect transient subdiffusion of inert tracers in flexible gel networks}, series = {New journal of physics : the open-access journal for physics}, volume = {16}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/16/9/092002}, pages = {13}, year = {2014}, abstract = {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.}, language = {en} } @article{GodecChechkinBarkaietal.2014, author = {Godec, Aljaz and Chechkin, Aleksei V. and Barkai, Eli and Kantz, Holger and Metzler, Ralf}, title = {Localisation and universal fluctuations in ultraslow diffusion processes}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {47}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {49}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/47/49/492002}, pages = {10}, year = {2014}, abstract = {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.}, language = {en} } @article{GodecMetzler2013, author = {Godec, Aljaz and Metzler, Ralf}, title = {Finite-Time effects and ultraweak ergodicity breaking in superdiffusive dynamics}, series = {Physical review letters}, volume = {110}, journal = {Physical review letters}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.110.020603}, pages = {5}, year = {2013}, abstract = {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.}, language = {en} } @article{GodecMetzler2015, author = {Godec, Aljaz and Metzler, Ralf}, title = {Optimization and universality of Brownian search in a basic model of quenched heterogeneous media}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {91}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {5}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.91.052134}, pages = {17}, year = {2015}, abstract = {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.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep20349}, pages = {11}, year = {2016}, abstract = {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.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {Universal Proximity Effect in Target Search Kinetics in the Few-Encounter Limit}, series = {Physical review : X, Expanding access}, volume = {6}, journal = {Physical review : X, Expanding access}, publisher = {American Physical Society}, address = {College Park}, issn = {2160-3308}, doi = {10.1103/PhysRevX.6.041037}, pages = {11}, year = {2016}, abstract = {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.}, language = {en} } @article{GodecMetzler2016, author = {Godec, Aljaz and Metzler, Ralf}, title = {Active transport improves the precision of linear long distance molecular signalling}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {49}, journal = {Journal of physics : A, Mathematical and theoretical}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8113/49/36/364001}, pages = {11}, year = {2016}, abstract = {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.}, language = {en} } @article{GodecMetzler2017, author = {Godec, Aljaž and Metzler, Ralf}, title = {First passage time statistics for two-channel diffusion}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {50}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {8}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aa5204}, pages = {17}, year = {2017}, abstract = {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.}, language = {en} } @article{GoychukKharchenkoMetzler2014, author = {Goychuk, Igor A. and Kharchenko, Vasyl O. and Metzler, Ralf}, title = {Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion}, series = {Physical Chemistry Chemical Physics}, journal = {Physical Chemistry Chemical Physics}, number = {16}, publisher = {the Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, pages = {16524 -- 16535}, year = {2014}, abstract = {The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.}, language = {en} } @article{GoychukKharchenkoMetzler2014, author = {Goychuk, Igor and Kharchenko, Vasyl O. and Metzler, Ralf}, title = {How Molecular Motors Work in the Crowded Environment of Living Cells: Coexistence and Efficiency of Normal and Anomalous Transport}, series = {PLoS one}, volume = {9}, journal = {PLoS one}, number = {3}, publisher = {PLoS}, address = {San Fransisco}, issn = {1932-6203}, doi = {10.1371/journal.pone.0091700}, pages = {7}, year = {2014}, abstract = {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.}, language = {en} } @article{GoychukKharchenkoMetzler2014, author = {Goychuk, Igor and Kharchenko, Vasyl O. and Metzler, Ralf}, title = {Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {16}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {31}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c4cp01234h}, pages = {16524 -- 16535}, year = {2014}, abstract = {The discovery of anomalous diffusion of larger biopolymers and submicron tracers such as endogenous granules, organelles, or virus capsids in living cells, attributed to the viscoelastic nature of the cytoplasm, provokes the question whether this complex environment equally impacts the active intracellular transport of submicron cargos by molecular motors such as kinesins: does the passive anomalous diffusion of free cargo always imply its anomalously slow active transport by motors, the mean transport distance along microtubule growing sublinearly rather than linearly in time? Here we analyze this question within the widely used two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubules driven by a flashing binding potential, where the cargo particle is elastically attached to the motor. Depending on the cargo size, the loading force, the amplitude of the binding potential, the turnover frequency of the molecular motor enzyme, and the linker stiffness we demonstrate that the motor transport may turn out either normal or anomalous, as indeed measured experimentally. We show how a highly efficient normal active transport mediated by motors may emerge despite the passive anomalous diffusion of the cargo, and study the intricate effects of the elastic linker. Under different, well specified conditions the microtubule-based motor transport becomes anomalously slow and thus significantly less efficient.}, language = {en} } @article{GranadoAbadMetzleretal.2020, author = {Granado, Felipe Le Vot and Abad, Enrique and Metzler, Ralf and Yuste, Santos B.}, title = {Continuous time random walk in a velocity field}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab9ae2}, pages = {27}, year = {2020}, abstract = {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{\´e}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{\´e}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.}, language = {en} } @article{GrebenkovMetzlerOshanin2021, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Distribution of first-reaction times with target regions on boundaries of shell-like domains}, series = {New Journal of Physics (NJP)}, volume = {2021}, journal = {New Journal of Physics (NJP)}, edition = {23}, publisher = {IOP Publishing}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac4282}, pages = {1 -- 23}, year = {2021}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2021, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {A molecular relay race: sequential first-passage events to the terminal reaction centre in a cascade of diffusion controlled processes}, series = {New Journal of Physics (NJP)}, volume = {23}, journal = {New Journal of Physics (NJP)}, publisher = {IOP - Institute of Physics Publishing}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac1e42}, pages = {18}, year = {2021}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2017, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains}, series = {New journal of physics : the open-access journal for physics}, volume = {19}, journal = {New journal of physics : the open-access journal for physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa8ed9}, pages = {11}, year = {2017}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2020, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {From single-particle stochastic kinetics to macroscopic reaction rates}, series = {New Journal of Physics}, volume = {22}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/abb1de}, pages = {28}, year = {2020}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2017, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains}, series = {New journal of physics}, volume = {19}, journal = {New journal of physics}, publisher = {IOP}, address = {London}, issn = {1367-2630}, doi = {10.1088/1367-2630/aa8ed9}, pages = {1 -- 11}, year = {2017}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2018, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Strong defocusing of molecular reaction times results from an interplay of geometry and reaction control}, series = {Communications Chemistry}, volume = {1}, journal = {Communications Chemistry}, publisher = {Macmillan Publishers Limited}, address = {London}, issn = {2399-3669}, doi = {10.1038/s42004-018-0096-x}, pages = {12}, year = {2018}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2019, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Full distribution of first exit times in the narrow escape problem}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab5de4}, pages = {23}, year = {2019}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2018, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Towards a full quantitative description of single-molecule reaction kinetics in biological cells}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {20}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {24}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp02043d}, pages = {16393 -- 16401}, year = {2018}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshanin2022, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Search efficiency in the Adam-Delbruck reduction-of-dimensionality scenario versus direct diffusive search}, series = {New journal of physics : the open-access journal for physics}, volume = {24}, journal = {New journal of physics : the open-access journal for physics}, number = {8}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac8824}, pages = {32}, year = {2022}, abstract = {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.}, language = {en} } @article{GrebenkovMetzlerOshaninetal.2019, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb and Dagdug, Leonardo and Berezhkovskii, Alexander M. and Skvortsov, Alexei T.}, title = {Trapping of diffusing particles by periodic absorbing rings on a cylindrical tube}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {150}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {20}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.5098390}, pages = {2}, year = {2019}, language = {en} } @article{GrebenkovSposiniMetzleretal.2020, author = {Grebenkov, Denis S. and Sposini, Vittoria and Metzler, Ralf and Oshanin, Gleb and Seno, Flavio}, title = {Exact distributions of the maximum and range of random diffusivity processes}, series = {New Journal of Physics}, volume = {23}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges.}, address = {Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/abd313}, pages = {23}, year = {2020}, abstract = {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.}, language = {en} } @article{GuggenbergerChechkinMetzler2022, author = {Guggenberger, Tobias and Chechkin, Aleksei and Metzler, Ralf}, title = {Absence of stationary states and non-Boltzmann distributions of fractional Brownian motion in shallow external potentials}, series = {New journal of physics : the open-access journal for physics}, volume = {24}, journal = {New journal of physics : the open-access journal for physics}, number = {7}, publisher = {Dt. Physikalische Ges.}, address = {[Bad Honnef]}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac7b3c}, pages = {18}, year = {2022}, abstract = {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.}, language = {en} } @article{GuggenbergerChechkinMetzler2021, author = {Guggenberger, Tobias and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {54}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {29}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac019b}, pages = {17}, year = {2021}, abstract = {We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.}, language = {en} } @article{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {New Journal of Physics}, volume = {21}, journal = {New Journal of Physics}, publisher = {Deutsche Physikalische Gesellschaft ; IOP, Institute of Physics}, address = {Bad Honnef und London}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab075f}, pages = {13}, year = {2019}, abstract = {Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.}, language = {en} } @article{HerrmannMetzler2017, author = {Herrmann, Carl J. J. and Metzler, Ralf}, title = {A self-avoiding walk with neural delays as a model of fixational eye movements}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Springer Nature}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-13489-8}, pages = {1 -- 17}, year = {2017}, abstract = {Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.}, language = {en} } @article{HerrmannMetzlerEngbert2017, author = {Herrmann, Carl J. J. and Metzler, Ralf and Engbert, Ralf}, title = {A self-avoiding walk with neural delays as a model of fixational eye movements}, series = {Scientific reports}, volume = {7}, journal = {Scientific reports}, publisher = {Nature Publ. Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/s41598-017-13489-8}, pages = {17}, year = {2017}, abstract = {Fixational eye movements show scaling behaviour of the positional mean-squared displacement with a characteristic transition from persistence to antipersistence for increasing time-lag. These statistical patterns were found to be mainly shaped by microsaccades (fast, small-amplitude movements). However, our re-analysis of fixational eye-movement data provides evidence that the slow component (physiological drift) of the eyes exhibits scaling behaviour of the mean-squared displacement that varies across human participants. These results suggest that drift is a correlated movement that interacts with microsaccades. Moreover, on the long time scale, the mean-squared displacement of the drift shows oscillations, which is also present in the displacement auto-correlation function. This finding lends support to the presence of time-delayed feedback in the control of drift movements. Based on an earlier non-linear delayed feedback model of fixational eye movements, we propose and discuss different versions of a new model that combines a self-avoiding walk with time delay. As a result, we identify a model that reproduces oscillatory correlation functions, the transition from persistence to antipersistence, and microsaccades.}, language = {en} } @article{HouCherstvyMetzleretal.2018, author = {Hou, Ru and Cherstvy, Andrey G. and Metzler, Ralf and Akimoto, Takuma}, title = {Biased continuous-time random walks for ordinary and equilibrium cases}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {20}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {32}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp01863d}, pages = {20827 -- 20848}, year = {2018}, abstract = {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.}, language = {en} } @article{JavanainenHammarenMonticellietal.2013, author = {Javanainen, Matti and Hammaren, Henrik and Monticelli, Luca and Jeon, Jae-Hyung and Miettinen, Markus S. and Martinez-Seara, Hector and Metzler, Ralf and Vattulainen, Ilpo}, title = {Anomalous and normal diffusion of proteins and lipids in crowded lipid membranes}, series = {Faraday discussions}, volume = {161}, journal = {Faraday discussions}, number = {1}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1359-6640}, doi = {10.1039/c2fd20085f}, pages = {397 -- 417}, year = {2013}, abstract = {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.}, language = {en} } @article{JavanainenMartinezSearaMetzleretal.2017, author = {Javanainen, Matti and Martinez-Seara, Hector and Metzler, Ralf and Vattulainen, Ilpo}, title = {Diffusion of Integral Membrane Proteins in Protein-Rich Membranes}, series = {The journal of physical chemistry letters}, volume = {8}, journal = {The journal of physical chemistry letters}, publisher = {American Chemical Society}, address = {Washington}, issn = {1948-7185}, doi = {10.1021/acs.jpclett.7b01758}, pages = {4308 -- 4313}, year = {2017}, abstract = {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.}, language = {en} } @article{JeonBarkaiMetzler2013, author = {Jeon, Jae-Hyung and Barkai, Eli and Metzler, Ralf}, title = {Noisy continuous time random walks}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {139}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {12}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4816635}, pages = {15}, year = {2013}, abstract = {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.}, language = {en} } @article{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, series = {Physical chemistry, chemical physics : PCCP}, volume = {30}, journal = {Physical chemistry, chemical physics : PCCP}, number = {16}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, doi = {10.1039/C4CP02019G}, pages = {15811 -- 15817}, year = {2014}, abstract = {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.}, language = {en} } @article{JeonChechkinMetzler2014, author = {Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Scaled Brownian motion: a paradoxical process with a time dependent diffusivity for the description of anomalous diffusion}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {16}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {30}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c4cp02019g}, pages = {15811 -- 15817}, year = {2014}, abstract = {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.}, language = {en} }