BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Eli Barkai

Eli Barkai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechBiological Physicsphysics.data-anquant-phcond-mat.dis-nnmath-phmath.MPcond-mat.softnlin.CDQuantitative Methodsastro-phcond-mat.mes-hallmath.DSmath.NAphysics.comp-ph

Catalog footprint

What is connected

52works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2022arXiv

Aging power spectrum of membrane protein transport and other subordinated random walks

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 scalefree 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discrete Sampling of Extreme Events Modifies Their Statistics

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 work shows that the EV statistics can be extremely sensitive to the sampling rate of the data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exponential Tails and Asymmetry Relations for the Spread of Biased Random Walks

Exponential, and not Gaussian, decay of probability density functions was studied by Laplace in the context of his analysis of errors. Such Laplace propagators for the diffusive motion of single particles in disordered media were recently observed in numerous experimental systems. What will happen to this universality when an external driving force is applied? Using the ubiquitous continuous time random walk with bias, and the Crooks relation in conjunction with large deviations theory, we derive two properties of the positional probability density function $P_F(x,t)$ that hold for a wide spectrum of random walk models: (I) Universal asymmetric exponential decay of $P_F(X,t)$ for large $|X|$, and (II) Existence of a time transformation that for large $|X|$ allows to express $P_F(X,t)$ in terms of the propagator of the unbiased process (measured at a shorter time). These findings allow us to establish how the symmetric exponential-like tails, measured in many unbiased processes, will transform into asymmetric Laplace tails when an external force is applied.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Infinite ergodic theory for three heterogeneous stochastic models with application to subrecoil laser cooling

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 devise a useful approximation for the momentum distribution in the HRW model. While the models are kinetically non-identical, 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 non-integrable observables, the models in general exhibit non-identical statistical laws. This implies that focusing on non-integrable observables, we discover non-universal features of the cooling process, that 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Large deviations for continuous time random walks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local equilibrium properties of ultraslow diffusion in the Sinai model

We perform numerical studies of a thermally driven, overdamped particle in a random quenched force field, known as the Sinai model. We compare the unbounded motion on an infinite 1-dimensional domain to the motion in bounded domains with reflecting boundaries and show that the unbounded motion is at every time close to the equilibrium state of a finite system of growing size. This is due to time scale separation: Inside wells of the random potential, there is relatively fast equilibration, while the motion across major potential barriers is ultraslow. Quantities studied by us are the time dependent mean squared displacement, the time dependent mean energy of an ensemble of particles, and the time dependent entropy of the probability distribution. Using a very fast numerical algorithm, we can explore times up top $10^{17}$ steps and thereby also study finite-time crossover phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Driving quantum systems with repeated conditional measurements

We investigate the effect of conditional null measurements on a quantum system and find a rich variety of behaviors. Specifically, quantum dynamics with a time independent $H$ in a finite dimensional Hilbert space are considered with repeated strong null measurements of a specified state. We discuss four generic behaviors that emerge in these monitored systems. The first arises in systems without symmetry, along with their associated degeneracies in the energy spectrum, and hence in the absence of dark states as well. In this case, a unique final state can be found which is determined by the largest eigenvalue of the survival operator, the non-unitary operator encoding both the unitary evolution between measurements and the measurement itself. For a three-level system, this is similar to the well known shelving effect. Secondly, for systems with built-in symmetry and correspondingly a degenerate energy spectrum, the null measurements dynamically select the degenerate energy levels, while the non-degenerate levels are effectively wiped out. Thirdly, in the absence of dark states, and for specific choices of parameters, two or more eigenvalues of the survival operator match in magnitude, and this leads to an oscillatory behavior controlled by the measurement rate and not solely by the energy levels. Finally, when the control parameters are tuned, such that the eigenvalues of the survival operator all coalesce to zero, one has exceptional points that corresponds to situations that violate the null measurement condition, making the conditional measurement process impossible.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Gas of sub-recoiled laser cooled atoms described by infinite ergodic theory

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 behaviour come into play. Here we show how tools from infinite ergodic theory describe the cool gas. Specifically, we derive the scaling function and the infinite invariant density of a stochastic model for the momentum of laser cooled atoms using two approaches. The first is a direct analysis of the master equation and the second following the analysis of Bertin and Bardou using the lifetime dynamics. The two methods are shown to be identical, but yield different insights into the problem. In the main part of the paper we focus on the case where the laser trapping is strong, namely the rate of escape from the velocity trap is $R(v) \propto |v|^α$ for $v \to 0$ and $α>1$. We construct a machinery to investigate the 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 generalisation 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 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 behaviour of the fluctuations. Since previous experimental work showed that both $α=2$ and $α=4$ are attainable we believe that both phases could be explored also experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Generalized equipartition for nonlinear multiplicative Langevin dynamics: application to laser-cooled atoms

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 the equipartition theorem for confined laser-cooled atoms in the semi-classical 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 semi-classical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-normalizable quasi-equilibrium solution of the Fokker-Planck equation for nonconfining fields

