Showing 1–50 of 161 results for author: Barkai, E

A tale of three approaches: dynamical phase transitions for weakly bound Brownian particles

Authors: Lucianno Defaveri, Eli Barkai, David A. Kessler

Abstract: We investigate a system of Brownian particles weakly bound by attractive parity-symmetric potentials that grow at large distances as $V(x) \sim |x|^α$, with $0 < α< 1$. The probability density function $P(x,t)$ at long times reaches the Boltzmann-Gibbs equilibrium state, with all moments finite. However, the system's relaxation is not exponential, as is usual for a confining system with a well-def… ▽ More We investigate a system of Brownian particles weakly bound by attractive parity-symmetric potentials that grow at large distances as $V(x) \sim |x|^α$, with $0 < α< 1$. The probability density function $P(x,t)$ at long times reaches the Boltzmann-Gibbs equilibrium state, with all moments finite. However, the system's relaxation is not exponential, as is usual for a confining system with a well-defined equilibrium, but instead follows a stretched exponential $e^{- \mathrm{const} \, t^ν}$ with exponent $ν=α/(2+α)$. This problem is studied from three perspectives. First, we propose a straightforward and general scaling rate-function solution for $P(x,t)$. This rate-function, which is an important tool from large deviation theory, also displays anomalous time scaling and a dynamical phase transition. Second, through the eigenfunctions of the Fokker-Planck operator, we obtain, using the WKB method, more complete solutions that reproduce the rate function approach. Finally, we show how the alternative path-integral formalism allows us to recover the same results, with the above rate-function being the solution of the classical Hamilton-Jacobi equation describing the most probable path. Properties such as parity, the role of initial conditions, and the dynamical phase transition are thoroughly studied in all three approaches. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Comments: 42 pages, 11 figures, 1 table

Directed Motion and Spatial Coherence in the Cell Nucleus

Authors: M. Hidalgo-Soria, Y. Haddad, E. Barkai, Y. Garini, S. Burov

Abstract: Investigating the dynamics of chromatin and the factors that are affecting it, has provided valuable insights into the organization and functionality of the genome in the cell nucleus. We control the expression of Lamin-A, an important organizer protein of the chromatin and nucleus structure. By simultaneously tracking tens of chromosomal loci (telomeres) in each nucleus, we find that the motion o… ▽ More Investigating the dynamics of chromatin and the factors that are affecting it, has provided valuable insights into the organization and functionality of the genome in the cell nucleus. We control the expression of Lamin-A, an important organizer protein of the chromatin and nucleus structure. By simultaneously tracking tens of chromosomal loci (telomeres) in each nucleus, we find that the motion of chromosomal loci in Lamin-A depleted cells is both faster and more directed on a scale of a few micrometers, which coincides with the size of chromosome territories. Moreover, in the absence of Lamin-A we reveal the existence of correlations among neighboring telomeres. We show how these pairwise correlations are linked with the intermittent and persistent character of telomere trajectories, underscoring the importance of Lamin-A protein in chromosomal organization. △ Less

Submitted 11 July, 2024; originally announced July 2024.

First Hitting Times on a Quantum Computer: Tracking vs. Local Monitoring, Topological Effects, and Dark States

Authors: Qingyuan Wang, Silin Ren, Ruoyu Yin, Klaus Ziegler, Eli Barkai, Sabine Tornow

Abstract: We investigate a quantum walk on a ring represented by a directed triangle graph with complex edge weights and monitored at a constant rate until the quantum walker is detected. To this end, the first hitting time statistics is recorded using unitary dynamics interspersed stroboscopically by measurements, which is implemented on IBM quantum computers with a midcircuit readout option. Unlike classi… ▽ More We investigate a quantum walk on a ring represented by a directed triangle graph with complex edge weights and monitored at a constant rate until the quantum walker is detected. To this end, the first hitting time statistics is recorded using unitary dynamics interspersed stroboscopically by measurements, which is implemented on IBM quantum computers with a midcircuit readout option. Unlike classical hitting times, the statistical aspect of the problem depends on the way we construct the measured path, an effect that we quantify experimentally. First, we experimentally verify the theoretical prediction that the mean return time to a target state is quantized, with abrupt discontinuities found for specific sampling times and other control parameters, which has a well-known topological interpretation. Second, depending on the initial state, system parameters, and measurement protocol, the detection probability can be less than one or even zero, which is related to dark-state physics. Both, return-time quantization and the appearance of the dark states are related to degeneracies in the eigenvalues of the unitary time evolution operator. We conclude that, for the IBM quantum computer under study, the first hitting times of monitored quantum walks are resilient to noise. Yet, a finite number of measurements leads to broadening effects, which modify the topological quantization and chiral effects of the asymptotic theory with an infinite number of measurements. Our results point the way for the development of novel quantum walk algorithms that exploit measurement-induced effects on quantum computers. △ Less

Submitted 10 April, 2024; v1 submitted 24 February, 2024; originally announced February 2024.

Laplace's first law of errors applied to diffusive motion

Authors: Omer Hamdi, Stanislav Burov, Eli Barkai

Abstract: In biological, glassy, and active systems, various tracers exhibit Laplace-like, i.e., exponential, spreading of the diffusing packet of particles. The limitations of the central limit theorem in fully capturing the behaviors of such diffusive processes, especially in the tails, have been studied using the continuous time random walk model. For cases when the jump length distribution is super-expo… ▽ More In biological, glassy, and active systems, various tracers exhibit Laplace-like, i.e., exponential, spreading of the diffusing packet of particles. The limitations of the central limit theorem in fully capturing the behaviors of such diffusive processes, especially in the tails, have been studied using the continuous time random walk model. For cases when the jump length distribution is super-exponential, e.g., a Gaussian, we use large deviations theory and relate it to the appearance of exponential tails. When the jump length distribution is sub-exponential the packet of spreading particles is described by the big jump principle. We demonstrate the applicability of our approach for finite time, indicating that rare events and the asymptotics of the large deviations rate function can be sampled for large length scales within a reasonably short measurement time. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Comments: 14 pages, 8 figures

