1
|
Dolgushev M, Mendes TV, Gorin B, Xie K, Levernier N, Bénichou O, Kellay H, Voituriez R, Guérin T. Evidence and quantification of memory effects in competitive first-passage events. SCIENCE ADVANCES 2025; 11:eadp2386. [PMID: 40117368 PMCID: PMC11927618 DOI: 10.1126/sciadv.adp2386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Accepted: 02/14/2025] [Indexed: 03/23/2025]
Abstract
Splitting probabilities quantify the likelihood of a given outcome out of competitive events. This key observable of random walk theory, historically introduced as the gambler's ruin problem, is well understood for memoryless (Markovian) processes. However, in complex systems such as polymer fluids, the motion of a particle should typically be described as a process with memory, for which splitting probabilities are much less characterized analytically. Here, we introduce an analytical approach that provides the splitting probabilities for one-dimensional isotropic non-Markovian Gaussian processes with stationary increments, in the case of two targets. This analysis shows that splitting probabilities are controlled by the out-of-equilibrium trajectories observed after the first passage. This is directly evidenced in a prototypical experimental reaction scheme in viscoelastic fluids. These results are extended to d-dimensional processes in large confining volumes, opening a path toward the study of competitive events in complex media.
Collapse
Affiliation(s)
- Maxim Dolgushev
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne University, 4 Place Jussieu, 75005 Paris, France
| | - Toni Vieira Mendes
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | - Benjamin Gorin
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | - Kaili Xie
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | | | - Olivier Bénichou
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne University, 4 Place Jussieu, 75005 Paris, France
| | - Hamid Kellay
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| | - Raphaël Voituriez
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/Sorbonne University, 4 Place Jussieu, 75005 Paris, France
- Laboratoire Jean Perrin, CNRS/Sorbonne University, 4 Place Jussieu, 75005 Paris, France
| | - Thomas Guérin
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400 Talence, France
| |
Collapse
|
2
|
Ye Y, Chaigneau A, Grebenkov DS. First-passage times to a fractal boundary: Local persistence exponent and its log-periodic oscillations. Phys Rev E 2025; 111:014153. [PMID: 39972924 DOI: 10.1103/physreve.111.014153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2024] [Accepted: 12/03/2024] [Indexed: 02/21/2025]
Abstract
We investigate the statistics of the first-passage time (FPT) to a fractal self-similar boundary of the Koch snowflake. When the starting position is fixed near the absorbing boundary, the FPT distribution exhibits an apparent power-law decay over a broad range of timescales, culminated by an exponential cutoff. By extensive Monte Carlo simulations, we compute the local persistence exponent of the survival probability and reveal its log-periodic oscillations in time due to self-similarity of the boundary. The effect of the starting point on this behavior is analyzed in depth. Theoretical bounds on the survival probability are derived from the analysis of diffusion in a circular sector. Physical rationales for the refined structure of the survival probability are presented.
Collapse
Affiliation(s)
- Yilin Ye
- Ecole Polytechnique, CNRS, Laboratoire de Physique de la Matière Condensée (UMR 7643), - , Institut Polytechnique de Paris, 91120 Palaiseau, France
| | - Adrien Chaigneau
- Ecole Polytechnique, CNRS, Laboratoire de Physique de la Matière Condensée (UMR 7643), - , Institut Polytechnique de Paris, 91120 Palaiseau, France
| | - Denis S Grebenkov
- Ecole Polytechnique, CNRS, Laboratoire de Physique de la Matière Condensée (UMR 7643), - , Institut Polytechnique de Paris, 91120 Palaiseau, France
- Université de Montréal, CNRS, - , CRM - CNRS, 6128 succ Centre-Ville, Montréal QC H3C 3J7, Canada
| |
Collapse
|
3
|
Barbier-Chebbah A, Bénichou O, Voituriez R, Guérin T. Long-term memory induced correction to Arrhenius law. Nat Commun 2024; 15:7408. [PMID: 39198409 PMCID: PMC11358423 DOI: 10.1038/s41467-024-50938-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Accepted: 07/25/2024] [Indexed: 09/01/2024] Open
Abstract
The Kramers escape problem is a paradigmatic model for the kinetics of rare events, which are usually characterized by Arrhenius law. So far, analytical approaches have failed to capture the kinetics of rare events in the important case of non-Markovian processes with long-term memory, as occurs in the context of reactions involving proteins, long polymers, or strongly viscoelastic fluids. Here, based on a minimal model of non-Markovian Gaussian process with long-term memory, we determine quantitatively the mean FPT to a rare configuration and provide its asymptotics in the limit of a large energy barrier E. Our analysis unveils a correction to Arrhenius law, induced by long-term memory, which we determine analytically. This correction, which we show can be quantitatively significant, takes the form of a second effective energy barrierE ' < E and captures the dependence of rare event kinetics on initial conditions, which is a hallmark of long-term memory. Altogether, our results quantify the impact of long-term memory on rare event kinetics, beyond Arrhenius law.
Collapse
Affiliation(s)
- A Barbier-Chebbah
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) and Neuroscience Department CNRS UMR 3751, Institut Pasteur, Université de Paris, CNRS, 75015, Paris, France
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
| | - O Bénichou
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France.
| | - R Voituriez
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
- Laboratoire Jean Perrin, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
| | - T Guérin
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400, Talence, France
| |
Collapse
|
4
|
Benkhadaj Z, Grebenkov DS. Encounter-based approach to diffusion with resetting. Phys Rev E 2022; 106:044121. [PMID: 36397494 DOI: 10.1103/physreve.106.044121] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2022] [Accepted: 09/27/2022] [Indexed: 06/16/2023]
Abstract
An encounter-based approach consists in using the boundary local time as a proxy for the number of encounters between a diffusing particle and a target to implement various surface reaction mechanisms on that target. In this paper, we investigate the effects of stochastic resetting onto diffusion-controlled reactions in bounded confining domains. We first discuss the effect of position resetting onto the propagator and related quantities; in this way, we retrieve a number of earlier results but also provide complementary insights into them. Second, we introduce boundary local time resetting and investigate its impact. Curiously, we find that this type of resetting does not alter the conventional propagator governing the diffusive dynamics in the presence of a partially reactive target with a constant reactivity. In turn, the generalized propagator for other surface reaction mechanisms can be significantly affected. Our general results are illustrated for diffusion on an interval with reactive end points. Further perspectives and some open problems are discussed.
Collapse
Affiliation(s)
| | - Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|
5
|
Everlasting impact of initial perturbations on first-passage times of non-Markovian random walks. Nat Commun 2022; 13:5319. [PMID: 36085151 PMCID: PMC9463153 DOI: 10.1038/s41467-022-32280-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Accepted: 07/25/2022] [Indexed: 11/09/2022] Open
Abstract
Persistence, defined as the probability that a signal has not reached a threshold up to a given observation time, plays a crucial role in the theory of random processes. Often, persistence decays algebraically with time with non trivial exponents. However, general analytical methods to calculate persistence exponents cannot be applied to the ubiquitous case of non-Markovian systems relaxing transiently after an imposed initial perturbation. Here, we introduce a theoretical framework that enables the non-perturbative determination of persistence exponents of Gaussian non-Markovian processes with non stationary dynamics relaxing to a steady state after an initial perturbation. Two situations are analyzed: either the system is subjected to a temperature quench at initial time, or its past trajectory is assumed to have been observed and thus known. Our theory covers the case of spatial dimension higher than one, opening the way to characterize non-trivial reaction kinetics for complex systems with non-equilibrium initial conditions.
Collapse
|
6
|
Chaigneau A, Grebenkov DS. First-passage times to anisotropic partially reactive targets. Phys Rev E 2022; 105:054146. [PMID: 35706315 DOI: 10.1103/physreve.105.054146] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Accepted: 05/05/2022] [Indexed: 06/15/2023]
Abstract
We investigate restricted diffusion in a bounded domain towards a small partially reactive target in three- and higher-dimensional spaces. We propose a simple explicit approximation for the principal eigenvalue of the Laplace operator with mixed Robin-Neumann boundary conditions. This approximation involves the harmonic capacity and the surface area of the target, the volume of the confining domain, the diffusion coefficient, and the reactivity. The accuracy of the approximation is checked by using a finite-elements method. The proposed approximation determines also the mean first-reaction time, the long-time decay of the survival probability, and the overall reaction rate on that target. We identify the relevant lengthscale of the target, which determines its trapping capacity, and we investigate its relation to the target shape. In particular, we study the effect of target anisotropy on the principal eigenvalue by computing the harmonic capacity of prolate and oblate spheroids in various space dimensions. Some implications of these results in chemical physics and biophysics are briefly discussed.
Collapse
Affiliation(s)
- Adrien Chaigneau
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| | - Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| |
Collapse
|
7
|
Grebenkov DS. Depletion of resources by a population of diffusing species. Phys Rev E 2022; 105:054402. [PMID: 35706291 DOI: 10.1103/physreve.105.054402] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2022] [Accepted: 04/19/2022] [Indexed: 06/15/2023]
Abstract
Depletion of natural and artificial resources is a fundamental problem and a potential cause of economic crises, ecological catastrophes, and death of living organisms. Understanding the depletion process is crucial for its further control and optimized replenishment of resources. In this paper, we investigate a stock depletion by a population of species that undergo an ordinary diffusion and consume resources upon each encounter with the stock. We derive the exact form of the probability density of the random depletion time, at which the stock is exhausted. The dependence of this distribution on the number of species, the initial amount of resources, and the geometric setting is analyzed. Future perspectives and related open problems are discussed.
Collapse
Affiliation(s)
- Denis S Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS-Ecole Polytechnique, IP Paris, 91128 Palaiseau, France
| |
Collapse
|
8
|
Grebenkov DS, Kumar A. Reversible target-binding kinetics of multiple impatient particles. J Chem Phys 2022; 156:084107. [DOI: 10.1063/5.0083849] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Certain biochemical reactions can only be triggered after binding a sufficient number of particles to a specific target region such as an enzyme or a protein sensor. We investigate the distribution of the reaction time, i.e., the first instance when all independently diffusing particles are bound to the target. When each particle binds irreversibly, this is equivalent to the first-passage time of the slowest (last) particle. In turn, reversible binding to the target renders the problem much more challenging and drastically changes the distribution of the reaction time. We derive the exact solution of this problem and investigate the short-time and long-time asymptotic behaviors of the reaction time probability density. We also analyze how the mean reaction time depends on the unbinding rate and the number of particles. Our exact and asymptotic solutions are compared to Monte Carlo simulations.
Collapse
Affiliation(s)
- Denis S. Grebenkov
- Laboratoire de Physique de la Matière Condensée (UMR 7643), CNRS–Ecole Polytechnique, IP Paris, 91120 Palaiseau, France
| | - Aanjaneya Kumar
- Department of Physics, Indian Institute of Science Education and Research, Dr. Homi Bhabha Road, Pune 411008, India
| |
Collapse
|
9
|
Ray A, Pal A, Ghosh D, Dana SK, Hens C. Mitigating long transient time in deterministic systems by resetting. CHAOS (WOODBURY, N.Y.) 2021; 31:011103. [PMID: 33754784 DOI: 10.1063/5.0038374] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Accepted: 12/24/2020] [Indexed: 06/12/2023]
Abstract
How long does a trajectory take to reach a stable equilibrium point in the basin of attraction of a dynamical system? This is a question of quite general interest and has stimulated a lot of activities in dynamical and stochastic systems where the metric of this estimation is often known as the transient or first passage time. In nonlinear systems, one often experiences long transients due to their underlying dynamics. We apply resetting or restart, an emerging concept in statistical physics and stochastic process, to mitigate the detrimental effects of prolonged transients in deterministic dynamical systems. We show that resetting the intrinsic dynamics intermittently to a spatial control line that passes through the equilibrium point can dramatically expedite its completion, resulting in a huge reduction in mean transient time and fluctuations around it. Moreover, our study reveals the emergence of an optimal restart time that globally minimizes the mean transient time. We corroborate the results with detailed numerical studies on two canonical setups in deterministic dynamical systems, namely, the Stuart-Landau oscillator and the Lorenz system. The key features-expedition of transient time-are found to be very generic under different resetting strategies. Our analysis opens up a door to control the mean and fluctuations in transient time by unifying the original dynamics with an external stochastic or periodic timer and poses open questions on the optimal way to harness transients in dynamical systems.
Collapse
Affiliation(s)
- Arnob Ray
- Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India
| | - Arnab Pal
- School of Chemistry, Faculty of Exact Sciences and The Center for Physics and Chemistry of Living Systems, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Dibakar Ghosh
- Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India
| | - Syamal K Dana
- Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata 700032, India
| | - Chittaranjan Hens
- Physics and Applied Mathematics Unit, Indian Statistical Institute, 203 B. T. Road, Kolkata 700108, India
| |
Collapse
|
10
|
Barbier-Chebbah A, Benichou O, Voituriez R. Anomalous persistence exponents for normal yet aging diffusion. Phys Rev E 2020; 102:062115. [PMID: 33466050 DOI: 10.1103/physreve.102.062115] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 11/13/2020] [Indexed: 11/07/2022]
Abstract
The persistence exponent θ, which characterizes the long-time decay of the survival probability of stochastic processes in the presence of an absorbing target, plays a key role in quantifying the dynamics of fluctuating systems. So far, anomalous values of the persistence exponent (θ≠1/2) were obtained, but only for anomalous processes (i.e., with Hurst exponent H≠1/2). Here we exhibit examples of ageing processes which, even if they display asymptotically a normal diffusive scaling (H=1/2), are characterized by anomalous persistent exponents that we determine analytically. Based on this analysis, we propose the following general criterion: The persistence exponent of asymptotically diffusive processes is anomalous if the increments display ageing and depend on the observation time T at all timescales.
Collapse
Affiliation(s)
- A Barbier-Chebbah
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS, UPMC, 75005 Paris, France
| | - O Benichou
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS, UPMC, 75005 Paris, France
| | - R Voituriez
- Laboratoire Jean Perrin, CNRS, UPMC, 75005 Paris, France
| |
Collapse
|