We investigate the overdamped Langevin motion for particles in a potential well that is asymptotically flat. When the potential well is deep compared to 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 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, valid for times which are long, but still short compared to the exponentially large escape time. The escaped particles follow a general free-particle statistics, where the solution is a an error function, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fractional advection-diffusion-asymmetry equation

Fractional kinetic equations employ non-integer calculus to model anomalous relaxation and diffusion in many systems. While this approach is well explored, it so far failed to describe an important class of transport in disordered systems. Motivated by work on contaminant spreading in geological formations we propose and investigate a fractional advection-diffusion equation describing the biased spreading packet. While usual transport is described by diffusion and drift, we find a third term describing symmetry breaking which is omnipresent for transport in disordered systems. Our work is based on continuous time random walks with a finite mean waiting time and a diverging variance, a case that on the one hand is very common and on the other was missing in the kaleidoscope literature of fractional equations. The fractional space derivatives stem from long trapping times while previously they were interpreted as a consequence of spatial L{é}vy flights.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infinite ergodic theory meets Boltzmann statistics

We investigate the overdamped stochastic dynamics of a particle in an asymptotically flat external potential field, in contact with a thermal bath. For an infinite system size, the particles may escape the force field and diffuse freely at large length scales. The partition function diverges and hence the standard canonical ensemble fails. This is replaced with tools stemming from infinite ergodic theory. Boltzmann-Gibbs statistics, even though not normalized, still describes integrable observables, like energy and occupation times. The Boltzmann infinite density is derived heuristically using an entropy maximization principle, as well as via a first-principles calculation using an eigenfunction expansion in the continuum of low-energy states. A generalized virial theorem is derived, showing how the virial coefficient describes the delay in the diffusive spreading of the particles, found at large distances. When the process is non-recurrent, e.g. diffusion in three dimensions with a Coulomb-like potential, we use weighted time averages to restore basic canonical relations between time and ensemble averages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Infinite invariant density in a semi-Markov process with continuous state variables

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 inter-event 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 inter-event-time distribution follows a fat-tailed distribution, which makes the state value more likely to be zero because long inter-event 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 non-stationary processes. We show that the infinite invariant density plays an important role in determining the distribution of time averages.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Packets of Diffusing Particles Exhibit Universal Exponential Tails

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 holds also for short times, a fact that makes experimental observations readily achievable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantization of the mean decay time for non-Hermitian quantum systems

We show that the mean time, which a quantum particle needs to escape from a system to the environment, is quantized and independent from most dynamical details of the system. In particular, we consider a quantum system with a general Hermitian Hamiltonian $\hat{H}$ and one decay channel, through which probability dissipates to the environment with rate $Γ$. When the system is initially prepared exactly in the decay state, the mean decay time $\langle T \rangle$ is quantized and equal to $w/(2Γ)$. $w$ is the number of distinct energy levels, i.e. eigenvalues of $\hat{H}$, that have overlap with the decay state, and is also the winding number of a transform of the resolvent in the complex plane. Apart from the integer $w$, $\langle T \rangle$ is completely independent of the system's dynamics. The complete decay time distribution can be obtained from an electrostatic analogy and features rare events of very large dissipation times for parameter choices close to critical points, where $w$ changes, e.g. when a degeneracy is lifted. Experiments of insufficient observation time may thus measure a too small value of $w$. We discuss our findings in a disordered tight-binding model and in the two-level atom in a continuous-wave field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rare events in generalized Lévy Walks and the Big Jump principle

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, which model 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The first detection time of a quantum state under random probing

We solve for the statistics of the first detection of a quantum system in a particular desired state, when the system is subject to a projective measurement at independent identically distributed random time intervals. We present formulas for the probability of detection in the $n$th attempt. We calculate as well the mean and mean square both of the number of the first successful detection attempt and the time till first detection. We present explicit results for a particle initially localized at a site on a ring of size $L$, probed at some arbitrary given site, in the case when the detection intervals are distributed exponentially. We prove that, for all interval distributions and finite-dimensional Hamiltonians, the mean detection time is equal to the mean attempt number times the mean time interval between attempts. We further prove that for the return problem when the initial and target state are identical, the total detection probability is unity and the mean attempts till detection is an integer, which is the size of the Hilbert space (symmetrized about the target state). We study an interpolation between the fixed time interval case to an exponential distribution of time intervals via the Gamma distribution with constant mean and varying width. The mean arrival time as a function of the mean interval changes qualitatively as we tune the inter-arrival time distribution from very narrow (delta peaked) to exponential, as resonances are wiped out by the randomness of the sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Uncertainty and symmetry bounds for the quantum total detection probability

