Refine
Year of publication
Document Type
- Article (40)
- Postprint (16)
- Doctoral Thesis (2)
- Review (1)
Language
- English (59)
Is part of the Bibliography
- yes (59) (remove)
Keywords
- anomalous diffusion (59) (remove)
Institute
Fractional Brownian motion in superharmonic potentials and non-Boltzmann stationary distributions
(2021)
We study the stochastic motion of particles driven by long-range correlated fractional Gaussian noise (FGN) in a superharmonic external potential of the form U(x) proportional to x(2n) (n is an element of N). When the noise is considered to be external, the resulting overdamped motion is described by the non-Markovian Langevin equation for fractional Brownian motion. For this case we show the existence of long time, stationary probability density functions (PDFs) the shape of which strongly deviates from the naively expected Boltzmann PDF in the confining potential U(x). We analyse in detail the temporal approach to stationarity as well as the shape of the non-Boltzmann stationary PDF. A typical characteristic is that subdiffusive, antipersistent (with negative autocorrelation) motion tends to effect an accumulation of probability close to the origin as compared to the corresponding Boltzmann distribution while the opposite trend occurs for superdiffusive (persistent) motion. For this latter case this leads to distinct bimodal shapes of the PDF. This property is compared to a similar phenomenon observed for Markovian Levy flights in superharmonic potentials. We also demonstrate that the motion encoded in the fractional Langevin equation driven by FGN always relaxes to the Boltzmann distribution, as in this case the fluctuation-dissipation theorem is fulfilled.
We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) < x(2)(t)> similar or equal to log(gamma)t. Comparison of annealed (renewal) continuous time random walks (CTRWs) with logarithmic waiting time distribution psi(tau) similar or equal to 1/(tau log(1+gamma)tau) and Sinai diffusion in quenched random landscapes reveals striking similarities, despite the great differences in their physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time-averaged and ensemble-averaged MSDs. Remarkably, for the CTRW we observe that the fluctuations of time averages become universal, with an exponential suppression of mobile trajectories. We discuss the fundamental connection between the Golosov localization effect and non-ergodicity in the sense of the disparity between ensemble-averaged MSD and time-averaged MSD.
Based on extensive Brownian dynamics simulations we study the thermal motion of a tracer bead in a cross-linked, flexible gel in the limit when the tracer particle size is comparable to or even larger than the equilibrium mesh size of the gel. The analysis of long individual trajectories of the tracer demonstrates the existence of pronounced transient anomalous diffusion. From the time averaged mean squared displacement and the time averaged van Hove correlation functions we elucidate the many-body origin of the non-Brownian tracer bead dynamics. Our results shed new light onto the ongoing debate over the physical origin of steric tracer interactions with structured environments.
We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer–obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer–obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer–crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids.
We study the diffusion of a tracer particle, which moves in continuum space between a lattice of excluded volume, immobile non-inert obstacles. In particular, we analyse how the strength of the tracer–obstacle interactions and the volume occupancy of the crowders alter the diffusive motion of the tracer. From the details of partitioning of the tracer diffusion modes between trapping states when bound to obstacles and bulk diffusion, we examine the degree of localisation of the tracer in the lattice of crowders. We study the properties of the tracer diffusion in terms of the ensemble and time averaged mean squared displacements, the trapping time distributions, the amplitude variation of the time averaged mean squared displacements, and the non-Gaussianity parameter of the diffusing tracer. We conclude that tracer–obstacle adsorption and binding triggers a transient anomalous diffusion. From a very narrow spread of recorded individual time averaged trajectories we exclude continuous type random walk processes as the underlying physical model of the tracer diffusion in our system. For moderate tracer–crowder attraction the motion is found to be fully ergodic, while at stronger attraction strength a transient disparity between ensemble and time averaged mean squared displacements occurs. We also put our results into perspective with findings from experimental single-particle tracking and simulations of the diffusion of tagged tracers in dense crowded suspensions. Our results have implications for the diffusion, transport, and spreading of chemical components in highly crowded environments inside living cells and other structured liquids.
A topic of intense current investigation pursues the question of how the highly crowded environment of biological cells affects the dynamic properties of passively diffusing particles. Motivated by recent experiments we report results of extensive simulations of the motion of a finite sized tracer particle in a heterogeneously crowded environment made up of quenched distributions of monodisperse crowders of varying sizes in finite circular two-dimensional domains. For given spatial distributions of monodisperse crowders we demonstrate how anomalous diffusion with strongly non-Gaussian features arises in this model system. We investigate both biologically relevant situations of particles released either at the surface of an inner domain or at the outer boundary, exhibiting distinctly different features of the observed anomalous diffusion for heterogeneous distributions of crowders. Specifically we reveal an asymmetric spreading of tracers even at moderate crowding. In addition to the mean squared displacement (MSD) and local diffusion exponent we investigate the magnitude and the amplitude scatter of the time averaged MSD of individual tracer trajectories, the non-Gaussianity parameter, and the van Hove correlation function. We also quantify how the average tracer diffusivity varies with the position in the domain with a heterogeneous radial distribution of crowders and examine the behaviour of the survival probability and the dynamics of the tracer survival probability. Inter alia, the systems we investigate are related to the passive transport of lipid molecules and proteins in two-dimensional crowded membranes or the motion in colloidal solutions or emulsions in effectively two-dimensional geometries, as well as inside supercrowded, surface adhered cells.
The mobile-immobile model (MIM) has been established in geoscience in the context of contaminant transport in groundwater. Here the tracer particles effectively immobilise, e.g., due to diffusion into dead-end pores or sorption. The main idea of the MIM is to split the total particle density into a mobile and an immobile density. Individual tracers switch between the mobile and immobile state following a two-state telegraph process, i.e., the residence times in each state are distributed exponentially. In geoscience the focus lies on the breakthrough curve (BTC), which is the concentration at a fixed location over time. We apply the MIM to biological experiments with a special focus on anomalous scaling regimes of the mean squared displacement (MSD) and non-Gaussian displacement distributions. As an exemplary system, we have analysed the motion of tau proteins, that diffuse freely inside axons of neurons. Their free diffusion thereby corresponds to the mobile state of the MIM. Tau proteins stochastically bind to microtubules, which effectively immobilises the tau proteins until they unbind and continue diffusing. Long immobilisation durations compared to the mobile durations give rise to distinct non-Gaussian Laplace shaped distributions. It is accompanied by a plateau in the MSD for initially mobile tracer particles at relevant intermediate timescales. An equilibrium fraction of initially mobile tracers gives rise to non-Gaussian displacements at intermediate timescales, while the MSD remains linear at all times. In another setting bio molecules diffuse in a biosensor and transiently bind to specific receptors, where advection becomes relevant in the mobile state. The plateau in the MSD observed for the advection-free setting and long immobilisation durations persists also for the case with advection. We find a new clear regime of anomalous diffusion with non-Gaussian distributions and a cubic scaling of the MSD. This regime emerges for initially mobile and for initially immobile tracers. For an equilibrium fraction of initially mobile tracers we observe an intermittent ballistic scaling of the MSD. The long-time effective diffusion coefficient is enhanced by advection, which we physically explain with the variance of mobile durations. Finally, we generalize the MIM to incorporate arbitrary immobilisation time distributions and focus on a Mittag-Leffler immobilisation time distribution with power-law tail ~ t^(-1-mu) with 0<mu<1 and diverging mean immobilisation durations. A fit of our model to the BTC of experimental data from tracer particles in aquifers matches the BTC including the power-law tail. We use the fit parameters for plotting the displacement distributions and the MSD. We find Gaussian normal diffusion at short times and long-time power-law decay of mobile mass accompanied by anomalous diffusion at long times. The long-time diffusion is subdiffusive in the advection-free setting, while it is either subdiffusive for 0<mu<1/2 or superdiffusive for 1/2<mu<1 when advection is present. In the long-time limit we show equivalence of our model to a bi-fractional diffusion equation.
Fluctuation relations for anomalous dynamics generated by time-fractional Fokker-Planck equations
(2015)
Anomalous dynamics characterized by non-Gaussian probability distributions (PDFs) and/or temporal long-range correlations can cause subtle modifications of conventional fluctuation relations (FRs). As prototypes we study three variants of a generic time-fractional Fokker-Planck equation with constant force. Type A generates superdiffusion, type B subdiffusion and type C both super-and subdiffusion depending on parameter variation. Furthermore type C obeys a fluctuation-dissipation relation whereas A and B do not. We calculate analytically the position PDFs for all three cases and explore numerically their strongly non-Gaussian shapes. While for type C we obtain the conventional transient work FR, type A and type B both yield deviations by featuring a coefficient that depends on time and by a nonlinear dependence on the work. We discuss possible applications of these types of dynamics and FRs to experiments.
We investigate a diffusion process with a time-dependent diffusion coefficient, both exponentially increasing and decreasing in time, D(t)=D-0(e +/- 2 alpha t). For this (hypothetical) nonstationary diffusion process we compute-both analytically and from extensive stochastic simulations-the behavior of the ensemble- and time-averaged mean-squared displacements (MSDs) of the particles, both in the over- and underdamped limits. Simple asymptotic relations derived for the short- and long-time behaviors are shown to be in excellent agreement with the results of simulations. The diffusive characteristics in the presence of ageing are also considered, with dramatic differences of the over- versus underdamped regime. Our results for D(t)=D-0(e +/- 2 alpha t) extend and generalize the class of diffusive systems obeying scaled Brownian motion featuring a power-law-like variation of the diffusivity with time, D(t) similar to t(alpha-1). We also examine the logarithmically increasing diffusivity, D(t)=D(0)log[t/tau(0)], as another fundamental functional dependence (in addition to the power-law and exponential) and as an example of diffusivity slowly varying in time. One of the main conclusions is that the behavior of the massive particles is predominantly ergodic, while weak ergodicity breaking is repeatedly found for the time-dependent diffusion of the massless particles at short times. The latter manifests itself in the nonequivalence of the (both nonaged and aged) MSD and the mean time-averaged MSD. The current findings are potentially applicable to a class of physical systems out of thermal equilibrium where a rapid increase or decrease of the particles' diffusivity is inherently realized. One biological system potentially featuring all three types of time-dependent diffusion (power-law-like, exponential, and logarithmic) is water diffusion in the brain tissues, as we thoroughly discuss in the end.
We consider diffusion processes with a spatially varying diffusivity giving rise to anomalous diffusion. Such heterogeneous diffusion processes are analysed for the cases of exponential, power-law, and logarithmic dependencies of the diffusion coefficient on the particle position. Combining analytical approaches with stochastic simulations, we show that the functional form of the space-dependent diffusion coefficient and the initial conditions of the diffusing particles are vital for their statistical and ergodic properties. In all three cases a weak ergodicity breaking between the time and ensemble averaged mean squared displacements is observed. We also demonstrate a population splitting of the time averaged traces into fast and slow diffusers for the case of exponential variation of the diffusivity as well as a particle trapping in the case of the logarithmic diffusivity. Our analysis is complemented by the quantitative study of the space coverage, the diffusive spreading of the probability density, as well as the survival probability.
We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.
We study the thermal Markovian diffusion of tracer particles in a 2D medium with spatially varying diffusivity D(r), mimicking recently measured, heterogeneous maps of the apparent diffusion coefficient in biological cells. For this heterogeneous diffusion process (HDP) we analyse the mean squared displacement (MSD) of the tracer particles, the time averaged MSD, the spatial probability density function, and the first passage time dynamics from the cell boundary to the nucleus. Moreover we examine the non-ergodic properties of this process which are important for the correct physical interpretation of time averages of observables obtained from single particle tracking experiments. From extensive computer simulations of the 2D stochastic Langevin equation we present an in-depth study of this HDP. In particular, we find that the MSDs along the radial and azimuthal directions in a circular domain obey anomalous and Brownian scaling, respectively. We demonstrate that the time averaged MSD stays linear as a function of the lag time and the system thus reveals a weak ergodicity breaking. Our results will enable one to rationalise the diffusive motion of larger tracer particles such as viruses or submicron beads in biological cells.
We study the effects of ageing-the time delay between initiation of the physical process at t = 0 and start of observation at some time t(a) > 0-and spatial confinement on the properties of heterogeneous diffusion processes (HDPs) with deterministic power-law space-dependent diffusivities, D(x) = D-0 vertical bar x vertical bar(alpha). From analysis of the ensemble and time averaged mean squared displacements and the ergodicity breaking parameter quantifying the inherent degree of irreproducibility of individual realizations of the HDP we obtain striking similarities to ageing subdiffusive continuous time random walks with scale-free waiting time distributions. We also explore how both processes can be distinguished. For confined HDPs we study the long-time saturation of the ensemble and time averaged particle displacements as well as the magnitude of the inherent scatter of time averaged displacements and contrast the outcomes to the results known for other anomalous diffusion processes under confinement.
In various biological systems and small scale technological applications particles transiently bind to a cylindrical surface. Upon unbinding the particles diffuse in the vicinal bulk before rebinding to the surface. Such bulk-mediated excursions give rise to an effective surface translation, for which we here derive and discuss the dynamic equations, including additional surface diffusion. We discuss the time evolution of the number of surface-bound particles, the effective surface mean squared displacement, and the surface propagator. In particular, we observe sub- and superdiffusive regimes. A plateau of the surface mean-squared displacement reflects a stalling of the surface diffusion at longer times. Finally, the corresponding first passage problem for the cylindrical geometry is analysed.
In various biological systems and small scale technological applications particles transiently bind to a cylindrical surface. Upon unbinding the particles diffuse in the vicinal bulk before rebinding to the surface. Such bulk-mediated excursions give rise to an effective surface translation, for which we here derive and discuss the dynamic equations, including additional surface diffusion. We discuss the time evolution of the number of surface-bound particles, the effective surface mean squared displacement, and the surface propagator. In particular, we observe sub- and superdiffusive regimes. A plateau of the surface mean-squared displacement reflects a stalling of the surface diffusion at longer times. Finally, the corresponding first passage problem for the cylindrical geometry is analysed.
We define and study in detail utraslow scaled Brownian motion (USBM) characterized by a time dependent diffusion coefficient of the form D(t) similar or equal to 1/t. For unconfined motion the mean squared displacement (MSD) of USBM exhibits an ultraslow, logarithmic growth as function of time, in contrast to the conventional scaled Brownian motion. In a harmonic potential the MSD of USBM does not saturate but asymptotically decays inverse-proportionally to time, reflecting the highly non-stationary character of the process. We show that the process is weakly non-ergodic in the sense that the time averaged MSD does not converge to the regular MSD even at long times, and for unconfined motion combines a linear lag time dependence with a logarithmic term. The weakly non-ergodic behaviour is quantified in terms of the ergodicity breaking parameter. The USBM process is also shown to be ageing: observables of the system depend on the time gap between initiation of the test particle and start of the measurement of its motion. Our analytical results are shown to agree excellently with extensive computer simulations.
Diffusion of finite-size particles in two-dimensional channels with random wall configurations
(2014)
Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick–Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].
Diffusion of finite-size particles in two-dimensional channels with random wall configurations
(2014)
Diffusion of chemicals or tracer molecules through complex systems containing irregularly shaped channels is important in many applications. Most theoretical studies based on the famed Fick–Jacobs equation focus on the idealised case of infinitely small particles and reflecting boundaries. In this study we use numerical simulations to consider the transport of finite-size particles through asymmetrical two-dimensional channels. Additionally, we examine transient binding of the molecules to the channel walls by applying sticky boundary conditions. We consider an ensemble of particles diffusing in independent channels, which are characterised by common structural parameters. We compare our results for the long-time effective diffusion coefficient with a recent theoretical formula obtained by Dagdug and Pineda [J. Chem. Phys., 2012, 137, 024107].
The Cattaneo or telegrapher's equation describes the crossover from initial ballistic to normal diffusion. Here we study and survey time-fractional generalisations of this equation that are shown to produce the crossover of the mean squared displacement from superdiffusion to subdiffusion. Conditional solutions are derived in terms of Fox H-functions and the dth-order moments as well as the diffusive flux of the different models are derived. Moreover, the concept of the distribution-like is proposed as an alternative to the probability density function.