Time-dependent density of diffusion with stochastic resetting is invariant to return speed

Optimal conditions for first passage of jump processes with resetting

We investigate the first passage time beyond a barrier located at b≥0 of a random walk with independent and identically distributed jumps, starting from x0=0. The walk is subject to stochastic resetting, meaning that after each step the evolution is restarted with fixed probability r. We consider a resetting protocol that is an intermediate situation between a random walk (r=0) and an uncorrelated sequence of jumps all starting from the origin (r=1) and derive a general condition for determining when restarting the process with 0 Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Mattia Radice Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
Giampaolo Cristadoro Dipartimento di Matematica e Applicazioni, Università degli Studi Milano-Bicocca, 20126 Milan, Italy
Samudrajit Thapa Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany

Collapse

Granular gases under resetting

We investigate the granular temperatures in force-free granular gases under exponential resetting. When a resetting event occurs, the granular temperature attains its initial value, whereas it decreases because of the inelastic collisions between the resetting events. We develop a theory and perform computer simulations for granular gas cooling in the presence of Poissonian resetting events. We also investigate the probability density function to quantify the distribution of granular temperatures. Our theory may help us to understand the behavior of nonperiodically driven granular systems.

Search with stochastic home returns can expedite classical first passage under resetting

Classical first passage under resetting is a paradigm in the search process. Despite its multitude of applications across interdisciplinary sciences, experimental realizations of such resetting processes posit practical challenges in calibrating these zero time irreversible transitions. Here, we consider a strategy in which resetting is performed using finite-time return protocols in lieu of instantaneous returns. These controls could also be accompanied with random fluctuations or errors allowing target detection even during the return phase. To better understand the phenomena, we develop a unified renewal approach that can encapsulate arbitrary search processes centered around home in a fairly general topography containing targets, various resetting times, and return mechanisms in arbitrary dimensions. While such finite-time protocols would apparently seem to prolong the overall search time in comparison to the instantaneous resetting process, we show on the contrary that a significant speed-up can be gained by leveraging the stochasticity in home returns. The formalism is then explored to reveal a universal criterion distilling the benefits of this strategy. We demonstrate how this general principle can be utilized to improve overall performance of a one-dimensional diffusive search process reinforced with experimentally feasible parameters. We believe that such strategies designed with inherent randomness can be made optimal with precise controllability in complex search processes.

Role of long jumps in Lévy noise-induced multimodality

Lévy noise is a paradigmatic noise used to describe out-of-equilibrium systems. Typically, properties of Lévy noise driven systems are very different from their Gaussian white noise driven counterparts. In particular, under action of Lévy noise, stationary states in single-well, super-harmonic, potentials are no longer unimodal. Typically, they are bimodal; however, for fine-tuned potentials, the number of modes can be further increased. The multimodality arises as a consequence of the competition between long displacements induced by the non-equilibrium stochastic driving and action of the deterministic force. Here, we explore robustness of bimodality in the quartic potential under action of the Lévy noise. We explore various scenarios of bounding long jumps and assess their ability to weaken and destroy multimodality. In general, we demonstrate that despite its robustness it is possible to destroy the bimodality, however it requires drastic reduction in the length of noise-induced jumps.

Lévy flights and Lévy walks under stochastic resetting

Stochastic resetting is a protocol of starting anew, which can be used to facilitate the escape kinetics. We demonstrate that restarting can accelerate the escape kinetics from a finite interval restricted by two absorbing boundaries also in the presence of heavy-tailed, Lévy-type, α-stable noise. However, the width of the domain where resetting is beneficial depends on the value of the stability index α determining the power-law decay of the jump length distribution. For heavier (smaller α) distributions, the domain becomes narrower in comparison to lighter tails. Additionally, we explore connections between Lévy flights (LFs) and Lévy walks (LWs) in the presence of stochastic resetting. First of all, we show that for Lévy walks, the stochastic resetting can also be beneficial in the domain where the coefficient of variation is smaller than 1. Moreover, we demonstrate that in the domain where LWs are characterized by a finite mean jump duration (length), with the increasing width of the interval, the LWs start to share similarities with LFs under stochastic resetting.