We investigate a generic discrete quantum system prepared in state $|ψ_\text{in}\rangle$, under repeated detection attempts aimed to find the particle in state $|d\rangle$, for example a quantum walker on a finite graph searching for a node. For the corresponding classical random walk, the total detection probability $P_\text{det}$ is unity. Due to destructive interference, one may find initial states $|ψ_\text{in}\rangle$ with $P_\text{det}<1$. We first obtain an uncertainty relation which yields insight on this deviation from classical behavior, showing the relation between $P_\text{det}$ and energy fluctuations: $ ΔP \,\mathrm{Var}[\hat{H}]_d \ge | \langle d| [\hat{H}, \hat{D}] | ψ_\text{in} \rangle |^2$ where $ΔP = P_\text{det} - |\langleψ_\text{in}|d\rangle |^2$, and $\hat{D} = |d\rangle\langle d|$ is the measurement projector. Secondly, exploiting symmetry we show that $P_\text{det}\le 1/ν$ where the integer $ν$ is the number of states equivalent to the initial state. These bounds are compared with the exact solution for small systems, obtained from an analysis of the dark and bright subspaces, showing the usefulness of the approach. The upper bounds works well even in large systems, and we show how to tighten the lower bound in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heavy-tailed phase-space distributions beyond Boltzmann-Gibbs and equipartition: Statistics of confined cold atoms

The Boltzmann-Gibbs density, a central result of equilibrium statistical mechanics, relates the energy of a system in contact with a thermal bath to its equilibrium statistics. This relation is lost for non-thermal systems such as cold atoms in optical lattices, where the heat bath is replaced by the laser beams of the lattice. We investigate in detail the stationary phase-space probability for Sisyphus cooling under harmonic confinement. In particular, we elucidate whether the total energy of the system still describes its stationary state statistics. We find that this is true for the center part of the phase-space density for deep lattices, where the Boltzmann-Gibbs density provides an approximate description. The relation between energy and statistics also persists for strong confinement and in the limit of high energies, where the system becomes underdamped. However, the phase-space density now exhibits heavy power-law tails. In all three cases we find expressions for the leading order phase-space density and corrections which break the equivalence of probability and energy and violate energy equipartition. The non-equilibrium nature of the steady state is confounded by explicit violations of detailed balance. We complement these analytical results with numerical simulations to map out the intricate structure of the phase-space density.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tempered Fractional Feynman-Kac Equation

Functionals of Brownian/non-Brownian motions have diverse applications and attracted a lot of interest of scientists. This paper focuses on deriving the forward and backward fractional Feynman-Kac equations describing the distribution of the functionals of the space and time tempered anomalous diffusion, belonging to the continuous time random walk class. Several examples of the functionals are explicitly treated, including the occupation time in half-space, the first passage time, the maximal displacement, the fluctuations of the occupation fraction, and the fluctuations of the time-averaged position.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal Fluctuations of Single-Particle Diffusivity in Quenched Environment

Local diffusion coefficients in disordered materials such as living cells are highly heterogeneous. Quenched disorder is utilized substantially to study such complex systems, whereas its analytical treatment is difficult to handle. We consider finite systems with quenched disorder in order to investigate the effects of sample disorder fluctuations and confinement on single-particle diffusivity. While the system is ergodic in a single disorder realization, the time-averaged mean squared displacement depends on the disorder, i.e., the system is ergodic but non-self-averaging. We find that the inverse Levy distribution is a universal distribution for diffusivity in the sense that it can be applied for arbitrary dimensions. Quantifying the degree of the non-self-averaging effect, we show that fluctuations of single-particle diffusivity far exceed the corresponding annealed theory and also find confinement effects. The relevance for experimental situations is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuations around equilibrium laws in ergodic continuous-time random walks

We study occupation time statistics in ergodic continuous-time random walks. Under thermal detailed balance conditions, the average occupation time is given by the Boltzmann-Gibbs canonical law. But close to the non-ergodic phase, the finite-time fluctuations around this mean are large and nontrivial. They exhibit dual time scaling and distribution laws: the infinite density of large fluctuations complements the Lévy-stable density of bulk fluctuations. Neither of the two should be interpreted as a stand-alone limiting law, as each has its own deficiency: the infinite density has an infinite norm (despite particle conservation), while the stable distribution has an infinite variance (although occupation times are bounded). These unphysical divergences are remedied by consistent use and interpretation of both formulas. Interestingly, while the system's canonical equilibrium laws naturally determine the mean occupation time of the ergodic motion, they also control the infinite and Lévy-stable densities of fluctuations. The duality of stable and infinite densities is in fact ubiquitous for these dynamics, as it concerns the time averages of general physical observables.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

