Diffusion under time-dependent resetting

Diffusion under time-dependent resetting

Abstract

We study a Brownian particle diffusing under a time-modulated stochastic resetting mechanism to a fixed position. The rate of resetting is a function of the time since the last reset event. We derive a sufficient condition on for a steady-state probability distribution of the position of the particle to exist. We derive the form of the steady-state distributions under some particular choices of and also consider the late time relaxation behavior of the probability distribution. We consider first passage time properties for the Brownian particle to reach the origin and derive a formula for the mean first passage time. Finally, we study optimal properties of the mean first passage time and show that a threshold function is at least locally optimal for the problem of minimizing the mean first passage time.

pacs:
02.50.-r,05.10.Gg,05.40.-a
\definecolor

dgreenrgb0,0.7,0

I Introduction

Over the last few decades there has been great interest in the study of search problems that appear in various contexts from animal foraging, protein binding on DNA, internet search algorithms to locate one’s misplaced keys. In such different situations, one important issue is to consider optimal search strategies. For example, in protein binding on DNA, a protein molecule binds to a specific target binding site on DNA by an appropriate mixture of 3D diffusion and 1D sliding motion on DNA – called facilitated diffusion (1).

Quite often efficient search strategies involve a mixture of local steps and long-range moves. Such search strategies are called intermittent search process and are observed in foraging animals such as humming birds or bumblebees (2); (3). The E. Coli bacteria alternatively uses ballistic moves, called “runs” with random “tumbles”, to change its direction in order to reach high food concentration regions (4).

Recently, an intermittent stochastic strategy containing such a mixture of local and long-range moves has been introduced (5); (6); (7) in which diffusion of a particle is interrupted by stochastically resetting it to a preferred position whereupon the diffusion process starts afresh. Examples of such stochastic resetting are found in a wide variety of situations. In daily life, while searching for some lost possession, after an unsuccessful search for some duration of time one often goes back to the starting place and recommences the search process again. In the ecological context, the movement of animals during foraging period often involves local diffusive search for food (8); (3). Interestingly, such local diffusive movements are interrupted by long range moves to relocate themselves in other areas and after that local diffusive motion restarts from the relocation region(9). Movement of free-ranging capuchin monkeys in the wild is described quite well by random walks with preferential relocations to places visited in the past (10). There are other examples in nature where similar notions of stochastic resetting can also be found. For example, in the biological context several living organisms use stochastic switching between different phenotypic states to adapt in a fluctuating environment (11); (12); (13); (15); (14). Stochastic restarts are often considered as a useful strategy to optimize computer search algorithms in hard combinatorial problems (16); (17); (18). Also in the context of population growth, random catastrophic events may cause a sudden reduction in the population size and reset it to some lower value. (15).

The mechanism of stochastic resetting fundamentally affects the properties of diffusion process. Consider a particle, starting from , diffusing in one dimension. Let be its position at time . Along with diffusion, the particle is subject to a resetting mechanism in which its motion is interrupted stochastically such that the position of the particle is reset to some fixed position at some rate and after every such event the particle recommences its diffusive motion. In the absence of resetting, the Gaussian distribution of the position of the particle never reaches a steady state. The width of the distribution keeps on growing with time as . On the other hand, in the presence of resetting the distribution of the position becomes a globally current-carrying non-equilibrium steady state with non-Gaussian fluctuations. Resetting also has another important consequence on the properties of diffusion process : the mean first passage time (MFPT) to a particular position of a diffusing particle is infinite whereas it becomes finite in presence of resetting. In fact there is an optimal resetting rate for which this finite MFPT is minimum.

In the last decades, several theoretical works have been dedicated to the study of first-passage properties of intermittent search processes, and in particular to the question of the minimization of the mean first-passage time i.e. the mean time to locate some target (19); (20); (21); (22); (23); (24). These works use different analytical and numerical approaches to minimize the MFPT in a broad range of contexts such as search processes mixing slow diffusive movement and fast ballistic motion (20), persistent random walks (21), finding targets in a domain (22) and Brownian search in spatial heterogeneous media (24). Also in (25), an intermittent search process with Levy flights interrupted by random resettings has been studied where a first order phase transition associated to a discontinuous change in the optimal parameters has been observed.