Intermittent random walks under stochastic resetting

We analyze a one-dimensional intermittent random walk on an unbounded domain in the presence of stochastic resetting. In this process, the walker alternates between local intensive search, diffusion, and rapid ballistic relocations in which it does not react to the target. We demonstrate that Poissonian resetting leads to the existence of a non-equilibrium steady state. We calculate the distribution of the first arrival time to a target along with its mean and show the existence of an optimal reset rate. In particular, we prove that the initial condition of the walker, i.e., either starting diffusely or relocating, can significantly affect the long-time properties of the search process. Moreover, we demonstrate the presence of distinct parameter regimes for the global optimization of the mean first arrival time when ballistic and diffusive movements are in direct competition.

Dynamically emergent correlations between particles in a switching harmonic trap

We study a one dimensional gas of N noninteracting diffusing particles in a harmonic trap, whose stiffness switches between two values μ_{1} and μ_{2} with constant rates r_{1} and r_{2}, respectively. Despite the absence of direct interaction between the particles, we show that strong correlations between them emerge in the stationary state at long times, induced purely by the dynamics itself. We compute exactly the joint distribution of the positions of the particles in the stationary state, which allows us to compute several physical observables analytically. In particular, we show that the extreme value statistics (EVS), i.e., the distribution of the position of the rightmost particle, has a nontrivial shape in the large N limit. The scaling function characterizing this EVS has a finite support with a tunable shape (by varying the parameters). Remarkably, this scaling function turns out to be universal. First, it also describes the distribution of the position of the kth rightmost particle in a 1d trap. Moreover, the distribution of the position of the particle farthest from the center of the harmonic trap in d dimensions is also described by the same scaling function for all d≥1. Numerical simulations are in excellent agreement with our analytical predictions.

Truncated stochastically switching processes

There are a large variety of hybrid stochastic systems that couple a continuous process with some form of stochastic switching mechanism. In many cases the system switches between different discrete internal states according to a finite-state Markov chain, and the continuous dynamics depends on the current internal state. The resulting hybrid stochastic differential equation (hSDE) could describe the evolution of a neuron's membrane potential, the concentration of proteins synthesized by a gene network, or the position of an active particle. Another major class of switching system is a search process with stochastic resetting, where the position of a diffusing or active particle is reset to a fixed position at a random sequence of times. In this case the system switches between a search phase and a reset phase, where the latter may be instantaneous. In this paper, we investigate how the behavior of a stochastically switching system is modified when the maximum number of switching (or reset) events in a given time interval is fixed. This is motivated by the idea that each time the system switches there is an additive energy cost. We first show that in the case of an hSDE, restricting the number of switching events is equivalent to truncating a Volterra series expansion of the particle propagator. Such a truncation significantly modifies the moments of the resulting renormalized propagator. We then investigate how restricting the number of reset events affects the diffusive search for an absorbing target. In particular, truncating a Volterra series expansion of the survival probability, we calculate the splitting probabilities and conditional MFPTs for the particle to be absorbed by the target or exceed a given number of resets, respectively.

Efficiency of a microscopic heat engine subjected to stochastic resetting

We explore the thermodynamics of stochastic heat engines in the presence of stochastic resetting. The setup comprises an engine whose working substance is a Brownian particle undergoing overdamped Langevin dynamics in a harmonic potential with a time-dependent stiffness, with the dynamics interrupted at random times with a resetting to a fixed location. The effect of resetting to the potential minimum is shown to enhance the efficiency of the engine, while the output work is shown to have a nonmonotonic dependence on the rate of resetting. The resetting events are found to drive the system out of the linear response regime, even for small differences in the bath temperatures. Shifting the reset point from the potential minimum is observed to reduce the engine efficiency. The experimental setup for the realization of such an engine is briefly discussed.