1/f noise for intermittent quantum dots exhibits non-stationarity and critical exponents

The power spectrum of quantum dot fluorescence exhibits $1/f^β$ noise, related to the intermittency of these nanosystems. As in other systems exhibiting $1/f$ noise, this power spectrum is not integrable at low frequencies, which appears to imply infinite total power. We report measurements of individual quantum dots that address this long-standing paradox. We find that the level of $1/f^β$ noise decays with the observation time. The change of the spectrum with time places a bound on the total power. These observations are in stark contrast with most measurements of noise in macroscopic systems which do not exhibit any evidence for non-stationarity. We show that the traditional description of the power spectrum with a single exponent $β$ is incomplete and three additional critical exponents characterize the dependence on experimental time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deviations from Boltzmann-Gibbs equilibrium in confined optical lattices

Cold atoms in dissipative optical lattices have long been known to exhibit anomalous kinetics due to an effective nonlinear friction force. Here we show that confining the spatial motion of the atoms will lead to an anomalous non-Boltzmann-Gibbs equilibrium state, with a power law tail at large energies. Only in the limit of deep optical lattices, do we regain the Boltzmann-Gibbs state. For strong confinement relative to the damping, we find an explicit expression for the equilibrium phase-space distribution, which generally differs from the canonical Boltzmann-Gibbs state at all energies. Both in the low and high energy limits, the equilibrium distribution is a function of the system's Hamiltonian. At intermediate energies, however, the distribution is not a function of energy only and equipartition is violated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fractional Edgeworth Expansion: Corrections to the Gaussian-Lévy Central Limit Theorem

In this article we generalize the classical Edgeworth expansion for the probability density function (PDF) of sums of a finite number of symmetric independent identically distributed random variables with a finite variance to sums of variables with an infinite variance which converge by the generalized central limit theorem to a Lévy $α$-stable density function. Our correction may be written by means of a series of fractional derivatives of the Lévy and the conjugate Lévy PDFs. This series expansion is general and applies also to the Gaussian regime. To describe the terms in the series expansion, we introduce a new family of special functions and briefly discuss their properties. We implement our generalization to the distribution of the momentum for atoms undergoing Sisyphus cooling, and show the improvement of our leading order approximation compared to previous approximations. In vicinity of the transition between Lévy and Gauss behaviors, convergence to asymptotic results slows down.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Infinite densities for Lévy walks

Motion of particles in many systems exhibits a mixture between periods of random diffusive like events and ballistic like motion. In many cases, such systems exhibit strong anomalous diffusion, where low order moments $< |x(t)|^q >$ with $q$ below a critical value $q_c$ exhibit diffusive scaling while for $q>q_c$ a ballistic scaling emerges. The mixed dynamics constitutes a theoretical challenge since it does not fall into a unique category of motion, e.g., the known diffusion equations and central limit theorems fail to describe both aspects. In this paper we resolve this problem by resorting to the concept of infinite density. Using the widely applicable Lévy walk model, we find a general expression for the corresponding non-normalized density which is fully determined by the particles velocity distribution, the anomalous diffusion exponent $α$ and the diffusion coefficient $K_α$. We explain how infinite densities play a central role in the description of dynamics of a large class of physical processes and discuss how they can be evaluated from experimental or numerical data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Localization and universal fluctuations in ultraslow diffusion processes

We study ultraslow diffusion processes with logarithmic mean squared displacement (MSD) $\langle x^2(t)\rangle\simeq\log^γt$. Comparison of annealed continuous time random walks (CTRWs) with logarithmic waiting time distribution $ψ(τ)\simeq1/(τ\log^{1+γ}τ)$ and Sinai diffusion in quenched random landscapes shows striking similarities, despite their very different physical nature. In particular, they exhibit a weakly non-ergodic disparity of the time and ensemble averaged MSDs. Remarkably, for the CTRW we observe that the fluctuations of time averages become universal with an exponential suppression of mobile trajectories. We discuss the fundamental connection between the Golosov localization effect and non-ergodicity in the sense of the disparity between ensemble and time averaged MSD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical algorithms for the forward and backward fractional Feynman-Kac equations

The Feynman-Kac equations are a type of partial differential equations describing the distribution of functionals of diffusive motion. The probability density function (PDF) of Brownian functionals satisfies the Feynman-Kac formula, being a Schrödinger equation in imaginary time. The functionals of no-Brownian motion, or anomalous diffusion, follow the fractional Feynman-Kac equation [J. Stat. Phys. 141, 1071-1092, 2010], where the fractional substantial derivative is involved. Based on recently developed discretized schemes for fractional substantial derivatives [arXiv:1310.3086], this paper focuses on providing algorithms for numerically solving the forward and backward fractional Feynman-Kac equations; since the fractional substantial derivative is non-local time-space coupled operator, new challenges are introduced comparing with the general fractional derivative. Two ways (finite difference and finite element) of discretizing the space derivative are considered. For the backward fractional Feynman-Kac equation, the numerical stability and convergence of the algorithms with first order accuracy are theoretically discussed; and the optimal estimates are obtained. For all the provided schemes, including the first order and high order ones, of both forward and backward Feynman-Kac equations, extensive numerical experiments are performed to show their effectiveness.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Occupation times on a comb with ramified teeth