Properties of Fractionally Quantized Recurrence Times for Interacting Spin Models

Authors: Quancheng Liu, David A. Kessler, Eli Barkai

Abstract: Recurrence time quantifies the duration required for a physical system to return to its initial state, playing a pivotal role in understanding the predictability of complex systems. In quantum systems with subspace measurements, recurrence times are governed by Anandan-Aharonov phases, yielding fractionally quantized recurrence times. However, the fractional quantization phenomenon in interacting… ▽ More Recurrence time quantifies the duration required for a physical system to return to its initial state, playing a pivotal role in understanding the predictability of complex systems. In quantum systems with subspace measurements, recurrence times are governed by Anandan-Aharonov phases, yielding fractionally quantized recurrence times. However, the fractional quantization phenomenon in interacting quantum systems remains poorly explored. Here, we address this gap by establishing universal lower and upper bounds for recurrence times in interacting spins. Notably, we investigate scenarios where these bounds are approached, shedding light on the speed of quantum processes under monitoring. In specific cases, our findings reveal that the complex many-body system can be effectively mapped onto a dynamical system with a single quasi-particle, leading to the discovery of integer quantized recurrence times. Our research yields a valuable link between recurrence times and the number of dark states in the system, thus providing a deeper understanding of the intricate interplay between quantum recurrence, measurements, and interaction effects. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 5+5 pages, 3+3 figures

Restart uncertainty relation for monitored quantum dynamics

Authors: Ruoyu Yin, Qingyuan Wang, Sabine Tornow, Eli Barkai

Abstract: We introduce a novel time-energy uncertainty relationship within the context of restarts in monitored quantum dynamics. Initially, we investigate the concept of ``first hitting time'' in quantum systems using an IBM quantum computer and a three-site ring graph as our starting point. Previous studies have established that the mean recurrence time, which represents the time taken to return to the in… ▽ More We introduce a novel time-energy uncertainty relationship within the context of restarts in monitored quantum dynamics. Initially, we investigate the concept of ``first hitting time'' in quantum systems using an IBM quantum computer and a three-site ring graph as our starting point. Previous studies have established that the mean recurrence time, which represents the time taken to return to the initial state, is quantized as an integer multiple of the sampling time, displaying pointwise discontinuous transitions at resonances. Our findings demonstrate that, the natural utilization of the restart mechanism in laboratory experiments, driven by finite data collection time spans, leads to a broadening effect on the transitions of the mean recurrence time. Our newly proposed uncertainty relation captures the underlying essence of these phenomena, by connecting the broadening of the mean hitting time near resonances, to the intrinsic energies of the quantum system and to the fluctuations of recurrence time. This work not only contributes to our understanding of fundamental aspects related to quantum measurements and dynamics, but also offers practical insights for the design of efficient quantum algorithms with mid-circuit measurements. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 10 pages, 6 figures

First passage times in compact domains exhibits bi-scaling

Authors: Talia Baravi, Eli Barkai

Abstract: The study of first passage times for diffusing particles reaching target states is foundational in various practical applications, including diffusion-controlled reactions. In this work, we present a bi-scaling theory for the probability density function of first passage times in confined compact processes, applicable to both Euclidean and Fractal domains, diverse geometries, and scenarios with or… ▽ More The study of first passage times for diffusing particles reaching target states is foundational in various practical applications, including diffusion-controlled reactions. In this work, we present a bi-scaling theory for the probability density function of first passage times in confined compact processes, applicable to both Euclidean and Fractal domains, diverse geometries, and scenarios with or without external force fields, accommodating Markovian and semi-Markovian random walks. In large systems, first passage time statistics exhibit a bi-scaling behavior, challenging the use of a single time scale. Our theory employs two distinct scaling functions: one for short times, capturing initial dynamics in unbounded systems, and the other for long times is sensitive to finite size effects. The combined framework provides a complete expression for first passage time statistics across all time scales. △ Less

Submitted 23 November, 2023; originally announced November 2023.

Comments: 7 pages (14 with supplementary material), 3 figures

Stretched-exponential relaxation in weakly-confined Brownian systems through large deviation theory

Authors: Lucianno Defaveri, Eli Barkai, David A. Kessler

Abstract: 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 a $r^{2/3}$ potential. Here, we study a Brownian particle under the influence of a general confining, albeit weak, potential field that gro… ▽ More 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 a $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 sub-linear 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. △ Less

Submitted 20 February, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 6 pages (12 with SM), 3 figures

Fractional Advection Diffusion Asymmetry Equation, derivation, solution and application

Authors: Wanli Wang, Eli Barkai

Abstract: The non-Markovian continuous-time random walk model, featuring fat-tailed waiting times and narrow distributed displacements with a non-zero mean, is a well studied model for anomalous diffusion. Using an analytical approach, we recently demonstrated how a fractional space advection diffusion asymmetry equation, usually associated with Markovian Lévy flights, describes the spreading of a packet of… ▽ More The non-Markovian continuous-time random walk model, featuring fat-tailed waiting times and narrow distributed displacements with a non-zero mean, is a well studied model for anomalous diffusion. Using an analytical approach, we recently demonstrated how a fractional space advection diffusion asymmetry equation, usually associated with Markovian Lévy flights, describes the spreading of a packet of particles. Since we use Gaussian statistics for jump lengths though fat-tailed distribution of waiting times, the appearance of fractional space derivatives in the kinetic equation demands explanations provided in this manuscript. As applications we analyse the spreading of tracers in two dimensions, breakthrough curves investigated in the field of contamination spreading in hydrology and first passage time statistics. We present a subordination scheme valid for the case when the mean waiting time is finite and the variance diverges, which is related to Lévy statistics for the number of renewals in the process. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 31pp