Thermodynamic trade-off relation for first passage time in resetting processes

Resetting is a strategy for boosting the speed of a target-searching process. Since its introduction over a decade ago, most studies have been carried out under the assumption that resetting takes place instantaneously. However, due to its irreversible nature, resetting processes incur a thermodynamic cost, which becomes infinite in the case of instantaneous resetting. Here, we take into consideration both the cost and the first passage time (FPT) required for a resetting process, in which the reset or return to the initial location is implemented using a trapping potential over a finite but random time period. An iterative generating function and a counting functional method à la Feynman and Kac are employed to calculate the FPT and the average work for this process. From these results, we obtain an explicit form of the time-cost trade-off relation, which provides the lower bound of the mean FPT for a given work input when the trapping potential is linear. This trade-off relation clearly shows that instantaneous resetting is achievable only when an infinite amount of work is provided. More surprisingly, the trade-off relation derived from the linear potential seems to be valid for a wide range of trapping potentials. In addition, we have also shown that the fixed-time or sharp resetting can further enhance the trade-off relation compared to that of the stochastic resetting.

Universal performance bounds of restart

As has long been known to computer scientists, the performance of probabilistic algorithms characterized by relatively large runtime fluctuations can be improved by applying a restart, i.e., episodic interruption of a randomized computational procedure followed by initialization of its new statistically independent realization. A similar effect of restart-induced process acceleration could potentially be possible in the context of enzymatic reactions, where dissociation of the enzyme-substrate intermediate corresponds to restarting the catalytic step of the reaction. To date, a significant number of analytical results have been obtained in physics and computer science regarding the effect of restart on the completion time statistics in various model problems, however, the fundamental limits of restart efficiency remain unknown. Here we derive a range of universal statistical inequalities that offer constraints on the effect that restart could impose on the completion time of a generic stochastic process. The corresponding bounds are expressed via simple statistical metrics of the original process such as harmonic mean h, median value m, and mode M, and, thus, are remarkably practical. We test our analytical predictions with multiple numerical examples, discuss implications arising from them and important avenues of future work.

Diffusion with partial resetting

Inspired by many examples in nature, stochastic resetting of random processes has been studied extensively in the past decade. In particular, various models of stochastic particle motion were considered where, upon resetting, the particle is returned to its initial position. Here we generalize the model of diffusion with resetting to account for situations where a particle is returned only a fraction of its distance to the origin, e.g., half way. We show that this model always attains a steady-state distribution which can be written as an infinite sum of independent, but not identical, Laplace random variables. As a result, we find that the steady-state transitions from the known Laplace form which is obtained in the limit of full resetting to a Gaussian form, which is obtained close to the limit of no resetting. A similar transition is shown to be displayed by drift diffusion whose steady state can also be expressed as an infinite sum of independent random variables. Finally, we extend our analysis to capture the temporal evolution of drift diffusion with partial resetting, providing a bottom-up probabilistic construction that yields a closed-form solution for the time-dependent distribution of this process in Fourier-Laplace space. Possible extensions and applications of diffusion with partial resetting are discussed.

First-passage time of run-and-tumble particles with noninstantaneous resetting

We study the statistics of the first-passage time of a single run-and-tumble particle (RTP) in one spatial dimension, with or without resetting, to a fixed target located at L>0. First, we compute the first-passage time distribution of a free RTP, without resetting or in a confining potential, but averaged over the initial position drawn from an arbitrary distribution p(x). Recent experiments used a noninstantaneous resetting protocol that motivated us to study in particular the case where p(x) corresponds to the stationary non-Boltzmann distribution of an RTP in the presence of a harmonic trap. This distribution p(x) is characterized by a parameter ν>0, which depends on the microscopic parameters of the RTP dynamics. We show that the first-passage time distribution of the free RTP, drawn from this initial distribution, develops interesting singular behaviors, depending on the value of ν. We then switch on resetting, mimicked by relaxation of the RTP in the presence of a harmonic trap. Resetting leads to a finite mean first-passage time and we study this as a function of the resetting rate for different values of the parameters ν and b=L/c, where c is the position of the right edge of the initial distribution p(x). In the diffusive limit of the RTP dynamics, we find a rich phase diagram in the (b,ν) plane, with an interesting reentrance phase transition. Away from the diffusive limit, qualitatively similar rich behaviors emerge for the full RTP dynamics.