Recently a number of generalizations of the simple diffusion with resetting have been made. The generalization from the one-dimensional case to higher dimensions has been considered in (26). In the context of a target search process, the effect of partial absorption has been considered in (27). Properties of non-equilibrium steady state for diffusion with resetting have been studied in the presence of a potential (28) or in a bounded domain (29). Other generalizations include resetting to the current maximum of the Brownian particle (30), resetting in continuous-time random walks (31), in Lévy flights (25) etc. The effect of resetting in the dynamics of interacting multi-particle systems such as fluctuating interfaces (32), coagulation-diffusion process (33) and in general chemical reaction schemes (34) have been studied.

In this paper we consider a different generalization of the resetting process: we consider a resetting rate that is time-dependent. This generalization is quite natural in the context of target search. While searching for a lost object, the searcher naturally would not like to reset in the beginning of the search, but as time progresses without success, the searcher would be more and more inclined to go back to the most likely location and start the search process anew. As a result it would be quite natural to consider a situation where the resetting rate grows from zero as time increases. Specifically, we consider a situation where the resetting rate depends on the time since the last reset. Thus the rate as well as the searcher is reset at the resetting event.

In the case of time-independent resetting rates, the renewal property of Brownian motion has been exploited to compute the propagator (the probability density of the particle being at at time given that it began at at time 0) as well as the steady-state distribution of the position of the particle (6); (7); (26); (37). More precisely, the propagator has been expressed as an integral over the time of the last resetting events. For time dependent resetting rates, we find that it is more convenient to express the propagator of the particle as an integral over the time of first resetting events. We refer to this formalism as the first renewal picture in contrast to the formalism used in (6); (7); (26); (37) which we refer to as the last renewal picture. We first recapitulate the constant resetting rate case in Sec. II.1 where we re-derive the known results using the first renewal picture. Then in Sec. II.2 we discuss the generalization of this approach to the time dependent case, where we compute the Laplace transform (w.r.t. time) of the propagator for a general resetting rate function .

Next we study the steady-state properties of the Brownian particle subject to such a resetting mechanism. In the absence of resetting, the Brownian particle keeps on diffusing freely over space with time; but the introduction of a resetting mechanism creates a current towards the reset position. This current may balance the outward current (flowing towards infinity) in the large time limit. As a result the system may reach a steady state eventually but not for all possible choices of the rate function . The question is then for what choices of will a steady state be attained? We address this question in Sec. III, where we study the steady-state behavior of the Brownian particle in presence of a time-dependent resetting rate. We find a sufficient condition on the late time behavior of for a steady state to be attained which is given by equation (21). Such a condition can be easily understood as follows : Let us consider the following two extreme cases, a) and b) . In the former case, the particle is always instantaneously reset to the resetting position and hence the probability distribution of the position is then given by irrespective of the initial position . On the other hand, when , the particle diffuses freely and never attains a steady state distribution (on an infinite system). Hence if decays fast enough to zero as time increases, then the particle may have finite probability of not being reset at all and as a result the distribution of the position of the particle may not reach a steady state. In fact we find that if decays slower than for large , then the particle reaches a steady state. The next question is: how does it relax to the steady state, if indeed one exists?

This question was first studied by Majumdar et al (38) in the context of constant resetting rate case. They found that for a given large time , the distribution in an inner core region around the reset position has already relaxed to the (time-independent) steady state whereas the distribution in the outer region has not yet relaxed. The front dividing the relaxed from the non-relaxed region moves linearly with time through the system. In Sec. IV we observe the same phenomena for a specific class of time dependent resetting rate function . However, depending on the value of , the motion of the front becomes linear, superlinear or sublinear with .