Statistics of Long-Range Force Fields in Random Environments: Beyond Holtsmark

Authors: Avraham Samama, Eli Barkai

Abstract: 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 desc… ▽ More 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 non-physical, 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 bi-scaling, 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 non-normalized similar to infinite invariant densities found in dynamical systems. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 9 pages 2 figures

Ergodic properties of Brownian motion under stochastic resetting

Authors: Eli Barkai, Rosa Flaquer-Galmes, Vicenç Méndez

Abstract: We study ergodic properties of one-dimensional Brownian motion with resetting. Using generic classes of statistics of times between resets, we find respectively for thin/fat tailed distributions, the normalized/non-normalised invariant density of this process. The former case corresponds to known results in the resetting literature and the latter to infinite ergodic theory. Two types of ergodic tr… ▽ More We study ergodic properties of one-dimensional Brownian motion with resetting. Using generic classes of statistics of times between resets, we find respectively for thin/fat tailed distributions, the normalized/non-normalised invariant density of this process. The former case corresponds to known results in the resetting literature and the latter to infinite ergodic theory. Two types of ergodic transitions are found in this system. The first is when the mean waiting time between resets diverges, when standard ergodic theory switches to infinite ergodic theory. The second is when the mean of the square root of time between resets diverges and the properties of the invariant density are drastically modified. We then find a fractional integral equation describing the density of particles. This finite time tool is particularly useful close to the ergodic transition where convergence to asymptotic limits is logarithmically slow. Our study implies rich ergodic behaviors for this non-equilibrium process which should hold far beyond the case of Brownian motion analyzed here. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Non-normalizable quasi-equilibrium states under fractional dynamics

Authors: Lucianno Defaveri, Maike A. F. dos Santos, David A. Kessler, Eli Barkai, Celia Anteneodo

Abstract: We study non-normalizable quasi-equilibrium states (NNQE) arising from anomalous diffusion. Initially, particles in contact with a thermal bath are released from an asymptotically flat potential well, with dynamics that is described by fractional calculus. For temperatures that are sufficiently low compared to the potential depth, the properties of the system remain almost constant in time. We use… ▽ More We study non-normalizable quasi-equilibrium states (NNQE) arising from anomalous diffusion. Initially, particles in contact with a thermal bath are released from an asymptotically flat potential well, with dynamics that is described by fractional calculus. For temperatures that are sufficiently low compared to the potential depth, the properties of the system remain almost constant in time. We use the fractional-time Fokker-Planck equation (FTFPE) and continuous-time random walk approaches to calculate the ensemble averages of observables. We obtain analytical estimates of the duration of NNQE, 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 time scales for stagnation depend exponentially on the activation energy in units of temperature multiplied by a function of the fractional exponent. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 9 pages, 6 figures

Journal ref: Phys. Rev. E 108, 024133 (2023)

Instability in the quantum restart problem

Authors: Ruoyu Yin, Qingyuan Wang, Eli Barkai

Abstract: Repeatedly-monitored quantum walks with a rate $1/τ$ yield discrete-time trajectories which are inherently random. With these paths the first-hitting time with sharp restart is studied. We find an instability in the optimal mean hitting time, which is not found in the corresponding classical random walk process. This instability implies that a small change in parameters can lead to a rather large… ▽ More Repeatedly-monitored quantum walks with a rate $1/τ$ yield discrete-time trajectories which are inherently random. With these paths the first-hitting time with sharp restart is studied. We find an instability in the optimal mean hitting time, which is not found in the corresponding classical random walk process. This instability implies that a small change in parameters can lead to a rather large change of the optimal restart time. We show that the optimal restart time versus $τ$, as a control parameter, exhibits sets of staircases and plunges. The plunges, are due to the mentioned instability, which in turn is related to the quantum oscillations of the first-hitting time probability, in the absence of restarts. Furthermore, we prove that there are only two patterns of staircase structures, dependent on the parity of the distance between the target and the source in units of lattice constant. The global minimum of the hitting time, is controlled not only by the restart time, as in classical problems, but also by the sampling time $τ$. We provide numerical evidence that this global minimum occurs for the $τ$ minimizing the mean hitting time, given restarts taking place after each measurement. Last but not least, we numerically show that the instability found in this work is relatively robust against stochastic perturbations in the sampling time $τ$. △ Less

Submitted 7 June, 2024; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: 17 pages (main text), 15 figures

Brownian particles in periodic potentials: coarse-graining versus fine structure

Authors: Lucianno Defaveri, Eli Barkai, David A. Kessler

Abstract: 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… ▽ More 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 $2 D^* 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 non-normalizable. 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. △ Less

Submitted 14 February, 2023; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: 17 pages, 8 figures

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

Authors: Stanislav Burov, Wanli Wang, Eli Barkai

Abstract: 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… ▽ More 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. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Controls that expedite first passage times in disordered systems

Authors: Marc Höll, Alon Nissan, Brian Berkowitz, Eli Barkai

Abstract: 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 employ the big jump principle that shows the dominance of th… ▽ More 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 employ the big jump principle that shows the dominance of the maximum trapping time on the first passage time. Inspired by the restart paradigm, 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. △ Less

Submitted 25 September, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Local equilibrium properties of ultraslow diffusion in the Sinai model

