Generalized Arcsine Laws for Fractional Brownian Motion

Occupation time of a system of Brownian particles on the line with steplike initial condition

We consider a system of noninteracting Brownian particles on the line with steplike initial condition and study the statistics of the occupation time on the positive half-line. We demonstrate that even at large times, the behavior of the occupation time exhibits long-lasting memory effects of the initialization. Specifically, we calculate the mean and the variance of the occupation time, demonstrating that the memory effects in the variance are determined by a generalized compressibility (or Fano factor), associated with the initial condition. In the particular case of the uncorrelated uniform initial condition we conduct a detailed study of two probability distributions of the occupation time: annealed (averaged over all possible initial configurations) and quenched (for a typical configuration). We show that at large times both the annealed and the quenched distributions admit large deviation form and we compute analytically the associated rate functions. We verify our analytical predictions via numerical simulations using importance sampling Monte Carlo strategy.

Arcsine Laws of Light

We demonstrate that the time-integrated light intensity transmitted by a coherently driven resonator obeys Lévy's arcsine laws-a cornerstone of extreme value statistics. We show that convergence to the arcsine distribution is algebraic, universal, and independent of nonequilibrium behavior due to nonconservative forces or nonadiabatic driving. We furthermore verify, numerically, that the arcsine laws hold in the presence of frequency noise and in Kerr-nonlinear resonators supporting non-Gaussian states. The arcsine laws imply a weak ergodicity breaking which can be leveraged to enhance the precision of resonant optical sensors with zero energy cost, as shown in our companion manuscript [V. G. Ramesh et al., companion paper, Phys. Rev. Res. (2024).PPRHAI2643-1564]. Finally, we discuss perspectives for probing the possible breakdown of the arcsine laws in systems with memory.

Extremal statistics for a resetting Brownian motion before its first-passage time

We study the extreme value statistics of a one-dimensional resetting Brownian motion (RBM) till its first passage through the origin starting from the position x_{0} (>0). By deriving the exit probability of RBM in an interval [0,M] from the origin, we obtain the distribution P_{r}(M|x_{0}) of the maximum displacement M and thus gives the expected value 〈M〉 of M as functions of the resetting rate r and x_{0}. We find that 〈M〉 decreases monotonically as r increases, and tends to 2x_{0} as r→∞. In the opposite limit, 〈M〉 diverges logarithmically as r→0. Moreover, we derive the propagator of RBM in the Laplace domain in the presence of both absorbing ends, and then leads to the joint distribution P_{r}(M,t_{m}|x_{0}) of M and the time t_{m} at which this maximum is achieved in the Laplace domain by using a path decomposition technique, from which the expected value 〈t_{m}〉 of t_{m} is obtained explicitly. Interestingly, 〈t_{m}〉 shows a nonmonotonic dependence on r, and attains its minimum at an optimal r^{*}≈2.71691D/x_{0}^{2}, where D is the diffusion coefficient. Finally, we perform extensive simulations to validate our theoretical results.

Occupation time statistics of the fractional Brownian motion in a finite domain

We study statistics of occupation times for a fractional Brownian motion (fBm), which is a typical model of a non-Markov process. Due to the non-Markovian nature, recurrence times to the origin depend on the history. Numerical simulations indicate that dependence on the sum of successive recurrence times becomes weak. As a result, the distribution of the occupation time in a finite domain follows the Mittag-Leffler distribution when the Hurst exponent of the fBm is close to 1/2. We show this distributional behavior of a time-averaged observable by renewal theory. This result is an extension of the distributional limit theorem known as the Darling-Kac theorem in general Markov processes to non-Markov processes.

Path integrals for fractional Brownian motion and fractional Gaussian noise

Wiener's path integral plays a central role in the study of Brownian motion. Here we derive exact path-integral representations for the more general fractional Brownian motion (FBM) and for its time derivative process, fractional Gaussian noise (FGN). These paradigmatic non-Markovian stochastic processes, introduced by Kolmogorov, Mandelbrot, and van Ness, found numerous applications across the disciplines, ranging from anomalous diffusion in cellular environments to mathematical finance. Their exact path-integral representations were previously unknown. Our formalism exploits the Gaussianity of the FBM and FGN, relies on the theory of singular integral equations, and overcomes some technical difficulties by representing the action functional for the FBM in terms of the FGN for the subdiffusive FBM and in terms of the derivative of the FGN for the super-diffusive FBM. We also extend the formalism to include external forcing. The exact and explicit path-integral representations make inroads in the study of the FBM and FGN.

