Distribution of the time between maximum and minimum of random walks

Universal to nonuniversal transition of the statistics of rare events during the spread of random walks

Through numerous experiments that analyzed rare event statistics in heterogeneous media, it was discovered that in many cases the probability density function for particle position, P(X,t), exhibits a slower decay rate than the Gaussian function. Typically, the decay behavior is exponential, referred to as Laplace tails. However, many systems exhibit an even slower decay rate, such as power-law, log-normal, or stretched exponential. In this study, we utilize the continuous-time random walk method to investigate the rare events in particle hopping dynamics and find that the properties of the hop size distribution induce a critical transition between the Laplace universality of rare events and a more specific, slower decay of P(X,t). Specifically, when the hop size distribution decays slower than exponential, such as e^{-|x|^{β}} (β>1), the Laplace universality no longer applies, and the decay is specific, influenced by a few large events, rather than by the accumulation of many smaller events that give rise to Laplace tails.

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.

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

Non-Newtonian Rheology in a Capillary Tube with Varying Radius

Fluid blobs in an immiscible Newtonian fluid flowing in a capillary tube with varying radius show highly nonlinear behavior. We consider here a generalization of previously obtained results to blobs of non-Newtonian fluids. We compute here the yield pressure drop and the mean flow rate in two cases: (i) When a single blob is injected, (ii) When many blobs are randomly injected into the tube. We find that the capillary effects emerge from the non-uniformity of the tube radius and contribute to the threshold pressure for flow to occur. Furthermore, in the presence of many blobs the threshold value depends on the number of blobs and their relative distances which are randomly distributed. For a capillary fiber bundle of identical parallel tubes, we calculate the probability distribution of the threshold pressure and the mean flow rate. We consider two geometries: tubes of sinusoidal shape, for which we derive explicit expressions, and triangular-shaped tubes, for which we find that essential singularities are developed. We perform numerical simulations confirming our analytical results.

Sampling rare trajectories using stochastic bridges

The numerical quantification of the statistics of rare events in stochastic processes is a challenging computational problem. We present a sampling method that constructs an ensemble of stochastic trajectories that are constrained to have fixed start and end points (so-called stochastic bridges). We then show that by carefully choosing a set of such bridges and assigning an appropriate statistical weight to each bridge, one can focus more processing power on the rare events of a target stochastic process while faithfully preserving the statistics of these rare trajectories. Further, we also compare the stochastic bridges we produce to the Wentzel-Kramers-Brillouin (WKB) optimal paths of the target process, derived in the limit of low noise. We see that the generated paths, encoding the full statistics of the process, collapse onto the WKB optimal path as the level of noise is reduced. We propose that the method can also be used to judge the accuracy of the WKB approximation at finite levels of noise.

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.

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.

Mean perimeter and area of the convex hull of a planar Brownian motion in the presence of resetting

We compute exactly the mean perimeter and the mean area of the convex hull of a two-dimensional isotropic Brownian motion of duration t and diffusion constant D, in the presence of resetting to the origin at a constant rate r. We show that for any t, the mean perimeter is given by 〈L(t)〉=2πsqrt[D/r]f_{1}(rt) and the mean area is given by 〈A(t)〉=2πD/rf_{2}(rt) where the scaling functions f_{1}(z) and f_{2}(z) are computed explicitly. For large t≫1/r, the mean perimeter grows extremely slowly as 〈L(t)〉∝ln(rt) with time. Likewise, the mean area also grows slowly as 〈A(t)〉∝ln^{2}(rt) for t≫1/r. Our exact results indicate that the convex hull, in the presence of resetting, approaches a circular shape at late times due to the isotropy of the Brownian motion. Numerical simulations are in perfect agreement with our analytical predictions.

Universal properties of a run-and-tumble particle in arbitrary dimension

We consider an active run-and-tumble particle (RTP) in d dimensions, starting from the origin and evolving over a time interval [0,t]. We examine three different models for the dynamics of the RTP: the standard RTP model with instantaneous tumblings, a variant with instantaneous runs and a general model in which both the tumblings and the runs are noninstantaneous. For each of these models, we use the Sparre Andersen theorem for discrete-time random walks to compute exactly the probability that the x component does not change sign up to time t, showing that it does not depend on d. As a consequence of this result, we compute exactly other x-component properties, namely, the distribution of the time of the maximum and the record statistics, showing that they are universal, i.e., they do not depend on d. Moreover, we show that these universal results hold also if the speed v of the particle after each tumbling is random, drawn from a generic probability distribution. Our findings are confirmed by numerical simulations. Some of these results have been announced in a recent Letter [Phys. Rev. Lett. 124, 090603 (2020)10.1103/PhysRevLett.124.090603].