One of the main motivations for studying Brownian motion with resetting is to improve our understanding of search paradigms, as discussed in the first and second paragraphs of the introduction. As a quantitative measure of the performance of a search process, one may consider the mean search time or the mean first passage time (MFPT) to a static target for a given starting position and diffusion constant . In the absence of resetting, the MFPT of a free Brownian particle to a static target at origin is infinite. On the other hand, in the presence of resetting to the initial position with constant rate , the MFPT becomes finite (6); (7). Moreover, it has also been shown (6); (7) that there exists an optimal choice of the constant rate for which MFPT becomes minimum. In this work we ask the following questions : can this scenario be improved further if one considers a time-dependent resetting rate ? What is the optimal rate function for which the MFPT becomes minimum ? We address these questions in Sec. V and Sec. VI respectively. We show that indeed there exists a few choices of for which it is possible to achieve a MFPT lower than the minimum MFPT obtained using constant resetting rate . In Sec. VI we study the optimal resetting rate function, where we conjecture that the optimal time-dependent rate function is given by a threshold function.

We note that several very recent works have also considered the problem of time-dependent resetting. Eule and Metzger (35) use a generalized (non-Markovian) master equation approach to consider the stationary distribution of the particle position and also some late time properties. In that work a Gamma-distribution of waiting times between successive reset events is considered and mean first passage times are numerically determined. Nagar and Gupta (36) consider the particular case of a power law distribution of waiting times between resets. They discuss the steady states and some first passage time properties. In contrast, in our work we consider the general scenario of a time-dependent resetting rate (with corresponding waiting time distribution given by (13)) for various different choices of . Furthermore we consider the late time relaxation of the probability distribution and we study analytically the problem of optimising the mean first passage time.

Ii The model

First let us define diffusion with resetting with time-dependent resetting rate. We consider a single particle (or searcher) in one dimension with initial position at and a resetting position . The position of the particle at time is updated by the following stochastic rule: in a small time interval the position becomes

(1)
(2)

where is the time of the last resetting event. Thus the resetting rate is a function of the time elapsed since the last resetting event. In (2) is a Gaussian white noise with mean and two-time correlator given by

(3)

where is the diffusion constant. Here, angular brackets denote averaging over noise realizations. The initial condition is . To simplify matters, from now on we shall take the initial position to coincide with the resetting position

(4)

unless otherwise specified. The dynamics thus consists of a stochastic mixture of resetting to the initial position with rate (long range move) and ordinary diffusion (local move) with diffusion constant .

We define as the probability of finding the particle at position at time , given that it was at at time in the presence of time-independent stochastic resetting, being the same quantity for the time-dependent stochastic resetting. In both cases the suffix ’’ indicates the presence of resetting.

ii.1 Recap of constant resetting rate

For completeness we review here the formalism and the results for the case of constant resetting rate (6); (26). In this case, one can write down a Master equation for from the dynamical rules for the evolution of the particle given in the preceding section (6)

(5)

with the initial condition . For convenience, we have omitted the arguments of in the above equation. Here, the second and third terms on the right hand side (RHS) account for the resetting events, denoting, respectively, the negative probability flux from each point and a corresponding positive probability flux into . The steady-state solution for the time-independent case satisfies

(6)

Alternatively, a renewal picture which we refer to as the last renewal picture may be used to write down an equation for in terms of the free propagator for a pure diffusive process (without resetting) as in (6); (26)

(7)

Here, we have divided the process into two contributions. The first term in the RHS signifies that there has been no reset at all between time with the probability of no resets given by . This probability is then simply multiplied by the free Brownian propagator

(8)

The second term in the RHS accounts for the fact that there can be multiple resets: the integral sums over contributions in which the last resetting event takes place between time and . Subsequently, the particle propagates freely until the observation time scale . It can be shown that  (7) satisfies the master equation  (5). The steady state can be obtained by taking the infinite time limit

(9)

which satisfies the steady state equation given by  (6). The integral (9) can be evaluated to yield (26)

(10)

The last resetting equation (7) also allows the long time relaxation to the steady state to be analyzed (38) as we shall review in Section IV.