Experimental verification of arcsine laws in mesoscopic nonequilibrium systems

A large number of processes in the mesoscopic world occur out of equilibrium, where the time evolution of a system becomes immensely important since it is driven principally by dissipative effects. Nonequilibrium steady states (NESS) represent a crucial category in such systems, where relaxation timescales are comparable to the operational timescales. In this study, we employ a model NESS stochastic system, which is comprised of a colloidal microparticle optically trapped in a viscous fluid, externally driven by a temporally correlated noise, and show that time-integrated observables such as the entropic current, the work done on the system or the work dissipated by it, follow the three Lévy arcsine laws [A. C. Barato et al., Phys. Rev. Lett. 121, 090601 (2018)0031-900710.1103/PhysRevLett.121.090601], in the large time limit. We discover that cumulative distributions converge faster to arcsine distributions when it is near equilibrium and the rate of entropy production is small, because in that case the entropic current has weaker temporal autocorrelation. We study this phenomenon by changing the strength of the added noise as well as by perturbing our system with a flow field produced by a microbubble at close proximity to the trapped particle. We confirm our experimental findings with theoretical simulations of the systems. Our work provides an interesting insight into the NESS statistics of the meso-regime, where stochastic fluctuations play a pivotal role.

Time to reach the maximum for a stationary stochastic process

We consider a one-dimensional stationary time series of fixed duration T. We investigate the time t_{m} at which the process reaches the global maximum within the time interval [0,T]. By using a path-decomposition technique, we compute the probability density function P(t_{m}|T) of t_{m} for several processes, that are either at equilibrium (such as the Ornstein-Uhlenbeck process) or out of equilibrium (such as Brownian motion with stochastic resetting). We show that for equilibrium processes the distribution of P(t_{m}|T) is always symmetric around the midpoint t_{m}=T/2, as a consequence of the time-reversal symmetry. This property can be used to detect nonequilibrium fluctuations in stationary time series. Moreover, for a diffusive particle in a confining potential, we show that the scaled distribution P(t_{m}|T) becomes universal, i.e., independent of the details of the potential, at late times. This distribution P(t_{m}|T) becomes uniform in the "bulk" 1≪t_{m}≪T and has a nontrivial universal shape in the "edge regimes" t_{m}→0 and t_{m}→T. Some of these results have been announced in a recent letter [Europhys. Lett. 135, 30003 (2021)0295-507510.1209/0295-5075/ac19ee].

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Extreme value statistics and arcsine laws for heterogeneous diffusion processes

Heterogeneous diffusion with a spatially changing diffusion coefficient arises in many experimental systems such as protein dynamics in the cell cytoplasm, mobility of cajal bodies, and confined hard-sphere fluids. Here, we showcase a simple model of heterogeneous diffusion where the diffusion coefficient D(x) varies in a power-law way, i.e., D(x)∼|x|^{-α} with the exponent α>-1. This model is known to exhibit anomalous scaling of the mean-squared displacement (MSD) of the form ∼t^{2/2+α} and weak ergodicity breaking in the sense that ensemble averaged and time averaged MSDs do not converge. In this paper, we look at the extreme value statistics of this model and derive, for all α, the exact probability distributions of the maximum spatial displacement M(t) and arg-maximum t_{m}(t) (i.e., the time at which this maximum is reached) till duration t. In the second part of our paper, we analyze the statistical properties of the residence time t_{r}(t) and the last-passage time t_{ℓ}(t) and compute their distributions exactly for all values of α. Our study unravels that the heterogeneous version (α≠0) displays many rich and contrasting features compared to that of the standard Brownian motion (BM). For example, while for BM (α=0), the distributions of t_{m}(t),t_{r}(t), and t_{ℓ}(t) are all identical (á la "arcsine laws" due to Lévy), they turn out to be significantly different for nonzero α. Another interesting property of t_{r}(t) is the existence of a critical α (which we denote by α_{c}=-0.3182) such that the distribution exhibits a local maximum at t_{r}=t/2 for α<α_{c} whereas it has minima at t_{r}=t/2 for α≥α_{c}. The underlying reasoning for this difference hints at the very contrasting natures of the process for α≥α_{c} and α<α_{c} which we thoroughly examine in our paper. All our analytical results are backed by extensive numerical simulations.