We investigate occupation time statistics for random walks on a comb with ramified teeth. This is achieved through the relation between the occupation time and the first passage times. Statistics of occupation times in half space follows Lamperti's distribution, i.e. the generalized arcsine law holds. Transitions between different behaviors are observed, which are controlled by the size of the backbone and teeth of the comb, as well as bias. Occupation time on a non-simply connected domain is analyzed with a mean-field theory and numerical simulations. In that case, the generalized arcsine law isn't valid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Distribution of the Area under a Bessel Excursion and its Moments

A Bessel excursion is a Bessel process that begins at the origin and first returns there at some given time $T$. We study the distribution of the area under such an excursion, which recently found application in the context of laser cooling. The area $A$ scales with the time as $A \sim T^{3/2}$, independent of the dimension, $d$, but the functional form of the distribution does depend on $d$. We demonstrate that for $d=1$, the distribution reduces as expected to the distribution for the area under a Brownian excursion, known as the Airy distribution, deriving a new expression for the Airy distribution in the process. We show that the distribution is symmetric in $d-2$, with nonanalytic behavior at $d=2$. We calculate the first and second moments of the distribution, as well as a particular fractional moment. We also analyze the analytic continuation from $d<2$ to $d>2$. In the limit where $d\to 4$ from below, this analytically continued distribution is described by a one-sided Lévy $α$-stable distribution with index $2/3$ and a scale factor proportional to $[(4-d)T]^{3/2}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Aging renewal theory and application to random walks

The versatility of renewal theory is owed to its abstract formulation. Renewals can be interpreted as steps of a random walk, switching events in two-state models, domain crossings of a random motion, etc. We here discuss a renewal process in which successive events are separated by scale-free waiting time periods. Among other ubiquitous long time properties, this process exhibits aging: events counted initially in a time interval [0,t] statistically strongly differ from those observed at later times [t_a,t_a+t]. In complex, disordered media, processes with scale-free waiting times play a particularly prominent role. We set up a unified analytical foundation for such anomalous dynamics by discussing in detail the distribution of the aging renewal process. We analyze its half-discrete, half-continuous nature and study its aging time evolution. These results are readily used to discuss a scale-free anomalous diffusion process, the continuous time random walk. By this we not only shed light on the profound origins of its characteristic features, such as weak ergodicity breaking. Along the way, we also add an extended discussion on aging effects. In particular, we find that the aging behavior of time and ensemble averages is conceptually very distinct, but their time scaling is identical at high ages. Finally, we show how more complex motion models are readily constructed on the basis of aging renewal dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Noisy continuous time random walks

Experimental studies of the diffusion of biomolecules in the environment of biological cells are routinely confronted with multiple sources of stochasticity, whose identification renders the detailed data analysis of single molecule trajectories quite intricate. Here we consider subdiffusive continuous time random walks, that represent a seminal model for the anomalous diffusion of tracer particles in complex environments. This motion is characterized by multiple trapping events with infinite mean sojourn time. In real physical situations, however, instead of the full immobilization predicted by the continuous time random walk model, the motion of the tracer particle shows additional jiggling, for instance, due to thermal agitation of the environment. We here present and analyze in detail an extension of the continuous time random walk model. Superimposing the multiple trapping behavior with additive Gaussian noise of variable strength, we demonstrate that the resulting process exhibits a rich variety of apparent dynamic regimes. In particular, such noisy continuous time random walks may appear ergodic while the naked continuous time random walk exhibits weak ergodicity breaking. Detailed knowledge of this behavior will be useful for the truthful physical analysis of experimentally observed subdiffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical estimate of infinite invariant densities: application to Pesin-type identity

Weakly chaotic maps with unstable fixed points are investigated in the regime where the invariant density is non-normalizable. We propose that the infinite invariant density of these maps can be estimated using as the long time limit of t^(1-alpha) rho(x, t), in agreement with earlier work of Thaler. Here rho(x, t) is the normalizable density of particles. This definition uniquely determines the infinite density and is a valuable tool for numerical estimations. We use this density to estimate the subexponential separation lambda_alpha of nearby trajectories. For a particular map introduced by Thaler we use an analytical expression for the infinite invariant density to calculate lambda_alpha exactly, which perfectly matches simulations without fitting. Misunderstanding which recently appeared in the literature is removed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Time averaged Einstein relation and fluctuating diffusivities for the Lévy walk

