@misc{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}, number = {168}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-74021}, 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{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} } @misc{LiMeiXuetal.2020, author = {Li, Yongge and Mei, Ruoxing and Xu, Yong and Kurths, J{\"u}rgen and Duan, Jinqiao and Metzler, Ralf}, title = {Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {974}, issn = {1866-8372}, doi = {10.25932/publishup-47454}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-474542}, pages = {28}, year = {2020}, abstract = {This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the It{\^o} calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.}, language = {en} } @article{LiMeiXuetal.2020, author = {Li, Yongge and Mei, Ruoxing and Xu, Yong and Kurths, J{\"u}rgen and Duan, Jinqiao and Metzler, Ralf}, title = {Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity}, 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/ab81b9}, pages = {27}, year = {2020}, abstract = {This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the It{\^o} calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.}, 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{NezhadhaghighiChechkinMetzler2014, author = {Nezhadhaghighi, M. Ghasemi and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Numerical approach to unbiased and driven generalized elastic model}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {140}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {2}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.4858425}, pages = {9}, year = {2014}, abstract = {From scaling arguments and numerical simulations, we investigate the properties of the generalized elastic model (GEM) that is used to describe various physical systems such as polymers, membranes, single-file systems, or rough interfaces. We compare analytical and numerical results for the subdiffusion exponent beta characterizing the growth of the mean squared displacement <(delta h)(2)> of the field h described by the GEM dynamic equation. We study the scaling properties of the qth order moments with time, finding that the interface fluctuations show no intermittent behavior. We also investigate the ergodic properties of the process h in terms of the ergodicity breaking parameter and the distribution of the time averaged mean squared displacement. Finally, we study numerically the driven GEM with a constant, localized perturbation and extract the characteristics of the average drift for a tagged probe.}, 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{VinodCherstvyWangetal.2022, author = {Vinod, Deepak and Cherstvy, Andrey G. and Wang, Wei and Metzler, Ralf and Sokolov, Igor M.}, title = {Nonergodicity of reset geometric Brownian motion}, 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 = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.105.L012106}, pages = {4}, year = {2022}, abstract = {We derive. the ensemble-and time-averaged mean-squared displacements (MSD, TAMSD) for Poisson-reset geometric Brownian motion (GBM), in agreement with simulations. We find MSD and TAMSD saturation for frequent resetting, quantify the spread of TAMSDs via the ergodicity-breaking parameter and compute distributions of prices. General MSD-TAMSD nonequivalence proves reset GBM nonergodic.}, 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} } @misc{GhoshCherstvyMetzler2014, author = {Ghosh, Surya K. and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Non-universal tracer diffusion in crowded media of non-inert obstacles}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-77128}, 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{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{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{KlettCherstvyShinetal.2021, author = {Klett, Kolja and Cherstvy, Andrey G. and Shin, Jaeoh and Sokolov, Igor M. and Metzler, Ralf}, title = {Non-Gaussian, transiently anomalous, and ergodic self-diffusion of flexible dumbbells in crowded two-dimensional environments}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {104}, 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.104.064603}, pages = {18}, year = {2021}, abstract = {We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed.}, 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} } @misc{MetzlerJeonCherstvy2016, author = {Metzler, Ralf and Jeon, J. -H. and Cherstvy, Andrey G.}, title = {Non-Brownian diffusion in lipid membranes: Experiments and simulations}, series = {Biochimica et biophysica acta : Biomembranes}, volume = {1858}, journal = {Biochimica et biophysica acta : Biomembranes}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0005-2736}, doi = {10.1016/j.bbamem.2016.01.022}, pages = {2451 -- 2467}, year = {2016}, abstract = {The dynamics of constituents and the surface response of cellular membranes also in connection to the binding of various particles and macromolecules to the membrane are still a matter of controversy in the membrane biophysics community, particularly with respect to crowded membranes of living biological cells. We here put into perspective recent single particle tracking experiments in the plasma membranes of living cells and supercomputing studies of lipid bilayer model membranes with and without protein crowding. Special emphasis is put on the observation of anomalous, non-Brownian diffusion of both lipid molecules and proteins embedded in the lipid bilayer. While single component, pure lipid bilayers in simulations exhibit only transient anomalous diffusion of lipid molecules on nanosecond time scales, the persistence of anomalous diffusion becomes significantly longer ranged on the addition of disorder through the addition of cholesterol or proteins and on passing of the membrane lipids to the gel phase. Concurrently, experiments demonstrate the anomalous diffusion of membrane embedded proteins up to macroscopic time scales in the minute time range. Particular emphasis will be put on the physical character of the anomalous diffusion, in particular, the occurrence of ageing observed in the experiments the effective diffusivity of the measured particles is a decreasing function of time. Moreover, we present results for the time dependent local scaling exponent of the mean squared displacement of the monitored particles. Recent results finding deviations from the commonly assumed Gaussian diffusion patterns in protein crowded membranes are reported. The properties of the displacement autocorrelation function of the lipid molecules are discussed in the light of their appropriate physical anomalous diffusion models, both for non-crowded and crowded membranes. In the last part of this review we address the upcoming field of membrane distortion by elongated membrane-binding particles. We discuss how membrane compartmentalisation and the particle-membrane binding energy may impact the dynamics and response of lipid membranes. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Rog. (C) 2016 The Authors. Published by Elsevier B.V.}, 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{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{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{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} } @misc{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}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76478}, 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{KosztolowiczMetzlerWąsiketal.2020, author = {Kosztolowicz, Tadeusz and Metzler, Ralf and Wąsik, Slawomir and Arabski, Michal}, title = {Modelling experimentally measured of ciprofloxacin antibiotic diffusion in Pseudomonas aeruginosa biofilm formed in artificial sputum medium}, series = {PLoS ONE}, volume = {15}, journal = {PLoS ONE}, publisher = {PLOS}, address = {San Francisco, California, US}, issn = {1932-6203}, doi = {10.1371/journal.pone.0243003}, pages = {14}, year = {2020}, abstract = {We study the experimentally measured ciprofloxacin antibiotic diffusion through a gel-like artificial sputum medium (ASM) mimicking physiological conditions typical for a cystic fibrosis layer, in which regions occupied by Pseudomonas aeruginosa bacteria are present. To quantify the antibiotic diffusion dynamics we employ a phenomenological model using a subdiffusion-absorption equation with a fractional time derivative. This effective equation describes molecular diffusion in a medium structured akin Thompson's plumpudding model; here the 'pudding' background represents the ASM and the 'plums' represent the bacterial biofilm. The pudding is a subdiffusion barrier for antibiotic molecules that can affect bacteria found in plums. For the experimental study we use an interferometric method to determine the time evolution of the amount of antibiotic that has diffused through the biofilm. The theoretical model shows that this function is qualitatively different depending on whether or not absorption of the antibiotic in the biofilm occurs. We show that the process can be divided into three successive stages: (1) only antibiotic subdiffusion with constant biofilm parameters, (2) subdiffusion and absorption of antibiotic molecules with variable biofilm transport parameters, (3) subdiffusion and absorption in the medium but the biofilm parameters are constant again. Stage 2 is interpreted as the appearance of an intensive defence build-up of bacteria against the action of the antibiotic, and in the stage 3 it is likely that the bacteria have been inactivated. Times at which stages change are determined from the experimentally obtained temporal evolution of the amount of antibiotic that has diffused through the ASM with bacteria. Our analysis shows good agreement between experimental and theoretical results and is consistent with the biologically expected biofilm response. We show that an experimental method to study the temporal evolution of the amount of a substance that has diffused through a biofilm is useful in studying the processes occurring in a biofilm. We also show that the complicated biological process of antibiotic diffusion in a biofilm can be described by a fractional subdiffusion-absorption equation with subdiffusion and absorption parameters that change over time.}, language = {en} } @misc{KosztolowiczMetzlerWąsiketal.2020, author = {Kosztolowicz, Tadeusz and Metzler, Ralf and Wąsik, Slawomir and Arabski, Michal}, title = {Modelling experimentally measured of ciprofloxacin antibiotic diffusion in Pseudomonas aeruginosa biofilm formed in artificial sputum medium}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1107}, issn = {1866-8372}, doi = {10.25932/publishup-49086}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-490866}, pages = {16}, year = {2020}, abstract = {We study the experimentally measured ciprofloxacin antibiotic diffusion through a gel-like artificial sputum medium (ASM) mimicking physiological conditions typical for a cystic fibrosis layer, in which regions occupied by Pseudomonas aeruginosa bacteria are present. To quantify the antibiotic diffusion dynamics we employ a phenomenological model using a subdiffusion-absorption equation with a fractional time derivative. This effective equation describes molecular diffusion in a medium structured akin Thompson's plumpudding model; here the 'pudding' background represents the ASM and the 'plums' represent the bacterial biofilm. The pudding is a subdiffusion barrier for antibiotic molecules that can affect bacteria found in plums. For the experimental study we use an interferometric method to determine the time evolution of the amount of antibiotic that has diffused through the biofilm. The theoretical model shows that this function is qualitatively different depending on whether or not absorption of the antibiotic in the biofilm occurs. We show that the process can be divided into three successive stages: (1) only antibiotic subdiffusion with constant biofilm parameters, (2) subdiffusion and absorption of antibiotic molecules with variable biofilm transport parameters, (3) subdiffusion and absorption in the medium but the biofilm parameters are constant again. Stage 2 is interpreted as the appearance of an intensive defence build-up of bacteria against the action of the antibiotic, and in the stage 3 it is likely that the bacteria have been inactivated. Times at which stages change are determined from the experimentally obtained temporal evolution of the amount of antibiotic that has diffused through the ASM with bacteria. Our analysis shows good agreement between experimental and theoretical results and is consistent with the biologically expected biofilm response. We show that an experimental method to study the temporal evolution of the amount of a substance that has diffused through a biofilm is useful in studying the processes occurring in a biofilm. We also show that the complicated biological process of antibiotic diffusion in a biofilm can be described by a fractional subdiffusion-absorption equation with subdiffusion and absorption parameters that change over time.}, language = {en} } @article{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Mixing and segregation of ring polymers: spatial confinement and molecular crowding effects}, 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/5/053047}, pages = {19}, year = {2014}, abstract = {During the life cycle of bacterial cells the non-mixing of the two ring-shaped daughter genomes is an important prerequisite for the cell division process. Mimicking the environments inside highly crowded biological cells, we study the dynamics and statistical behavior of two flexible ring polymers in the presence of cylindrical confinement and crowding molecules. From extensive computer simulations we determine the degree of ring-ring overlap and the number of inter-monomer contacts for varying volume fractions phi of crowders. We also examine the entropic demixing of polymer rings in the presence of mobile crowders and determine the characteristic times of the internal polymer dynamics. Effects of the ring length on ring-ring overlap are also analyzed. In particular, on systematic variation of the fraction of crowding molecules, a (1 - phi)-scaling is found for the ring-ring overlap length along the cylinder axis, and a non-monotonic dependence of the 3D ring-ring contact number with a maximum at phi approximate to 0.2 is obtained. Our results demonstrate that polymer rings are demixed and separated by particular entropy-favourable partitioning of crowders along the axis of the cylindrical simulation box. These findings help to rationalize the implications of macromolecular crowding for circular DNA molecules in confined spaces inside bacteria as well as in localized cellular compartments inside eukaryotic cells.}, language = {en} } @article{LomholtLizanaMetzleretal.2013, author = {Lomholt, Michael A. and Lizana, Ludvig and Metzler, Ralf and Ambjoernsson, Tobias}, title = {Microscopic origin of the logarithmic time evolution of aging processes in complex systems}, series = {Physical review letters}, volume = {110}, journal = {Physical review letters}, number = {20}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.110.208301}, pages = {5}, year = {2013}, abstract = {There exists compelling experimental evidence in numerous systems for logarithmically slow time evolution, yet its full theoretical understanding remains elusive. We here introduce and study a generic transition process in complex systems, based on nonrenewal, aging waiting times. Each state n of the system follows a local clock initiated at t = 0. The random time tau between clock ticks follows the waiting time density psi (tau). Transitions between states occur only at local clock ticks and are hence triggered by the local forward waiting time, rather than by psi (tau). For power-law forms psi (tau) similar or equal to tau(-1-alpha) (0 < alpha < 1) we obtain a logarithmic time evolution of the state number < n(t)> similar or equal to log(t/t(0)), while for alpha > 2 the process becomes normal in the sense that < n(t)> similar or equal to t. In the intermediate range 1 < alpha < 2 we find the power-law growth < n(t)> similar or equal to t(alpha-1). Our model provides a universal description for transition dynamics between aging and nonaging states.}, language = {en} } @misc{NorregaardMetzlerRitteretal.2017, author = {Norregaard, Kamilla and Metzler, Ralf and Ritter, Christine M. and Berg-Sorensen, Kirstine and Oddershede, Lene Broeng}, title = {Manipulation and Motion of Organelles and Single Molecules in Living Cells}, series = {Chemical reviews}, volume = {117}, journal = {Chemical reviews}, number = {5}, publisher = {American Chemical Society}, address = {Washington}, issn = {0009-2665}, doi = {10.1021/acs.chemrev.6b00638}, pages = {4342 -- 4375}, year = {2017}, abstract = {The biomolecule is among the most important building blocks of biological systems, and a full understanding of its function forms the scaffold for describing the mechanisms of higher order structures as organelles and cells. Force is a fundamental regulatory mechanism of biomolecular interactions driving many cellular processes. The forces on a molecular scale are exactly in the range that can be manipulated and probed with single molecule force spectroscopy. The natural environment of a biomolecule is inside a living cell, hence, this is the most relevant environment for probing their function. In vivo studies are, however, challenged by the complexity of the cell. In this review, we start with presenting relevant theoretical tools for analyzing single molecule data obtained in intracellular environments followed by a description of state-of-the art visualization techniques. The most commonly used force spectroscopy techniques, namely optical tweezers, magnetic tweezers, and atomic force microscopy, are described in detail, and their strength and limitations related to in vivo experiments are discussed. Finally, recent exciting discoveries within the field of in vivo manipulation and dynamics of single molecule and organelles are reviewed.}, language = {en} } @article{XuZhouMetzleretal.2020, author = {Xu, Pengbo and Zhou, Tian and Metzler, Ralf and Deng, Weihua}, title = {L{\´e}vy walk dynamics in an external harmonic potential}, series = {Physical review : E, Statistical, nonlinear, and soft matter physics}, volume = {101}, 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.101.062127}, pages = {12}, year = {2020}, abstract = {Levy walks (LWs) are spatiotemporally coupled random-walk processes describing superdiffusive heat conduction in solids, propagation of light in disordered optical materials, motion of molecular motors in living cells, or motion of animals, humans, robots, and viruses. We here investigate a key feature of LWs-their response to an external harmonic potential. In this generic setting for confined motion we demonstrate that LWs equilibrate exponentially and may assume a bimodal stationary distribution. We also show that the stationary distribution has a horizontal slope next to a reflecting boundary placed at the origin, in contrast to correlated superdiffusive processes. Our results generalize LWs to confining forces and settle some longstanding puzzles around LWs.}, 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{PadashAghionSchulzetal.2022, author = {Padash, Amin and Aghion, Erez and Schulz, Alexander and Barkai, Eli and Chechkin, Aleksei V. and Metzler, Ralf and Kantz, Holger}, title = {Local equilibrium properties of ultraslow diffusion in the Sinai model}, series = {New journal of physics}, volume = {24}, journal = {New journal of physics}, number = {7}, publisher = {IOP Publishing}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac7df8}, pages = {14}, year = {2022}, abstract = {We perform numerical studies of a thermally driven, overdamped particle in a random quenched force field, known as the Sinai model. We compare the unbounded motion on an infinite 1-dimensional domain to the motion in bounded domains with reflecting boundaries and show that the unbounded motion is at every time close to the equilibrium state of a finite system of growing size. This is due to time scale separation: inside wells of the random potential, there is relatively fast equilibration, while the motion across major potential barriers is ultraslow. Quantities studied by us are the time dependent mean squared displacement, the time dependent mean energy of an ensemble of particles, and the time dependent entropy of the probability distribution. Using a very fast numerical algorithm, we can explore times up top 10(17) steps and thereby also study finite-time crossover phenomena.}, 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{PalyulinChechkinMetzler2014, author = {Palyulin, Vladimir V. and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Levy flights do not always optimize random blind search for sparse targets}, series = {Proceedings of the National Academy of Sciences of the United States of America}, volume = {111}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {8}, publisher = {National Acad. of Sciences}, address = {Washington}, issn = {0027-8424}, doi = {10.1073/pnas.1320424111}, pages = {2931 -- 2936}, year = {2014}, abstract = {It is generally believed that random search processes based on scale-free, Levy stable jump length distributions (Levy flights) optimize the search for sparse targets. Here we show that this popular search advantage is less universal than commonly assumed. We study the efficiency of a minimalist search model based on Levy flights in the absence and presence of an external drift (underwater current, atmospheric wind, a preference of the walker owing to prior experience, or a general bias in an abstract search space) based on two different optimization criteria with respect to minimal search time and search reliability (cumulative arrival probability). Although Levy flights turn out to be efficient search processes when the target is far from the starting point, or when relative to the starting point the target is upstream, we show that for close targets and for downstream target positioning regular Brownian motion turns out to be the advantageous search strategy. Contrary to claims that Levy flights with a critical exponent alpha = 1 are optimal for the search of sparse targets in different settings, based on our optimization parameters the optimal a may range in the entire interval (1, 2) and especially include Brownian motion as the overall most efficient search strategy.}, language = {en} } @article{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {New Journal of Physics}, volume = {23}, journal = {New Journal of Physics}, publisher = {Dt. Physikalische Ges. ; IOP}, address = {Bad Honnef ; London}, issn = {1367-2630}, doi = {10.1088/1367-2630/abd50e}, pages = {22}, year = {2021}, abstract = {Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.}, language = {en} } @misc{ThapaWyłomańskaSikoraetal.2021, author = {Thapa, Samudrajit and Wyłomańska, Agnieszka and Sikora, Grzegorz and Wagner, Caroline E. and Krapf, Diego and Kantz, Holger and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1118}, issn = {1866-8372}, doi = {10.25932/publishup-49349}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-493494}, pages = {24}, year = {2021}, abstract = {Extensive time-series encoding the position of particles such as viruses, vesicles, or individualproteins are routinely garnered insingle-particle tracking experiments or supercomputing studies.They contain vital clues on how viruses spread or drugs may be delivered in biological cells.Similar time-series are being recorded of stock values in financial markets and of climate data.Such time-series are most typically evaluated in terms of time-averaged mean-squareddisplacements (TAMSDs), which remain random variables for finite measurement times. Theirstatistical properties are different for differentphysical stochastic processes, thus allowing us toextract valuable information on the stochastic process itself. To exploit the full potential of thestatistical information encoded in measured time-series we here propose an easy-to-implementand computationally inexpensive new methodology, based on deviations of the TAMSD from itsensemble average counterpart. Specifically, we use the upper bound of these deviations forBrownian motion (BM) to check the applicability of this approach to simulated and real data sets.By comparing the probability of deviations fordifferent data sets, we demonstrate how thetheoretical bound for BM reveals additional information about observed stochastic processes. Weapply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracerbeads measured in mucin hydrogels, and of geographic surface temperature anomalies. Ouranalysis shows how the large-deviation properties can be efficiently used as a simple yet effectiveroutine test to reject the BM hypothesis and unveil relevant information on statistical propertiessuch as ergodicity breaking and short-time correlations.}, language = {en} } @article{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, series = {Soft Matter}, journal = {Soft Matter}, editor = {Metzler, Ralf}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, pages = {472 -- 488}, year = {2014}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, language = {en} } @misc{ShinCherstvyMetzler2014, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, publisher = {The Royal Society of Chemistry}, address = {Cambridge}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76961}, pages = {472 -- 488}, year = {2014}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, language = {en} } @article{ShinCherstvyMetzler2015, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size}, series = {Soft matter}, volume = {11}, journal = {Soft matter}, number = {3}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1744-683X}, doi = {10.1039/c4sm02007c}, pages = {472 -- 488}, year = {2015}, abstract = {The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.}, 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} } @misc{deCarvalhoMetzlerCherstvy2015, author = {de Carvalho, Sidney J. and Metzler, Ralf and Cherstvy, Andrey G.}, title = {Inverted critical adsorption of polyelectrolytes in confinement}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-89562}, 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{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} } @misc{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}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-95882}, 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{XuMetzlerWang2022, author = {Xu, Pengbo and Metzler, Ralf and Wang, Wanli}, title = {Infinite density and relaxation for Levy walks in an external potential}, series = {Physical review}, volume = {105}, journal = {Physical review}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.105.044118}, pages = {15}, year = {2022}, abstract = {Levy walks are continuous-time random-walk processes with a spatiotemporal coupling of jump lengths and waiting times. We here apply the Hermite polynomial method to study the behavior of LWs with power-law walking time density for four different cases. First we show that the known result for the infinite density of an unconfined, unbiased LW is consistently recovered. We then derive the asymptotic behavior of the probability density function (PDF) for LWs in a constant force field, and we obtain the corresponding qth-order moments. In a harmonic external potential we derive the relaxation dynamic of the LW. For the case of a Poissonian walking time an exponential relaxation behavior is shown to emerge. Conversely, a power-law decay is obtained when the mean walking time diverges. Finally, we consider the case of an unconfined, unbiased LW with decaying speed v(r ) = v0/./r. When the mean walking time is finite, a universal Gaussian law for the position-PDF of the walker is obtained explicitly.}, 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{JeonMetzler2012, author = {Jeon, Jae-Hyung and Metzler, Ralf}, title = {Inequivalence of time and ensemble averages in ergodic systems: exponential versus power-law relaxation in confinement}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {85}, 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.85.021147}, pages = {8}, year = {2012}, abstract = {Single-particle tracking has become a standard tool for the investigation of diffusive properties, especially in small systems such as biological cells. Usually the resulting time series are analyzed in terms of time averages over individual trajectories. Here we study confined normal as well as anomalous diffusion, modeled by fractional Brownian motion and the fractional Langevin equation, and show that even for such ergodic systems time-averaged quantities behave differently from their ensemble-averaged counterparts, irrespective of how long the measurement time becomes. Knowledge of the exact behavior of time averages is therefore fundamental for the proper physical interpretation of measured time series, in particular, for extraction of the relaxation time scale from data.}, language = {en} } @article{StojkoskiJolakoskiPaletal.2022, author = {Stojkoski, Viktor and Jolakoski, Petar and Pal, Arnab and Sandev, Trifce and Kocarev, Ljupco and Metzler, Ralf}, title = {Income inequality and mobility in geometric Brownian motion with stochastic resetting: theoretical results and empirical evidence of non-ergodicity}, series = {Philosophical transactions of the Royal Society A: Mathematical, physical and engineering sciences}, volume = {380}, journal = {Philosophical transactions of the Royal Society A: Mathematical, physical and engineering sciences}, number = {2224}, publisher = {Royal Society}, address = {London}, issn = {1364-503X}, doi = {10.1098/rsta.2021.0157}, pages = {17}, year = {2022}, abstract = {We explore the role of non-ergodicity in the relationship between income inequality, the extent of concentration in the income distribution, and income mobility, the feasibility of an individual to change their position in the income rankings. For this purpose, we use the properties of an established model for income growth that includes 'resetting' as a stabilizing force to ensure stationary dynamics. We find that the dynamics of inequality is regime-dependent: it may range from a strictly non-ergodic state where this phenomenon has an increasing trend, up to a stable regime where inequality is steady and the system efficiently mimics ergodicity. Mobility measures, conversely, are always stable over time, but suggest that economies become less mobile in non-ergodic regimes. By fitting the model to empirical data for the income share of the top earners in the USA, we provide evidence that the income dynamics in this country is consistently in a regime in which non-ergodicity characterizes inequality and immobility. Our results can serve as a simple rationale for the observed real-world income dynamics and as such aid in addressing non-ergodicity in various empirical settings across the globe.This article is part of the theme issue 'Kinetic exchange models of societies and economies'.}, 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{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{PalyulinMetzler2012, author = {Palyulin, Vladimir V. and Metzler, Ralf}, title = {How a finite potential barrier decreases the mean first-passage time}, series = {Journal of statistical mechanics: theory and experiment}, journal = {Journal of statistical mechanics: theory and experiment}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1742-5468}, doi = {10.1088/1742-5468/2012/03/L03001}, pages = {10}, year = {2012}, abstract = {We consider the mean first-passage time of a random walker moving in a potential landscape on a finite interval, the starting and end points being at different potentials. From analytical calculations and Monte Carlo simulations we demonstrate that the mean first-passage time for a piecewise linear curve between these two points is minimized by the introduction of a potential barrier. Due to thermal fluctuations, this barrier may be crossed. It turns out that the corresponding expense for this activation is less severe than the gain from an increased slope towards the end point. In particular, the resulting mean first-passage time is shorter than for a linear potential drop between the two points.}, language = {en} } @article{SandevDomazetoskiKocarevetal.2022, author = {Sandev, Trifce and Domazetoski, Viktor and Kocarev, Ljupco and Metzler, Ralf and Chechkin, Aleksei}, title = {Heterogeneous diffusion with stochastic resetting}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {7}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac491c}, pages = {26}, year = {2022}, abstract = {We study a heterogeneous diffusion process (HDP) with position-dependent diffusion coefficient and Poissonian stochastic resetting. We find exact results for the mean squared displacement and the probability density function. The nonequilibrium steady state reached in the long time limit is studied. We also analyse the transition to the non-equilibrium steady state by finding the large deviation function. We found that similarly to the case of the normal diffusion process where the diffusion length grows like t (1/2) while the length scale xi(t) of the inner core region of the nonequilibrium steady state grows linearly with time t, in the HDP with diffusion length increasing like t ( p/2) the length scale xi(t) grows like t ( p ). The obtained results are verified by numerical solutions of the corresponding Langevin equation.}, language = {en} } @article{XuLiuLietal.2020, author = {Xu, Yong and Liu, Xuemei and Li, Yongge and Metzler, Ralf}, title = {Heterogeneous diffusion processes and nonergodicity with Gaussian colored noise in layered diffusivity landscapes}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {102}, 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.102.062106}, pages = {16}, year = {2020}, abstract = {Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent a. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent alpha is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.}, 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{MardoukhiJeonMetzler2015, author = {Mardoukhi, Yousof and Jeon, Jae-Hyung and Metzler, Ralf}, title = {Geometry controlled anomalous diffusion in random fractal geometries: looking beyond the infinite cluster}, 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 = {44}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c5cp03548a}, pages = {30134 -- 30147}, year = {2015}, abstract = {We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law similar to T-h with h < 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.}, language = {en} } @article{MardoukhiJeonMetzler2015, author = {Mardoukhi, Yousof and Jeon, Jae-Hyung and Metzler, Ralf}, title = {Geometry controlled anomalous diffusion in random fractal geometries}, 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}, publisher = {Wiley-VCH Verl.}, address = {Weinheim}, issn = {1439-7641}, doi = {10.1039/c5cp03548a}, pages = {30134 -- 30147}, year = {2015}, abstract = {We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law BT� h with h o 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.}, language = {en} } @misc{MardoukhiJeonMetzler2015, author = {Mardoukhi, Yousof and Jeon, Jae-Hyung and Metzler, Ralf}, title = {Geometry controlled anomalous diffusion in random fractal geometries}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-85247}, year = {2015}, abstract = {We investigate the ergodic properties of a random walker performing (anomalous) diffusion on a random fractal geometry. Extensive Monte Carlo simulations of the motion of tracer particles on an ensemble of realisations of percolation clusters are performed for a wide range of percolation densities. Single trajectories of the tracer motion are analysed to quantify the time averaged mean squared displacement (MSD) and to compare this with the ensemble averaged MSD of the particle motion. Other complementary physical observables associated with ergodicity are studied, as well. It turns out that the time averaged MSD of individual realisations exhibits non-vanishing fluctuations even in the limit of very long observation times as the percolation density approaches the critical value. This apparent non-ergodic behaviour concurs with the ergodic behaviour on the ensemble averaged level. We demonstrate how the non-vanishing fluctuations in single particle trajectories are analytically expressed in terms of the fractal dimension and the cluster size distribution of the random geometry, thus being of purely geometrical origin. Moreover, we reveal that the convergence scaling law to ergodicity, which is known to be inversely proportional to the observation time T for ergodic diffusion processes, follows a power-law BT� h with h o 1 due to the fractal structure of the accessible space. These results provide useful measures for differentiating the subdiffusion on random fractals from an otherwise closely related process, namely, fractional Brownian motion. Implications of our results on the analysis of single particle tracking experiments are provided.}, language = {en} } @article{TomovskiSandevMetzleretal.2012, author = {Tomovski, Zivorad and Sandev, Trifce and Metzler, Ralf and Dubbeldam, Johan}, title = {Generalized space-time fractional diffusion equation with composite fractional time derivative}, series = {Physica : europhysics journal ; A, Statistical mechanics and its applications}, volume = {391}, journal = {Physica : europhysics journal ; A, Statistical mechanics and its applications}, number = {8}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0378-4371}, doi = {10.1016/j.physa.2011.12.035}, pages = {2527 -- 2542}, year = {2012}, abstract = {We investigate the solution of space-time fractional diffusion equations with a generalized Riemann-Liouville time fractional derivative and Riesz-Feller space fractional derivative. The Laplace and Fourier transform methods are applied to solve the proposed fractional diffusion equation. The results are represented by using the Mittag-Leffler functions and the Fox H-function. Special cases of the initial and boundary conditions are considered. Numerical scheme and Grunwald-Letnikov approximation are also used to solve the space-time fractional diffusion equation. The fractional moments of the fundamental solution of the considered space-time fractional diffusion equation are obtained. Many known results are special cases of those obtained in this paper. We investigate also the solution of a space-time fractional diffusion equations with a singular term of the form delta(x). t-beta/Gamma(1-beta) (beta > 0).}, 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{StojkoskiSandevBasnarkovetal.2020, author = {Stojkoski, Viktor and Sandev, Trifce and Basnarkov, Lasko and Kocarev, Ljupco and Metzler, Ralf}, title = {Generalised geometric Brownian motion}, series = {Entropy}, volume = {22}, journal = {Entropy}, number = {12}, publisher = {MDPI}, address = {Basel}, issn = {1099-4300}, doi = {10.3390/e22121432}, pages = {34}, year = {2020}, abstract = {Classical option pricing schemes assume that the value of a financial asset follows a geometric Brownian motion (GBM). However, a growing body of studies suggest that a simple GBM trajectory is not an adequate representation for asset dynamics, due to irregularities found when comparing its properties with empirical distributions. As a solution, we investigate a generalisation of GBM where the introduction of a memory kernel critically determines the behaviour of the stochastic process. We find the general expressions for the moments, log-moments, and the expectation of the periodic log returns, and then obtain the corresponding probability density functions using the subordination approach. Particularly, we consider subdiffusive GBM (sGBM), tempered sGBM, a mix of GBM and sGBM, and a mix of sGBMs. We utilise the resulting generalised GBM (gGBM) in order to examine the empirical performance of a selected group of kernels in the pricing of European call options. Our results indicate that the performance of a kernel ultimately depends on the maturity of the option and its moneyness.}, language = {en} } @misc{Metzler2017, author = {Metzler, Ralf}, title = {Gaussianity Fair}, series = {Biophysical journal}, volume = {112}, journal = {Biophysical journal}, number = {3}, publisher = {Cell Press}, address = {Cambridge}, issn = {0006-3495}, doi = {10.1016/j.bpj.2016.12.019}, pages = {413 -- 415}, year = {2017}, language = {en} } @misc{GrebenkovMetzlerOshanin2019, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Full distribution of first exit times in the narrow escape problem}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {810}, issn = {1866-8372}, doi = {10.25932/publishup-44288}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-442883}, pages = {24}, 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{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} } @misc{GrebenkovMetzlerOshanin2020, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {From single-particle stochastic kinetics to macroscopic reaction rates}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1018}, issn = {1866-8372}, doi = {10.25932/publishup-48405}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-484059}, pages = {29}, 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{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{SandevMetzlerChechkin2018, author = {Sandev, Trifce and Metzler, Ralf and Chechkin, Aleksei V.}, title = {From continuous time random walks to the generalized diffusion equation}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {21}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {1}, publisher = {De Gruyter}, address = {Berlin}, issn = {1311-0454}, doi = {10.1515/fca-2018-0002}, pages = {10 -- 28}, year = {2018}, abstract = {We obtain a generalized diffusion equation in modified or Riemann-Liouville form from continuous time random walk theory. The waiting time probability density function and mean squared displacement for different forms of the equation are explicitly calculated. We show examples of generalized diffusion equations in normal or Caputo form that encode the same probability distribution functions as those obtained from the generalized diffusion equation in modified form. The obtained equations are general and many known fractional diffusion equations are included as special cases.}, language = {en} } @article{TomovskiMetzlerGerhold2022, author = {Tomovski, Živorad and Metzler, Ralf and Gerhold, Stefan}, title = {Fractional characteristic functions, and a fractional calculus approach for moments of random variables}, series = {Fractional calculus and applied analysis : an international journal for theory and applications}, volume = {25}, journal = {Fractional calculus and applied analysis : an international journal for theory and applications}, number = {4}, publisher = {De Gruyter}, address = {Berlin ; Boston}, issn = {1314-2224}, doi = {10.1007/s13540-022-00047-x}, pages = {1307 -- 1323}, year = {2022}, abstract = {In this paper we introduce a fractional variant of the characteristic function of a random variable. It exists on the whole real line, and is uniformly continuous. We show that fractional moments can be expressed in terms of Riemann-Liouville integrals and derivatives of the fractional characteristic function. The fractional moments are of interest in particular for distributions whose integer moments do not exist. Some illustrative examples for particular distributions are also presented.}, language = {en} } @article{WangCherstvyChechkinetal.2020, author = {Wang, Wei and Cherstvy, Andrey G. and Chechkin, Aleksei V. and Thapa, Samudrajit and Seno, Flavio and Liu, Xianbin and Metzler, Ralf}, title = {Fractional Brownian motion with random diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {53}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {47}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aba467}, pages = {34}, year = {2020}, abstract = {Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.}, 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} } @misc{GuggenbergerPagniniVojtaetal.2019, author = {Guggenberger, Tobias and Pagnini, Gianni and Vojta, Thomas and Metzler, Ralf}, title = {Fractional Brownian motion in a finite interval}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {755}, issn = {1866-8372}, doi = {10.25932/publishup-43666}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-436665}, 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{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{MardoukhiJeonChechkinetal.2018, author = {Mardoukhi, Yousof and Jeon, Jae-Hyung and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Fluctuations of random walks in critical random environments}, 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 = {31}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c8cp03212b}, pages = {20427 -- 20438}, year = {2018}, abstract = {Percolation networks have been widely used in the description of porous media but are now found to be relevant to understand the motion of particles in cellular membranes or the nucleus of biological cells. Random walks on the infinite cluster at criticality of a percolation network are asymptotically ergodic. On any finite size cluster of the network stationarity is reached at finite times, depending on the cluster's size. Despite of this we here demonstrate by combination of analytical calculations and simulations that at criticality the disorder and cluster size average of the ensemble of clusters leads to a non-vanishing variance of the time averaged mean squared displacement, regardless of the measurement time. Fluctuations of this relevant experimental quantity due to the disorder average of such ensembles are thus persistent and non-negligible. The relevance of our results for single particle tracking analysis in complex and biological systems is discussed.}, language = {en} } @article{KruesemannGodecMetzler2014, author = {Kr{\"u}semann, Henning and Godec, Aljaz and Metzler, Ralf}, title = {First-passage statistics for aging diffusion in systems with annealed and quenched disorder}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {89}, 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.89.040101}, pages = {5}, year = {2014}, abstract = {Aging, the dependence of the dynamics of a physical process on the time t(a) since its original preparation, is observed in systems ranging from the motion of charge carriers in amorphous semiconductors over the blinking dynamics of quantum dots to the tracer dispersion in living biological cells. Here we study the effects of aging on one of the most fundamental properties of a stochastic process, the first-passage dynamics. We find that for an aging continuous time random walk process, the scaling exponent of the density of first-passage times changes twice as the aging progresses and reveals an intermediate scaling regime. The first-passage dynamics depends on t(a) differently for intermediate and strong aging. Similar crossovers are obtained for the first-passage dynamics for a confined and driven particle. Comparison to the motion of an aged particle in the quenched trap model with a bias shows excellent agreement with our analytical findings. Our results demonstrate how first-passage measurements can be used to unravel the age t(a) of a physical system.}, language = {en} } @article{PadashChechkinDybiecetal.2019, author = {Padash, Amin and Chechkin, Aleksei V. and Dybiec, Bartlomiej and Pavlyukevich, Ilya and Shokri, Babak and Metzler, Ralf}, title = {First-passage properties of asymmetric Levy flights}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {45}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ab493e}, pages = {48}, year = {2019}, abstract = {L{\´e}vy flights are paradigmatic generalised random walk processes, in which the independent stationary increments—the 'jump lengths'—are drawn from an -stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of L{\´e}vy flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering L{\´e}vy flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of L{\´e}vy flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index and the skewness (asymmetry) parameter . The other approach is based on the stochastic Langevin equation with -stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times.}, language = {en} } @article{MattosMejiaMonasterioMetzleretal.2012, author = {Mattos, Thiago G. and Mejia-Monasterio, Carlos and Metzler, Ralf and Oshanin, Gleb}, title = {First passages in bounded domains When is the mean first passage time meaningful?}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {86}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {3}, publisher = {American Physical Society}, address = {College Park}, issn = {1539-3755}, doi = {10.1103/PhysRevE.86.031143}, pages = {8}, year = {2012}, abstract = {We study the first passage statistics to adsorbing boundaries of a Brownian motion in bounded two-dimensional domains of different shapes and configurations of the adsorbing and reflecting boundaries. From extensive numerical analysis we obtain the probability P(omega) distribution of the random variable omega = tau(1)/(tau(1) + tau(2)), which is a measure for how similar the first passage times tau(1) and tau(2) are of two independent realizations of a Brownian walk starting at the same location. We construct a chart for each domain, determining whether P(omega) represents a unimodal, bell-shaped form, or a bimodal, M-shaped behavior. While in the former case the mean first passage time (MFPT) is a valid characteristic of the first passage behavior, in the latter case it is an insufficient measure for the process. Strikingly we find a distinct turnover between the two modes of P(omega), characteristic for the domain shape and the respective location of absorbing and reflective boundaries. Our results demonstrate that large fluctuations of the first passage times may occur frequently in two-dimensional domains, rendering quite vague the general use of the MFPT as a robust measure of the actual behavior even in bounded domains, in which all moments of the first passage distribution exist.}, 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{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{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{SposiniChechkinMetzler2018, author = {Sposini, Vittoria and Chechkin, Aleksei V. and Metzler, Ralf}, title = {First passage statistics for diffusing diffusivity}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {52}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {4}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/aaf6ff}, pages = {11}, year = {2018}, abstract = {A rapidly increasing number of systems is identified in which the stochastic motion of tracer particles follows the Brownian law < r(2)(t)> similar or equal to Dt yet the distribution of particle displacements is strongly non-Gaussian. A central approach to describe this effect is the diffusing diffusivity (DD) model in which the diffusion coefficient itself is a stochastic quantity, mimicking heterogeneities of the environment encountered by the tracer particle on its path. We here quantify in terms of analytical and numerical approaches the first passage behaviour of the DD model. We observe significant modifications compared to Brownian-Gaussian diffusion, in particular that the DD model may have a faster first passage dynamics. Moreover we find a universal crossover point of the survival probability independent of the initial condition.}, language = {en} } @article{MutothyaXuLietal.2021, author = {Mutothya, Nicholas Mwilu and Xu, Yong and Li, Yongge and Metzler, Ralf and Mutua, Nicholas Muthama}, title = {First passage dynamics of stochastic motion in heterogeneous media driven by correlated white Gaussian and coloured non-Gaussian noises}, series = {Journal of physics. Complexity}, volume = {2}, journal = {Journal of physics. Complexity}, publisher = {IOP Publishing}, address = {Bristol}, issn = {2632-072X}, doi = {10.1088/2632-072X/ac35b5}, pages = {24}, year = {2021}, abstract = {We study the first passage dynamics for a diffusing particle experiencing a spatially varying diffusion coefficient while driven by correlated additive Gaussian white noise and multiplicative coloured non-Gaussian noise. We consider three functional forms for position dependence of the diffusion coefficient: power-law, exponential, and logarithmic. The coloured non-Gaussian noise is distributed according to Tsallis' q-distribution. Tracks of the non-Markovian systems are numerically simulated by using the fourth-order Runge-Kutta algorithm and the first passage times (FPTs) are recorded. The FPT density is determined along with the mean FPT (MFPT). Effects of the noise intensity and self-correlation of the multiplicative noise, the intensity of the additive noise, the cross-correlation strength, and the non-extensivity parameter on the MFPT are discussed.}, language = {en} } @misc{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V and Blackburn, George and Lomholt, Michael A and Watkins, Nicholas W and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of L{\´e}vy flights and L{\´e}vy walks}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {785}, issn = {1866-8372}, doi = {10.25932/publishup-43983}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-439832}, pages = {25}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}vy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For L{\´e}vy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the L{\´e}vy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.}, language = {en} } @article{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V and Blackburn, George and Lomholt, Michael A and Watkins, Nicholas W and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of L{\´e}vy flights and L{\´e}vy walks}, 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/ab41bb}, pages = {24}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}vy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For L{\´e}vy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the L{\´e}vy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.}, language = {en} } @article{PalyulinBlackburnLomholtetal.2019, author = {Palyulin, Vladimir V. and Blackburn, George and Lomholt, Michael A. and Watkins, Nicholas W. and Metzler, Ralf and Klages, Rainer and Chechkin, Aleksei V.}, title = {First passage and first hitting times of Levy flights and Levy walks}, series = {New journal of physics : the open-access journal for physics}, volume = {21}, journal = {New journal of physics : the open-access journal for physics}, number = {10}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ab41bb}, pages = {23}, year = {2019}, abstract = {For both L{\´e}vy flight and L{\´e}vy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For L{\´e}vy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the L{\´e}vy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.}, 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{ShinCherstvyKimetal.2015, author = {Shin, Jaeoh and Cherstvy, Andrey G. and Kim, Won Kyu and Metzler, Ralf}, title = {Facilitation of polymer looping and giant polymer diffusivity in crowded solutions of active particles}, 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/11/113008}, pages = {12}, year = {2015}, abstract = {We study the dynamics of polymer chains in a bath of self-propelled particles (SPP) by extensive Langevin dynamics simulations in a two-dimensional model system. Specifically, we analyse the polymer looping properties versus the SPP activity and investigate how the presence of the active particles alters the chain conformational statistics. We find that SPPs tend to extend flexible polymer chains, while they rather compactify stiffer semiflexible polymers, in agreement with previous results. Here we show that higher activities of SPPs yield a higher effective temperature of the bath and thus facilitate the looping kinetics of a passive polymer chain. We explicitly compute the looping probability and looping time in a wide range of the model parameters. We also analyse the motion of a monomeric tracer particle and the polymer's centre of mass in the presence of the active particles in terms of the time averaged mean squared displacement, revealing a giant diffusivity enhancement for the polymer chain via SPP pooling. Our results are applicable to rationalising the dimensions and looping kinetics of biopolymers at constantly fluctuating and often actively driven conditions inside biological cells or in suspensions of active colloidal particles or bacteria cells.}, language = {en} } @article{LiuCherstvyMetzler2017, author = {Liu, Lin and Cherstvy, Andrey G. and Metzler, Ralf}, title = {Facilitated Diffusion of Transcription Factor Proteins with Anomalous Bulk Diffusion}, series = {The journal of physical chemistry : B, Condensed matter, materials, surfaces, interfaces \& biophysical chemistry}, volume = {121}, journal = {The journal of physical chemistry : B, Condensed matter, materials, surfaces, interfaces \& biophysical chemistry}, publisher = {American Chemical Society}, address = {Washington}, issn = {1520-6106}, doi = {10.1021/acs.jpcb.6b12413}, pages = {1284 -- 1289}, year = {2017}, abstract = {What are the physical laws of the diffusive search of proteins for their specific binding sites on DNA in the presence of the macromolecular crowding in cells? We performed extensive computer simulations to elucidate the protein target search on DNA. The novel feature is the viscoelastic non-Brownian protein bulk diffusion recently observed experimentally. We examine the influence of the protein-DNA binding affinity and the anomalous diffusion exponent on the target search time. In all cases an optimal search time is found. The relative contribution of intermittent three-dimensional bulk diffusion and one-dimensional sliding of proteins along the DNA is quantified. Our results are discussed in the light of recent single molecule tracking experiments, aiming at a better understanding of the influence of anomalous kinetics of proteins on the facilitated diffusion mechanism.}, language = {en} } @article{ScottWeissSelhuberUnkeletal.2022, author = {Scott, Shane and Weiss, Matthias and Selhuber-Unkel, Christine and Barooji, Younes F. and Sabri, Adal and Erler, Janine T. and Metzler, Ralf and Oddershede, Lene B.}, title = {Extracting, quantifying, and comparing dynamical and biomechanical properties of living matter through single particle tracking}, series = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, volume = {25}, journal = {Physical chemistry, chemical physics : a journal of European Chemical Societies}, number = {3}, publisher = {RSC Publ.}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/d2cp01384c}, pages = {1513 -- 1537}, year = {2022}, abstract = {A panoply of new tools for tracking single particles and molecules has led to an explosion of experimental data, leading to novel insights into physical properties of living matter governing cellular development and function, health and disease. In this Perspective, we present tools to investigate the dynamics and mechanics of living systems from the molecular to cellular scale via single-particle techniques. In particular, we focus on methods to measure, interpret, and analyse complex data sets that are associated with forces, materials properties, transport, and emergent organisation phenomena within biological and soft-matter systems. Current approaches, challenges, and existing solutions in the associated fields are outlined in order to support the growing community of researchers at the interface of physics and the life sciences. Each section focuses not only on the general physical principles and the potential for understanding living matter, but also on details of practical data extraction and analysis, discussing limitations, interpretation, and comparison across different experimental realisations and theoretical frameworks. Particularly relevant results are introduced as examples. While this Perspective describes living matter from a physical perspective, highlighting experimental and theoretical physics techniques relevant for such systems, it is also meant to serve as a solid starting point for researchers in the life sciences interested in the implementation of biophysical methods.}, 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} } @misc{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 = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1142}, issn = {1866-8372}, doi = {10.25932/publishup-50397}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-503976}, pages = {24}, 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{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{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} } @misc{BarkaiGariniMetzler2013, author = {Barkai, Eli and Garini, Yuval and Metzler, Ralf}, title = {Electrostatic effects in living cells Reply}, series = {PHYSICS TODAY}, volume = {66}, journal = {PHYSICS TODAY}, number = {7}, publisher = {AMER INST PHYSICS}, address = {MELVILLE}, issn = {0031-9228}, pages = {11 -- 11}, year = {2013}, 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} } @misc{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}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-403726}, 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. 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{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{MejiaMonasterioMetzlerVollmer2020, author = {Mejia-Monasterio, Carlos and Metzler, Ralf and Vollmer, J{\"u}rgen}, title = {Editorial: anomalous transport}, series = {Frontiers in Physics}, volume = {8}, journal = {Frontiers in Physics}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-424X}, doi = {10.3389/fphy.2020.622417}, pages = {4}, year = {2020}, language = {en} } @article{SarabadaniMetzlerAlaNissila2022, author = {Sarabadani, Jalal and Metzler, Ralf and Ala-Nissila, Tapio}, title = {Driven polymer translocation into a channel: Isoflux tension propagation theory and Langevin dynamics simulations}, series = {Physical Review Research}, volume = {4}, journal = {Physical Review Research}, edition = {3}, publisher = {American Physical Society}, address = {College Park, Maryland, USA}, issn = {2643-1564}, doi = {10.1103/PhysRevResearch.4.033003}, pages = {033003-1 -- 033003-14}, year = {2022}, abstract = {Isoflux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f βc for the translocation time, where β = -1 is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics and the force exponent β = -1 regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant and the resulting entropic force causes the force exponent to drop below unity.}, language = {en} } @misc{SarabadaniMetzlerAlaNissila2022, author = {Sarabadani, Jalal and Metzler, Ralf and Ala-Nissila, Tapio}, title = {Driven polymer translocation into a channel: Isoflux tension propagation theory and Langevin dynamics simulations}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1292}, issn = {1866-8372}, doi = {10.25932/publishup-57438}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-574387}, pages = {033003-1 -- 033003-14}, year = {2022}, abstract = {Isoflux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f βc for the translocation time, where β = -1 is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics and the force exponent β = -1 regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant and the resulting entropic force causes the force exponent to drop below unity.}, language = {en} } @misc{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 = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, issn = {1866-8372}, doi = {10.25932/publishup-55754}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-557542}, 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 = {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{SandevChechkinKorabeletal.2015, author = {Sandev, Trifce and Chechkin, Aleksei V. and Korabel, Nickolay and Kantz, Holger and Sokolov, Igor M. and Metzler, Ralf}, title = {Distributed-order diffusion equations and multifractality: Models and solutions}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {92}, 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.92.042117}, pages = {19}, year = {2015}, abstract = {We study distributed-order time fractional diffusion equations characterized by multifractal memory kernels, in contrast to the simple power-law kernel of common time fractional diffusion equations. Based on the physical approach to anomalous diffusion provided by the seminal Scher-Montroll-Weiss continuous time random walk, we analyze both natural and modified-form distributed-order time fractional diffusion equations and compare the two approaches. The mean squared displacement is obtained and its limiting behavior analyzed. We derive the connection between the Wiener process, described by the conventional Langevin equation and the dynamics encoded by the distributed-order time fractional diffusion equation in terms of a generalized subordination of time. A detailed analysis of the multifractal properties of distributed-order diffusion equations is provided.}, language = {en} } @article{PulkkinenMetzler2013, author = {Pulkkinen, Otto and Metzler, Ralf}, title = {Distance matters the impact of gene proximity in bacterial gene regulation}, series = {Physical review letters}, volume = {110}, journal = {Physical review letters}, number = {19}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.110.198101}, pages = {5}, year = {2013}, abstract = {Following recent discoveries of colocalization of downstream-regulating genes in living cells, the impact of the spatial distance between such genes on the kinetics of gene product formation is increasingly recognized. We here show from analytical and numerical analysis that the distance between a transcription factor (TF) gene and its target gene drastically affects the speed and reliability of transcriptional regulation in bacterial cells. For an explicit model system, we develop a general theory for the interactions between a TF and a transcription unit. The observed variations in regulation efficiency are linked to the magnitude of the variation of the TF concentration peaks as a function of the binding site distance from the signal source. Our results support the role of rapid binding site search for gene colocalization and emphasize the role of local concentration differences.}, language = {en} }