In this work we introduce a first renewal picture where instead of the last resetting, we consider the first resetting between time and having started from . Subsequently, the particle diffuses from until time , in the presence of multiple resets. It is again straightforward to write down an equation for the probability

(11)

where the first term in the RHS corresponds to trajectories in which there are no resets at all. The integral in the second term sums over trajectories in which there has been a first reset between time and and then there can be multiple resets which is taken care of by the reset propagator inside the integral. The equivalence between  (7) and  (11) is easy to show by taking Laplace transforms of both equations. We find that they result in identical expression in the Laplace space

(12)

where is the Laplace transform of and is the Laplace transform of . Taking the inverse Laplace transform with respect to , one can easily obtain  (7).

ii.2 Time-dependent resetting rate

We now turn to the main subject of this paper, that of time-dependent resetting rates , as defined above. In this case one cannot simply write a Master equation for in the presence of time-dependent resetting rate. This is because one must in addition keep track of the time since the last reset. The renewal pictures are more useful than the Master equation formalism to describe the time-dependent rate process. To this end, we define the following time-integrated quantity

(13)

Then the probability of no resets subsequent to an initial reset at is given by and is the probability density for a first reset to occur in the interval .

In the case of time-dependent resetting the last renewal equation (7) is modified to

(14)

where is the probability density for a reset (which turns out to be the last) to occur in . We shall return to this picture in section IV where we consider the late time relaxation behavior.

However, to study the steady-state behavior it is most convenient to use the first renewal framework. In the presence of the time-dependent reset process the first renewal equation becomes

(15)

By taking the Laplace transform, we obtain

(16)

where we have defined

(17)
(18)
(19)

Equation (16) is the main result of this section. In principle, the Laplace transform can be inverted although in practice this is difficult for arbitrary .

Iii Steady state behaviour

In the previous section we have studied the propagator of a Brownian particle in presence of a time-dependent resetting mechanism using the first renewal picture. We have computed the Laplace transform (16) of the propagator in terms of the Laplace transforms of the free propagator weighted by the probability of no reset in duration . In this section we are interested in the large time behavior of . More precisely, we are interested in whether takes a time independent form in limit. If so then, in the Laplace transform language, this simply means that should be expressible as multiplied by a quantity that does not depend on in limit and the steady state distribution would be given by

(20)

provided the limit exists and is not zero. Now it is easy to check from the definitions (18) and (19) that if then also . Thus a sufficient condition on the choice of to achieve a steady state is that which implies that

(21)

This condition implies that must grow sufficiently fast for large so that sufficiently quickly.

Recalling the definition (13), condition (21) will hold if is an increasing function of time since then . Thus for an increasing resetting rate, there always exists a steady state. Also if tends to some constant value greater than zero for large , then grows linearly with time and again (21) is satisfied. On the other hand if is a function decreasing to zero at large time since the last reset, to achieve a steady state, it must decrease sufficiently slowly that diverges with time. Thus (21) implies that if the decreasing rate function is bounded one fold i.e. if decays more slowly than , then there exists a unique steady state.

In the following, we explore few plausible choices of the time-dependent rates.

Case I: linear rate

As mentioned in the introduction, it is natural in the context of search processes to consider a rate which is an increasing function of time. The simplest such function is a linear function. For the linear resetting process , we find the following relations

(22)

One can then write the steady state, in terms of a simple integral depending on one variable by inserting the relations (22) in  (20)

(23)
(24)

This integral can be computed using Mathematica and its full expression is given in terms of hypergeometric functions in Appendix VIII. In Fig. 1a, we compare this result against the same obtained from direct numerical simulation of the dynamics and observe an excellent agreement.

Figure 1: (Color online) Stationary distributions of a single Brownian particle moving on a one-dimensional line in the presence of time-dependent resetting. In this plot we consider three choices of the resetting rate functions : (a) linearly increasing and (b) increasing but bounded and (c) decreasing rate , where is some parameter in appropriate units. In all the figures the circles represent data obtained from direct numerical simulations whereas the solid lines correspond to the respective theoretical expressions (23), (29) and (32) in the main text. The intrinsic parameters of the system are : , and ( in respective units ) for all the three cases.