The Lévy walk model is a stochastic framework of enhanced diffusion with many applications in physics and biology. Here we investigate the time averaged mean squared displacement $\bar{δ^2}$ often used to analyze single particle tracking experiments. The ballistic phase of the motion is non-ergodic and we obtain analytical expressions for the fluctuations of $\bar{δ^2}$. For enhanced sub-ballistic diffusion we observe numerically apparent ergodicity breaking on long time scales. As observed by Akimoto \textit{Phys. Rev. Lett.} \textbf{108}, 164101 (2012) deviations of temporal averages $\bar{δ^2}$ from the ensemble average $< x^2 >$ depend on the initial preparation of the system, and here we quantify this discrepancy from normal diffusive behavior. Time averaged response to a bias is considered and the resultant generalized Einstein relations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transport and the first passage time problem with application to cold atoms in optical traps

Measurements of spatial diffusion of cold atoms in optical lattices have revealed anomalous super-diffusion, which is controlled by the depth of the optical lattice. We use first passage time statistics to derive the diffusion front of the atoms. In particular, the distributions of areas swept under the first passage curve till its first arrival, and of areas under the Bessel excursion are shown to be powerful tools in the analysis of the atomic cloud. A rather general relation between first passage time statistics and diffusivity is discussed, showing that first passage time analysis is a useful tool in the calculation of transport coefficients. A brief introduction to the semi-classical description of Sisyphus cooling is provided which yields a rich phase diagram for the dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ageing effects in single particle trajectory averages

We study time averages of single particle trajectories in scale free anomalous diffusion processes, in which the measurement starts at some time t_a>0 after initiation of the process at the time origin, t=0. Using ageing renewal theory we show that for such non-stationary processes a large class of observables are affected by a unique ageing function, which is independent of boundary conditions or the external forces. We quantify the weakly non-ergodic nature of this process in terms of the distribution of time averages and the ergodicity breaking parameter which both explicitly depend on the ageing time t_a. Consequences for the interpretation of single particle tracking data are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Aging generates regular motions in weakly chaotic systems

Using intermittent maps with infinite invariant measures, we investigate the universality of time-averaged observables under aging conditions. According to Aaronson-Darling-Kac theorem, in non-aged dynamical systems with infinite invariant measures, the distribution of the normalized time averages of integrable functions converge to the Mittag-Leffler distribution. This well known theorem holds when the start of observations coincides with the start of the dynamical processes. Introducing a concept of the aging limit where the aging time $t_a$ and the total measurement time $t$ goes to infinity while the aging ratio $t_a/t$ is a constant, we obtain a novel distributional limit theorem of time-averaged observables integrable with respect to the infinite invariant density. Applying the theorem to the Lyapunov exponent in intermittent maps, we find that regular motions and a weakly chaotic behavior coexist in the aging limit. This mixed type of dynamics is controlled by the aging ratio and hence is very different from the usual scenario of regular and chaotic motions in Hamiltonian systems. The probability of finding regular motions in non-aged processes is zero, while in the aging regime it is finite and it increases when system ages.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuations of $1/f$ noise and the low frequency cutoff paradox

Recent experiments on blinking quantum dots and weak turbulence in liquid crystals reveal the fundamental connection between $1/f$ noise and power law intermittency. The non-stationarity of the process implies that the power spectrum is random -- a manifestation of weak ergodicity breaking. Here we obtain the universal distribution of the power spectrum, which can be used to identify intermittency as the source of the noise. We solve an outstanding paradox on the non integrability of $1/f$ noise and the violation of Parseval's theorem. We explain why there is no physical low frequency cutoff and therefore cannot be found in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theory of fractional-Lévy kinetics for cold atoms diffusing in optical lattices

Recently, anomalous superdiffusion of ultra cold 87Rb atoms in an optical lattice has been observed along with a fat-tailed, Lévy type, spatial distribution. The anomalous exponents were found to depend on the depth of the optical potential. We find, within the framework of the semiclassical theory of Sisyphus cooling, three distinct phases of the dynamics, as the optical potential depth is lowered: normal diffusion; Lévy diffusion; and x ~ t^3/2 scaling, the latter related to Obukhov's model (1959) of turbulence. The process can be formulated as a Lévy walk, with strong correlations between the length and duration of the excursions. We derive a fractional diffusion equation describing the atomic cloud, and the corresponding anomalous diffusion coefficient.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak subordination breaking for the quenched trap model

