Researcher profile

Eli Barkai

Eli Barkai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.