Case II: increasing, but bounded rate

One unrealistic feature of an increasing rate, such as the linear rate considered in the previous paragraph, is that it may increase indefinitely. Here we consider a resetting processes in which the rate is an increasing function but bounded by an upper limit:

(25)

Thus and as . In this case, the integrated rate (13) is given by implying for . Hence we expect to have a steady state which can be computed from  (20). The function can be easily computed and its limit as is given by

(26)

where is an incomplete Gamma-function. To complete the evaluation of , we now need to compute from  (18) in the limit. After some straightforward manipulations one can show that is given by

(27)
(28)

where the symbol is given by . Hence, the steady state distribution is given by

(29)

In Fig. 1b, we compare the analytical results against the same obtained from direct numerical simulation of the dynamics. The results are found to be in excellent agreement with the numerical results.

Case III: Rates decreasing as a power law

We now explore the case where the rate decreases as a function of time. If the rate decreases too quickly then we expect from (21) that a steady state will not be attained. Therefore it is of natural interest to inquire more details of the condition for which the steady states exist and consequently their characteristic forms if they do.

We first consider the case for which we may take

(30)

(it turns out that the singularity at does not affect the steady state). Then we find

(31)

In this case, criterion (21) is satisfied and we have steady states given by

(32)

where,

(33)

We now consider the case for which . In order to ensure convergence of we take

(34)

The time-integrated rates can be found as

(35)

and therefore criterion (21) is not satisfied. Moreover, it can be checked that for small , whereas so that (20) as and there is no steady state.

Finally, we notice that, is a marginal case. Taking as in (34)

(36)

we obtain

(37)

For , diverges as whereas converges as . Thus according to (20) there is no steady state. However, if , converges as , thus condition (21) is satisfied and there is a steady state given by (20) with

(38)

In Fig. 1c, we plot the steady state distribution corresponding to the power law choice of with and compare with the same obtained from numerical simulations. We again observe an excellent agreement.

Iv Late time relaxation of the probability distribution

For the case of time-independent constant resetting rate , the late time relaxation of the probability distribution has been studied in (38). It was shown that at large time an inner region of the distribution has relaxed to the (time-independent) non-equilibrium steady state, whereas the outer region has not yet relaxed and the time-dependent probability is dominated by trajectories that have not yet undergone any resetting. Thus a front dividing the ‘equilibrated’ region from the ‘unequilibrated’ one moves linearly with time through the system with speed (38).

Here we extend this analysis to the case of time-dependent resetting rates. We begin with the last renewal equation (14)

(39)

and assume that the second term in the right hand side (RHS) dominates at large

(40)

where in the second line, a new integration variable is defined through . We now wish to evaluate the integral by the saddle-point method for large .

Consider a class of rate processes such as the following

(41)

This ensures us of both increasing and decreasing resetting rates along with the existence of a steady state in each case. Now, we have the integrated quantity

(42)

We observe in the scale where i.e. we consider and with held fixed. The integral (40) becomes

Now we expect as and therefore we may evaluate the integral in  (LABEL:ltP-2) by the saddle-point method. The saddle-point equation simply reads

(44)

which defines the value of that dominates the integral

(45)

However this value can only dominate the integral if it is within the integration range which is the case if or equivalently

(46)

On the other hand, if is out of the integration range in which case (46) no longer holds and it turns out that the first term in the right hand side (RHS) of the last renewal equation (39) dominates. This contribution essentially represents those trajectories of the dynamics which have not participated in the resetting process yet.

Figure 2: (Color online) Dynamical relaxation of a single Brownian particle under the time-dependent stochastic resetting to the fixed point at the origin is shown. The resetting rate has the time dependence, according to  (41), with fixed for all the dynamical processes observed at different time scales respectively. The intrinsic parameters for the system are : . The distribution is characterized by two regimes: inside one where the non-equilibrium steady state has already set in and the outside one where the system is still in a transient state. The front given by  (46) separates these two regimes and have been depicted by the dashed vertical lines in the figures. The simulation results, indicated by magenta circles, agree well with the theoretical results obtained from  (47), for both the steady state (dashed red curves) and the transient (the blue triangles) regimes.