We map the problem of diffusion in the quenched trap model onto a new stochastic process: Brownian motion which is terminated at the coverage "time" ${\cal S}_α=\sum_{x=-\infty} ^\infty (n_x)^α$ with $n_x$ being the number of visits to site $x$. Here $0<α=T/T_g<1$ is a measure of the disorder in the original model. This mapping allows us to treat the intricate correlations in the underlying random walk in the random environment. The operational "time" ${\cal S}_α$ is changed to laboratory time $t$ with a Lévy time transformation. Investigation of Brownian motion stopped at "time" ${\cal S}_α$ yields the diffusion front of the quenched trap model which is favorably compared with numerical simulations. In the zero temperature limit of $α\to 0$ we recover the renormalization group solution obtained by C. Monthus. Our theory surmounts critical slowing down which is found when $α\to 1$. Above the critical dimension two mapping the problem to a continuous time random walk becomes feasible though still not trivial.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A fractional Feynman-Kac equation for weak ergodicity breaking

Continuous-time random walk (CTRW) is a model of anomalous sub-diffusion in which particles are immobilized for random times between successive jumps. A power-law distribution of the waiting times, $ψ(τ) τ^{-(1+α)}$, leads to sub-diffusion ($<x^2>~t^α$) for 0<α<1. In closed systems, the long stagnation periods cause time-averages to divert from the corresponding ensemble averages, which is a manifestation of weak ergodicity breaking. The time-average of a general observable $\bar{U} = \int_0^t U[x(τ)]dτ/ t$ is a functional of the path and is described by the well known Feynman-Kac equation if the motion is Brownian. Here, we derive forward and backward fractional Feynman-Kac equations for functionals of CTRW in a binding potential. We use our equations to study two specific time-averages: the fraction of time spent by a particle in half box, and the time-average of the particle's position in a harmonic field. In both cases, we obtain the probability density function of the time-averages for $t \rightarrow \infty$ and the first two moments. Our results show that both the occupation fraction and the time-averaged position are random variables even for long-times, except for α=1 when they are identical to their ensemble averages. Using the fractional Feynman-Kac equation, we also study the dynamics leading to weak ergodicity breaking, namely the convergence of the fluctuations to their asymptotic values.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anomalous Infiltration

Infiltration of anomalously diffusing particles from one material to another through a biased interface is studied using continuous time random walk and Levy walk approaches. Subdiffusion in both systems may lead to a net drift from one material to another (e.g. <x(t)> > 0) even if particles eventually flow in the opposite direction (e.g. number of particles in x>0 approaches zero). A weaker paradox is found for a symmetric interface: a flow of particles is observed while the net drift is zero. For a subdiffusive sample coupled to a superdiffusive system we calculate the average occupation fractions and the scaling of the particles distribution. We find a net drift in this system, which is always directed to the superdiffusive material, while the particles flow to the material with smaller sub or superdiffusion exponent. We report the exponents of the first passage times distribution of Levy walks, which are needed for the calculation of anomalous infiltration

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

In vivo anomalous diffusion and weak ergodicity breaking of lipid granules

Combining extensive single particle tracking microscopy data of endogenous lipid granules in living fission yeast cells with analytical results we show evidence for anomalous diffusion and weak ergodicity breaking. Namely we demonstrate that at short times the granules perform subdiffusion according to the laws of continuous time random walk theory. The associated violation of ergodicity leads to a characteristic turnover between two scaling regimes of the time averaged mean squared displacement. At longer times the granule motion is consistent with fractional Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Foundation of Fractional Langevin Equation: Harmonization of a Many Body Problem

In this study we derive a single-particle equation of motion, from first-principles, starting out with a microscopic description of a tracer particle in a one-dimensional many-particle system with a general two-body interaction potential. Using a new harmonization technique, we show that the resulting dynamical equation belongs to the class of fractional Langevin equations, a stochastic framework which has been proposed in a large body of works as a means of describing anomalous dynamics. Our work sheds light on the fundamental assumptions of these phenomenological models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Infinite invariant densities for anomalous diffusion in optical lattices and other logarithmic potentials

We solve the Fokker-Planck equation for Brownian motion in a logarithmic potential. When the diffusion constant is below a critical value the solution approaches a non-normalizable scaling state, reminiscent of an infinite invariant density. With this non-normalizable density we obtain the phase diagram of anomalous diffusion for this important process. We briefly discuss the consequence for a range of physical systems including atoms in optical lattices and charges in vicinity of long polyelectrolytes. Our work explains in what sense the infinite invariant density and not Boltzmann's equilibrium describes the long time limit of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On distributions of functionals of anomalous diffusion paths

Functionals of Brownian motion have diverse applications in physics, mathematics, and other fields. The probability density function (PDF) of Brownian functionals satisfies the Feynman-Kac formula, which is a Schrodinger equation in imaginary time. In recent years there is a growing interest in particular functionals of non-Brownian motion, or anomalous diffusion, but no equation existed for their PDF. Here, we derive a fractional generalization of the Feynman-Kac equation for functionals of anomalous paths based on sub-diffusive continuous-time random walk. We also derive a backward equation and a generalization to Levy flights. Solutions are presented for a wide number of applications including the occupation time in half space and in an interval, the first passage time, the maximal displacement, and the hitting probability. We briefly discuss other fractional Schrodinger equations that recently appeared in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Paradoxes of Subdiffusive Infiltration in Disordered Systems