Authors: Amin Padash, Erez Aghion, Alexander Schulz, Eli Barkai, Aleksei V Chechkin, Ralf Metzler, Holger Kantz

Abstract: 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 s… ▽ More 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. △ Less

Submitted 21 July, 2022; v1 submitted 6 July, 2022; originally announced July 2022.

Journal ref: New J. Phys. 24 (2022) 073026

Restart expedites quantum walk hitting times

Authors: Ruoyu Yin, Eli Barkai

Abstract: Classical first-passage times under restart are used in a wide variety of models, yet the quantum version of the problem still misses key concepts. We study the quantum hitting time with restart using a monitored quantum walk. The restart strategy eliminates the problem of dark states, i.e. cases where the particle evades detection, while maintaining the ballistic propagation which is important fo… ▽ More Classical first-passage times under restart are used in a wide variety of models, yet the quantum version of the problem still misses key concepts. We study the quantum hitting time with restart using a monitored quantum walk. The restart strategy eliminates the problem of dark states, i.e. cases where the particle evades detection, while maintaining the ballistic propagation which is important for fast search. We find profound effects of quantum oscillations on the restart problem, namely a type of instability of the mean detection time, and optimal restart times that form staircases, with sudden drops as the rate of sampling is modified. In the absence of restart and in the Zeno limit, the detection of the walker is not possible and we examine how restart overcomes this well-known problem, showing that the optimal restart time becomes insensitive to the sampling period. △ Less

Submitted 27 December, 2022; v1 submitted 4 May, 2022; originally announced May 2022.

Comments: 13 pages, 11 figures

Designing exceptional-point-based graphs yielding topologically guaranteed quantum search

Authors: Quancheng Liu, David A. Kessler, Eli Barkai

Abstract: Quantum walks underlie an important class of quantum computing algorithms, and represent promising approaches in various simulations and practical applications. Here we design stroboscopically monitored quantum walks and their subsequent graphs that can naturally boost target searches. We show how to construct walks with the property that all the eigenvalues of the non-Hermitian survival operator,… ▽ More Quantum walks underlie an important class of quantum computing algorithms, and represent promising approaches in various simulations and practical applications. Here we design stroboscopically monitored quantum walks and their subsequent graphs that can naturally boost target searches. We show how to construct walks with the property that all the eigenvalues of the non-Hermitian survival operator, describing the mixed effects of unitary dynamics and the back-action of measurement, coalesce to zero, corresponding to an exceptional point whose degree is the size of the system. Generally, the resulting search is guaranteed to succeed in a bounded time for any initial condition, which is faster than classical random walks or quantum walks on typical graphs. We then show how this efficient quantum search is related to a quantized topological winding number and further discuss the connection of the problem to an effective massless Dirac particle. △ Less

Submitted 24 September, 2022; v1 submitted 7 February, 2022; originally announced February 2022.

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

Authors: Takuma Akimoto, Eli Barkai, Günter Radons

Abstract: 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 t… ▽ More 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. △ Less

Submitted 4 June, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures

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

Authors: Zachary R Fox, Eli Barkai, Diego Krapf

Abstract: 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 stat… ▽ More 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. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Journal ref: Nat Commun 12, 6162 (2021)

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

Authors: Gianmaria Falasco, Eli Barkai, Marco Baiesi

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: 7 pages, 3 figures

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

Authors: Eli Barkai, Günter Radons, Takuma Akimoto

Abstract: 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 densit… ▽ More 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. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 32 pages, 11 figures

Discrete Sampling of Extreme Events Modifies Their Statistics

Authors: Lior Zarfaty, Eli Barkai, David A. Kessler

Abstract: 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… ▽ More 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. △ Less

Submitted 26 August, 2022; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: Main text: 7 pages, 7 figures. Supplemental Material: 11 pages. This work is a short version of a still-in-preparation study (arXiv:2108.06778)

Journal ref: Phys. Rev. Lett. 129, 094101 (2022)

Measurement induced quantum walks

Authors: A. Didi, E. Barkai

Abstract: We investigate a tight binding quantum walk on a graph. Repeated stroboscopic measurements of the position of the particle yield a measured "trajectory", and a combination of classical and quantum mechanical properties for the walk are observed. We explore the effects of the measurements on the spreading of the packet on a one dimensional line, showing that except for the Zeno limit, the system co… ▽ More We investigate a tight binding quantum walk on a graph. Repeated stroboscopic measurements of the position of the particle yield a measured "trajectory", and a combination of classical and quantum mechanical properties for the walk are observed. We explore the effects of the measurements on the spreading of the packet on a one dimensional line, showing that except for the Zeno limit, the system converges to Gaussian statistics similarly to a classical random walk. A large deviation analysis and an Edgeworth expansion yield quantum corrections to this normal behavior. We then explore the first passage time to a target state using a generating function method, yielding properties like the quantization of the mean first return time. In particular, we study the effects of certain sampling rates which cause remarkable change in the behavior in the system, like divergence of the mean detection time in finite systems and a decomposition of the phase space into mutually exclusive regions, an effect that mimics ergodicity breaking, whose origin here is the destructive interference in quantum mechanics. For a quantum walk on a line we show that in our system the first detection probability decays classically like $(\text{time})^{-3/2}$, this is dramatically different compared to local measurements which yield a decay rate of $(\text{time})^{-3}$, indicating that the exponents of the first passage time depends on the type of measurements used. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 22 pages, 18 figures

Discrete sampling of correlated random variables modifies the long-time behavior of their extreme value statistics

Authors: Lior Zarfaty, Eli Barkai, David A. Kessler