One-dimensional telegraphic process with noninstantaneous stochastic resetting

In this paper, we consider the one-dimensional dynamical evolution of a particle traveling at constant speed and performing, at a given rate, random reversals of the velocity direction. The particle is subject to stochastic resetting, meaning that at random times it is forced to return to the starting point. Here we consider a return mechanism governed by a deterministic law of motion, so that the time cost required to return is correlated to the position occupied at the time of the reset. We show that in such conditions the process reaches a stationary state which, for some kinds of deterministic return dynamics, is independent of the return phase. Furthermore, we investigate the first-passage properties of the system and provide explicit formulas for the mean first-hitting time. Our findings are supported by numerical simulations.

Time averaging and emerging nonergodicity upon resetting of fractional Brownian motion and heterogeneous diffusion processes

How different are the results of constant-rate resetting of anomalous-diffusion processes in terms of their ensemble-averaged versus time-averaged mean-squared displacements (MSDs versus TAMSDs) and how does stochastic resetting impact nonergodicity? We examine, both analytically and by simulations, the implications of resetting on the MSD- and TAMSD-based spreading dynamics of particles executing fractional Brownian motion (FBM) with a long-time memory, heterogeneous diffusion processes (HDPs) with a power-law space-dependent diffusivity D(x)=D_{0}|x|^{γ} and their "combined" process of HDP-FBM. We find, inter alia, that the resetting dynamics of originally ergodic FBM for superdiffusive Hurst exponents develops disparities in scaling and magnitudes of the MSDs and mean TAMSDs indicating weak ergodicity breaking. For subdiffusive HDPs we also quantify the nonequivalence of the MSD and TAMSD and observe a new trimodal form of the probability density function. For reset FBM, HDPs and HDP-FBM we compute analytically and verify by simulations the short-time MSD and TAMSD asymptotes and long-time plateaus reminiscent of those for processes under confinement. We show that certain characteristics of these reset processes are functionally similar despite a different stochastic nature of their nonreset variants. Importantly, we discover nonmonotonicity of the ergodicity-breaking parameter EB as a function of the resetting rate r. For all reset processes studied we unveil a pronounced resetting-induced nonergodicity with a maximum of EB at intermediate r and EB∼(1/r)-decay at large r. Alongside the emerging MSD-versus-TAMSD disparity, this r-dependence of EB can be an experimentally testable prediction. We conclude by discussing some implications to experimental systems featuring resetting dynamics.

Stochastic resetting by a random amplitude

Stochastic resetting, a diffusive process whose amplitude is reset to the origin at random times, is a vividly studied strategy to optimize encounter dynamics, e.g., in chemical reactions. Here we generalize the resetting step by introducing a random resetting amplitude such that the diffusing particle may be only partially reset towards the trajectory origin or even overshoot the origin in a resetting step. We introduce different scenarios for the random-amplitude stochastic resetting process and discuss the resulting dynamics. Direct applications are geophysical layering (stratigraphy) and population dynamics or financial markets, as well as generic search processes.

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.

Resetting transition is governed by an interplay between thermal and potential energy

