BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ralf Metzler

Ralf Metzler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechBiological Physicscond-mat.softQuantitative Methodsphysics.data-anSubcellular ProcessesBiomoleculesmath-phmath.MPMolecular Networksphysics.chem-phq-bioCell Behaviorcond-mat.othercond-mat.supr-conecon.GNMethodologynlin.AOnlin.PSphysics.comp-phphysics.soc-phPopulations and Evolutionq-fin.ECq-fin.PR

Catalog footprint

What is connected

95works
24topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

95 published item(s)

preprint2026arXiv

Conformable Scaling and Critical Dynamics: A Unified Framework for Phase Transitions

We investigate the application of conformable derivatives to model critical phenomena near continuous phase transitions. By incorporating a deformation parameter into the differential structure, we derive unified expressions for thermodynamic observables such as heat capacity, magnetization, susceptibility, and coherence length, each exhibiting power-law behavior near the critical temperature. The conformable derivative framework naturally embeds scale invariance and critical slowing down into the dynamics without resorting to fully nonlocal fractional calculus. Modified Ginzburg-Landau equations are constructed to model superconducting transitions, leading to analytical expressions for the order parameter and London penetration depth. Experimental data from niobium confirm the model's applicability, showing excellent fits and capturing asymmetric scaling behavior around Tc. This work offers a bridge between classical mean-field theory and generalized scaling frameworks, with implications for both theoretical modeling and experimental analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Absence of confinement and non-Boltzmann stationary states of fractional Brownian motion in shallow external potentials

We study the diffusive motion of a particle in a subharmonic potential of the form $U(x)=|x|^c$ ($0<c<2$) driven by long-range correlated, stationary fractional Gaussian noise $ξ_α(t)$ with $0<α\le2$. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent $α$. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation $c>2(1-1/α)$ holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function (PDF). Moreover we discuss analogies of non-stationarity of L{é}vy flights in shallow external potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bayesian inference of scaled versus fractional Brownian motion

We present a Bayesian inference scheme for scaled Brownian motion, and investigate its performance on synthetic data for parameter estimation and model selection in a combined inference with fractional Brownian motion. We include the possibility of measurement noise in both models. We find that for trajectories of a few hundred time points the procedure is able to resolve well the true model and parameters. Using the prior of the synthetic data generation process also for the inference, the approach is optimal based on decision theory. We include a comparison with inference using a prior different from the data generating one.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Classification of anomalous diffusion in animal movement data using power spectral analysis

The field of movement ecology has seen a rapid increase in high-resolution data in recent years, leading to the development of numerous statistical and numerical methods to analyse relocation trajectories. Data are often collected at the level of the individual and for long periods that may encompass a range of behaviours. Here, we use the power spectral density (PSD) to characterise the random movement patterns of a black-winged kite (Elanus caeruleus) and a white stork (Ciconia ciconia). The tracks are first segmented and clustered into different behaviours (movement modes), and for each mode we measure the PSD and the ageing properties of the process. For the foraging kite we find $1/f$ noise, previously reported in ecological systems mainly in the context of population dynamics, but not for movement data. We further suggest plausible models for each of the behavioural modes by comparing both the measured PSD exponents and the distribution of the single-trajectory PSD to known theoretical results and simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Driven polymer translocation into a channel: Iso-flux tension propagation theory and Langevin dynamics simulations