Abstract: We consider the extreme value statistics of correlated random variables that arise from a Langevin equation. Recently, it was shown that the extreme values of the Ornstein-Uhlenbeck process follow a different distribution than those originating from its equilibrium measure, composed of independent and identically distributed Gaussian random variables. Here, we first focus on the discretely sampled… ▽ More We consider the extreme value statistics of correlated random variables that arise from a Langevin equation. Recently, it was shown that the extreme values of the Ornstein-Uhlenbeck process follow a different distribution than those originating from its equilibrium measure, composed of independent and identically distributed Gaussian random variables. Here, we first focus on the discretely sampled Ornstein-Uhlenbeck process, which interpolates between these two limits. We show that in the limit of large times, its extreme values converge to those of the equilibrium distribution, instead of those of the continuously sampled process. This finding folds for any positive sampling interval, with an abrupt transition at zero. We then analyze the Langevin equation for any force that gives rise to a stable equilibrium distribution. For forces which asymptotically grow with the distance from the equilibrium point, the above conclusion continues to hold, and the extreme values for large times correspond to those of independent variables drawn from the equilibrium distribution. However, for forces which asymptotically decay to zero with the distance, the discretely sampled extreme value statistics at large times approach those of the continuously sampled process. △ Less

Submitted 15 August, 2021; originally announced August 2021.

Comments: 21 pages, 10 figures

Anomalous statistics of laser-cooled atoms in dissipative optical lattices

Authors: Gadi Afek, Nir Davidson, David A. Kessler, Eli Barkai

Abstract: Diffusion occurs in numerous physical systems throughout nature, drawing its generality from the universality of the central limit theorem. Around a century ago it was realized that an extension to this type of dynamics can be obtained in the form of ``anomalous" diffusion where distributions are allowed to have heavy, power-law tails. Due to a unique feature of its momentum-dependant dissipative… ▽ More Diffusion occurs in numerous physical systems throughout nature, drawing its generality from the universality of the central limit theorem. Around a century ago it was realized that an extension to this type of dynamics can be obtained in the form of ``anomalous" diffusion where distributions are allowed to have heavy, power-law tails. Due to a unique feature of its momentum-dependant dissipative friction force, laser-cooled atomic ensembles can be used as a test bed for such dynamics. The interplay between laser cooling and anomalous dynamics bears deep, predictive implications for fundamental concepts in both equilibrium and non-equilibrium statistical physics. The high degree of control available in cold-atom experiments allows for tuning of the parameters of the friction force, revealing transitions in the dynamical properties of the system. Rare events, in both the momentum and spatial distributions, are described by non-normalized states using tools adapted from infinite ergodic theory. This leads to new experimental and theoretical results, illuminating the various features of the system. △ Less

Submitted 25 January, 2023; v1 submitted 20 July, 2021; originally announced July 2021.

Big jump principle for heavy-tailed random walks with correlated increments

Authors: Marc Höll, Eli Barkai

Abstract: The big jump principle explains the emergence of extreme events for physical quantities modelled by a sum of independent and identically distributed random variables which are heavy-tailed. Extreme events are large values of the sum and they are solely dominated by the largest summand called the big jump. Recently, the principle was introduced into physical sciences where systems usually exhibit c… ▽ More The big jump principle explains the emergence of extreme events for physical quantities modelled by a sum of independent and identically distributed random variables which are heavy-tailed. Extreme events are large values of the sum and they are solely dominated by the largest summand called the big jump. Recently, the principle was introduced into physical sciences where systems usually exhibit correlations. Here, we study the principle for a random walk with correlated increments. Examples are the autoregressive model of first order and the discretized Ornstein-Uhlenbeck process both with heavy-tailed noise. The correlation leads to the dependence of large values of the sum not only on the big jump but also on the following increments. We describe this behaviour by two big jump principles, namely unconditioned and conditioned on the step number when the big jump occurs. The unconditional big jump principle is described by a correlation dependent shift between the sum and maximum distribution tails. For the conditional big jump principle, the shift depends also on the step number of the big jump. △ Less

Submitted 27 June, 2021; originally announced June 2021.

Driving quantum systems with repeated conditional measurements

Authors: Quancheng Liu, Klaus Ziegler, David A. Kessler, Eli Barkai

Abstract: 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… ▽ More 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. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Transitions in the ergodicity of subrecoil-laser-cooled gases

Authors: Eli Barkai, Günter Radons, Takuma Akimoto

Abstract: With subrecoil-laser-cooled atoms one may reach nano-Kelvin 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… ▽ More With subrecoil-laser-cooled atoms one may reach nano-Kelvin 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 non-integrable observable, with respect to the infinite (non-normalised) invariant density. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 5 pages, 1 figure, 4 pages of supplementary material

Journal ref: Phys. Rev. Lett. 127, 140605 (2021)

Uncertainty relation between detection probability and energy fluctuations

Authors: Felix Thiel, Itay Mualem, David Kessler, Eli Barkai

Abstract: 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 is fully exploring 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 detec… ▽ More 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 is fully exploring 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. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 12 pages, special issue Axiomatic Approaches to Quantum Mechanics, 4 figures

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

Authors: Celia Anteneodo, Lucianno Defaveri, Eli Barkai, David A. Kessler

Abstract: 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… ▽ More 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. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Comments: 22 pages, 8 figures

Journal ref: Entropy 2021, 23(2), 131

Randomly repeated measurements on quantum systems: Correlations and topological invariants of the quantum evolution

Authors: K. Ziegler, E. Barkai, D. Kessler