A dynamical process that takes a random time to complete, e.g., a chemical reaction, may either be accelerated or hindered due to resetting. Tuning system parameters, such as temperature, viscosity, or concentration, can invert the effect of resetting on the mean completion time of the process, which leads to a resetting transition. Although the resetting transition has been recently studied for diffusion in a handful of model potentials, it is yet unknown whether the results follow any universality in terms of well-defined physical parameters. To bridge this gap, we propose a general framework that reveals that the resetting transition is governed by an interplay between the thermal and potential energy. This result is illustrated for different classes of potentials that are used to model a wide variety of stochastic processes with numerous applications.

Stochastic resetting of active Brownian particles with Lorentz force

Equilibrium properties of a system of passive diffusing particles in an external magnetic field are unaffected by Lorentz force. In contrast, active Brownian particles exhibit steady-state phenomena that depend on both the strength and the polarity of the applied magnetic field. The intriguing effects of the Lorentz force, however, can only be observed when out-of-equilibrium density gradients are maintained in the system. To this end, we use the method of stochastic resetting on active Brownian particles in two dimensions by resetting them to the line x = 0 at a constant rate and periodicity in the y direction. Under stochastic resetting, an active system settles into a nontrivial stationary state which is characterized by an inhomogeneous density distribution, polarization and bulk fluxes perpendicular to the density gradients. We show that whereas for a uniform magnetic field the properties of the stationary state of the active system can be obtained from its passive counterpart, novel features emerge in the case of an inhomogeneous magnetic field which have no counterpart in passive systems. In particular, there exists an activity-dependent threshold rate such that for smaller resetting rates, the density distribution of active particles becomes non-monotonic. We also study the mean first-passage time to the x axis and find a surprising result: it takes an active particle more time to reach the target from any given point for the case when the magnetic field increases away from the axis. The theoretical predictions are validated using Brownian dynamics simulations.

Continuous time random walks under Markovian resetting

We investigate the effects of Markovian resetting events on continuous time random walks where the waiting times and the jump lengths are random variables distributed according to power-law probability density functions. We prove the existence of a nonequilibrium stationary state and finite mean first arrival time. However, the existence of an optimum reset rate is conditioned to a specific relationship between the exponents of both power-law tails. We also investigate the search efficiency by finding the optimal random walk which minimizes the mean first arrival time in terms of the reset rate, the distance of the initial position to the target, and the characteristic transport exponents.

Asymmetric stochastic resetting: Modeling catastrophic events

In the classical stochastic resetting problem, a particle, moving according to some stochastic dynamics, undergoes random interruptions that bring it to a selected domain, and then the process recommences. Hitherto, the resetting mechanism has been introduced as a symmetric reset about the preferred location. However, in nature, there are several instances where a system can only reset from certain directions, e.g., catastrophic events. Motivated by this, we consider a continuous stochastic process on the positive real line. The process is interrupted at random times occurring at a constant rate, and then the former relocates to a value only if the current one exceeds a threshold; otherwise, it follows the trajectory defined by the underlying process without resetting. An approach to obtain the exact nonequilibrium steady state of such systems and the mean first passage time to reach the origin is presented. Furthermore, we obtain the explicit solutions for two different model systems. Some of the classical results found in symmetric resetting, such as the existence of an optimal resetting, are strongly modified. Finally, numerical simulations have been performed to verify the analytical findings, showing an excellent agreement.

Experimental Realization of Diffusion with Stochastic Resetting

Stochastic resetting is prevalent in natural and man-made systems, giving rise to a long series of nonequilibrium phenomena. Diffusion with stochastic resetting serves as a paradigmatic model to study these phenomena, but the lack of a well-controlled platform by which this process can be studied experimentally has been a major impediment to research in the field. Here, we report the experimental realization of colloidal particle diffusion and resetting via holographic optical tweezers. We provide the first experimental corroboration of central theoretical results and go on to measure the energetic cost of resetting in steady-state and first-passage scenarios. In both cases, we show that this cost cannot be made arbitrarily small because of fundamental constraints on realistic resetting protocols. The methods developed herein open the door to future experimental study of resetting phenomena beyond diffusion.

