Refine
Has Fulltext
- no (228)
Year of publication
Document Type
- Article (228) (remove)
Language
- English (228) (remove)
Is part of the Bibliography
- yes (228)
Keywords
- anomalous diffusion (34)
- diffusion (28)
- stochastic processes (10)
- ageing (5)
- first passage (5)
- geometric Brownian motion (5)
- Brownian motion (4)
- Levy flights (4)
- first passage time (4)
- fractional Brownian motion (4)
Institute
Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous–infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.
Polymer looping is controlled by macromolecular crowding, spatial confinement, and chain stiffness
(2015)
We study by extensive computer simulations the looping characteristics of linear polymers with varying persistence length inside a spherical cavity in the presence of macromolecular crowding. For stiff chains, the looping probability and looping time reveal wildly oscillating patterns as functions of the chain length. The effects of crowding differ dramatically for flexible versus stiff polymers. While for flexible chains the looping kinetics is slowed down by the crowders, for stiffer chains the kinetics turns out to be either decreased or facilitated, depending on the polymer length. For severe confinement, the looping kinetics may become strongly facilitated by crowding. Our findings are of broad impact for DNA looping in the crowded and compartmentalized interior of living biological cells.
“A chain is only as strong as its weakest link” says the proverb. But what about a collection of statistically identical chains: How long till all chains fail? The answer to this question is given by the max-min of a matrix whose (i,j)entry is the failure time of link j of chain i: take the minimum of each row, and then the maximum of the rows' minima. The corresponding min-max is obtained by taking the maximum of each column, and then the minimum of the columns' maxima. The min-max applies to the storage of critical data. Indeed, consider multiple backup copies of a set of critical data items, and consider the (i,j) matrix entry to be the time at which item j on copy i is lost; then, the min-max is the time at which the first critical data item is lost. In this paper we address random matrices whose entries are independent and identically distributed random variables. We establish Poisson-process limit laws for the row's minima and for the columns' maxima. Then, we further establish Gumbel limit laws for the max-min and for the min-max. The limit laws hold whenever the entries' distribution has a density, and yield highly applicable approximation tools and design tools for the max-min and min-max of large random matrices. A brief of the results presented herein is given in: Gumbel central limit theorem for max-min and min-max
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.
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.
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ô 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.
Optimization and universality of Brownian search in a basic model of quenched heterogeneous media
(2015)
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.
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 <vertical bar delta h vertical bar(q)> 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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
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.
Molecular motors pulling cargos in the viscoelastic cytosol: how power strokes beat subdiffusion
(2014)
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.
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.
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.
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.
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.
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.
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.
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.
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.
Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories
(2021)
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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'.
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.
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.
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.
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.
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.
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.
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.
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 BTh 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.
Generalized space-time fractional diffusion equation with composite fractional time derivative
(2012)
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).
Generalized facilitated diffusion model for DNA-binding proteins with search and recognition states
(2012)
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.
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.
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.