Abstract: Randomly repeated measurements during the evolution of a closed quantum system create a sequence of probabilities for the first detection of a certain quantum state. The related discrete monitored evolution for the return of the quantum system to its initial state is investigated. We found that the mean number of measurements until the first detection is an integer, namely the dimensionality of th… ▽ More Randomly repeated measurements during the evolution of a closed quantum system create a sequence of probabilities for the first detection of a certain quantum state. The related discrete monitored evolution for the return of the quantum system to its initial state is investigated. We found that the mean number of measurements until the first detection is an integer, namely the dimensionality of the accessible Hilbert space. Moreover, the mean first detected return time is equal to the average time step between successive measurements times the mean number of measurements. Thus, the mean first detected return time scales linearly with the dimensionality of the accessible Hilbert space. The main goal of this work is to explain the quantization of the mean return time in terms of a quantized Berry phase. △ Less

Submitted 10 June, 2021; v1 submitted 30 December, 2020; originally announced December 2020.

Comments: 21 pages, 3 figures

Journal ref: J. Phys. A: Math. Theor. 54 395302 (2021)

Cusp of non-Gaussian density of particles for a diffusing diffusivity model

Authors: M. Hidalgo-Soria, E. Barkai, S. Burov

Abstract: We study a two state ``jumping diffusivity'' model for a Brownian process alternating between two different diffusion constants, $D_{+}>D_{-}$, with random waiting times in both states whose distribution is rather general. In the limit of long measurement times Gaussian behavior with an effective diffusion coefficient is recovered. We show that for equilibrium initial conditions and when the limit… ▽ More We study a two state ``jumping diffusivity'' model for a Brownian process alternating between two different diffusion constants, $D_{+}>D_{-}$, with random waiting times in both states whose distribution is rather general. In the limit of long measurement times Gaussian behavior with an effective diffusion coefficient is recovered. We show that for equilibrium initial conditions and when the limit of the diffusion coefficient $D_-\to0$ is taken, the short time behavior leads to a cusp, namely a non - analytical behavior, in the distribution of the displacements $P(x,t)$ for $x\longrightarrow 0$. Visually this cusp, or tent-like shape, resembles similar behavior found in many experiments of diffusing particles in disordered environments, such as glassy systems and intracellular media. This general result depends only on the existence of finite mean values of the waiting times at the different states of the model. Gaussian statistics in the long time limit is achieved due to ergodicity and convergence of the distribution of the temporal occupation fraction in state $D_{+}$ to a $δ$-function. The short time behavior of the same quantity converges to a uniform distribution, which leads to the non - analyticity in $P(x,t)$. We demonstrate how super - statistical framework is a zeroth order short time expansion of $P(x,t)$, in the number of transitions, that does not yield the cusp like shape. The latter, considered as the key feature of experiments in the field, is found with the first correction in perturbation theory. △ Less

Submitted 8 February, 2021; v1 submitted 21 December, 2020; originally announced December 2020.

The first detection time of a quantum state under random probing

Authors: David A. Kessler, Eli Barkai, Klaus Ziegler

Abstract: 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… ▽ More 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. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. A 103, 022222 (2021)

Large deviations of the ballistic Lévy walk model

Authors: Wanli Wang, Marc Höll, Eli Barkai

Abstract: 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<x<v_0 t$. In particular we study this density close to its maximum in the vicinity of the `light' cone. The spreading density follows the Lamperti-arcsine law de… ▽ More 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<x<v_0 t$. In particular we study this density close to its maximum in the vicinity of the `light' cone. The spreading density follows the Lamperti-arcsine law describing typical fluctuations far from the `light' cone. However this law blows up in the vicinity of the `light' cone horizon which is nonphysical, in the sense that any finite time observation will never diverge. We claim that one can find two laws for the spatial density, the first one is the mentioned Lamperti-arcsine law describing the central part of the distribution and the second is an infinite density illustrating the dynamics for large $x$. We identify the relationship between a large position and the longest traveling time describing the single big jump principle. From the renewal theory we find that the distribution of rare events of the position is related to the derivative of the average of the number of renewals at a short `time' using a rate formalism. △ Less

Submitted 15 October, 2020; v1 submitted 5 July, 2020; originally announced July 2020.

Comments: 12 pages

Journal ref: Phys. Rev. E 102, 052115 (2020)

Accurately approximating extreme value statistics

Authors: Lior Zarfaty, Eli Barkai, David A. Kessler

Abstract: We consider the extreme value statistics of $N$ independent and identically distributed random variables, which is a classic problem in probability theory. When $N\to\infty$, fluctuations around the maximum of the variables are described by the Fisher-Tippett-Gnedenko theorem, which states that the distribution of maxima converges to one out of three limiting forms. Among these is the Gumbel distr… ▽ More We consider the extreme value statistics of $N$ independent and identically distributed random variables, which is a classic problem in probability theory. When $N\to\infty$, fluctuations around the maximum of the variables are described by the Fisher-Tippett-Gnedenko theorem, which states that the distribution of maxima converges to one out of three limiting forms. Among these is the Gumbel distribution, for which the convergence rate with $N$ is of a logarithmic nature. Here, we present a theory that allows one to use the Gumbel limit to accurately approximate the exact extreme value distribution. We do so by representing the scale and width parameters as power series, and by a transformation of the underlying distribution. We consider functional corrections to the Gumbel limit as well, showing they are obtainable via Taylor expansion. Our method also improves the description of large deviations from the mean extreme value. Additionally, it helps to characterize the extreme value statistics when the underlying distribution is unknown, for example when fitting experimental data. △ Less

Submitted 6 July, 2021; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: 17 pages, 10 figures

Journal ref: J. Phys. A: Math. Theor. 54 315205 (2021)

Extreme value theory for constrained physical systems

Authors: Marc Höll, Wanli Wang, Eli Barkai