Infiltration of diffusing particles from one material to another where the diffusion mechanism is either normal or anomalous is a widely observed phenomena. When the diffusion is anomalous we find interesting behaviors: diffusion may lead to an averaged net drift <x> from one material to another even if all particles eventually flow in the opposite direction, or may lead to a flow without drift. Starting with an underlying continuous time random walk model we solve diffusion equations describing this problem. Similar drift against flow is found in the quenched trap model. We argue that such a behavior is a general feature of diffusion in disordered systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Separation of trajectories and its Relation to Entropy for Intermittent Systems with a Zero Lyapunov exponent

One dimensional intermittent maps with stretched exponential separation of nearby trajectories are considered. When time goes infinity the standard Lyapunov exponent is zero. We investigate the distribution of $λ_α= \sum_{i=0}^{t-1} \ln \left| M'(x_i) \right|/t^α$, where $α$ is determined by the nonlinearity of the map in the vicinity of marginally unstable fixed points. The mean of $λ_α$ is determined by the infinite invariant density. Using semi analytical arguments we calculate the infinite invariant density for the Pomeau-Manneville map, and with it obtain excellent agreement between numerical simulation and theory. We show that $α\left< λ_α\right>$ is equal to Krengel's entropy and to the complexity calculated by the Lempel-Ziv compression algorithm. This generalized Pesin's identity shows that $\left< λ_α\right>$ and Krengel's entropy are the natural generalizations of usual Lyapunov exponent and entropy for these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

From non-Brownian Functionals to a Fractional Schrödinger Equation

We derive backward and forward fractional Schrödinger type of equations for the distribution of functionals of the path of a particle undergoing anomalous diffusion. Fractional substantial derivatives introduced by Friedrich and co-workers [PRL {\bf 96}, 230601 (2006)] provide the correct fractional framework for the problem at hand. In the limit of normal diffusion we recover the Feynman-Kac treatment of Brownian functionals. For applications, we calculate the distribution of occupation times in half space and show how statistics of anomalous functionals is related to weak ergodicity breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Ergodic Properties of Fractional Brownian-Langevin Motion

We investigate the time average mean square displacement $\overline{δ^2}(x(t))=\int_0^{t-Δ}[x(t^\prime+Δ)-x(t^\prime)]^2 dt^\prime/(t-Δ)$ for fractional Brownian and Langevin motion. Unlike the previously investigated continuous time random walk model $\overline{δ^2}$ converges to the ensemble average $<x^2 > \sim t^{2 H}$ in the long measurement time limit. The convergence to ergodic behavior is however slow, and surprisingly the Hurst exponent $H=3/4$ marks the critical point of the speed of convergence. When $H<3/4$, the ergodicity breaking parameter ${EB} = {Var} (\overline{δ^2}) / < \overline{δ^2} >^2\sim k(H) \cdotΔ\cdot t^{-1}$, when $H=3/4$, ${EB} \sim (9/16)(\ln t) \cdotΔ\cdot t^{-1}$, and when $3/4<H <1, {EB} \sim k(H)Δ^{4-4H} t^{4H-4}$. In the ballistic limit $H \to 1$ ergodicity is broken and ${EB} \sim 2$. The critical point $H=3/4$ is marked by the divergence of the coefficient $k(H)$. Fractional Brownian motion as a model for recent experiments of sub-diffusion of mRNA in the cell is briefly discussed and comparison with the continuous time random walk model is made.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Residence Time Statistics for Normal and Fractional Diffusion in a Force Field

We investigate statistics of occupation times for an over-damped Brownian particle in an external force field. A backward Fokker-Planck equation introduced by Majumdar and Comtet describing the distribution of occupation times is solved. The solution gives a general relation between occupation time statistics and probability currents which are found from solutions of the corresponding problem of first passage time. This general relationship between occupation times and first passage times, is valid for normal Markovian diffusion and for non-Markovian sub-diffusion, the latter modeled using the fractional Fokker-Planck equation. For binding potential fields we find in the long time limit ergodic behavior for normal diffusion, while for the fractional framework weak ergodicity breaking is found, in agreement with previous results of Bel and Barkai on the continuous time random walk on a lattice. For non-binding potential rich physical behaviors are obtained, and classification of occupation time statistics is made possible according to whether or not the underlying random walk is recurrent and the averaged first return time to the origin is finite. Our work establishes a link between fractional calculus and ergodicity breaking.

0 citations0 reviewsNo code linkedOpen access