Brownian motion under noninstantaneous resetting in higher dimensions

We consider Brownian motion under resetting in higher dimensions for the case when the return of the particle to the origin occurs at a constant speed. We investigate the behavior of the probability density function (PDF) and of the mean-squared displacement (MSD) in this process. We study two different resetting protocols: exponentially distributed time intervals between the resetting events (Poissonian resetting) and resetting at fixed time intervals (deterministic resetting). We moreover discuss a general problem of the invariance of the PDF with respect to the return speed, as observed in the one-dimensional system for Poissonian resetting, and show that this one-dimensional situation is the only one in which such an invariance can be found. However, the invariance of the MSD can still be observed in higher dimensions.

Continuous-time random walks under power-law resetting

We study continuous-time random walks (CTRW) with power-law distribution of waiting times under resetting which brings the walker back to the origin, with a power-law distribution of times between the resetting events. Two situations are considered. Under complete resetting, the CTRW after the resetting event starts anew, with a new waiting time, independent of the prehistory. Under incomplete resetting, the resetting of the coordinate does not influence the waiting time until the next jump. We focus on the behavior of the mean-squared displacement (MSD) of the walker from its initial position, on the conditions under which the probability density functions of the walker's displacement show universal behavior, and on this universal behavior itself. We show, that the behavior of the MSD is the same as in the scaled Brownian motion (SBM), being the mean-field model of the CTRW. The intermediate asymptotics of the probability density functions (PDF) for CTRW under complete resetting (provided they exist) are also the same as in the corresponding case for SBM. For incomplete resetting, however, the behavior of the PDF for CTRW and SBM is vastly different.

Resetting processes with noninstantaneous return

We consider a random two-phase process which we call a reset-return one. The particle starts its motion at the origin. The first, displacement, phase corresponds to a stochastic motion of a particle and is finished at a resetting event. The second, return, phase corresponds to the particle's motion toward the origin from the position it attained at the end of the displacement phase. This motion toward the origin takes place according to a given equation of motion. The whole process is a renewal one. We provide general expressions for the stationary probability density function of the particle's position and for the mean hitting time in one dimension. We perform explicit analysis for the Brownian motion during the displacement phase and three different types of the return motion: return at a constant speed, return at a constant acceleration with zero initial speed, and return under the action of a harmonic force. We assume that the waiting times for resetting events follow an exponential distribution or that resetting takes place after a fixed waiting period. For the first two types of return motion and the exponential resetting, the stationary probability density function of the particle's position is invariant under return speed (acceleration), while no such invariance is found for deterministic resetting, and for exponential resetting with return under the action of the harmonic force. We discuss necessary conditions for such invariance of the stationary PDF of the positions with respect to the properties of the return process, and we demonstrate some additional examples when this invariance does or does not take place.

Work Fluctuations and Jarzynski Equality in Stochastic Resetting

We consider the paradigm of an overdamped Brownian particle in a potential well, which is modulated through an external protocol, in the presence of stochastic resetting. Thus, in addition to the short range diffusive motion, the particle also experiences intermittent long jumps that reset the particle back at a preferred location. Due to the modulation of the trap, work is done on the system and we investigate the statistical properties of the work fluctuations. We find that the distribution function of the work typically, in asymptotic times, converges to a universal Gaussian form for any protocol as long as that is also renewed after each resetting event. When observed for a finite time, we show that the system does not generically obey the Jarzynski equality that connects the finite time work fluctuations to the difference in free energy. Nonetheless, we identify herein a restricted set of protocols which embraces the relation. In stark contrast, the Jarzynski equality is always fulfilled when the protocols continue to evolve without being reset. We present a set of exactly solvable models, demonstrate the validation of our theory and carry out numerical simulations to illustrate these findings. Finally, we have pointed out possible realistic implementations for resetting in experiments using the so-called engineered swift equilibration.