Functionals of fractional Brownian motion and the three arcsine laws

Fractional Brownian motion is a non-Markovian Gaussian process indexed by the Hurst exponent H∈(0,1), generalizing standard Brownian motion to account for anomalous diffusion. Functionals of this process are important for practical applications as a standard reference point for nonequilibrium dynamics. We describe a perturbation expansion allowing us to evaluate many nontrivial observables analytically: We generalize the celebrated three arcsine laws of standard Brownian motion. The functionals are: (i) the fraction of time the process remains positive, (ii) the time when the process last visits the origin, and (iii) the time when it achieves its maximum (or minimum). We derive expressions for the probability of these three functionals as an expansion in ɛ=H-1/2, up to second order. We find that the three probabilities are different, except for H=1/2, where they coincide. Our results are confirmed to high precision by numerical simulations.

Arcsine law and multistable Brownian dynamics in a tilted periodic potential

Multistability is one of the most important phenomena in dynamical systems, e.g., bistability enables the implementation of logic gates and therefore computation. Recently multistability has attracted a greatly renewed interest related to memristors and graphene structures, to name only a few. We investigate tristability in velocity dynamics of a Brownian particle subjected to a tilted periodic potential. It is demonstrated that the origin of this effect is attributed to the arcsine law for the velocity dynamics at the zero temperature limit. We analyze the impact of thermal fluctuations and construct the phase diagram for the stability of the velocity dynamics. It suggests an efficient strategy to control the multistability by changing solely the force acting on the particle or temperature of the system. Our findings for the paradigmatic model of nonequilibrium statistical physics apply to, inter alia, Brownian motors, Josephson junctions, cold atoms dwelling in optical lattices, and colloidal systems.

Extremal statistics for stochastic resetting systems

While averages and typical fluctuations often play a major role in understanding the behavior of a nonequilibrium system, this nonetheless is not always true. Rare events and large fluctuations are also pivotal when a thorough analysis of the system is being done. In this context, the statistics of extreme fluctuations in contrast to the average plays an important role, as has been discussed in fields ranging from statistical and mathematical physics to climate, finance, and ecology. Herein, we study extreme value statistics (EVS) of stochastic resetting systems, which have recently gained significant interest due to its ubiquitous and enriching applications in physics, chemistry, queuing theory, search processes, and computer science. We present a detailed analysis for the finite and large time statistics of extremals (maximum and arg-maximum, i.e., the time when the maximum is reached) of the spatial displacement in such system. In particular, we derive an exact renewal formula that relates the joint distribution of maximum and arg-maximum of the reset process to the statistical measures of the underlying process. Benchmarking our results for the maximum of a reset trajectory that pertain to the Gumbel class for large sample size, we show that the arg-maximum density attains a uniform distribution independent of the underlying process at a large observation time. This emerges as a manifestation of the renewal property of the resetting mechanism. The results are augmented with a wide spectrum of Markov and non-Markov stochastic processes under resetting, namely, simple diffusion, diffusion with drift, Ornstein-Uhlenbeck process, and random acceleration process in one dimension. Rigorous results are presented for the first two setups, while the latter two are supported with heuristic and numerical analysis.

Aging arcsine law in Brownian motion and its generalization

Classical arcsine law states that the fraction of occupation time on the positive or the negative side in Brownian motion does not converge to a constant but converges in distribution to the arcsine distribution. Here we consider how a preparation of the system affects the arcsine law, i.e., aging of the arcsine law. We derive an aging distributional theorem for occupation time statistics in Brownian motion, where the ratio of time when measurements start to the measurement time plays an important role in determining the shape of the distribution. Furthermore, we show that this result can be generalized as an aging distributional limit theorem in renewal processes.

Extreme events for fractional Brownian motion with drift: Theory and numerical validation

