Optimal mean first-passage time of a Brownian searcher with resetting in one and two dimensions: Experiments, theory and numerical tests F. Faisant, B. Besga, A. Petrosyan, S. Ciliberto∗ Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS, Laboratoire de Physique, UMR 5672, F-69342 Lyon, France Satya N. Majumdar LPTMS, CNRS, Univ. Paris-Sud, Université Paris-Saclay, UMR 8626, 91405 Orsay, France (Dated: June 18, 2021) We study experimentally, numerically and theoretically the optimal mean time needed by a Brownian particle, freely diffusing either in one or two dimensions, to reach, within a tolerance radius Rtol, a target at a distance L from an initial position in the presence of resetting. The reset position is Gaussian distributed with width s. We derived and tested two resetting protocols, one with a periodic and one with random (Poissonian) resetting times. We computed and measured the full first-passage probability distribution that displays spectacular spikes immediately after each resetting time for close targets. We study the optimal mean first-passage time as a function of the resetting period/rate for different target distances (values of the ratios b = L=s) and target size (a = Rtol=L). We find an interesting phase transition at a critical value of b, both in one and two dimensions. The details of the calculations as well as experimental setup and limitations are discussed. Contents E. Two dimensions 18 1. Poissonian Resetting 19 I. Introduction 2 a. The limit s 0: resetting to the origin 20 b. Finite s > 0! 21 II. Experimental set-up. 2 2. Periodic Resetting 22 a. The limit s 0: resetting to the origin 23 III. The optimal resetting protocols in one dimension 3 b. Finite s > 0! 25 A. Periodic resetting protocol.3 B. Poissonian resetting protocol.5 F. Anisotropy of the initial position distribution 27 IV. Optimal resetting in 2 dimensions 6 G. Numerical integration of the Langevin equation 28 A. Periodic resetting in 2 dimensions6 1. Numerical simulation of Langevin dynamics7 2. Experimental results7 B. Poissonian resetting in 2 dimensions7 V. Discussion and conclusions 8 Acknowledgments 9 References 9 A. Preliminaries 10 B. Protocol-1: Periodic Resetting to a random initial position 10 1. Reset to a Gaussian initial distribution 12 a. Mean first-passage time 13 b. Full first-passage probability density F(t) 13 arXiv:2106.09113v1 [cond-mat.stat-mech] 16 Jun 2021 C. Protocol-2: Random resetting to a random initial position 15 D. General Setting: the survival probability of a Brownian walker in arbitrary dimensions 16 ∗E-mail me at: [email protected] 2 I. INTRODUCTION The time needed by a freely diffusing Brownian particle to reach a fixed target in space is a widely studied problem not 2.0 only for its fundamental aspects but also for its numerous ap- 1.5 plications from chemical reactions, astrophysics, all the way to animal foraging in ecology. This time is a random variable 1.0 t (s) and its probability distribution is called the first-passage time 0.5 distribution [1,2]. The first moment of this distribution, when it exists, is called the mean first-passage time (MFPT) which 0.0 is often used as a measure of the efficiency of a search pro- cess. The smaller the MFPT, the more efficient is the search 0.5 0.4 process. The MFPT can be optimized using various stochastic − 0.0 0.0 µm) search algorithms (such as simulated annealing) which speed 0.4 y( x(µ 0.5 up the search process [3–6]. More recently it has been demon- m) − strated that the MFPT for a single particle (searcher) to a fixed 1.0 target can be minimized using the reset protocols (for a re- 0.5 cent review see [7]). Searching a target via resetting is an m) µ 0.0 example of the so called intermittent search strategy [8] con- x( 0.5 sisting of a mixture of short-range localised moves (where the − searcher actually performs the search) with intermittent long- 0.0 0.5 1.0 1.5 2.0 range moves where the searcher moves to a new location and time (s) starts a local search in the new place. Consider, for example, a fixed target at some point in space FIG. 1: Top panel : Brownian trajectory in the 2D plane (x;y) over and a particle (searcher) starts its dynamics from a fixed ini- time t with resetting. Bottom panel : Corresponding 1D Brown- tial position in space. The dynamics of the particle is rep- ian trajectory with periodic resetting after T = 0:5s (red colored resented by a stochastic process, e.g. it may just be simple area). The equilibrium standard deviation of the trapped particle diffusion [9, 10]. Under resetting protocol, the motion of sx = 43nm is shown (dotted lines). The target (red square 1µm away from the center of the trap) is reached when x crosses the red line for the searcher is interrupted either randomly at a constant rate the first time. r [9, 10] or periodically with a period T [11, 12], and then the particle is sent back to its initial position (usually instanta- neously) and the search restarts again from the initial position In the next section we describe the experimental set-up and after each resetting event. For a given fixed distance between the difficulties that one encounters in performing such exper- the target and the initial position of the searcher, if one plots iments. In section III we recall the main experimental results the MFPT as a function of the resetting rate r (or T for periodi- in 1D, whereas the theoretical predictions for 1D are summa- cal resetting), the MFPT typically displays a unique minimum rized in the AppendicesA-C. Although the 1D results have at some optimal value r∗. This led to the paradigm that reset- been recently published in a shortened form without details in ting typically makes the search process more efficient if the Ref. [32], we think that it is useful to recall them with details resetting rate is chosen to have its optimal value r∗ [9, 10]. in this article, in order to have a full comparisons with the 2D While this ‘optimal resetting’ paradigm has been tested and results described in SectionIV. The details of the 2D theo- verified in a large number of recent theoretical and numerical retical computations, which are rather nontrivial, are provided studies [11, 13–30] (see also the review [7]), it has been been in the AppendicesD andE. We finally conclude in SectionV verified only in a few experiments limited to the one dimen- and some further details are relegated to the rest of the Appen- sional (1D) case [31, 32]. In Ref. [32] we have demonstrated dices. the emergence of new physical effects when the resetting posi- tion is not fixed to be exactly the initial position, but has a dis- tribution of finite width s around an average initial position. A finite value of s led to interesting metastability and phase II. EXPERIMENTAL SET-UP. transition in the MFPT as a function of resetting rate/period. Very recently [33], we have shown that a finite s in the ini- The optimal search protocol is implemented in an optical tial condition also affects, rather profoundly, the first-passage tweezer experiment [34]. The Brownian searcher is a silica probability distribution of the searcher, even in the absence of micro-sphere of radius R = 1µm ( 5%) immersed in pure resetting. water. The fluid chamber is made± of two parts : The col- The purpose of this article is to extend these results to the loids are inserted in the thicker part of the chamber (reservoir) two dimensional (2D) case, showing both theoretically and and the measurements are performed in a second part (thick- experimentally that the new physical effects induced by a ness around 50µm) with very few particles which allows us to nonzero s discussed above are not restricted just to 1D, but take long duration measurements (tenths of hours). The opti- are also observed in 2D. Our article is organized as follows. cal trap is implemented with a near infrared laser with wave- 3 x2=2s 2 2 length l = 1064nm focused into the chamber through an oil in a harmonic trap: P(x0) = e− 0 =p2ps with width immersion objective (Leica 63 and 1:4 NA). The particle s = k =k, where is the temperature. We also set a × BT T is trapped in the transverse plane (x;y) by a harmonic poten- a fixed target at location L. p tial with stiffness kx;y. The stiffness is controlled by changing At time t = 0, the trap is switched off for an interval T and the optical power (typically hundreds of mW) in the chamber the particle undergoes free diffusion (overdamped) with diffu- directly by laser current modulation (up to 10 kHz) or faster sion constant D = kBT =G, where G is the friction coefficient. using an electro-optical modulator (EOM). The position of the At the end of the period T, the particle’s position is reset (see particle in the (x;y) plane can be tracked (see Fig. (1)) either Fig. (1)). Performing the resetting of the position poses the by a white light imaging on a camera or from the deviation of real experimental challenge. a second laser on a quadrant photodiode (QPD). The camera In standard models of resetting one usually assumes instan- allows us to extract the position of the bead at maximum speed taneous resetting [7], which is however impossible to achieve around 1000 frames per second over a broad range (> 100µm) experimentally.