Refine
Has Fulltext
- no (22)
Year of publication
- 2020 (22) (remove)
Document Type
- Article (22) (remove)
Language
- English (22)
Is part of the Bibliography
- yes (22)
Keywords
- diffusion (7)
- anomalous diffusion (3)
- first-passage (2)
- fractional Brownian motion (2)
- random diffusivity (2)
- Absorption (1)
- Antibiotics (1)
- Bacterial biofilms (1)
- Biofilms (1)
- Biological defense mechanisms (1)
Institute
Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent a. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent alpha is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.
The escape from a potential well is an archetypal problem in the study of stochastic dynamical systems, representing real-world situations from chemical reactions to leaving an established home range in movement ecology. Concurrently, Levy noise is a well-established approach to model systems characterized by statistical outliers and diverging higher order moments, ranging from gene expression control to the movement patterns of animals and humans. Here, we study the problem of Levy noise-driven escape from an almost rectangular, arctangent potential well restricted by two absorbing boundaries, mostly under the action of the Cauchy noise. We unveil analogies of the observed transient dynamics to the general properties of stationary states of Levy processes in single-well potentials. The first-escape dynamics is shown to exhibit exponential tails. We examine the dependence of the escape on the shape parameters, steepness, and height of the arctangent potential. Finally, we explore in detail the behavior of the probability densities of the first-escape time and the last-hitting point.
This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D ₀), as well as a low (D ₘ) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the Itô calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D ₘ will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D ₘ, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed.
Diffusion of antibiotics through a biofilm in the presence of diffusion and absorption barriers
(2020)
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 subdiffusion and absorption parameters of the biofilm change over time. We propose a four-stage model of antibiotic diffusion in biofilm based on the following physical characteristics: whether there is absorption of the antibiotic in the biofilm and whether all biofilm parameters remain unchanged over time. The biological interpretation of the stages, in particular their relation with the bacterial defense 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.
Stochastic models based on random diffusivities, such as the diffusing-diffusivity approach, are popular concepts for the description of non-Gaussian diffusion in heterogeneous media. Studies of these models typically focus on the moments and the displacement probability density function. Here we develop the complementary power spectral description for a broad class of random-diffusivity processes. In our approach we cater for typical single particle tracking data in which a small number of trajectories with finite duration are garnered. Apart from the diffusing-diffusivity model we study a range of previously unconsidered random-diffusivity processes, for which we obtain exact forms of the probability density function. These new processes are different versions of jump processes as well as functionals of Brownian motion. The resulting behaviour subtly depends on the specific model details. Thus, the central part of the probability density function may be Gaussian or non-Gaussian, and the tails may assume Gaussian, exponential, log-normal, or even power-law forms. For all these models we derive analytically the moment-generating function for the single-trajectory power spectral density. We establish the generic 1/f²-scaling of the power spectral density as function of frequency in all cases. Moreover, we establish the probability density for the amplitudes of the random power spectral density of individual trajectories. The latter functions reflect the very specific properties of the different random-diffusivity models considered here. Our exact results are in excellent agreement with extensive numerical simulations.
We consider the first-passage problem for N identical independent particles that are initially released uniformly in a finite domain Ω and then diffuse toward a reactive area Γ, which can be part of the outer boundary of Ω or a reaction centre in the interior of Ω. For both cases of perfect and partial reactions, we obtain the explicit formulas for the first two moments of the fastest first-passage time (fFPT), i.e., the time when the first out of the N particles reacts with Γ. Moreover, we investigate the full probability density of the fFPT. We discuss a significant role of the initial condition in the scaling of the average fFPT with the particle number N, namely, a much stronger dependence (1/N and 1/N² for partially and perfectly reactive targets, respectively), in contrast to the well known inverse-logarithmic behaviour found when all particles are released from the same fixed point. We combine analytic solutions with scaling arguments and stochastic simulations to rationalise our results, which open new perspectives for studying the relevance of multiple searchers in various situations of molecular reactions, in particular, in living cells.
We study the extremal properties of a stochastic process xt defined by the Langevin equation ẋₜ =√2Dₜ ξₜ, in which ξt is a Gaussian white noise with zero mean and Dₜ is a stochastic‘diffusivity’, defined as a functional of independent Brownian motion Bₜ.We focus on threechoices for the random diffusivity Dₜ: cut-off Brownian motion, Dₜt ∼ Θ(Bₜ), where Θ(x) is the Heaviside step function; geometric Brownian motion, Dₜ ∼ exp(−Bₜ); and a superdiffusive process based on squared Brownian motion, Dₜ ∼ B²ₜ. For these cases we derive exact expressions for the probability density functions of the maximal positive displacement and of the range of the process xₜ on the time interval ₜ ∈ (0, T).We discuss the asymptotic behaviours of the associated probability density functions, compare these against the behaviour of the corresponding properties of standard Brownian motion with constant diffusivity (Dₜ = D0) and also analyse the typical behaviour of the probability density functions which is observed for a majority of realisations of the stochastic diffusivity process.
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 a priori unexpected dependencies of the effective diffusion coefficient on the correlation exponent. Our observations of the 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.
What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.
Transition path dynamics have been widely studied in chemical, physical, and technological systems. Mostly, the transition path dynamics is obtained for smooth barrier potentials, for instance, generic inverse-parabolic shapes. We here present analytical results for the mean transition path time, the distribution of transition path times, the mean transition path velocity, and the mean transition path shape in a rough inverted parabolic potential function under the driving of Gaussian white noise. These are validated against extensive simulations using the forward flux sampling scheme in parallel computations. We observe how precisely the potential roughness, the barrier height, and the noise intensity contribute to the particle transition in the rough inverted barrier potential.