Refine
Year of publication
Language
- English (31)
Is part of the Bibliography
- yes (31)
Keywords
- anomalous diffusion (5)
- Levy flights (2)
- diffusion (2)
- fractional Brownian motion (2)
- non-Gaussianity (2)
- stochastic processes (2)
- Anomalous diffusion (1)
- Bulk-mediated diffusion (1)
- Bulk-mediated diffusion; (1)
- Comb model (1)
Institute
We consider the diffusion-advection problem in two simple cellular flow models ( often invoked as examples of subdiffusive tracer motion) and concentrate on the intermediate time range, in which the tracer motion indeed may show subdiffusion. We perform extensive numerical simulations of the systems under different initial conditions and show that the pure intermediate-time subdiffusion regime is only evident when the particles start at the border between different cells, i.e., at the separatrix, and is less pronounced or absent for other initial conditions. The motion moreover shows quite peculiar aging properties, which are also mirrored in the behavior of the time-averaged mean squared displacement for single trajectories. This kind of behavior is due to the complex motion of tracers trapped inside the cell and is absent in classical models based on continuous-time random walks with no dynamics in the trapped state.
Spot variation fluorescence correlation spectroscopy (SV-FCS) is a variant of the FCS techniques which may give useful information about the structural organisation of the medium in which the diffusion takes place. We show that the same results can be obtained by post-processing the photon count data from ordinary FCS measurements. By using this method, one obtains the fluorescence autocorrelation functions for sizes of confocal volume, which are effectively smaller than that of the initial FCS measurement. The photon counts of the initial experiment are first transformed into smooth intensity trace using kernel smoothing method or to a piecewise-continuous intensity trace using binning and then a non-linear transformation is applied to this trace. The result of this transformation mimics the photon count rate in an experiment performed with a smaller confocal volume. The applicability of the method is established in extensive numerical simulations and directly supported in in-vitro experiments. The procedure is then applied to the diffusion of AlexaFluor647-labeled streptavidin in living cells.
The diffusion in the comb structures is a popular model of geometrically induced anomalous diffusion. In the present work we concentrate on the diffusion along the backbone in a system where sidebranches are planes, and the diffusion thereon is anomalous and described by continuous time random walks (CTRW). We show that the mean squared displacement (MSD) in the backbone of the comb behaves differently depending on whether the waiting time periods in the sidebranches are reset after the step in the backbone is done (a rejuvenating junction model), or not (a non-rejuvenating junction model). In the rejuvenating case the subdiffusion in the sidebranches only changes the prefactor in the ultra-slow (logarithmic) diffusion along the backbone, while in the non-rejuvenating case the ultraslow, logarithmic subdiffusion is changed to a much faster power-law subdiffusion (with a logarithmic correction) as it was found earlier by Iomin and Mendez [25]. Moreover, in the first case the result does not change if the diffusion in the backbone is itself anomalous, while in the second case it does. Two of the special cases of the considered models (the non-rejuvenating junction under normal diffusion in the backbone, and rejuvenating junction for the same waiting time distribution in the sidebranches and in junction points) were also investigated within the approach based on the corresponding generalized Fokker-Planck equations. (c) 2017 Elsevier Ltd. All rights reserved.
How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x) = D0|x|gamma and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicitybreaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB similar to(1/r )-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics.
Lennard-Jones mixtures represent one of the popular systems for the study of glass-forming liquids.
Spatio/temporal heterogeneity and rare (activated) events are at the heart of the slow dynamics typical of these systems. Such slow dynamics is characterised by the development of a plateau in the mean-squared displacement (MSD) at intermediate times, accompanied by a non-Gaussianity in the displacement distribution identified by exponential tails.
As pointed out by some recent works, the non-Gaussianity persists at times beyond the MSD plateau, leading to a Brownian yet non-Gaussian regime and thus highlighting once again the relevance of rare events in such systems.
Single-particle motion of glass-forming liquids is usually interpreted as an alternation of rattling within the local cage and cage-escape motion and therefore can be described as a sequence of waiting times and jumps. In this work, by using a simple yet robust algorithm, we extract jumps and waiting times from single-particle trajectories obtained via molecular dynamics simulations.
We investigate the presence of correlations between waiting times and find negative correlations, which becomes more and more pronounced when lowering the temperature.
We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed.
How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM.
The relaxation of a dissipative system to its equilibrium state often shows a multiexponential pattern with relaxation rates, which are typically considered to be independent of the initial condition. The rates follow from the spectrum of a Hermitian operator obtained by a similarity transformation of the initial Fokker-Planck operator. However, some initial conditions are mapped by this similarity transformation to functions which growat infinity. These cannot be expanded in terms of the eigenfunctions of a Hermitian operator, and show different relaxation patterns. Considering the exactly solvable examples of Gaussian and generalized Levy Ornstein-Uhlenbeck processes (OUPs) we show that the relaxation rates belong to the Hermitian spectrum only if the initial condition belongs to the domain of attraction of the stable distribution defining the noise. While for an ordinary OUP initial conditions leading to nonspectral relaxation can be considered exotic, for generalized OUPs driven by Levy noise, such initial conditions are the rule. DOI: 10.1103/PhysRevLett.110.150602
We present a detailed analysis of the eigenfunctions of the Fokker-Planck operator for the LevyOrnstein- Uhlenbeck process, their asymptotic behavior and recurrence relations, explicit expressions in coordinate space for the special cases of the Ornstein-Uhlenbeck process with Gaussian and with Cauchy white noise and for the transformation kernel, which maps the fractional Fokker-Planck operator of the Cauchy-Ornstein-Uhlenbeck process to the non-fractional Fokker-Planck operator of the usual Gaussian Ornstein-Uhlenbeck process. We also describe how non-spectral relaxation can be observed in bounded random variables of the Levy-Ornstein-Uhlenbeck process and their correlation functions.