Stretched-exponential relaxation is a widely observed phenomenon found in ordered ferromagnets as well as glassy systems. One modeling approach connects this behavior to a droplet dynamics described by an effective Langevin equation for the droplet radius with an r^{2/3} potential. Here, we study a Brownian particle under the influence of a general confining, albeit weak, potential field that grows with distance as a sublinear power law. We find that for this memoryless model, observables display stretched-exponential relaxation. The probability density function of the system is studied using a rate-function ansatz. We obtain analytically the stretched-exponential exponent along with an anomalous power-law scaling of length with time. The rate function exhibits a point of nonanalyticity, indicating a dynamical phase transition. In particular, the rate function is double valued both to the left and right of this point, leading to four different rate functions, depending on the choice of initial conditions and symmetry.
Since the times of Holtsmark (1911), statistics of fields in random environments have been widely studied, for example in astrophysics, active matter, and line-shape broadening. The power-law decay of the two-body interaction of the form 1/|r|^{δ}, and assuming spatial uniformity of the medium particles exerting the forces, imply that the fields are fat-tailed distributed, and in general are described by stable Lévy distributions. With this widely used framework, the variance of the field diverges, which is nonphysical, due to finite size cutoffs. We find a complementary statistical law to the Lévy-Holtsmark distribution describing the large fields in the problem, which is related to the finite size of the tracer particle. We discover biscaling with a sharp statistical transition of the force moments taking place when the order of the moment is d/δ, where d is the dimension. The high-order moments, including the variance, are described by the framework presented in this paper, which is expected to hold for many systems. The new scaling solution found here is nonnormalized similar to infinite invariant densities found in dynamical systems.
In many disordered systems, the diffusion of classical particles is described by a displacement distribution P(x, t) that displays exponential tails instead of Gaussian statistics expected for Brownian motion. However, the experimental demonstration of control of this behavior by increasing the disorder strength has remained challenging. In this work, we explore the Gaussian-to-exponential transition by using diffusion of poly(ethylene glycol) (PEG) in attractive nanoparticle-polymer mixtures and controlling the volume fraction of the nanoparticles. In this work, we find "knobs", namely nanoparticle concentration and interaction, which enable the change in the shape of P(x,t) in a well-defined way. The Gaussian-to-exponential transition is consistent with a modified large deviation approach for a continuous time random walk and also with Monte Carlo simulations involving a microscopic model of polymer trapping via reversible adsorption to the nanoparticle surface. Our work bears significance in unraveling the fundamental physics behind the exponential decay of the displacement distribution at the tails, which is commonly observed in soft materials and nanomaterials.
First-passage time statistics in disordered systems exhibiting scale invariance are studied widely. In particular, long trapping times in energy or entropic traps are fat-tailed distributed, which slow the overall transport process. We study the statistical properties of the first-passage time of biased processes in different models, and we employ the big-jump principle that shows the dominance of the maximum trapping time on the first-passage time. We demonstrate that the removal of this maximum significantly expedites transport. As the disorder increases, the system enters a phase where the removal shows a dramatic effect. Our results show how we may speed up transport in strongly disordered systems exploiting scale invariance. In contrast to the disordered systems studied here, the removal principle has essentially no effect in homogeneous systems; this indicates that improving the conductance of a poorly conducting system is, theoretically, relatively easy as compared to a homogeneous system.
Particles anomalously diffusing in contact with a thermal bath are initially released from an asymptotically flat potential well. For temperatures that are sufficiently low compared to the potential depth, the dynamical and thermodynamical observables of the system remain almost constant for long times. We show how these stagnated states are characterized as non-normalizable quasiequilibrium (NNQE) states. We use the fractional-time Fokker-Planck equation (FTFPE) and continuous-time random walk approaches to calculate ensemble averages. We obtain analytical estimates of the durations of NNQE states, depending on the fractional order, from approximate theoretical solutions of the FTFPE. We study and compare two types of observables, the mean square displacement typically used to characterize diffusion, and the thermodynamic energy. We show that the typical timescales for transient stagnation depend exponentially on the value of the depth of the potential well, in units of temperature, multiplied by a function of the fractional exponent.
We study the motion of an overdamped particle connected to a thermal heat bath in the presence of an external periodic potential in one dimension. When we coarse-grain, i.e., bin the particle positions using bin sizes that are larger than the periodicity of the potential, the packet of spreading particles, all starting from a common origin, converges to a normal distribution centered at the origin with a mean-squared displacement that grows as 2D^{*}t, with an effective diffusion constant that is smaller than that of a freely diffusing particle. We examine the interplay between this coarse-grained description and the fine structure of the density, which is given by the Boltzmann-Gibbs (BG) factor e^{-V(x)/k_{B}T}, the latter being nonnormalizable. We explain this result and construct a theory of observables using the Fokker-Planck equation. These observables are classified as those that are related to the BG fine structure, like the energy or occupation times, while others, like the positional moments, for long times, converge to those of the large-scale description. Entropy falls into a special category as it has a coarse-grained and a fine structure description. The basic thermodynamic formula F=TS-E is extended to this far-from-equilibrium system. The ergodic properties are also studied using tools from infinite ergodic theory.
Extreme value (EV) statistics of correlated systems are widely investigated in many fields, spanning the spectrum from weather forecasting to earthquake prediction. Does the unavoidable discrete sampling of a continuous correlated stochastic process change its EV distribution? We explore this question for correlated random variables modeled via Langevin dynamics for a particle in a potential field. For potentials growing at infinity faster than linearly and for long measurement times, we find that the EV distribution of the discretely sampled process diverges from that of the full continuous dataset and converges to that of independent and identically distributed random variables drawn from the process's equilibrium measure. However, for processes with sublinear potentials, the long-time limit is the EV statistics of the continuously sampled data. We treat processes whose equilibrium measures belong to the three EV attractors: Gumbel, Fréchet, and Weibull. Our Letter shows that the EV statistics can be extremely sensitive to the sampling rate of the data.
We compare ergodic properties of the kinetic energy for three stochastic models of subrecoil-laser-cooled gases. One model is based on a heterogeneous random walk (HRW), another is an HRW with long-range jumps (the exponential model), and the other is a mean-field-like approximation of the exponential model (the deterministic model). All the models show an accumulation of the momentum at zero in the long-time limit, and a formal steady state cannot be normalized, i.e., there exists an infinite invariant density. We obtain the exact form of the infinite invariant density and the scaling function for the exponential and deterministic models, and we devise a useful approximation for the momentum distribution in the HRW model. While the models are kinetically nonidentical, it is natural to wonder whether their ergodic properties share common traits, given that they are all described by an infinite invariant density. We show that the answer to this question depends on the type of observable under study. If the observable is integrable, the ergodic properties, such as the statistical behavior of the time averages, are universal as they are described by the Darling-Kac theorem. In contrast, for nonintegrable observables, the models in general exhibit nonidentical statistical laws. This implies that focusing on nonintegrable observables, we discover nonuniversal features of the cooling process, which hopefully can lead to a better understanding of the particular model most suitable for a statistical description of the process. This result is expected to hold true for many other systems, beyond laser cooling.
The virial theorem, and the equipartition theorem in the case of quadratic degrees of freedom, are handy constraints on the statistics of equilibrium systems. Their violation is instrumental in determining how far from equilibrium a driven system might be. We extend the virial theorem to nonequilibrium conditions for Langevin dynamics with nonlinear friction and multiplicative noise. In particular, we generalize it for confined laser-cooled atoms in the semiclassical regime. The resulting relation between the lowest moments of the atom position and velocity allows to measure in experiments how dissipative the cooling mechanism is. Moreover, its violation can reveal the departure from a strictly harmonic confinement or from the semiclassical regime.
The velocity distribution of a classical gas of atoms in thermal equilibrium is the normal Maxwell distribution. It is well known that for sub-recoiled laser cooled atoms, Lévy statistics and deviations from usual ergodic behavior come into play. In a recent letter, we showed how tools from infinite ergodic theory describe the cool gas. Here, using the master equation, we derive the scaling function and the infinite invariant density of a stochastic model for the momentum of laser cooled atoms, recapitulating results obtained by Bertin and Bardou [Am. J. Phys. 76, 630 (2008)] using life-time statistics. We focus on the case where the laser trapping is strong, namely, the rate of escape from the velocity trap is R(v) â |v|α for v â 0 and α > 1. We construct a machinery to investigate time averages of physical observables and their relation to ensemble averages. The time averages are given in terms of functionals of the individual stochastic paths, and here we use a generalization of Lévy walks to investigate the ergodic properties of the system. Exploring the energy of the system, we show that when α = 3, it exhibits a transition between phases where it is either an integrable or a non-integrable observable with respect to the infinite invariant measure. This transition corresponds to very different properties of the mean energy and to a discontinuous behavior of fluctuations. While the integrable phase is described by universal statistics and the Darling-Kac law, the more challenging case is the exploration of statistical properties of non-integrable observables. Since previous experimental work showed that both α = 2 and α = 4 are attainable, we believe that both phases could also be explored experimentally.
With subrecoil-laser-cooled atoms, one may reach nanokelvin temperatures while the ergodic properties of these systems do not follow usual statistical laws. Instead, due to an ingenious trapping mechanism in momentum space, power-law-distributed sojourn times are found for the cooled particles. Here, we show how this gives rise to a statistical-mechanical framework based on infinite ergodic theory, which replaces ordinary ergodic statistical physics of a thermal gas of atoms. In particular, the energy of the system exhibits a sharp discontinuous transition in its ergodic properties. Physically, this is controlled by the fluorescence rate, but, more profoundly, it is a manifestation of a transition for any observable, from being an integrable to becoming a nonintegrable observable, with respect to the infinite (non-normalized) invariant density.
Single-particle tracking offers detailed information about the motion of molecules in complex environments such as those encountered in live cells, but the interpretation of experimental data is challenging. One of the most powerful tools in the characterization of random processes is the power spectral density. However, because anomalous diffusion processes in complex systems are usually not stationary, the traditional Wiener-Khinchin theorem for the analysis of power spectral densities is invalid. Here, we employ a recently developed tool named aging Wiener-Khinchin theorem to derive the power spectral density of fractional Brownian motion coexisting with a scale-free continuous time random walk, the two most typical anomalous diffusion processes. Using this analysis, we characterize the motion of voltage-gated sodium channels on the surface of hippocampal neurons. Our results show aging where the power spectral density can either increase or decrease with observation time depending on the specific parameters of both underlying processes.
Membrane Transport Proteins/physiology , Computer Simulation , Diffusion , Models, Biological , Motion , Reproducibility of Results , Single Molecule Imaging , Time Factors
A classical random walker starting on a node of a finite graph will always reach any other node since the search is ergodic, namely it fully explores space, hence the arrival probability is unity. For quantum walks, destructive interference may induce effectively non-ergodic features in such search processes. Under repeated projective local measurements, made on a target state, the final detection of the system is not guaranteed since the Hilbert space is split into a bright subspace and an orthogonal dark one. Using this we find an uncertainty relation for the deviations of the detection probability from its classical counterpart, in terms of the energy fluctuations.
We investigate the overdamped Langevin motion for particles in a potential well that is asymptotically flat. When the potential well is deep as compared to the temperature, physical observables, like the mean square displacement, are essentially time-independent over a long time interval, the stagnation epoch. However, the standard Boltzmann-Gibbs (BG) distribution is non-normalizable, given that the usual partition function is divergent. For this regime, we have previously shown that a regularization of BG statistics allows for the prediction of the values of dynamical and thermodynamical observables in the non-normalizable quasi-equilibrium state. In this work, based on the eigenfunction expansion of the time-dependent solution of the associated Fokker-Planck equation with free boundary conditions, we obtain an approximate time-independent solution of the BG form, being valid for times that are long, but still short as compared to the exponentially large escape time. The escaped particles follow a general free-particle statistics, where the solution is an error function, which is shifted due to the initial struggle to overcome the potential well. With the eigenfunction solution of the Fokker-Planck equation in hand, we show the validity of the regularized BG statistics and how it perfectly describes the time-independent regime though the quasi-stationary state is non-normalizable.
We study the ballistic Lévy walk stemming from an infinite mean traveling time between collision events. Our study focuses on the density of spreading particles all starting from a common origin, which is limited by a "light" cone -v_{0}t
Recently observation of random walks in complex environments like the cell and other glassy systems revealed that the spreading of particles, at its tails, follows a spatial exponential decay instead of the canonical Gaussian. We use the widely applicable continuous time random walk model and obtain the large deviation description of the propagator. Under mild conditions that the microscopic jump lengths distribution is decaying exponentially or faster i.e., Lévy like power law distributed jump lengths are excluded, and that the distribution of the waiting times is analytical for short waiting times, the spreading of particles follows an exponential decay at large distances, with a logarithmic correction. Here we show how anti-bunching of jump events reduces the effect, while bunching and intermittency enhances it. We employ exact solutions of the continuous time random walk model to test the large deviation theory.
We investigate extreme value theory for physical systems with a global conservation law which describes renewal processes, mass transport models, and long-range interacting spin models. As shown previously, a special feature is that the distribution of the extreme value exhibits a nonanalytical point in the middle of the support. We expose exact relationships between constrained extreme value theory and well-known quantities of the underlying stochastic dynamics, all valid beyond the midpoint in general, i.e., even far from the thermodynamic limit. For example, for renewal processes the distribution of the maximum time between two renewal events is exactly related to the mean number of these events. In the thermodynamic limit, we show how our theory is suitable to describe typical and rare events which deviate from classical extreme value theory. For example, for the renewal process we unravel dual scaling of the extreme value distribution, pointing out two types of limiting laws: a normalizable scaling function for the typical statistics and a non-normalized state describing the rare events.
We report on a fundamental role of a non-normalized formal steady state, i.e., an infinite invariant density, in a semi-Markov process where the state is determined by the interevent time of successive renewals. The state describes certain observables found in models of anomalous diffusion, e.g., the velocity in the generalized Lévy walk model and the energy of a particle in the trap model. In our model, the interevent-time distribution follows a fat-tailed distribution, which makes the state value more likely to be zero because long interevent times imply small state values. We find two scaling laws describing the density for the state value, which accumulates in the vicinity of zero in the long-time limit. These laws provide universal behaviors in the accumulation process and give the exact expression of the infinite invariant density. Moreover, we provide two distributional limit theorems for time-averaged observables in these nonstationary processes. We show that the infinite invariant density plays an important role in determining the distribution of time averages.
Brownian motion is a Gaussian process described by the central limit theorem. However, exponential decays of the positional probability density function P(X,t) of packets of spreading random walkers, were observed in numerous situations that include glasses, live cells, and bacteria suspensions. We show that such exponential behavior is generally valid in a large class of problems of transport in random media. By extending the large deviations approach for a continuous time random walk, we uncover a general universal behavior for the decay of the density. It is found that fluctuations in the number of steps of the random walker, performed at finite time, lead to exponential decay (with logarithmic corrections) of P(X,t). This universal behavior also holds for short times, a fact that makes experimental observations readily achievable.
The prediction and control of rare events is an important task in disciplines that range from physics and biology, to economics and social science. The Big Jump principle deals with a peculiar aspect of the mechanism that drives rare events. According to the principle, in heavy-tailed processes a rare huge fluctuation is caused by a single event and not by the usual coherent accumulation of small deviations. We consider generalized Lévy walks, a class of stochastic processes with power law distributed step durations and with complex microscopic dynamics in the single stretch. We derive the bulk of the probability distribution and using the big jump principle, the exact form of the tails that describes rare events. We show that the tails of the distribution present non-universal and non-analytic behaviors, which depend crucially on the dynamics of the single step. The big jump estimate also provides a physical explanation of the processes driving the rare events, opening new possibilities for their correct prediction.