We study the first-passage time, the distribution of the maximum, and the absorption probability of fractional Brownian motion of Hurst parameter H with both a linear and a nonlinear drift. The latter appears naturally when applying nonlinear variable transformations. Via a perturbative expansion in ɛ=H-1/2, we give the first-order corrections to the classical result for Brownian motion analytically. Using a recently introduced adaptive-bisection algorithm, which is much more efficient than the standard Davies-Harte algorithm, we test our predictions for the first-passage time on grids of effective sizes up to N_{eff}=2^{28}≈2.7×10^{8} points. The agreement between theory and simulations is excellent, and by far exceeds in precision what can be obtained by scaling alone.

Sampling first-passage times of fractional Brownian motion using adaptive bisections

We present an algorithm to efficiently sample first-passage times for fractional Brownian motion. To increase the resolution, an initial coarse lattice is successively refined close to the target, by adding exactly sampled midpoints, where the probability that they reach the target is non-negligible. Compared to a path of N equally spaced points, the algorithm achieves the same numerical accuracy N_{eff}, while sampling only a small fraction of all points. Though this induces a statistical error, the latter is bounded for each bridge, allowing us to bound the total error rate by a number of our choice, say P_{error}^{tot}=10^{-6}. This leads to significant improvements in both memory and speed. For H=0.33 and N_{eff}=2^{32}, we need 5000 times less CPU time and 10000 times less memory than the classical Davies-Harte algorithm. The gain grows for H=0.25 and N_{eff}=2^{42} to 3×10^{5} for CPU and 10^{6} for memory. We estimate our algorithmic complexity as C^{ABSec}(N_{eff})=O[(lnN_{eff})^{3}], to be compared to Davies-Harte, which has complexity C^{DH}(N)=O(NlnN). Decreasing P_{error}^{tot} results in a small increase in complexity, proportional to ln(1/P_{error}^{tot}). Our current implementation is limited to the values of N_{eff} given above, due to a loss of floating-point precision. Our algorithm can be adapted to other extreme events and arbitrary Gaussian processes. It enables one to numerically validate theoretical predictions that were hitherto inaccessible.

Distribution of the time between maximum and minimum of random walks

We consider a one-dimensional Brownian motion of fixed duration T. Using a path-integral technique, we compute exactly the probability distribution of the difference τ=t_{min}-t_{max} between the time t_{min} of the global minimum and the time t_{max} of the global maximum. We extend this result to a Brownian bridge, i.e., a periodic Brownian motion of period T. In both cases, we compute analytically the first few moments of τ, as well as the covariance of t_{max} and t_{min}, showing that these times are anticorrelated. We demonstrate that the distribution of τ for Brownian motion is valid for discrete-time random walks with n steps and with a finite jump variance, in the limit n→∞. In the case of Lévy flights, which have a divergent jump variance, we numerically verify that the distribution of τ differs from the Brownian case. For random walks with continuous and symmetric jumps we numerically verify that the probability of the event "τ=n" is exactly 1/(2n) for any finite n, independently of the jump distribution. Our results can be also applied to describe the distance between the maximal and minimal height of (1+1)-dimensional stationary-state Kardar-Parisi-Zhang interfaces growing over a substrate of finite size L. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 123, 200201 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.200201].

Analytical representation of Gaussian processes in the A-T plane

Closed-form expressions, parametrized by the Hurst exponent H and the length n of a time series, are derived for paths of fractional Brownian motion (fBm) and fractional Gaussian noise (fGn) in the A-T plane, composed of the fraction of turning points T and the Abbe value A. The exact formula for A_{fBm} is expressed via Riemann ζ and Hurwitz ζ functions. A very accurate approximation, yielding a simple exponential form, is obtained. Finite-size effects, introduced by the deviation of fGn's variance from unity, and asymptotic cases are discussed. Expressions for T for fBm, fGn, and differentiated fGn are also presented. The same methodology, valid for any Gaussian process, is applied to autoregressive moving average processes, for which regions of availability of the A-T plane are derived and given in analytic form. Locations in the A-T plane of some real-world examples as well as generated data are discussed for illustration.

Time Between the Maximum and the Minimum of a Stochastic Process