Iso-flux 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 $f_{\rm c}$. In the high driving force limit, regardless of the channel width, IFTP theory predicts $τ\propto f_{\textrm{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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First-passage times in conical varying-width channels biased by a transverse gravitational force: Comparison of analytical and numerical results

We study the crossing time statistic of diffusing point particles between the two ends of expanding and narrowing two-dimensional conical channels under a transverse external gravitational field. The theoretical expression for the mean first-passage time for such a system is derived under the assumption that the axial diffusion in a two-dimensional channel of smoothly varying geometry can be approximately described as a one-dimensional diffusion in an entropic potential with position-dependent effective diffusivity in terms of the modified Fick-Jacobs equation. We analyze the channel crossing dynamics in terms of the mean first-passage time, combining our analytical results with extensive two-dimensional Brownian dynamics simulations, allowing us to find the range of applicability of the one-dimensional approximation. We find that the effective particle diffusivity decreases with increasing amplitude of the external potential. Remarkably, the mean first-passage time for crossing the channel is shown to assume a minimum at finite values of the potential amplitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local equilibrium properties of ultraslow diffusion in the Sinai model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Steering self-organisation through confinement

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement occurs through interactions with boundaries, and can function as either a catalyst or inhibitor of self-organisation. It can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework for future research, we examine the role of confinement in self-organisation and identify overarching scientific challenges across disciplines that need to be addressed to harness its full scientific and technological potential. This framework will not only accelerate the generation of a common deeper understanding of self-organisation but also trigger the development of innovative strategies to steer it through confinement, with impact, e.g., on the design of smarter materials, tissue engineering for biomedicine and crowd management.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Unravelling the origins of anomalous diffusion: from molecules to migrating storks

Anomalous diffusion or, more generally, anomalous transport, with nonlinear dependence of the mean-squared displacement on the measurement time, is ubiquitous in nature. It has been observed in processes ranging from microscopic movement of molecules to macroscopic, large-scale paths of migrating birds. Using data from multiple empirical systems, spanning 12 orders of magnitude in length and 8 orders of magnitude in time, we employ a method to detect the individual underlying origins of anomalous diffusion and transport in the data. This method decomposes anomalous transport into three primary effects: long-range correlations ("Joseph effect"), fat-tailed probability density of increments ("Noah effect"), and non-stationarity ("Moses effect"). We show that such a decomposition of real-life data allows to infer nontrivial behavioral predictions, and to resolve open questions in the fields of single particle tracking in living cells and movement ecology.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Income inequality and mobility in geometric Brownian motion with stochastic resetting: theoretical results and empirical evidence of non-ergodicity

We explore the role of non-ergodicity in the relationship between income inequality, the extent of concentration in the income distribution, and mobility, the feasibility of an individual to change their position in the income distribution. For this purpose, we explore the properties of an established model for income growth that includes "resetting" as a stabilising force which ensures stationary dynamics. We find that the dynamics of inequality is regime-dependent and 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 ergodic behaviour. Mobility measures, conversely, are always stable over time, but the stationary value is dependent on the regime, suggesting that economies become less mobile in non-ergodic regimes. By fitting the model to empirical data for the dynamics of income share of the top earners in the United States, we provide evidence that the income dynamics in this country is consistently in a regime in which non-ergodicity characterises inequality and immobility dynamics. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

AnDi: The Anomalous Diffusion Challenge

The deviation from pure Brownian motion generally referred to as anomalous diffusion has received large attention in the scientific literature to describe many physical scenarios. Several methods, based on classical statistics and machine learning approaches, have been developed to characterize anomalous diffusion from experimental data, which are usually acquired as particle trajectories. With the aim to assess and compare the available methods to characterize anomalous diffusion, we have organized the Anomalous Diffusion (AnDi) Challenge (\url{http://www.andi-challenge.org/}). Specifically, the AnDi Challenge will address three different aspects of anomalous diffusion characterization, namely: (i) Inference of the anomalous diffusion exponent. (ii) Identification of the underlying diffusion model. (iii) Segmentation of trajectories. Each problem includes sub-tasks for different number of dimensions (1D, 2D and 3D). In order to compare the various methods, we have developed a dedicated open-source framework for the simulation of the anomalous diffusion trajectories that are used for the training and test datasets. The challenge was launched on March 1, 2020, and consists of three phases. Currently, the participation to the first phase is open. Submissions will be automatically evaluated and the performance of the top-scoring methods will be thoroughly analyzed and compared in an upcoming article.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diffusion of antibiotics through a biofilm in the presence of diffusion and absorption barriers

We propose a model of antibiotic diffusion through a bacterial biofilm when diffusion and/or absorption barriers develop in the biofilm. The idea of this model is: We deduce details of the diffusion process in a medium in which direct experimental study is difficult, based on probing diffusion in external regions. Since a biofilm has a gel-like consistency, we suppose that subdiffusion of particles in the biofilm may occur. To describe this process we use a fractional subdiffusion-absorption equation with an adjustable anomalous diffusion exponent. The boundary conditions at the boundaries of the biofilm are derived by means of a particle random walk model on a discrete lattice leading to an expression involving a fractional time derivative. We show that the temporal evolution of the total amount of substance that has diffused through the biofilm explicitly depends on whether there is antibiotic absorption in the biofilm. This fact is used to experimentally check for antibiotic absorption in the biofilm and if the biofilm parameters change over time. We propose a four-stage model of antibiotic diffusion in biofilm based on the mentioned above physical characteristics. The biological interpretation of the stages, in particular their relation with the bacterial defence mechanisms, is discussed. Theoretical results are compared with empirical results of ciprofloxacin diffusion through Pseudomonas aeruginosa biofilm, and ciprofloxacin and gentamicin diffusion through Proteus mirabilis biofilm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generalised geometric Brownian motion: Theory and applications to option pricing

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 develop a generalisation of GBM where the introduction of a memory kernel critically determines the behavior of the stochastic process. We find the general expressions for the moments, log-moments, and the expectation of the periodic log returns, and obtain the corresponding probability density functions by 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) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Leveraging large-deviation statistics to decipher the stochastic properties of measured trajectories

Extensive time-series encoding the position of particles such as viruses, vesicles, or individual proteins are routinely garnered in single-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-average mean-squared displacements, which remain random variables for finite measurement times. Their statistical properties are different for different physical stochastic processes, thus allowing us to extract valuable information on the stochastic process itself. To exploit the full potential of the statistical information encoded in measured time-series we here propose an easy-to-implement and computationally inexpensive new methodology, based on deviations of the time-averaged mean-squared displacement from its ensemble average counterpart. Specifically, we use the upper bound of these deviations for Brownian motion to check the applicability of this approach to simulated and real data sets. By comparing the probability of deviations for different data sets, we demonstrate how the theoretical bound for Brownian motion reveals additional information about observed stochastic processes. We apply the large-deviation method to data sets of tracer beads tracked in aqueous solution, tracer beads measured in mucin hydrogels, and of geographic surface temperature anomalies. Our analysis shows how the large-deviation properties can be efficiently used as a simple yet effective routine test to reject the Brownian motion hypothesis and unveil crucial information on statistical properties such as ergodicity breaking and short-time correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lévy noise-driven escape from arctan potential wells

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, L{é}vy 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 Lévy noise-driven escape from an almost rectangular, arctan potential well restricted by two absorbing boundaries. We unveil analogies of the observed transient dynamics to the general properties of stationary states of L{é}vy 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 arctan potential. Finally, we explore in detail the behavior of the probability densities of the first-escape time and the last-hitting point.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lévy walk dynamics in an external harmonic potential

Lévy 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 long-standing puzzles around LWs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reflected fractional Brownian motion in one and higher dimensions

Fractional Brownian motion (FBM), a non-Markovian self-similar Gaussian stochastic process with long-ranged correlations, represents a widely applied, paradigmatic mathematical model of anomalous diffusion. We report the results of large-scale computer simulations of FBM in one, two, and three dimensions in the presence of reflecting boundaries that confine the motion to finite regions in space. Generalizing earlier results for finite and semi-infinite one-dimensional intervals, we observe that the interplay between the long-time correlations of FBM and the reflecting boundaries leads to striking deviations of the stationary probability density from the uniform density found for normal diffusion. Particles accumulate at the boundaries for superdiffusive FBM while their density is depleted at the boundaries for subdiffusion. Specifically, the probability density $P$ develops a power-law singularity, $P\sim r^κ$, as function of the distance $r$ from the wall. We determine the exponent $κ$ as function of the dimensionality, the confining geometry, and the anomalous diffusion exponent $α$ of the FBM. We also discuss implications of our results, including an application to modeling serotonergic fiber density patterns in vertebrate brains.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unexpected crossovers in correlated random-diffusivity processes

The passive and active motion of micron-sized tracer particles in crowded liquids and inside living biological cells is ubiquitously characterised by "viscoelastic" anomalous diffusion, in which the increments of the motion feature long-ranged negative and positive correlations. While viscoelastic anomalous diffusion is typically modelled by a Gaussian process with correlated increments, so-called fractional Gaussian noise, an increasing number of systems are reported, in which viscoelastic anomalous diffusion is paired with non-Gaussian displacement distributions. Following recent advances in Brownian yet non-Gaussian diffusion we here introduce and discuss several possible versions of random-diffusivity models with long-ranged correlations. While all these models show a crossover from non-Gaussian to Gaussian distributions beyond some correlation time, their mean squared displacements exhibit strikingly different behaviours: depending on the model crossovers from anomalous to normal diffusion are observed, as well as unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the strong non-universality of random-diffusivity viscoelastic anomalous diffusion are important for the analysis of experiments and a better understanding of the physical origins of "viscoelastic yet non-Gaussian" diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Brownian motion and beyond: first-passage, power spectrum, non-Gaussianity, and anomalous diffusion

Brownian motion is a ubiquitous physical phenomenon across the sciences. After its discovery by Brown and intensive study since the first half of the 20th century, many different aspects of Brownian motion and stochastic processes in general have been addressed in Statistical Physics. In particular, there now exist a very large range of applications of stochastic processes in various disciplines. Here, we highlight some of the advances in stochastic processes prompted by novel experimental methods such as superresolution microscopy. Here we provide a summary of some of the recent developments highlighting both the experimental findings and theoretical frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Correlations and transport in exclusion processes with general finite memory

We consider the correlations and the hydrodynamic description of random walkers with a general finite memory moving on a $d$ dimensional hypercubic lattice. We derive a drift-diffusion equation and identify a memory-dependent critical density. Above the critical density, the effective diffusion coefficient decreases with the particles' propensity to move forward and below the critical density it increases with their propensity to move forward. If the correlations are neglected the critical density is exactly $1/2$. We also derive a low-density approximation for the same time correlations between different sites. We perform simulations on a one-dimensional system with one-step memory and find good agreement between our analytical derivation and the numerical results. We also consider the previously unexplored special case of totally anti-persistent particles. Generally, the correlation length converges to a finite value. However in the special case of totally anti-persistent particles and density $1/2$, the correlation length diverges with time. Furthermore, connecting a system of totally anti-persistent particles to external particle reservoirs creates a new phenomenon: In almost all systems, regardless of the precise details of the microscopic dynamics, when a system is connected to a reservoir, the mean density of particle at the edge is the same as the reservoir following the zeroth law of thermodynamics. In a totally anti-persistent system, however, the density at the edge is always higher than in the reservoir. We find a qualitative description of this phenomenon which agrees reasonably well with the numerics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Active transport improves the precision of linear long distance molecular signalling

Molecular signalling in living cells occurs at low copy numbers and is thereby inherently limited by the noise imposed by thermal diffusion. The precision at which biochemical receptors can count signalling molecules is intimately related to the noise correlation time. In addition to passive thermal diffusion, messenger RNA and vesicle-engulfed signalling molecules can transiently bind to molecular motors and are actively transported across biological cells. Active transport is most beneficial when trafficking occurs over large distances, for instance up to the order of 1 metre in neurons. Here we explain how intermittent active transport allows for faster equilibration upon a change in concentration triggered by biochemical stimuli. Moreover, we show how intermittent active excursions induce qualitative changes in the noise in effectively one-dimensional systems such as dendrites. Thereby they allow for significantly improved signalling precision in the sense of a smaller relative deviation in the concentration read-out by the receptor. On the basis of linear response theory we derive the exact mean field precision limit for counting actively transported molecules. We explain how intermittent active excursions disrupt the recurrence in the molecular motion, thereby facilitating improved signalling accuracy. Our results provide a deeper understanding of how recurrence affects molecular signalling precision in biological cells and novel diagnostic devices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comb model with slow and ultraslow diffusion

We consider a generalised diffusion equation in two dimensions for modeling diffusion on a comb-like structures. We analyse the probability distribution functions and we derive the mean squared displacement in $x$ and $y$ directions. Different forms of the memory kernels (Dirac delta, power-law, and distributed order) are considered. It is shown that anomalous diffusion may occur along both $x$ and $y$ directions. Ultraslow diffusion and some more general diffusive processes are observed as well. We give the corresponding continuous time random walk model for the considered two dimensional diffusion-like equation on a comb, and we derive the probability distribution functions which subordinate the process governed by this equation to the Wiener process.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ageing first passage time density in continuous time random walks and quenched energy landscapes

We study the first passage dynamics of an ageing stochastic process in the continuous time random walk (CTRW) framework. In such CTRW processes the test particle performs a random walk, in which successive steps are separated by random waiting times distributed in terms of the waiting time probability density function $ψ(t) \simeq t^{-1-α}$) ($0\le α\le 2$). An ageing stochastic process is defined by the explicit dependence of its dynamic quantities on the ageing time $t_a$, the time elapsed between its preparation and the start of the observation. Subdiffusive ageing CTRWs describe systems such as charge carriers in amorphous semiconductors, tracer dispersion in geological and biological systems, or the dynamics of blinking quantum dots. We derive the exact forms of the first passage time density for an ageing subdiffusive CTRW in the semi-infinite, confined, and biased case, finding different scaling regimes for weakly, intermediately, and strongly aged systems: these regimes, with different scaling laws, are also found when the scaling exponent is in the range $1<α<2$, for sufficiently long $t_a$. We compare our results with the ageing motion of a test particle in a quenched energy landscape. We test our theoretical results in the quenched landscape against simulations: only when the bias is strong enough, the correlations from returning to previously visited sites become insignificant and the results approach the aging CTRW results. With small or without bias, the ageing effects disappear and a change in the exponent compared to the case of a completely annealed landscape can be found, reflecting the build-up of correlations in the quenched landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Anomalous, non-Gaussian tracer diffusion in heterogeneously crowded environments

A topic of intense current investigation pursues the question how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain (nucleus), or at the outer boundary (cell membrane), exhibiting distinctly different behaviour of the observed anomalous diffusion for heterogeneous crowder distributions. Specifically we reveal an extremely asymmetric spreading of the tracer even at moderate crowding fractions. In addition to the standard mean squared displacement and the local diffusion exponent of the tracer particles we investigate the magnitude and the amplitude scatter of the time averaged mean squared displacement of individual trajectories, the non-Gaussianity parameter, and the van Hove correlation function of the particle displacements. We also quantify how the average tracer diffusivity varies with the position in the domain with heterogeneous radial distribution of the crowders and examine the behaviour of the survival probability and the dynamics of first passage events of the tracer. Finally, we discuss the relevance of our results to single particle tracking measurements in biological cells.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Conformational properties of complex polymers: rosette versus star-like structures

Multiple loop formation in polymer macromolecules is an important feature of the chromatin organization and DNA compactification in the nuclei. We analyze the size and shape characteristics of complex polymer structures, containing in general $f_1$ loops (petals) and $f_2$ linear chains (branches). Within the frames of continuous model of Gaussian macromolecule, we apply the path integration method and obtain the estimates for gyration radius $R_g$ and asphericity $\hat{A}$ of typical conformation as functions of parameters $f_1$, $f_2$. In particular, our results qualitatively reveal the extent of anisotropy of star-like topologies as compared to the rosette structures of the same total molecular weight.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ergodicity breaking and particle spreading in noisy heterogeneous diffusion processes

We study noisy heterogeneous diffusion processes with a position dependent diffusivity of the form $D(x)\sim D_0|x|^α$ in the presence of annealed and quenched disorder of the environment, corresponding to an effective variation of the exponent $α$ in time and space. In the case of annealed disorder, for which effectively $α=α(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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ergodicity breaking, ageing, and confinement in generalised diffusion processes with position and time dependent diffusivity

We study generalised anomalous diffusion processes whose diffusion coefficient $D(x,t)\sim D_0|x|^αt^β$ 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 generalised 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 generalised 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

First passage time distribution in heterogeneity controlled kinetics: going beyond the mean first passage time

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 a 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverted critical adsorption of polyelectrolytes in confinement

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 $σ_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 $σ_c$ for the concave interfaces versus the Debye screening length $1/κ$ and the extent of confinement $a$ for these three interfaces for small $κa$ values. For large $κa$ the critical adsorption condition approaches the 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/κ$. We also rationalize how $σ_c(κ)$ gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length onto critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonergodic dynamics of force-free granular gases

We study analytically and by event-driven molecular dynamics simulations the nonergodic and aging properties of force-free cooling granular gases with both constant and velocity-dependent (viscoelastic) restitution coefficient $\varepsilon$ for particle pair collisions. We compare the granular gas dynamics with an effective single particle stochastic model based on an underdamped Langevin equation with time dependent diffusivity. We find that both models share the same behavior of the ensemble mean squared displacement (MSD) and the velocity correlations in the small dissipation limit. However, we reveal that the time averaged MSD of granular gas particles significantly differs from this effective model due to ballistic correlations for systems with constant $\varepsilon$. For velocity-dependent $\varepsilon$ these corrections become weaker at longer times. Qualitatively the reported non-ergodic behavior is generic for granular gases with any realistic dependence of $\varepsilon$ on the impact velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimization and universality of Brownian search in quenched heterogeneous media

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 piece-wise 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 {\em 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantifying the non-ergodicity of scaled Brownian motion

We examine the non-ergodic properties of scaled Brownian motion, a non-stationary stochastic process with a time dependent diffusivity of the form $D(t)\simeq t^{α-1}$. We compute the ergodicity breaking parameter EB in the entire range of scaling exponents $α$, both analytically and via extensive computer simulations of the stochastic Langevin equation. We demonstrate that in the limit of long trajectory lengths $T$ and short lag times $Δ$ the EB parameter as function of the scaling exponent $α$ has no divergence at $α=1/2$ and present the asymptotes for EB in different limits. We generalise the analytical and simulations results for the time averaged and ergodic properties of scaled Brownian motion in the presence of ageing, that is, when the observation of the system starts only a finite time span after its initiation. The approach developed here for the calculation of the higher time averaged moments of the particle displacement can be applied to derive the ergodic properties of other stochastic processes such as fractional Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real sequence effects on the search dynamics of transcription factors on DNA

Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Self-subdiffusion in solutions of star-shaped crowders: non-monotonic effects of inter-particle interactions

We examine by extensive computer simulations the self-diffusion of anisotropic star like particles in crowded two-dimensional solutions. We investigate the implications of the area coverage fraction $ϕ$ of the crowders and the crowder-crowder adhesion properties on the regime of transient anomalous diffusion. We systematically compute the mean squared displacement (MSD) of the particles, their time averaged MSD, as well as the effective diffusion coefficient. The diffusion appears ergodic in the limit of long traces, such that the time averaged MSD converges towards the ensemble averaged MSD and features a small residual amplitude spread of the time averaged MSD from individual trajectories. At intermediate time scales we quantify the anomalous diffusion in the system. Also, we show that the translational---but not rotational---diffusivity of the particles $D$ is a non-monotonic function of the attraction strength between them. Both diffusion coefficients decrease as $D(ϕ)\sim (1-ϕ/ϕ^*)^2$ with the area fraction $ϕ$ occupied by the crowders. Our results might be applicable to rationalising the experimental observations of non-Brownian diffusion for a number of standard macromolecular crowders used in vitro to mimic the cytoplasmic conditions of living cells.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ultraslow scaled Brownian motion

We define and study in detail \emph{utraslow scaled Brownian motion (USBM)\/} characterised by a time dependent diffusion coefficient of the form $D(t)\simeq 1/t$. For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In an harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Breathing dynamics based parameter sensitivity analysis of hetero-polymeric DNA

We study the parameter sensitivity of hetero-polymeric DNA within the purview of DNA breathing dynamics. The degree of correlation between the mean bubble size and the model parameters are estimated for this purpose for three different DNA sequences. The analysis leads us to a better understanding of the sequence dependent nature of the breathing dynamics of hetero-polymeric DNA. Out of the fourteen model parameters for DNA stability in the statistical Poland-Scheraga approach, the hydrogen bond interaction $ε_{hb}(\mathtt{AT})$ for an $\mathtt{AT}$ base pair and the ring factor $ξ$ turn out to be the most sensitive parameters. In addition, the stacking interaction $ε_{st}(\mathtt{TA}-\mathtt{TA})$ for an $\mathtt{TA}-\mathtt{TA}$ nearest neighbor pair of base-pairs is found to be the most sensitive one among all stacking interactions. Moreover, we also establish that the nature of stacking interaction has a deciding effect on the DNA breathing dynamics, not the number of times a particular stacking interaction appears in a sequence. We show that the sensitivity analysis can be used as an effective measure to guide a stochastic optimization technique to find the kinetic rate constants related to the dynamics as opposed to the case where the rate constants are measured using the conventional unbiased way of optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collective dynamics effect transient subdiffusion of inert tracers in gel networks

Based on extensive Brownian dynamics simulations we study the thermally driven motion of a tracer bead in a cross-linked, dynamic gel network in the limit when the tracer bead's size is of the same size or even larger than the equilibrium mesh size of the gel. The analysis of long individual trajectories of the tracer bead demonstrates the existence of pronounced transient anomalous diffusion, accompanied by a drastic slow-down of the gel-bead relaxation dynamics. From the time averaged mean squared displacement and the van Hove cross-correlation function we elucidate the many-body origin of the non-Brownian tracer bead dynamics. Our results shed new light on the ongoing debate over the physical origin of sterical tracer interactions with structured environments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deformation propagation in responsive polymer network films

We study the elastic deformations in a cross-linked polymer network film triggered by the binding of submicron particles with a sticky surface, mimicking the interactions of viral pathogens with thin films of stimulus-responsive polymeric materials such as hydrogels. From extensive Langevin Dynamics simulations we quantify how far the network deformations propagate depending on the elasticity parameters of the network and the adhesion strength of the particles. We examine the dynamics of the collective area shrinkage of the network and obtain some simple relations for the associated characteristic decay lengths. A detailed analysis elucidates how the elastic energy of the network is distributed between stretching and compression modes in response to the particle binding. We also examine the force-distance curves of the repulsion or attraction interactions for a pair of sticky particles in the polymer network film as a function of the particle-particle separation. The results of this computational study provide new insight into collective phenomena in soft polymer network films and may, in particular, be applied to applications for visual detection of pathogens such as viruses via a macroscopic response of thin films of cross-linked hydrogels.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First passage statistics for aging diffusion in annealed and quenched disorder

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

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 slow 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 onto the looping kinetics in terms of an effective solution viscosity and standard excluded volume effects. Thus 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Localization and universal fluctuations in ultraslow diffusion processes

We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) $\langle x^2(t)\rangle\simeq\log^γt$. Comparison of annealed continuous time random walks (CTRWs) with logarithmic waiting time distribution $ψ(τ)\simeq1/(τ\log^{1+γ}τ)$ and Sinai diffusion in quenched random landscapes shows striking similarities, despite their very different physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time 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 and time averaged MSD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mixing and segregation of ring polymers: spatial confinement and molecular crowding effects

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 behaviour 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 $ϕ$ of crowders. We also examine the entropic de-mixing 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 analysed. In particular, on systematic variation of the fraction of crowding molecules a $(1-ϕ)$-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 is predicted. Our results help to rationalise the implications of macromolecular crowding for circular DNA molecules in confined spaces inside bacteria as well as in localised cellular compartments inside eukaryotic cells.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-ergodicity, fluctuations, and criticality in heterogeneous diffusion processes

We study the stochastic behavior of heterogeneous diffusion processes with the power-law dependence $D(x)\sim|x|^α$ 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 non-ergodic character of the heterogeneous diffusion process and, particularly, the degree of irreproducibility of individual realization. As we show, the fluctuations between individual realizations increase with growing modulus $|α|$ of the scaling exponent. The fluctuations appear to diverge when the critical value $α=2$ is approached, while for even larger $α$ 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 seem to 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-universal tracer diffusion in crowded media of non-inert obstacles

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 the 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimization of random search processes in the presence of an external bias

We study the efficiency of random search processes based on L{é}vy flights with power-law distributed jump lengths in the presence of an external drift, for instance, an underwater current, an airflow, or simply the bias of the searcher based on prior experience. While Lévy flights turn out to be efficient search processes when relative to the starting point the target is upstream, in the downstream scenario regular Brownian motion turns out to be advantageous. This is caused by the occurrence of leapovers of L{é}vy flights, due to which L{é}vy flights typically overshoot a point or small interval. Extending our recent work on biased LF search [V. V. Palyulin, A. V. Chechkin, and R. Metzler, Proc. Natl. Acad. Sci. USA, DOI:10.1073/pnas.1320424111] we establish criteria when the combination of the external stream and the initial distance between the starting point and the target favors L{é}vy flights over regular Brownian search. Contrary to the common belief that L{é}vy flights with a L{é}vy index $α=1$ (i.e., Cauchy flights) are optimal for sparse targets, we find that the optimal value for $α$ may range in the entire interval $(1,2)$ and include Brownian motion as the overall most efficient search strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Aging dynamics in interacting many-body systems

Low-dimensional, complex systems are often characterized by logarithmically slow dynamics. We study the generic motion of a labeled particle in an ensemble of identical diffusing particles with hardcore interactions in a strongly disordered, one-dimensional environment. Each particle in this single file is trapped for a random waiting time $τ$ with power law distribution $ψ(τ)\simeqτ^{-1- α}$, such that the $τ$ values are independent, local quantities for all particles. From scaling arguments and simulations, we find that for the scale-free waiting time case $0<α<1$, the tracer particle dynamics is ultra-slow with a logarithmic mean square displacement (MSD) $\langle x^2(t)\rangle\simeq(\log t)^{1/2}$. This extreme slowing down compared to regular single file motion $\langle x^2(t)\rangle\simeq t^{1/2}$ is due to the high likelihood that the labeled particle keeps encountering strongly immobilized neighbors. For the case $1<α<2$ we observe the MSD scaling $\langle x^2(t)\rangle\simeq t^γ$, where $γ<1/2$, while for $α>2$ we recover Harris law $\simeq t^{1/2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Aging renewal theory and application to random walks

The versatility of renewal theory is owed to its abstract formulation. Renewals can be interpreted as steps of a random walk, switching events in two-state models, domain crossings of a random motion, etc. We here discuss a renewal process in which successive events are separated by scale-free waiting time periods. Among other ubiquitous long time properties, this process exhibits aging: events counted initially in a time interval [0,t] statistically strongly differ from those observed at later times [t_a,t_a+t]. In complex, disordered media, processes with scale-free waiting times play a particularly prominent role. We set up a unified analytical foundation for such anomalous dynamics by discussing in detail the distribution of the aging renewal process. We analyze its half-discrete, half-continuous nature and study its aging time evolution. These results are readily used to discuss a scale-free anomalous diffusion process, the continuous time random walk. By this we not only shed light on the profound origins of its characteristic features, such as weak ergodicity breaking. Along the way, we also add an extended discussion on aging effects. In particular, we find that the aging behavior of time and ensemble averages is conceptually very distinct, but their time scaling is identical at high ages. Finally, we show how more complex motion models are readily constructed on the basis of aging renewal dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes

We demonstrate the non-ergodicity of a simple Markovian stochastic processes with space-dependent diffusion coefficient $D(x)$. For power-law forms $D(x) \simeq|x|^α$, this process yield anomalous diffusion of the form $\ < x^2(t)\ > \simeq t^{2/(2-α)}$. Interestingly, in both the sub- and superdiffusive regimes we observe weak ergodicity breaking: the scaling of the time averaged mean squared displacement $\{δ^2}$ remains \emph{linear} and thus differs from the corresponding ensemble average $\ <x^2(t)\ >$. We analyze the non-ergodic behavior of this process in terms of the ergodicity breaking parameters and the distribution of amplitude scatter of $\{δ^2}$. This model represents an alternative approach to non-ergodic, anomalous diffusion that might be particularly relevant for diffusion in heterogeneous media.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Correlated continuous-time random walks: combining scale-invariance with long-range memory for spatial and temporal dynamics

Standard continuous time random walk (CTRW) models are renewal processes in the sense that at each jump a new, independent pair of jump length and waiting time are chosen. Globally, anomalous diffusion emerges through action of the generalized central limit theorem leading to scale-free forms of the jump length or waiting time distributions. Here we present a modified version of recently proposed correlated CTRW processes, where we incorporate a power-law correlated noise on the level of both jump length and waiting time dynamics. We obtain a very general stochastic model, that encompasses key features of several paradigmatic models of anomalous diffusion: discontinuous, scale-free displacements as in Levy flights, scale-free waiting times as in subdiffusive CTRWs, and the long-range temporal correlations of fractional Brownian motion (FBM). We derive the exact solutions for the single-time probability density functions and extract the scaling behaviours. Interestingly, we find that different combinations of the model parameters lead to indistinguishable shapes of the emerging probability density functions and identical scaling laws. Our model will be useful to describe recent experimental single particle tracking data, that feature a combination of CTRW and FBM properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Diffusion of finite-size particles in channels with random walls

Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick-Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-sized particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. With the application of diffusing pathogens in hydrogels in mind, we consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Distance matters: the impact of gene proximity in bacterial gene regulation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First passage behaviour of multi-dimensional fractional Brownian motion and application to reaction phenomena

Fractional Brownian motion is a generalised Gaussian diffusive process that is found to describe numerous stochastic phenomena in physics and biology. Here we introduce a multi-dimensional fractional Brownian motion (FBM) defined as a superposition of conventional FBM for each coordinate in analogy to multi-dimensional Brownian motion, and study its first passage properties. Starting from the well-established first passage time statistics of one-dimensional FBM and the associated approximation schemes, we explore the first passage time behaviour of multi-dimensional FBM and compare these results with simulations. The asymptotic kinetic behaviour of diffusion-limited reactions of reactant particles performing FBM in a one- and multi-dimensional space is studied based on the corresponding first passage time statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

In vivo facilitated diffusion model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular motors transporting cargos in viscoelastic cytosol: how to beat subdiffusion with a power stroke?

Anomalously slow passive diffusion, $\langle δx^2(t)\rangle\simeq t^α$, with $0<α<1$, of larger tracers such as messenger RNA and endogenous submicron granules in the cytoplasm of living biological cells has been demonstrated in a number of experiments and has been attributed to the viscoelastic physical nature of the cellular cytoplasm. This finding provokes the question to which extent active intracellular transport is affected by this viscoelastic environment: does the subdiffusion of free submicron cargo such as vesicles and organelles always imply anomalously slow transport by molecular motors such as kinesins, that is, directed transport characterized by a sublinear growth of the mean distance, $\langle x(t)\rangle\simeq t^{α_{\rm eff}}$, with $0<α_{\rm eff}<1$? Here we study a generic model approach combining the commonly accepted two-state Brownian ratchet model of kinesin motors based on the continuous-state diffusion along microtubule driven by a flashing binding potential. The motor is elastically coupled to a cargo particle, which in turn is subject to the viscoelastic cytoplasmic environment. Depending on the physical parameters of cargo size, loading force, amplitude of the binding potential, and the turnover frequency of the molecular motor, the transport can be both normal ($α_{\rm eff}=1$) and anomalous ($α\leq α_{\rm eff}<1$). In particular, we demonstrate in detail how highly efficient normal motor transport can emerge despite the anomalously slow passive diffusion of cargo particles, and how the active motion of the same motor in the same cell may turn anomalously slow when the parameters are changed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical approach to unbiased and driven generalized elastic model

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 <|delta h|^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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Particle invasion, survival, and non-ergodicity in 2D diffusion processes with space-dependent diffusivity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantification of noise in the bifunctionality-induced post-translational modification

We present a generic analytical scheme for the quantification of fluctuations due to bifunctionality-induced signal transduction within the members of bacterial two-component system. The proposed model takes into account post-translational modifications in terms of elementary phosphotransfer kinetics. Sources of fluctuations due to autophosphorylation, kinase and phosphatase activity of the sensor kinase have been considered in the model via Langevin equations, which are then solved within the framework of linear noise approximation. The resultant analytical expression of phosphorylated response regulators are then used to quantify the noise profile of biologically motivated single and branched pathways. Enhancement and reduction of noise in terms of extra phosphate outflux and influx, respectively, have been analyzed for the branched system. Furthermore, role of fluctuations of the network output in the regulation of a promoter with random activation/deactivation dynamics has been analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Searching with and against the stream: Levy or Brown?

We study the efficiency of search processes based on Levy flights (LFs) with power-law distributed jump lengths in the presence of an external drift. While LFs turn out to be efficient search processes when relative to the starting point the target is upstream, in the downstream scenario regular Brownian motion turns out to be advantageous. This is caused by the occurrence of leapovers of LFs, due to which LFs typically overshoot a point in space. We establish criteria when the combination of the external stream and the initial distance between the starting point and the target favors LFs over regular Brownian search. Contrary to the common belief that LFs with a stable index alpha=1 are optimal, we find that the optimal alpha may range in the entire interval (1,2) and even include Brownian motion as the overall most efficient search strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sensing viruses by mechanical tension of DNA in responsive hydrogels

The rapid worldwide spread of severe viral infections, often involving novel modifications of viruses, poses major challenges to our health care systems. This means that tools that can efficiently and specifically diagnose viruses are much needed. To be relevant for a broad application in local health care centers, such tools should be relatively cheap and easy to use. Here we discuss the biophysical potential for the macroscopic detection of viruses based on the induction of a mechanical stress in a bundle of pre-stretched DNA molecules upon binding of viruses to the DNA. We show that the affinity of the DNA to the charged virus surface induces a local melting of the double-helix into two single-stranded DNA. This process effects a mechanical stress along the DNA chains leading to an overall contraction of the DNA. Our results suggest that when such DNA bundles are incorporated in a supporting matrix such as a responsive hydrogel, the presence of viruses may indeed lead to a significant, macroscopic mechanical deformation of the matrix. We discuss the biophysical basis for this effect and characterize the physical properties of the associated DNA melting transition. In particular, we reveal several scaling relations between the relevant physical parameters of the system. We promote this DNA-based assay for efficient and specific virus screening.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Speeding up the first-passage for subdiffusion by introducing a finite potential barrier

We show that for a subdiffusive continuous time random walk with scale-free waiting time distribution the first-passage dynamics on a finite interval can be optimised by introduction of a piecewise linear potential barrier. Analytical results for the survival probability and first-passage density based on the fractional Fokker-Planck equation are shown to agree well with Monte Carlo simulations results. As an application we discus an improved deign for efficient translocation of gradient copolymers compared to homopolymer translocation in a quasi-equilibrium approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic Optimization Based Study of Dimerization Kinetics

We investigate the potential of numerical algorithms to decipher the kinetic parameters involved in multi-step chemical reactions. To this end we study a dimerization kinetics of protein as a model system. We follow the dimerization kinetics using a stochastic simulation algorithm and combine it with three different optimization techniques (Genetic Algorithm, Simulated Annealing and Parallel Tempering) to obtain the rate constants involved in each reaction step. We find good convergence of the numerical scheme to the rate constants of the process. We also perform a sensitivity test on the reaction kinetic parameters to see the relative effects of the parameters for the associated profile of the monomer/dimer distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transient ageing in fractional Brownian and Langevin equation motion

Stochastic processes driven by stationary fractional Gaussian noise, that is, fractional Brownian motion and fractional Langevin equation motion, are usually considered to be ergodic in the sense that, after an algebraic relaxation, time and ensemble averages of physical observables coincide. Recently it was demonstrated that fractional Brownian motion and fractional Langevin equation motion under external confinement are transiently non-ergodic---time and ensemble averages behave differently---from the moment when the particle starts to sense the confinement. Here we show that these processes also exhibit transient ageing, that is, physical observables such as the time averaged mean squared displacement depend on the time lag between the initiation of the system at time $t=0$ and the start of the measurement at the ageing time $t_a$. In particular, it turns out that for fractional Langevin equation motion the ageing dependence on $t_a$ is different between the cases of free and confined motion. We obtain explicit analytical expressions for the aged moments of the particle position as well as the time averaged mean squared displacement and present a numerical analysis of this transient ageing phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transition dynamics in aging systems: microscopic origin of logarithmic time evolution

There exists compelling experimental evidence in numerous systems for logarithmically slow time evolution, yet its theoretical understanding remains elusive. We here introduce and study a generic transition process in complex systems, based on non-renewal, aging waiting times. Each state n of the system follows a local clock initiated at t=0. The random time τbetween clock ticks follows the waiting time density ψ(τ). Transitions between states occur only at local clock ticks and are hence triggered by the local forward waiting time, rather than by ψ(τ). For power-law forms ψ(τ) ~ τ^{-1-α} (0<α<1) we obtain a logarithmic time evolution of the state number <n(t)> ~ log(t/t_0), while for α>2 the process becomes normal in the sense that <n(t)> ~ t. In the intermediate range 1<α<2 we find the power-law growth <n(t)> ~ t^{α-1}. Our model provides a universal description for transition dynamics between aging and non-aging states.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A solution to the subdiffusion-efficiency paradox: Inactive states enhance reaction efficiency at subdiffusion conditions in living cells

Macromolecular crowding in living biological cells effects subdiffusion of larger biomolecules such as proteins and enzymes. Mimicking this subdiffusion in terms of random walks on a critical percolation cluster, we here present a case study of EcoRV restriction enzymes involved in vital cellular defence. We show that due to its so far elusive propensity to an inactive state the enzyme avoids non-specific binding and remains well-distributed in the bulk cytoplasm of the cell. Despite the reduced volume exploration capability of subdiffusion processes, this mechanism guarantees a high efficiency of the enzyme. By variation of the non-specific binding constant and the bond occupation probability on the percolation network, we demonstrate that reduced non-specific binding are beneficial for efficient subdiffusive enzyme activity even in relatively small bacteria cells. Our results corroborate a more local picture of cellular regulation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ageing effects in single particle trajectory averages

We study time averages of single particle trajectories in scale free anomalous diffusion processes, in which the measurement starts at some time t_a>0 after initiation of the process at the time origin, t=0. Using ageing renewal theory we show that for such non-stationary processes a large class of observables are affected by a unique ageing function, which is independent of boundary conditions or the external forces. We quantify the weakly non-ergodic nature of this process in terms of the distribution of time averages and the ergodicity breaking parameter which both explicitly depend on the ageing time t_a. Consequences for the interpretation of single particle tracking data are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins

Combining extensive molecular dynamics simulations of lipid bilayer systems of varying chemical composition with single-trajectory analyses we systematically elucidate the stochastic nature of the lipid motion. We observe subdiffusion over more than four orders of magnitude in time, clearly stretching into the sub-microsecond domain. The lipid motion delicately depends on the lipid chemistry, the lipid phase, and especially on the presence of cholesterol. We demonstrate that fractional Langevin equation motion universally describes the lipid motion in all phases including the gel phase, and in the presence of cholesterol. The results underline the relevance of anomalous diffusion in lipid bilayers and the strong effects of the membrane composition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Area coverage of radial Levy flights with periodic boundary conditions

We consider the time evolution of two-dimensional Levy flights in a finite area with periodic boundary conditions. From simulations we show that the fractal path dimension d_f and thus the degree of area coverage grows in time until it reaches the saturation value d_f=2 at sufficiently long times. We also investigate the time evolution of the probability density function and associated moments in these boundary conditions. Finally we consider the mean first passage time as function of the stable index. Our findings are of interest to assess the ergodic behavior of Levy flights, to estimate their efficiency as stochastic search mechanisms and to discriminate them from other types of search processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bulk-mediated diffusion on a planar surface: full solution

We consider the effective surface motion of a particle that intermittently unbinds from a planar surface and performs bulk excursions. Based on a random walk approach we derive the diffusion equations for surface and bulk diffusion including the surface-bulk coupling. From these exact dynamic equations we analytically obtain the propagator of the effective surface motion. This approach allows us to deduce a superdiffusive, Cauchy-type behavior on the surface, together with exact cutoffs limiting the Cauchy form. Moreover we study the long-time dynamics for the surface motion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-time effects and ultraweak ergodicity breaking in superdiffusive dynamics

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**2 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 trajectory-to-trajectory average of delta**2 as 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

First passages in bounded domains: When is the mean first passage time meaningful?

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(ω) distribution of the random variable ω=τ_1/(τ_1+τ_2), which is a measure for how similar the first passage times τ_1 and τ_2 are of two independent realisations of a Brownian walk starting at the same location. We construct a chart for each domain, determining whether P(ω) represents a unimodal, bell-shaped form, or a bimodal, M-shaped behaviour. While in the former case the mean first passage time (MFPT) is a valid characteristic of the first passage behaviour, in the latter case it is an insufficient measure for the process. Strikingly we find a distinct turnover between the two modes of P(ω), 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 behaviour even in bounded domains, in which all moments of the first passage distribution exist.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How a finite potential barrier decreases the mean first passage time

We consider the mean first passage time of a random walker moving in a potential landscape on a finite interval, 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 minimised by 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Inequivalence of time and ensemble averages in ergodic systems: exponential versus power-law relaxation in confinement

Single particle tracking has become a standard tool to investigate 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The RARE model: a generalized approach to random relaxation processes in disordered systems

This paper introduces and analyses a general statistical model, termed the RARE model, of random relaxation processes in disordered systems. The model considers excitations, that are randomly scattered around a reaction center in a general embedding space. The model's input quantities are the spatial scattering statistics of the excitations around the reaction center, and the chemical reaction rates between the excitations and the reaction center as a function of their mutual distance. The framework of the RARE model is robust, and a detailed stochastic analysis of the random relaxation processes is established. Analytic results regarding the duration and the range of the random relaxation processes, as well as the model's thermodynamic limit, are obtained in closed form. In particular, the case of power-law inputs, which turn out to yield stretched exponential relaxation patterns and asymptotically Paretian relaxation ranges, is addressed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Determining the DNA stability parameters for the breathing dynamics of heterogeneous DNA by stochastic optimization

We suggest that the thermodynamic stability parameters (nearest neighbor stacking and hydrogen bonding free energies) of double-stranded DNA molecules can be inferred reliably from time series of the size fluctuations (breathing) of local denaturation zones (bubbles). On the basis of the reconstructed bubble size distribution, this is achieved through stochastic optimization of the free energies in terms of Simulated Annealing. In particular, it is shown that even noisy time series allow the identification of the stability parameters at remarkable accuracy. This method will be useful to obtain the DNA stacking and hydrogen bonding free energies from single bubble breathing assays rather than equilibrium data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effective surface motion on a reactive cylinder of particles that perform intermittent bulk diffusion

In many biological and small scale technological applications particles may transiently bind to a cylindrical surface. In between two binding events the particles diffuse in the bulk, thus producing an effective translation on the cylinder surface. We here derive the effective motion on the surface, allowing for additional diffusion on the cylinder surface itself. We find explicit solutions for the number of adsorbed particles at one given instant, the effective surface displacement, as well as the surface propagator. In particular sub- and superdiffusive regimes are found, as well as an effective stalling of diffusion visible as a plateau in the mean squared displacement. We also investigate the corresponding first passage and first return problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

In vivo anomalous diffusion and weak ergodicity breaking of lipid granules

Combining extensive single particle tracking microscopy data of endogenous lipid granules in living fission yeast cells with analytical results we show evidence for anomalous diffusion and weak ergodicity breaking. Namely we demonstrate that at short times the granules perform subdiffusion according to the laws of continuous time random walk theory. The associated violation of ergodicity leads to a characteristic turnover between two scaling regimes of the time averaged mean squared displacement. At longer times the granule motion is consistent with fractional Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Residual mean first-passage time for jump processes: theory and applications to Lévy flights and fractional Brownian motion

We derive a functional equation for the mean first-passage time (MFPT) of a generic self-similar Markovian continuous process to a target in a one-dimensional domain and obtain its exact solution. We show that the obtained expression of the MFPT for continuous processes is actually different from the large system size limit of the MFPT for discrete jump processes allowing leapovers. In the case considered here, the asymptotic MFPT admits non-vanishing corrections, which we call residual MFPT. The case of L/'evy flights with diverging variance of jump lengths is investigated in detail, in particular, with respect to the associated leapover behaviour. We also show numerically that our results apply with good accuracy to fractional Brownian motion, despite its non-Markovian nature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The chain sucker: translocation dynamics of a polymer chain into a long narrow channel driven by longitudinal flow

Using analytical techniques and Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a narrow channel of width $R$ embedded in two dimensions, driven by a force proportional to the number of monomers in the channel. Such a setup mimics typical experimental situations in nano/micro-fluidics. During the the translocation process if the monomers in the channel can sufficiently quickly assume steady state motion, we observe the scaling $τ\sim N/F$ of the translocation time $τ$ with the driving force $F$ per bead and the number $N$ of monomers per chain. With smaller channel width $R$, steady state motion cannot be achieved, effecting a non-universal dependence of $τ$ on $N$ and $F$. From the simulations we also deduce the waiting time distributions under various conditions for the single segment passage through the channel entrance. For different chain lengths but the same driving force, the curves of the waiting time as a function of the translocation coordinate $s$ feature a maximum located at identical $s_{\mathrm{max}}$, while with increasing the driving force or the channel width the value of $s_{\mathrm{max}}$ decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anomalous diffusion in correlated continuous time random walks

We demonstrate that continuous time random walks in which successive waiting times are correlated by Gaussian statistics lead to anomalous diffusion with mean squared displacement <r^2(t)>~t^{2/3}. Long-ranged correlations of the waiting times with power-law exponent alpha (0<alpha<=2) give rise to subdiffusion of the form <r^2(t)>~t^{alpha/(1+alpha)}. In contrast correlations in the jump lengths are shown to produce superdiffusion. We show that in both cases weak ergodicity breaking occurs. Our results are in excellent agreement with simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fractional Brownian motion and generalized Langevin equation motion in confined geometries

Motivated by subdiffusive motion of bio-molecules observed in living cells we study the stochastic properties of a non-Brownian particle whose motion is governed by either fractional Brownian motion or the fractional Langevin equation and restricted to a finite domain. We investigate by analytic calculations and simulations how time-averaged observables (e.g., the time averaged mean squared displacement and displacement correlation) are affected by spatial confinement and dimensionality. In particular we study the degree of weak ergodicity breaking and scatter between different single trajectories for this confined motion in the subdiffusive domain. The general trend is that deviations from ergodicity are decreased with decreasing size of the movement volume, and with increasing dimensionality. We define the displacement correlation function and find that this quantity shows distinct features for fractional Brownian motion, fractional Langevin equation, and continuous time subdiffusion, such that it appears an efficient measure to distinguish these different processes based on single particle trajectory data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kramers escape driven by fractional Brownian motion

We investigate the Kramers escape from a potential well of a test particle driven by fractional Gaussian noise with Hurst exponent 0<H<1. From a numerical analysis we demonstrate the exponential distribution of escape times from the well and analyze in detail the dependence of the mean escape time as function of H and the particle diffusivity D. We observe different behavior for the subdiffusive (antipersistent) and superdiffusive (persistent) domains. In particular we find that the escape becomes increasingly faster for decreasing values of H, consistent with previous findings on the first passage behavior. Approximate analytical calculations are shown to support the numerically observed dependencies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polymer translocation into a fluidic channel through a nanopore

Using two dimensional Langevin dynamics simulations, we investigate the dynamics of polymer translocation into a fluidic channel with diameter $R$ through a nanopore under a driving force $F$. Due to the crowding effect induced by the partially translocated monomers, the translocation dynamics is significantly altered in comparison to an unconfined environment, namely, we observe a nonuniversal dependence of the translocation time $τ$ on the chain length $N$. $τ$ initially decreases rapidly and then saturates with increasing $R$, and a dependence of the scaling exponent $α$ of $τ$ with $N$ on the channel width $R$ is observed. The otherwise inverse linear scaling of $τ$ with $F$ breaks down and we observe a minimum of $α$ as a function of $F$. These behaviors are interpreted in terms of the waiting time of an individual segment passing through the pore during translocation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polymer translocation into laterally unbounded confined environments

Using Langevin dynamics simulations in three dimensions (3D), we investigate the dynamics of polymer translocation into the regions between two parallel plane walls with separation $R$ under a driving force $F$, respectively. Compared with an unconfined environment, the translocation dynamics is greatly changed due to the crowding effect of the partially translocated monomers. Translocation time $τ$ initially decreases rapidly with increasing $R$ and then saturates for larger $R$, and the confined environment leads to a nonuniversal dependence of $τ$ on $F$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantitative analysis of single particle trajectories: mean maximal excursion method

An increasing number of experimental studies employ single particle tracking to probe the physical environment in complex systems. We here propose and discuss new methods to analyze the time series of the particle traces, in particular, for subdiffusion phenomena. We discuss the statistical properties of mean maximal excursions, i.e., the maximal distance covered by a test particle up to time t. Compared to traditional methods focusing on the mean squared displacement we show that the mean maximal excursion analysis performs better in the determination of the anomalous diffusion exponent. We also demonstrate that combination of regular moments with moments of the mean maximal excursion method provides additional criteria to determine the exact physical nature of the underlying stochastic subdiffusion processes. We put the methods to test using experimental data as well as simulated time series from different models for normal and anomalous dynamics, such as diffusion on fractals, continuous time random walks, and fractional Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Driven polymer translocation through nanopores: slow versus fast dynamics

We investigate the dynamics of polymer translocation through nanopores under external driving by 3D Langevin Dynamics simulations, focusing on the scaling of the average translocation time $τ$ versus the length of the polymer, $τ\sim N^α$. For slow translocation, i.e., under low driving force and/or high friction, we find $α\approx 1+ν\approx 1.588$ where $ν$ denotes the Flory exponent. In contrast, $α\approx 1.37$ is observed for fast translocation due to the highly deformed chain conformation on the trans side, reflecting a pronounced non-equilibrium situation. The dependence of the translocation time on the driving force is given by $τ\sim F^{-1}$ and $τ\sim F^{-0.80}$ for slow and fast translocation, respectively. These results clarify the controversy on the magnitude of the scaling exponent $α$ for driven translocation.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Bulk-mediated surface diffusion on a cylinder: propagators and crossovers

We consider the effective surface motion of a particle that freely diffuses in the bulk and intermittently binds to that surface. From an exact approach we derive various regimes of the effective surface motion characterized by physical rates for binding/unbinding and the bulk diffusivity. We obtain a transient regime of superdiffusion and, in particular, a saturation regime characteristic for the cylindrical geometry. This saturation, however, in a finite system is not terminal but eventually turns over to normal surface diffusion. The first passage behavior of particles to the cylinder surface is derived. Consequences for actual systems are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

DNA bubble dynamics as a quantum Coulomb problem

We study the dynamics of denaturation bubbles in double-stranded DNA on the basis of the Poland-Scheraga model. We demonstrate that the associated Fokker-Planck equation is equivalent to a Coulomb problem. Below the melting temperature the bubble lifetime is associated with the continuum of scattering states of the repulsive Coulomb potential, at the melting temperature the Coulomb potential vanishes and the underlying first exit dynamics exhibits a long time power law tail, above the melting temperature, corresponding to an attractive Coulomb potential, the long time dynamics is controlled by the lowest bound state. Correlations and finite size effects are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Master equation approach to DNA-breathing in heteropolymer DNA

After crossing an initial barrier to break the first base-pair (bp) in double-stranded DNA, the disruption of further bps is characterized by free energies between less than one to a few kT. This causes the opening of intermittent single-stranded bubbles. Their unzipping and zipping dynamics can be monitored by single molecule fluorescence or NMR methods. We here establish a dynamic description of this DNA-breathing in a heteropolymer DNA in terms of a master equation that governs the time evolution of the joint probability distribution for the bubble size and position along the sequence. The transfer coefficients are based on the Poland-Scheraga free energy model. We derive the autocorrelation function for the bubble dynamics and the associated relaxation time spectrum. In particular, we show how one can obtain the probability densities of individual bubble lifetimes and of the waiting times between successive bubble events from the master equation. A comparison to results of a stochastic Gillespie simulation shows excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Chaperone assisted translocation

We investigate the translocation of a stiff polymer through a nanopore in a membrane, in the presence of binding particles (chaperones) that bind reversibly to the polymer on both sides of the membrane. A bound chaperone covers one (univalent binding) or many (multivalent binding) binding sites. Assuming that the diffusion of the chaperones is fast compared to the rate of translocation we describe the process by a one-dimensional master equation. We expand previous models by a detailed study of the effective force in the master equation, which is obtained by the appropriate statistical mechanical average over the chaperone states. The dependence of the force on the degree of valency (the number of binding sites occupied by a chaperone) is studied in detail. We obtain finite size corrections (to the thermodynamical expression for the force), which, for univalent binding, can be expressed analytically. We finally investigate the mean velocity for translocation as a function of chaperone binding strength and size. For both univalent and multivalent binding simple results are obtained for the case of a sufficiently long translocating polymer.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Comment on "Why is the DNA denaturation transition first order?"

In this comment we argue that while the conclusions in the original paper (Y. Kafri, D. Mukamel and L. Peliti, Phys. Rev. Lett. 85, 4988 (2000)) are correct for asymptotically long DNA chains, they do not apply to the chains used in typical experiments. In the added last paragraph, we point out that for real DNA the average distance between denatured loops is not of the order of the persistence length of a single-stranded chain but much larger. This corroborates our reasoning that the double helix between loops is quite rigid, and thereby our conclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Exact solution of a linear molecular motor model driven by two-step fluctuations and subject to protein friction

We investigate by analytical means the stochastic equations of motion of a linear molecular motor model based on the concept of protein friction. Solving the coupled Langevin equations originally proposed by Mogilner et al. (A. Mogilner et al., Phys. Lett. {\bf 237}, 297 (1998)), and averaging over both the two-step internal conformational fluctuations and the thermal noise, we present explicit, analytical expressions for the average motion and the velocity-force relationship. Our results allow for a direct interpretation of details of this motor model which are not readily accessible from numerical solutions. In particular, we find that the model is able to predict physiologically reasonable values for the load-free motor velocity and the motor mobility.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Physical Pictures of Transport in Heterogeneous Media: Advection-Dispersion, Random Walk and Fractional Derivative Formulations

The basic conceptual picture and theoretical basis for development of transport equations in porous media are examined. The general form of the governing equations is derived for conservative chemical transport in heterogeneous geological formations, for single realizations and for ensemble averages of the domain. The application of these transport equations is focused on accounting for the appearance of non-Fickian (anomalous) transport behavior. The general ensemble-averaged transport equation is shown to be equivalent to a continuous time random walk (CTRW) and reduces to the conventional forms of the advection-dispersion equation (ADE) under highly restrictive conditions. Fractional derivative formulations of the transport equations, both temporal and spatial, emerge as special cases of the CTRW. In particular, the use in this context of L{é}vy flights is critically examined. In order to determine chemical transport in field-scale situations, the CTRW approach is generalized to non-stationary systems. We outline a practical numerical scheme, similar to those used with extended geological models, to account for the often important effects of unresolved heterogeneities.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Equilibrium shapes of flat knots

We study the equilibrium shapes of prime and composite knots confined to two dimensions. Using rigorous scaling arguments we show that, due to self-avoiding effects, the topological details of prime knots are localised on a small portion of the larger ring polymer. Within this region, the original knot configuration can assume a hierarchy of contracted shapes, the dominating one given by just one small loop. This hierarchy is investigated in detail for the flat trefoil knot, and corroborated by Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Fractional Klein-Kramers equation for superdiffusive transport: normal versus anomalous time evolution in a differential L{é}vy walk model

We introduce a fractional Klein-Kramers equation which describes sub-ballistic superdiffusion in phase space in the presence of a space-dependent external force field. This equation defines the differential L{é}vy walk model whose solution is shown to be non-negative. In the velocity coordinate, the probability density relaxes in Mittag-Leffler fashion towards the Maxwell distribution whereas in the space coordinate, no stationary solution exists and the temporal evolution of moments exhibits a competition between Brownian and anomalous contributions.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The future is noisy: on the role of spatial fluctuations in genetic switching

A genetic switch may be realised by a certain operator sector on the DNA strand from which either genetic code, to the left or to the right of this operator sector, can be transcribed and the corresponding information processed. This switch is controlled by messenger molecules, i.e., they determine to which side the switch is flipped. Recently, it has been realised that noise plays an elementary role in genetic switching, and the effect of number fluctuations of the messenger molecules have been explored. Here we argue that the assumption of well-stirredness taken in the previous models may not be sufficient to characterise the influence of noise: {\em spatial} fluctuations play a non-negligible part in cellular genetic switching processes.

0 citations0 reviewsNo code linkedOpen access