530 Physik
Refine
Has Fulltext
- no (106) (remove)
Year of publication
Document Type
- Article (106) (remove)
Language
- English (106)
Is part of the Bibliography
- yes (106)
Keywords
- diffusion (21)
- anomalous diffusion (18)
- fractional Brownian motion (4)
- Fokker-Planck equation (3)
- first passage time (3)
- first-passage (3)
- first-passage time (3)
- geometric Brownian motion (3)
- nonergodicity (3)
- stochastic processes (3)
Institute
What are the physical laws of the mutual interactions of objects bound to cell membranes, such as various membrane proteins or elongated virus particles? To rationalise this, we here investigate by extensive computer simulations mutual interactions of rod-like particles adsorbed on the surface of responsive elastic two-dimensional sheets. Specifically, we quantify sheet deformations as a response to adhesion of such filamentous particles. We demonstrate that tip-to-tip contacts of rods are favoured for relatively soft sheets, while side-by-side contacts are preferred for stiffer elastic substrates. These attractive orientation-dependent substrate-mediated interactions between the rod-like particles on responsive sheets can drive their aggregation and self-assembly. The optimal orientation of the membrane-bound rods is established via responding to the elastic energy profiles created around the particles. We unveil the phase diagramme of attractive–repulsive rod–rod interactions in the plane of their separation and mutual orientation. Applications of our results to other systems featuring membrane-associated particles are also discussed.
Levy walks are continuous-time random-walk processes with a spatiotemporal coupling of jump lengths and waiting times. We here apply the Hermite polynomial method to study the behavior of LWs with power-law walking time density for four different cases. First we show that the known result for the infinite density of an unconfined, unbiased LW is consistently recovered. We then derive the asymptotic behavior of the probability density function (PDF) for LWs in a constant force field, and we obtain the corresponding qth-order moments. In a harmonic external potential we derive the relaxation dynamic of the LW. For the case of a Poissonian walking time an exponential relaxation behavior is shown to emerge. Conversely, a power-law decay is obtained when the mean walking time diverges. Finally, we consider the case of an unconfined, unbiased LW with decaying speed v(r ) = v0/./r. When the mean walking time is finite, a universal Gaussian law for the position-PDF of the walker is obtained explicitly.
How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM.
We study a heterogeneous diffusion process (HDP) with position-dependent diffusion coefficient and Poissonian stochastic resetting.
We find exact results for the mean squared displacement and the probability density function. The nonequilibrium steady state reached in the long time limit is studied.
We also analyse the transition to the non-equilibrium steady state by finding the large deviation function.
We found that similarly to the case of the normal diffusion process where the diffusion length grows like t (1/2) while the length scale xi(t) of the inner core region of the nonequilibrium steady state grows linearly with time t, in the HDP with diffusion length increasing like t ( p/2) the length scale xi(t) grows like t ( p ).
The obtained results are verified by numerical solutions of the corresponding Langevin equation.
Heterogeneous diffusion processes (HDPs) with space-dependent diffusion coefficients D(x) are found in a number of real-world systems, such as for diffusion of macromolecules or submicron tracers in biological cells. Here, we examine HDPs in quenched-disorder systems with Gaussian colored noise (GCN) characterized by a diffusion coefficient with a power-law dependence on the particle position and with a spatially random scaling exponent. Typically, D(x) is considered to be centerd at the origin and the entire x axis is characterized by a single scaling exponent a. In this work we consider a spatially random scenario: in periodic intervals ("layers") in space D(x) is centerd to the midpoint of each interval. In each interval the scaling exponent alpha is randomly chosen from a Gaussian distribution. The effects of the variation of the scaling exponents, the periodicity of the domains ("layer thickness") of the diffusion coefficient in this stratified system, and the correlation time of the GCN are analyzed numerically in detail. We discuss the regimes of superdiffusion, subdiffusion, and normal diffusion realisable in this system. We observe and quantify the domains where nonergodic and non-Gaussian behaviors emerge in this system. Our results provide new insights into the understanding of weak ergodicity breaking for HDPs driven by colored noise, with potential applications in quenched layered systems, typical model systems for diffusion in biological cells and tissues, as well as for diffusion in geophysical systems.
The max-min and min-max of matrices arise prevalently in science and engineering. However, in many real-world situations the computation of the max-min and min-max is challenging as matrices are large and full information about their entries is lacking. Here we take a statistical-physics approach and establish limit laws—akin to the central limit theorem—for the max-min and min-max of large random matrices. The limit laws intertwine random-matrix theory and extreme-value theory, couple the matrix dimensions geometrically, and assert that Gumbel statistics emerge irrespective of the matrix entries' distribution. Due to their generality and universality, as well as their practicality, these results are expected to have a host of applications in the physical sciences and beyond.
Classical option pricing schemes assume that the value of a financial asset follows a geometric Brownian motion (GBM). However, a growing body of studies suggest that a simple GBM trajectory is not an adequate representation for asset dynamics, due to irregularities found when comparing its properties with empirical distributions. As a solution, we investigate a generalisation of GBM where the introduction of a memory kernel critically determines the behaviour of the stochastic process. We find the general expressions for the moments, log-moments, and the expectation of the periodic log returns, and then obtain the corresponding probability density functions using the subordination approach. Particularly, we consider subdiffusive GBM (sGBM), tempered sGBM, a mix of GBM and sGBM, and a mix of sGBMs. We utilise the resulting generalised GBM (gGBM) in order to examine the empirical performance of a selected group of kernels in the pricing of European call options. Our results indicate that the performance of a kernel ultimately depends on the maturity of the option and its moneyness.
In the scenario of the narrow escape problem (NEP) a particle diffuses in a finite container and eventually leaves it through a small 'escape window' in the otherwise impermeable boundary, once it arrives to this window and crosses an entropic barrier at the entrance to it. This generic problem is mathematically identical to that of a diffusion-mediated reaction with a partially-reactive site on the container's boundary. Considerable knowledge is available on the dependence of the mean first-reaction time (FRT) on the pertinent parameters. We here go a distinct step further and derive the full FRT distribution for the NEP. We demonstrate that typical FRTs may be orders of magnitude shorter than the mean one, thus resulting in a strong defocusing of characteristic temporal scales. We unveil the geometry-control of the typical times, emphasising the role of the initial distance to the target as a decisive parameter. A crucial finding is the further FRT defocusing due to the barrier, necessitating repeated escape or reaction attempts interspersed with bulk excursions. These results add new perspectives and offer a broad comprehension of various features of the by-now classical NEP that are relevant for numerous biological and technological systems.
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 obtain a generalized diffusion equation in modified or Riemann-Liouville form from continuous time random walk theory. The waiting time probability density function and mean squared displacement for different forms of the equation are explicitly calculated. We show examples of generalized diffusion equations in normal or Caputo form that encode the same probability distribution functions as those obtained from the generalized diffusion equation in modified form. The obtained equations are general and many known fractional diffusion equations are included as special cases.
Numerous examples for a priori unexpected non-Gaussian behaviour for normal and anomalous diffusion have recently been reported in single-particle tracking experiments. Here, we address the case of non-Gaussian anomalous diffusion in terms of a random-diffusivity mechanism in the presence of power-law correlated fractional Gaussian noise. We study the ergodic properties of this model via examining the ensemble- and time-averaged mean-squared displacements as well as the ergodicity breaking parameter EB quantifying the trajectory-to-trajectory fluctuations of the latter. For long measurement times, interesting crossover behaviour is found as function of the correlation time tau characterising the diffusivity dynamics. We unveil that at short lag times the EB parameter reaches a universal plateau. The corresponding residual value of EB is shown to depend only on tau and the trajectory length. The EB parameter at long lag times, however, follows the same power-law scaling as for fractional Brownian motion. We also determine a corresponding plateau at short lag times for the discrete representation of fractional Brownian motion, absent in the continuous-time formulation. These analytical predictions are in excellent agreement with results of computer simulations of the underlying stochastic processes. Our findings can help distinguishing and categorising certain nonergodic and non-Gaussian features of particle displacements, as observed in recent single-particle tracking experiments.
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.
Fractional Brownian motion (FBM) is a Gaussian stochastic process with stationary, long-time correlated increments and is frequently used to model anomalous diffusion processes. We study numerically FBM confined to a finite interval with reflecting boundary conditions. The probability density function of this reflected FBM at long times converges to a stationary distribution showing distinct deviations from the fully flat distribution of amplitude 1/L in an interval of length L found for reflected normal Brownian motion. While for superdiffusion, corresponding to a mean squared displacement (MSD) 〈X² (t)〉 ⋍ tᵅ with 1 < α < 2, the probability density function is lowered in the centre of the interval and rises towards the boundaries, for subdiffusion (0 < α < 1) this behaviour is reversed and the particle density is depleted close to the boundaries. The MSD in these cases at long times converges to a stationary value, which is, remarkably, monotonically increasing with the anomalous diffusion exponent α. Our a priori surprising results may have interesting consequences for the application of FBM for processes such as molecule or tracer diffusion in the confines of living biological cells or organelles, or other viscoelastic environments such as dense liquids in microfluidic chambers.
Percolation networks have been widely used in the description of porous media but are now found to be relevant to understand the motion of particles in cellular membranes or the nucleus of biological cells. Random walks on the infinite cluster at criticality of a percolation network are asymptotically ergodic. On any finite size cluster of the network stationarity is reached at finite times, depending on the cluster's size. Despite of this we here demonstrate by combination of analytical calculations and simulations that at criticality the disorder and cluster size average of the ensemble of clusters leads to a non-vanishing variance of the time averaged mean squared displacement, regardless of the measurement time. Fluctuations of this relevant experimental quantity due to the disorder average of such ensembles are thus persistent and non-negligible. The relevance of our results for single particle tracking analysis in complex and biological systems is discussed.
Lévy flights are paradigmatic generalised random walk processes, in which the independent stationary increments—the 'jump lengths'—are drawn from an -stable jump length distribution with long-tailed, power-law asymptote. As a result, the variance of Lévy flights diverges and the trajectory is characterised by occasional extremely long jumps. Such long jumps significantly decrease the probability to revisit previous points of visitation, rendering Lévy flights efficient search processes in one and two dimensions. To further quantify their precise property as random search strategies we here study the first-passage time properties of Lévy flights in one-dimensional semi-infinite and bounded domains for symmetric and asymmetric jump length distributions. To obtain the full probability density function of first-passage times for these cases we employ two complementary methods. One approach is based on the space-fractional diffusion equation for the probability density function, from which the survival probability is obtained for different values of the stable index and the skewness (asymmetry) parameter . The other approach is based on the stochastic Langevin equation with -stable driving noise. Both methods have their advantages and disadvantages for explicit calculations and numerical evaluation, and the complementary approach involving both methods will be profitable for concrete applications. We also make use of the Skorokhod theorem for processes with independent increments and demonstrate that the numerical results are in good agreement with the analytical expressions for the probability density function of the first-passage times.
We present rigorous results for the mean first passage time and first passage time statistics for two-channel Markov additive diffusion in a 3-dimensional spherical domain. Inspired by biophysical examples we assume that the particle can only recognise the target in one of the modes, which is shown to effect a non-trivial first passage behaviour. We also address the scenario of intermittent immobilisation. In both cases we prove that despite the perfectly non-recurrent motion of two-channel Markov additive diffusion in 3 dimensions the first passage statistics at long times do not display Poisson-like behaviour if none of the phases has a vanishing diffusion coefficient. This stands in stark contrast to the standard (one-channel) Markov diffusion counterpart. We also discuss the relevance of our results in the context of cellular signalling.
A rapidly increasing number of systems is identified in which the stochastic motion of tracer particles follows the Brownian law < r(2)(t)> similar or equal to Dt yet the distribution of particle displacements is strongly non-Gaussian. A central approach to describe this effect is the diffusing diffusivity (DD) model in which the diffusion coefficient itself is a stochastic quantity, mimicking heterogeneities of the environment encountered by the tracer particle on its path. We here quantify in terms of analytical and numerical approaches the first passage behaviour of the DD model. We observe significant modifications compared to Brownian-Gaussian diffusion, in particular that the DD model may have a faster first passage dynamics. Moreover we find a universal crossover point of the survival probability independent of the initial condition.
We study the first passage dynamics for a diffusing particle experiencing a spatially varying diffusion coefficient while driven by correlated additive Gaussian white noise and multiplicative coloured non-Gaussian noise. We consider three functional forms for position dependence of the diffusion coefficient: power-law, exponential, and logarithmic. The coloured non-Gaussian noise is distributed according to Tsallis' q-distribution. Tracks of the non-Markovian systems are numerically simulated by using the fourth-order Runge-Kutta algorithm and the first passage times (FPTs) are recorded. The FPT density is determined along with the mean FPT (MFPT). Effects of the noise intensity and self-correlation of the multiplicative noise, the intensity of the additive noise, the cross-correlation strength, and the non-extensivity parameter on the MFPT are discussed.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
For both Lévy flight and Lévy walk search processes we analyse the full distribution of first-passage and first-hitting (or first-arrival) times. These are, respectively, the times when the particle moves across a point at some given distance from its initial position for the first time, or when it lands at a given point for the first time. For Lévy motions with their propensity for long relocation events and thus the possibility to jump across a given point in space without actually hitting it ('leapovers'), these two definitions lead to significantly different results. We study the first-passage and first-hitting time distributions as functions of the Lévy stable index, highlighting the different behaviour for the cases when the first absolute moment of the jump length distribution is finite or infinite. In particular we examine the limits of short and long times. Our results will find their application in the mathematical modelling of random search processes as well as computer algorithms.
What are the physical laws of the diffusive search of proteins for their specific binding sites on DNA in the presence of the macromolecular crowding in cells? We performed extensive computer simulations to elucidate the protein target search on DNA. The novel feature is the viscoelastic non-Brownian protein bulk diffusion recently observed experimentally. We examine the influence of the protein-DNA binding affinity and the anomalous diffusion exponent on the target search time. In all cases an optimal search time is found. The relative contribution of intermittent three-dimensional bulk diffusion and one-dimensional sliding of proteins along the DNA is quantified. Our results are discussed in the light of recent single molecule tracking experiments, aiming at a better understanding of the influence of anomalous kinetics of proteins on the facilitated diffusion mechanism.
A panoply of new tools for tracking single particles and molecules has led to an explosion of experimental data, leading to novel insights into physical properties of living matter governing cellular development and function, health and disease. In this Perspective, we present tools to investigate the dynamics and mechanics of living systems from the molecular to cellular scale via single-particle techniques. In particular, we focus on methods to measure, interpret, and analyse complex data sets that are associated with forces, materials properties, transport, and emergent organisation phenomena within biological and soft-matter systems. Current approaches, challenges, and existing solutions in the associated fields are outlined in order to support the growing community of researchers at the interface of physics and the life sciences. Each section focuses not only on the general physical principles and the potential for understanding living matter, but also on details of practical data extraction and analysis, discussing limitations, interpretation, and comparison across different experimental realisations and theoretical frameworks. Particularly relevant results are introduced as examples. While this Perspective describes living matter from a physical perspective, highlighting experimental and theoretical physics techniques relevant for such systems, it is also meant to serve as a solid starting point for researchers in the life sciences interested in the implementation of biophysical methods.
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.
Bias plays an important role in the enhancement of diffusion in periodic potentials. Using the continuous-time random walk in the presence of a bias, we report on an interesting phenomenon for the enhancement of diffusion by the start of the measurement in a random energy landscape. When the variance of the waiting time diverges, in contrast to the bias-free case, the dynamics with bias becomes superdiffusive. In the superdiffusive regime, we find a distinct initial ensemble dependence of the diffusivity. Moreover, the diffusivity can be increased by the aging time when the initial ensemble is not in equilibrium. We show that the time-averaged variance converges to the corresponding ensemble-averaged variance; i.e., ergodicity is preserved. However, trajectory-to-trajectory fluctuations of the time-averaged variance decay unexpectedly slowly. Our findings provide a rejuvenation phenomenon in the superdiffusive regime, that is, the diffusivity for a nonequilibrium initial ensemble gradually increases to that for an equilibrium ensemble when the start of the measurement is delayed.
Effects of the target aspect ratio and intrinsic reactivity onto diffusive search in bounded domains
(2017)
Westudy the mean first passage time (MFPT) to a reaction event on a specific site in a cylindrical geometry—characteristic, for instance, for bacterial cells, with a concentric inner cylinder representing the nuclear region of the bacterial cell. Asimilar problem emerges in the description of a diffusive search by a transcription factor protein for a specific binding region on a single strand of DNA.We develop a unified theoretical approach to study the underlying boundary value problem which is based on a self-consistent approximation of the mixed boundary condition. Our approach permits us to derive explicit, novel, closed-form expressions for the MFPT valid for a generic setting with an arbitrary relation between the system parameters.Weanalyse this general result in the asymptotic limits appropriate for the above-mentioned biophysical problems. Our investigation reveals the crucial role of the target aspect ratio and of the intrinsic reactivity of the binding region, which were disregarded in previous studies. Theoretical predictions are confirmed by numerical simulations.
Isoflux 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 fc. In the high driving force limit, regardless of the channel width, IFTP theory predicts τ ∝ f β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.
We study the probability density function (PDF) of the first-reaction times between a diffusive ligand and a membrane-bound, immobile imperfect target region in a restricted 'onion-shell' geometry bounded by two nested membranes of arbitrary shapes. For such a setting, encountered in diverse molecular signal transduction pathways or in the narrow escape problem with additional steric constraints, we derive an exact spectral form of the PDF, as well as present its approximate form calculated by help of the so-called self-consistent approximation. For a particular case when the nested domains are concentric spheres, we get a fully explicit form of the approximated PDF, assess the accuracy of this approximation, and discuss various facets of the obtained distributions. Our results can be straightforwardly applied to describe the PDF of the terminal reaction event in multi-stage signal transduction processes.
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 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.
Abstract
The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive–diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive–diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion–diffusion and subdiffusion–subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling.
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.
We study the adsorption–desorption transition of polyelectrolyte chains onto planar, cylindrical and spherical surfaces with arbitrarily high surface charge densities by massive Monte Carlo computer simulations. We examine in detail how the well known scaling relations for the threshold transition—demarcating the adsorbed and desorbed domains of a polyelectrolyte near weakly charged surfaces—are altered for highly charged interfaces. In virtue of high surface potentials and large surface charge densities, the Debye–Hückel approximation is often not feasible and the nonlinear Poisson–Boltzmann approach should be implemented. At low salt conditions, for instance, the electrostatic potential from the nonlinear Poisson–Boltzmann equation is smaller than the Debye–Hückel result, such that the required critical surface charge density for polyelectrolyte adsorption σc increases. The nonlinear relation between the surface charge density and electrostatic potential leads to a sharply increasing critical surface charge density with growing ionic strength, imposing an additional limit to the critical salt concentration above which no polyelectrolyte adsorption occurs at all. We contrast our simulations findings with the known scaling results for weak critical polyelectrolyte adsorption onto oppositely charged surfaces for the three standard geometries. Finally, we discuss some applications of our results for some physical–chemical and biophysical systems.
Employing extensive Monte Carlo computer simulations, we investigate in detail the properties of multichain adsorption of charged flexible polyelectrolytes (PEs) onto oppositely charged spherical nanoparticles (SNPs). We quantify the conditions of critical adsorption-the phase-separation curve between the adsorbed and desorbed states of the PEs-as a function of the SNP surface-charge density and the concentration of added salt. We study the degree of fluctuations of the PE-SNP electrostatic binding energy, which we use to quantify the emergence of the phase subtransitions, including a series of partially adsorbed PE configurations. We demonstrate how the phase-separation adsorption-desorption boundary shifts and splits into multiple subtransitions at low-salt conditions, thereby generalizing and extending the results for critical adsorption of a single PE onto the SNP. The current findings are relevant for finite concentrations of PEs around the attracting SNP, such as the conditions for PE adsorption onto globular proteins carrying opposite electric charges.
We consider the emerging dynamics of a separable continuous time random walk (CTRW) in the case when the random walker is biased by a velocity field in a uniformly growing domain. Concrete examples for such domains include growing biological cells or lipid vesicles, biofilms and tissues, but also macroscopic systems such as expanding aquifers during rainy periods, or the expanding Universe. The CTRW in this study can be subdiffusive, normal diffusive or superdiffusive, including the particular case of a Lévy flight. We first consider the case when the velocity field is absent. In the subdiffusive case, we reveal an interesting time dependence of the kurtosis of the particle probability density function. In particular, for a suitable parameter choice, we find that the propagator, which is fat tailed at short times, may cross over to a Gaussian-like propagator. We subsequently incorporate the effect of the velocity field and derive a bi-fractional diffusion-advection equation encoding the time evolution of the particle distribution. We apply this equation to study the mixing kinetics of two diffusing pulses, whose peaks move towards each other under the action of velocity fields acting in opposite directions. This deterministic motion of the peaks, together with the diffusive spreading of each pulse, tends to increase particle mixing, thereby counteracting the peak separation induced by the domain growth. As a result of this competition, different regimes of mixing arise. In the case of Lévy flights, apart from the non-mixing regime, one has two different mixing regimes in the long-time limit, depending on the exact parameter choice: in one of these regimes, mixing is mainly driven by diffusive spreading, while in the other mixing is controlled by the velocity fields acting on each pulse. Possible implications for encounter–controlled reactions in real systems are discussed.
Levy walks are continuous time random walks with spatio-temporal coupling of jump lengths and waiting times, often used to model superdiffusive spreading processes such as animals searching for food, tracer motion in weakly chaotic systems, or even the dynamics in quantum systems such as cold atoms. In the simplest version Levy walks move with a finite speed. Here, we present an extension of the Levy walk scenario for the case when external force fields influence the motion. The resulting motion is a combination of the response to the deterministic force acting on the particle, changing its velocity according to the principle of total energy conservation, and random velocity reversals governed by the distribution of waiting times. For the fact that the motion stays conservative, that is, on a constant energy surface, our scenario is fundamentally different from thermal motion in the same external potentials. In particular, we present results for the velocity and position distributions for single well potentials of different steepness. The observed dynamics with its continuous velocity changes enriches the theory of Levy walk processes and will be of use in a variety of systems, for which the particles are externally confined.
We show that the codifference is a useful tool in studying the ergodicity breaking and non-Gaussianity properties of stochastic time series. While the codifference is a measure of dependence that was previously studied mainly in the context of stable processes, we here extend its range of applicability to random-parameter and diffusing-diffusivity models which are important in contemporary physics, biology and financial engineering. We prove that the codifference detects forms of dependence and ergodicity breaking which are not visible from analysing the covariance and correlation functions. We also discuss a related measure of dispersion, which is a nonlinear analogue of the mean squared displacement.
Anomalous diffusion with a power-law time dependence vertical bar R vertical bar(2)(t) similar or equal to t(alpha i) of the mean squared displacement occurs quite ubiquitously in numerous complex systems. Often, this anomalous diffusion is characterised by crossovers between regimes with different anomalous diffusion exponents alpha(i). Here we consider the case when such a crossover occurs from a first regime with alpha(1) to a second regime with alpha(2) such that alpha(2) > alpha(1), i.e., accelerating anomalous diffusion. A widely used framework to describe such crossovers in a one-dimensional setting is the bi-fractional diffusion equation of the so-called modified type, involving two time-fractional derivatives defined in the Riemann-Liouville sense. We here generalise this bi-fractional diffusion equation to higher dimensions and derive its multidimensional propagator (Green's function) for the general case when also a space fractional derivative is present, taking into consideration long-ranged jumps (Levy flights). We derive the asymptotic behaviours for this propagator in both the short- and long-time as well the short- and long-distance regimes. Finally, we also calculate the mean squared displacement, skewness and kurtosis in all dimensions, demonstrating that in the general case the non-Gaussian shape of the probability density function changes.
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.
Characterising stochastic motion in heterogeneous media driven by coloured non-Gaussian noise
(2021)
We study the stochastic motion of a test particle in a heterogeneous medium in terms of a position dependent diffusion coefficient mimicking measured deterministic diffusivity gradients in biological cells or the inherent heterogeneity of geophysical systems. Compared to previous studies we here investigate the effect of the interplay of anomalous diffusion effected by position dependent diffusion coefficients and coloured non-Gaussian noise. The latter is chosen to be distributed according to Tsallis' q-distribution, representing a popular example for a non-extensive statistic. We obtain the ensemble and time averaged mean squared displacements for this generalised process and establish its non-ergodic properties as well as analyse the non-Gaussian nature of the associated displacement distribution. We consider both non-stratified and stratified environments.
The stable operation of a turbulent combustor is not completely silent; instead, there is a background of small amplitude aperiodic acoustic fluctuations known as combustion noise. Pressure fluctuations during this state of combustion noise are multifractal due to the presence of multiple temporal scales that contribute to its dynamics. However, existing models are unable to capture the multifractality in the pressure fluctuations. We conjecture an underlying fractional dynamics for the thermoacoustic system and obtain a fractional-order model for pressure fluctuations. The data from this model has remarkable visual similarity to the experimental data and also has a wide multifractal spectrum during the state of combustion noise. Quantitative similarity with the experimental data in terms of the Hurst exponent and the multifractal spectrum is observed during the state of combustion noise. This model is also able to produce pressure fluctuations that are qualitatively similar to the experimental data acquired during intermittency and thermoacoustic instability. Furthermore, we argue that the fractional dynamics vanish as we approach the state of thermoacoustic instability.
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.
Brownian motion and beyond: first-passage, power spectrum, non-Gaussianity, and anomalous diffusion
(2019)
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 exists a very large range of applications of stochastic processes in various disciplines. Here we provide a summary of some of the recent developments in the field of stochastic processes, highlighting both the experimental findings and theoretical frameworks.
We examine renewal processes with power-law waiting time distributions (WTDs) and non-zero drift via computing analytically and by computer simulations their ensemble and time averaged spreading characteristics. All possible values of the scaling exponent alpha are considered for the WTD psi(t) similar to 1/t(1+alpha). We treat continuous-time random walks (CTRWs) with 0 < alpha < 1 for which the mean waiting time diverges, and investigate the behaviour of the process for both ordinary and equilibrium CTRWs for 1 < alpha < 2 and alpha > 2. We demonstrate that in the presence of a drift CTRWs with alpha < 1 are ageing and non-ergodic in the sense of the non-equivalence of their ensemble and time averaged displacement characteristics in the limit of lag times much shorter than the trajectory length. In the sense of the equivalence of ensemble and time averages, CTRW processes with 1 < alpha < 2 are ergodic for the equilibrium and non-ergodic for the ordinary situation. Lastly, CTRW renewal processes with alpha > 2-both for the equilibrium and ordinary situation-are always ergodic. For the situations 1 < alpha < 2 and alpha > 2 the variance of the diffusion process, however, depends on the initial ensemble. For biased CTRWs with alpha > 1 we also investigate the behaviour of the ergodicity breaking parameter. In addition, we demonstrate that for biased CTRWs the Einstein relation is valid on the level of the ensemble and time averaged displacements, in the entire range of the WTD exponent alpha.
Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusionmodel and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a wellcalibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output.
Modern single-particle-tracking techniques produce extensive time-series of diffusive motion in a wide variety of systems, from single-molecule motion in living-cells to movement ecology. The quest is to decipher the physical mechanisms encoded in the data and thus to better understand the probed systems. We here augment recently proposed machine-learning techniques for decoding anomalous-diffusion data to include an uncertainty estimate in addition to the predicted output. To avoid the Black-Box-Problem a Bayesian-Deep-Learning technique named Stochastic-Weight-Averaging-Gaussian is used to train models for both the classification of the diffusion model and the regression of the anomalous diffusion exponent of single-particle-trajectories. Evaluating their performance, we find that these models can achieve a well-calibrated error estimate while maintaining high prediction accuracies. In the analysis of the output uncertainty predictions we relate these to properties of the underlying diffusion models, thus providing insights into the learning process of the machine and the relevance of the output. <br /> Diffusive motions in complex environments such as living biological cells or soft matter systems can be analyzed with single-particle-tracking approaches, where accuracy of output may vary. The authors involve a machine-learning technique for decoding anomalous-diffusion data and provide an uncertainty estimate together with predicted output.
We employ Bayesian statistics using the nested-sampling algorithm to compare and rank multiple models of ergodic diffusion (including anomalous diffusion) as well as to assess their optimal parameters for in silico-generated and real time-series. We focus on the recently-introduced model of Brownian motion with "diffusing diffusivity'-giving rise to widely-observed non-Gaussian displacement statistics-and its comparison to Brownian and fractional Brownian motion, also for the time-series with some measurement noise. We conduct this model-assessment analysis using Bayesian statistics and the nested-sampling algorithm on the level of individual particle trajectories. We evaluate relative model probabilities and compute best-parameter sets for each diffusion model, comparing the estimated parameters to the true ones. We test the performance of the nested-sampling algorithm and its predictive power both for computer-generated (idealised) trajectories as well as for real single-particle-tracking trajectories. Our approach delivers new important insight into the objective selection of the most suitable stochastic model for a given time-series. We also present first model-ranking results in application to experimental data of tracer diffusion in polymer-based hydrogels.
We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta.
We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times.
We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient.
When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity.
We analyse mobile-immobile transport of particles that switch between the mobile and immobile phases with finite rates. Despite this seemingly simple assumption of Poissonian switching, we unveil a rich transport dynamics including significant transient anomalous diffusion and non-Gaussian displacement distributions. Our discussion is based on experimental parameters for tau proteins in neuronal cells, but the results obtained here are expected to be of relevance for a broad class of processes in complex systems. Specifically, we obtain that, when the mean binding time is significantly longer than the mean mobile time, transient anomalous diffusion is observed at short and intermediate time scales, with a strong dependence on the fraction of initially mobile and immobile particles. We unveil a Laplace distribution of particle displacements at relevant intermediate time scales. For any initial fraction of mobile particles, the respective mean squared displacement (MSD) displays a plateau. Moreover, we demonstrate a short-time cubic time dependence of the MSD for immobile tracers when initially all particles are immobile.