We present an exact solution for the probability density function P(τ=t_{min}-t_{max}|T) of the time difference between the minimum and the maximum of a one-dimensional Brownian motion of duration T. We then generalize our results to a Brownian bridge, i.e., a periodic Brownian motion of period T. We demonstrate that these results can be directly applied to study the position difference between the minimal and the maximal heights of a fluctuating (1+1)-dimensional Kardar-Parisi-Zhang interface on a substrate of size L, in its stationary state. We show that the Brownian motion result is universal and, asymptotically, holds for any discrete-time random walk with a finite jump variance. We also compute this distribution numerically for Lévy flights and find that it differs from the Brownian motion result.

First passage in an interval for fractional Brownian motion

Let X_{t} be a random process starting at x∈[0,1] with absorbing boundary conditions at both ends of the interval. Denote by P_{1}(x) the probability to first exit at the upper boundary. For Brownian motion, P_{1}(x)=x, which is equivalent to P_{1}^{'}(x)=1. For fractional Brownian motion with Hurst exponent H, we establish that P_{1}^{'}(x)=N[x(1-x)]^{1/H-2}e^{εF(x)+O(ε^{2})}, where ε=H-1/2. The function F(x) is analytic and well approximated by its Taylor expansion F(x)≃16(C-1)(x-1/2)^{2}+O(x-1/2)^{4}, where C=0.915... is the Catalan constant. A similar result holds for moments of the exit time starting at x. We then consider the span of X_{t}, i.e., the size of the (compact) domain visited up to time t. For Brownian motion, we derive an analytic expression for the probability that the span reaches 1 for the first time and then generalize it to fractional Brownian motion. Using large-scale numerical simulations with system sizes up to N=2^{24} and a broad range of H, we confirm our analytic results. There are important finite-discretization corrections which we quantify. They are most severe for small H, necessitating going to the large systems mentioned above.

Arcsine Laws in Stochastic Thermodynamics

We show that the fraction of time that a thermodynamic current spends above its average value follows the arcsine law, a prominent result obtained by Lévy for Brownian motion. Stochastic currents with long streaks above or below their average are much more likely than those that spend similar fractions of time above and below their average. Our result is confirmed with experimental data from a Brownian Carnot engine. We also conjecture that two other random times associated with currents obey the arcsine law: the time a current reaches its maximum value and the last time a current crosses its average value. These results apply to, inter alia, molecular motors, quantum dots, and colloidal systems.

Statistics of zero crossings in rough interfaces with fractional elasticity

We study numerically the distribution of zero crossings in one-dimensional elastic interfaces described by an overdamped Langevin dynamics with periodic boundary conditions. We model the elastic forces with a Riesz-Feller fractional Laplacian of order z=1+2ζ, such that the interfaces spontaneously relax, with a dynamical exponent z, to a self-affine geometry with roughness exponent ζ. By continuously increasing from ζ=-1/2 (macroscopically flat interface described by independent Ornstein-Uhlenbeck processes [Phys. Rev. 36, 823 (1930)PHRVAO0031-899X10.1103/PhysRev.36.823]) to ζ=3/2 (super-rough Mullins-Herring interface), three different regimes are identified: (I) -1/2<ζ<0, (II) 0<ζ<1, and (III) 1<ζ<3/2. Starting from a flat initial condition, the mean number of zeros of the discretized interface (I) decays exponentially in time and reaches an extensive value in the system size, or decays as a power-law towards (II) a subextensive or (III) an intensive value. In the steady state, the distribution of intervals between zeros changes from an exponential decay in (I) to a power-law decay P(ℓ)∼ℓ^{-γ} in (II) and (III). While in (II) γ=1-θ with θ=1-ζ the steady-state persistence exponent, in (III) we obtain γ=3-2ζ, different from the exponent γ=1 expected from the prediction θ=0 for infinite super-rough interfaces with ζ>1. The effect on P(ℓ) of short-scale smoothening is also analyzed numerically and analytically. A tight relation between the mean interval, the mean width of the interface, and the density of zeros is also reported. The results drawn from our analysis of rough interfaces subject to particular boundary conditions or constraints, along with discretization effects, are relevant for the practical analysis of zeros in interface imaging experiments or in numerical analysis.