Abstract: We investigate extreme value theory for physical systems with a global conservation law which describe 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 non-analytical point in the middle of the support. We expose exact relationships between constrained extreme value theor… ▽ More We investigate extreme value theory for physical systems with a global conservation law which describe 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 non-analytical 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 generality, 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 normalisable scaling function for the typical statistics and a non-normalised state describing the rare events. △ Less

Submitted 13 October, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. E 102, 042141 (2020)

Large deviations for continuous time random walks

Authors: Wanli Wang, Eli Barkai, Stanislav Burov

Abstract: 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… ▽ More 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. △ Less

Submitted 21 March, 2022; v1 submitted 27 May, 2020; originally announced May 2020.

Comments: 19 pages and 14 figures

Fractional advection-diffusion-asymmetry equation

Authors: Wanli Wang, Eli Barkai

Abstract: 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 s… ▽ More 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. △ Less

Submitted 8 November, 2020; v1 submitted 12 May, 2020; originally announced May 2020.

Comments: 17 pages

Journal ref: Phys. Rev. Lett. 125, 240606 (2020)

Regularized Boltzmann-Gibbs statistics for a Brownian particle in a non-confining field

Authors: Lucianno Defaveri, Celia Anteneodo, David A. Kessler, Eli Barkai

Abstract: We consider an overdamped Brownian particle subject to an asymptotically flat potential with a trap of depth $U_0$ around the origin. When the temperature is small compared to the trap depth ($ξ=k_B T/U_0 \ll 1$), there exists a range of timescales over which physical observables remain practically constant. This range can be very long, of the order of the Arrhenius factor ${\rm e}^{1/ξ}$. For the… ▽ More We consider an overdamped Brownian particle subject to an asymptotically flat potential with a trap of depth $U_0$ around the origin. When the temperature is small compared to the trap depth ($ξ=k_B T/U_0 \ll 1$), there exists a range of timescales over which physical observables remain practically constant. This range can be very long, of the order of the Arrhenius factor ${\rm e}^{1/ξ}$. For these quasi-equilibrium states, the usual Boltzmann-Gibbs recipe does not work, since the partition function is divergent due to the flatness of the potential at long distances. However, we show that the standard Boltzmann-Gibbs (BG) statistical framework and thermodynamic relations can still be applied through proper regularization. This can be a valuable tool for the analysis of metastability in the non-confining potential fields that characterize a vast number of systems. △ Less

Submitted 15 June, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 9 pages, 6 figures

Journal ref: Phys. Rev. Research 2, 043088 (2020)

Quantum walks: the first detected transition time

Authors: Q. Liu, R. Yin, K. Ziegler, E. Barkai

Abstract: We consider the quantum first detection problem for a particle evolving on a graph under repeated projective measurements with fixed rate $1/τ$. A general formula for the mean first detected transition time is obtained for a quantum walk in a finite-dimensional Hilbert space where the initial state $|ψ_{\rm in}\rangle$ of the walker is orthogonal to the detected state $|ψ_{\rm d}\rangle$. We focus… ▽ More We consider the quantum first detection problem for a particle evolving on a graph under repeated projective measurements with fixed rate $1/τ$. A general formula for the mean first detected transition time is obtained for a quantum walk in a finite-dimensional Hilbert space where the initial state $|ψ_{\rm in}\rangle$ of the walker is orthogonal to the detected state $|ψ_{\rm d}\rangle$. We focus on diverging mean transition times, where the total detection probability exhibits a discontinuous drop of its value, by mapping the problem onto a theory of fields of classical charges located on the unit disk. Close to the critical parameter of the model, which exhibits a blow-up of the mean transition time, we get simple expressions for the mean transition time. Using previous results on the fluctuations of the return time, corresponding to $|ψ_{\rm in}\rangle = |ψ_{\rm d}\rangle$, we find close to these critical parameters that the mean transition time is proportional to the fluctuations of the return time, an expression reminiscent of the Einstein relation. △ Less

Submitted 1 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. Research 2, 033113 (2020)

Infinite ergodic theory meets Boltzmann statistics

Authors: Erez Aghion, David A. Kessler, Eli Barkai

Abstract: 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… ▽ More 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. △ Less

Submitted 29 May, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 21 pages (including Refs), 18 figures

Journal ref: https://doi.org/10.1016/j.chaos.2020.109890

The Hitchhiker model for Laplace diffusion processes in the cell environment

Authors: M. Hidalgo-Soria, E. Barkai

Abstract: Aggregation and fragmentation of single molecules in the cell environment lead to a spectrum of diffusivities and to statistical laws of movement very different from typical Brownian motion. Current models of intracellular transport do not explain at a microscopical level the emergence of theses deviations. Employing a many body approach, which we call the Hitchhiker model, we elucidate how the wi… ▽ More Aggregation and fragmentation of single molecules in the cell environment lead to a spectrum of diffusivities and to statistical laws of movement very different from typical Brownian motion. Current models of intracellular transport do not explain at a microscopical level the emergence of theses deviations. Employing a many body approach, which we call the Hitchhiker model, we elucidate how the widely observed exponential tails in the particle spreading, i.e. the Laplace distribution and the modulations of the diffusivities, are controlled by size fluctuations of single molecules. By means of numerical simulations Laplace distributions are obtained whether we track one molecule or many molecules in parallel. However, we show that the diffusivity varies significantly depending on which tracking protocol is applied. Using a renewal process in the space of sizes, we quantify to what extent the average diffusivity in the single molecule technique is decreased compared with the ensemble average. △ Less

Submitted 6 January, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 8 pages, 7 figures

Journal ref: Phys. Rev. E 102, 012109 (2020)

Quantum total detection probability from repeated measurements II. Exploiting symmetry

Authors: Felix Thiel, Itay Mualem, David A. Kessler, Eli Barkai