The dominant asymptotic large time, and large behavior is therefore given by the following large deviation form:

(47)

where . The interpretation of this result is that the inner region or equivalently, has a time-independent behavior and therefore already has reached the non-equilibrium steady state form. This is the relaxed regime. On the other hand, the outer region or equivalently, is yet to equilibrate and still in a transient state. Thus an equilibration front moves through the system: the motion is superlinear if ; linear if with speed , and sublinear but superdiffusive if . The case recovers the results of (38).

In Fig. 2, we observe this dynamical transition from numerical simulations, where we plot the distribution as a function of for a given value of at three different observation time scales. For comparison we plot the theoretical large deviation form of the distribution given by  (47) along with the simulation results. We find very good agreement between them in both the transient and the steady state regimes.

V Survival Probability under time-dependent resetting

We now look into the first passage probabilities of a diffusing particle under the influence of time-dependent stochastic resetting. First passage properties are generically of importance since they characterise the performance of search processes in various contexts. The most intuitive and important observable is the first passage time (FPT), for the particle to reach the origin say, which is itself a stochastic quantity and the PDF of FPT has its own rich characteristics. In particular, one of the main goals is to optimize the mean first passage time (MFPT) with respect to the system parameters. A review of the topic can be found in (39).

The problem of computing the survival probability in presence of a static target has been well studied for free Brownian motion (see e.g. (40) for a review). The MFPT was well investigated for diffusion under stochastic resetting at a constant rate in (6); (7). It was found that the MFPT attains a minima with respect to the rate in contrast to the case of a single diffusive searcher in the absence of resetting. Thus there in fact exists an optimal resetting strategy which can be useful for all kinds of persistence problems. In this section, we address the question whether this strategy can be made even more efficient by introducing time-dependent resetting.

To this end, we compute the survival probability in presence of resetting at a time-dependent rate. To do so it is useful to introduce that defines the survival probability of a free Brownian particle at until time , starting from . This result is well known (see e.g. (41)) and is given by

(48)

Let be the probability that the Brownian particle starting from survives until time without being absorbed at in the presence of stochastic resetting to . The suffix ’’ indicates the presence of resetting, as before. Without any loss of generality, once again we will assume .

Using the first renewal formalism, we now write the survival probability of the Brownian particle subject to time-dependent resetting

(49)

The first term on the RHS represents trajectories in which there has been neither reset to nor absorption at the origin until time . The integral in the second term represents a sum over first reset times and implies that there has been no reset and no absorption until time then a reset between and . The factor inside the integral gives the probability of no resetting and no absorption up to time and a resetting event between and . This is multiplied with the term which simply implies that no absorption occurs from to .

Equation (49) belongs to a Wiener-Hopf class of integrals which can be solved in the Laplace space

(50)

where we have defined the following quantities

(51)
(52)

The mean first passage time to the origin starting from in the presence of stochastic resetting to is then given by

(53)

which upon using  (50) yields

(54)

Inserting the expressions of and from eqs. (51) and (52) and performing some simple manipulations, one can write in a compact form :

(55)

The integral is given by

(56)

For time-independent resetting rate , one can easily compute to find where , is the inverse distance diffused between two resetting events. Now using this expression for in  (55) we recover the result

(57)

obtained in (6). In (6), it was also shown that there is an optimal choice of the constant value of the rate for given initial/resetting position and diffusion constant , such that the value of MFPT becomes minimum : .

The question now is whether one can make the search process more efficient by introducing time-dependent resetting rates. To answer this question, we compute the MFPT for two choices of rate functions in the following. Then varying the respective parameters in the two rate functions we find that in both the cases one can get smaller than the minimum possible using a constant resetting rate .

v.1

For this case, grows quadratically with time. Inserting this explicit form of in  (56) and performing some simple variable changes one finds that