Abstract: A quantum walker on a graph, prepared in the state $| ψ_{\rm in} \rangle$, e.g. initially localized at node $r_{\rm in}$, is repeatedly probed, with fixed frequency $1/τ$, to test its presence at some target node $r_{\rm d}$ until the first successful detection. This is a quantum version of the first-passage problem. We investigate the total detection probability $P_{\rm det}$, i.e. the probabilit… ▽ More A quantum walker on a graph, prepared in the state $| ψ_{\rm in} \rangle$, e.g. initially localized at node $r_{\rm in}$, is repeatedly probed, with fixed frequency $1/τ$, to test its presence at some target node $r_{\rm d}$ until the first successful detection. This is a quantum version of the first-passage problem. We investigate the total detection probability $P_{\rm det}$, i.e. the probability to eventually detect the particle after an arbitrary number of detection attempts. It is demonstrated that this total detection probability is less than unity in symmetric systems, where it is possible to find initial states which are shielded from the detector by destructive interference, so-called dark states. The identification of physically equivalent initial states yields an upper bound for $P_{\rm det}$ in terms of the reciprocal of the number $ν$ of physically equivalent states. The relevant subgroup of the system's symmetry operations is found to be the stabilizer of the detection state. Using this, we prove that all bright, i.e. surely detectable, states are symmetric with respect to the stabilizer. This implies that $P_{\rm det}$ can be obtained from a diagonalization of the "symmetrized" Hamiltonian, instead of having to find all eigenstates of the Hamiltonian. △ Less

Submitted 4 September, 2019; originally announced September 2019.

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

Authors: Alessandro Vezzani, Eli Barkai, Raffaella Burioni

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2020; v1 submitted 28 August, 2019; originally announced August 2019.

Comments: 14 pages, 6 figures

Journal ref: Sci. Rep. 10, 2732 (2020)

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

Authors: Takuma Akimoto, Eli Barkai, Günter Radons

Abstract: 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 ou… ▽ More 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. △ Less

Submitted 10 July, 2020; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 16 pages, 7 figures

Journal ref: Phys. Rev. E 101, 052112 (2020)

Infinite horizon billiards: Transport at the border between Gauss and Lévy universality classes

Authors: Lior Zarfaty, Alexander Peletskyi, Eli Barkai, Sergey Denisov

Abstract: We consider transport in two billiard models, the infinite horizon Lorentz gas and the stadium channel, presenting analytical results for the spreading packet of particles. We first obtain the cumulative distribution function of traveling times between collisions, which exhibits non-analytical behavior. Using a renewal assumption and the Lévy walk model, we obtain the particles' probability densit… ▽ More We consider transport in two billiard models, the infinite horizon Lorentz gas and the stadium channel, presenting analytical results for the spreading packet of particles. We first obtain the cumulative distribution function of traveling times between collisions, which exhibits non-analytical behavior. Using a renewal assumption and the Lévy walk model, we obtain the particles' probability density. For the Lorentz gas, it shows a distinguished difference when compared with the known Gaussian propagator, as the latter is valid only for extremely long times. In particular, we show plumes of particles spreading along the infinite corridors, creating power-law tails of the density. We demonstrate the slow convergence rate via summation of independent identically distributed random variables on the border between Lévy and Gauss laws. The renewal assumption works well for the Lorentz gas with intermediately sized scattering centers, but fails for the stadium channel due to strong temporal correlations. Our analytical results are supported with numerical samplings. △ Less

Submitted 8 October, 2019; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 14 pages, 13 figures, 5 appendixes (25 pages in total)

Journal ref: Phys. Rev. E 100, 042140 (2019)

Packets of Diffusing Particles Exhibit Universal Exponential Tails

Authors: Eli Barkai, Stanislav Burov

Abstract: 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… ▽ More 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. △ Less

Submitted 15 February, 2020; v1 submitted 23 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. Lett. 124, 060603 (2020)

Dark states of quantum search cause imperfect detection

Authors: Felix Thiel, Itay Mualem, Dror Meidan, Eli Barkai, David A. Kessler

Abstract: We consider a quantum walk where a detector repeatedly probes the system with fixed rate $1/τ$ until the walker is detected. This is a quantum version of the first-passage problem. We focus on the total probability, $P_{\mathrm{det}}$, that the particle is eventually detected in some target state, for example on a node $r_{\mathrm{d}}$ on a graph, after an arbitrary number of detection attempts. A… ▽ More We consider a quantum walk where a detector repeatedly probes the system with fixed rate $1/τ$ until the walker is detected. This is a quantum version of the first-passage problem. We focus on the total probability, $P_{\mathrm{det}}$, that the particle is eventually detected in some target state, for example on a node $r_{\mathrm{d}}$ on a graph, after an arbitrary number of detection attempts. Analyzing the dark and bright states for finite graphs, and more generally for systems with a discrete spectrum, we provide an explicit formula for $P_{\mathrm{det}}$ in terms of the energy eigenstates which is generically $τ$ independent. We find that disorder in the underlying Hamiltonian renders perfect detection: $P_{\mathrm{det}}=1$, and then expose the role of symmetry with respect to sub-optimal detection. Specifically, we give a simple upper bound for $P_{\mathrm{det}}$ that is controlled by the number of equivalent (with respect to the detection) states in the system. We also extend our results to infinite systems, for example the detection probability of a quantum walk on a line, which is $τ$-dependent and less than half, well below Polya's optimal detection for a classical random walk. △ Less

Submitted 7 October, 2020; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: 18 pages, 8 figures, revised and enlarged version, partly merged with arxiv:1909.02114v1, corrected a mistake in Eq. (17)

Journal ref: Phys. Rev. Research 2, 043107 (2020)