Single-particle survival in gated trapping

Controlling Uncertainty of Empirical First-Passage Times in the Small-Sample Regime

We derive general bounds on the probability that the empirical first-passage time τ[over ¯]_{n}≡∑_{i=1}^{n}τ_{i}/n of a reversible ergodic Markov process inferred from a sample of n independent realizations deviates from the true mean first-passage time by more than any given amount in either direction. We construct nonasymptotic confidence intervals that hold in the elusive small-sample regime and thus fill the gap between asymptotic methods and the Bayesian approach that is known to be sensitive to prior belief and tends to underestimate uncertainty in the small-sample setting. We prove sharp bounds on extreme first-passage times that control uncertainty even in cases where the mean alone does not sufficiently characterize the statistics. Our concentration-of-measure-based results allow for model-free error control and reliable error estimation in kinetic inference, and are thus important for the analysis of experimental and simulation data in the presence of limited sampling.

Gated reactions in discrete time and space

How much time does it take for two molecules to react? If a reaction occurs upon contact, the answer to this question boils down to the classic first-passage time problem: find the time it takes for the two molecules to meet. However, this is not always the case as molecules switch stochastically between reactive and non-reactive states. The reaction is then said to be "gated" by the internal states of the molecules involved, which could have a dramatic influence on kinetics. A unified, continuous-time, approach to gated reactions on networks was presented in a recent paper [Scher and Reuveni, Phys. Rev. Lett. 127, 018301 (2021)]. Here, we build on this recent advancement and develop an analogous discrete-time version of the theory. Similar to continuous-time, we employ a renewal approach to show that the gated reaction time can always be expressed in terms of the corresponding ungated first-passage and return times, which yields formulas for the generating function of the gated reaction-time distribution and its corresponding mean and variance. In cases where the mean reaction time diverges, we show that the long-time asymptotics of the gated problem is inherited from its ungated counterpart. However, when molecules spend most of their time non-reactive, an interim regime of slower power-law decay emerges prior to the terminal asymptotics. The discretization of time also gives rise to resonances and anti-resonances, which were absent from the continuous-time picture. These features are illustrated using two case studies that also demonstrate how the general approach presented herein greatly simplifies the analysis of gated reactions.

Unified Approach to Gated Reactions on Networks

For two molecules to react they first have to meet. Yet, reaction times are rarely on par with the first-passage times that govern such molecular encounters. A prime reason for this discrepancy is stochastic transitions between reactive and nonreactive molecular states, which results in effective gating of product formation and altered reaction kinetics. To better understand this phenomenon we develop a unifying approach to gated reactions on networks. We first show that the mean and distribution of the gated reaction time can always be expressed in terms of ungated first-passage and return times. This relation between gated and ungated kinetics is then explored to reveal universal features of gated reactions. The latter are exemplified using a diverse set of case studies which are also used to expose the exotic kinetics that arises due to molecular gating.

Diffusive transport in the presence of stochastically gated absorption

We analyze a population of Brownian particles moving in a spatially uniform environment with stochastically gated absorption. The state of the environment at time t is represented by a discrete stochastic variable k(t)∈{0,1} such that the rate of absorption is γ[1-k(t)], with γ a positive constant. The variable k(t) evolves according to a two-state Markov chain. We focus on how stochastic gating affects the attenuation of particle absorption with distance from a localized source in a one-dimensional domain. In the static case (no gating), the steady-state attenuation is given by an exponential with length constant sqrt[D/γ], where D is the diffusivity. We show that gating leads to slower, nonexponential attenuation. We also explore statistical correlations between particles due to the fact that they all diffuse in the same switching environment. Such correlations can be determined in terms of moments of the solution to a corresponding stochastic Fokker-Planck equation.

Effect of stochastic gating on channel-facilitated transport of non-interacting and strongly repelling solutes

Ligand- or voltage-driven stochastic gating-the structural rearrangements by which the channel switches between its open and closed states-is a fundamental property of biological membrane channels. Gating underlies the channel's ability to respond to different stimuli and, therefore, to be functionally regulated by the changing environment. The accepted understanding of the gating effect on the solute flux through the channel is that the mean flux is the product of the flux through the open channel and the probability of finding the channel in the open state. Here, using a diffusion model of channel-facilitated transport, we show that this is true only when the gating is much slower than the dynamics of solute translocation through the channel. If this condition breaks, the mean flux could differ from this simple estimate by orders of magnitude.

Diffusive flux in a model of stochastically gated oxygen transport in insect respiration

Oxygen delivery to insect tissues is controlled by transport through a branched tubular network that is connected to the atmosphere by valve-like gates, known as spiracles. In certain physiological regimes, the spiracles appear to be randomly switching between open and closed states. Quantitative analysis of this regime leads a reaction-diffusion problem with stochastically switching boundary condition. We derive an expression for the diffusive flux at long times in this problem. Our approach starts with the derivation of the passage probability for a single particle that diffuses between a stochastically gated boundary, which models the opening and closing spiracle, and the perfectly absorbing boundary, which models oxygen absorption by the tissue. This passage probability is then used to derive an expression giving the diffusive flux as a function of the geometric parameters of the tube and characteristic time scales of diffusion and gate dynamics.

Stochastically gated diffusion-limited reactions for a small target in a bounded domain

We calculate the reaction rate for stochastically gated ligands diffusing in a two-dimensional and a three-dimensional bounded domain with a single small target. Each ligand independently switches between an open and a closed state according to a two-state Markov process; a reaction between ligand and target can only occur when the former is an open state. In the large-time limit the reaction-rate is an exponentially decaying function of time, whose rate of decay is given by the principal eigenvalue of the Laplacian. We calculate the principal eigenvalue using matched asymptotics and determine the leading-order reduction in the reaction rate due to stochastic gating. We also develop a probabilistic interpretation of the reaction rate in terms of the first-passage time density to the target.

Narrow-escape-time problem: the imperfect trapping case

We present a master equation approach to the narrow escape time (NET) problem, i.e., the time needed for a particle contained in a confining domain with a single narrow opening to exit the domain for the first time. We introduce a finite transition probability, ν, at the narrow escape window, allowing the study of the imperfect trapping case. Ranging from 0 to ∞, ν allowed the study of both extremes of the trapping process: that of a highly deficient capture and situations where escape is certain ("perfect trapping" case). We have obtained analytic results for the basic quantity studied in the NET problem, the mean escape time, and we have studied its dependence in terms of the transition (desorption) probability over (from) the surface boundary, the confining domain dimensions, and the finite transition probability at the escape window. Particularly we show that the existence of a global minimum in the NET depends on the "imperfection" of the trapping process. In addition to our analytical approach, we have implemented Monte Carlo simulations, finding excellent agreement between the theoretical results and simulations.

Intermittent pathways towards a dynamical target

In this paper, we investigate the quest for a single target, which remains fixed in a lattice, by a set of independent walkers. The target exhibits fluctuating behavior between a trap and an ordinary site of the lattice, whereas the walkers perform an intermittent kind of search strategy. Our searchers carry out their movements in one of two states, between which they switch randomly. One of these states (the exploratory phase) is a symmetric nearest-neighbor random walk and the other state (the relocating phase) is a symmetric next-nearest-neighbor random walk. By using the multistate continuous-time random-walk approach we are able to show that for dynamical targets, the intermittent strategy (despite the simplicity of the kinetics chosen for searching) improves detection, in comparison to displacements in a single state. We have obtained analytic results, which can be numerically evaluated, for the survival probability and for the lifetime of the target. Thus, we have studied the dependence of these quantities both in terms of the transition probability that describes the dynamics of the target and in terms of the parameter that characterizes the walkers' intermittency. In addition to our analytical approach, we have implemented Monte Carlo simulations, finding excellent agreement between the theoretical-numerical results and simulations.

Kinetics of diffusion-assisted reactions in microheterogeneous systems

This review is focused on the basic theory of diffusion-assisted reactions in microheterogeneous systems, from porous solids to self-organized colloids and biomolecules. Rich kinetic behaviors observed experimentally are explained in a unified fashion using simple concepts of competing distance and time scales of the reaction and the embedding structure. We mainly consider pseudo-first-order reactions, such as luminescence quenching, described by the Smoluchowski type of equation for the reactant pair distribution function with a sink term defined by the reaction mechanism. Microheterogeneity can affect the microscopic rate constant. It also enters the evolution equation through various spatial constraints leading to complicated boundary conditions and, possibly, to the reduction of dimensionality of the diffusion space. The reaction coordinate and diffusive motion along this coordinate are understood in a general way, depending on the problem at hand. Thus, the evolution operator can describe translational and rotational diffusion of molecules in a usual sense, it can be a discrete random walk operator when dealing with hopping of adsorbates in solids, or it can correspond to conformational fluctuations in proteins. Mathematical formulation is universal but physical consequences can be different. Understanding the principal features of reaction kinetics in microheterogeneous systems enables one to extract important structural and dynamical information about the host environments by analyzing suitably designed experiments, it helps building effective strategies for computer simulations, and ultimately opens possibilities for designing systems with controllable reactivity properties.

Kinetics of stochastically gated diffusion-limited reactions and geometry of random walk trajectories

In this paper we study the kinetics of diffusion-limited, pseudo-first-order A+B-->B reactions in situations in which the particles' intrinsic reactivities are not constant but vary randomly in time. That is, we suppose that the particles are bearing "gates" which fluctuate in time, randomly and independently of each other, between two states-an active state, when the reaction may take place between A and B particles appearing in close contact; and a blocked state, when the reaction is completely inhibited. We focus here on two customary limiting cases of pseudo-first-order reactions-the so-called target annihilation and the Rosenstock trapping model-and consider four different particular models, such that the A particle can be either mobile or immobile or gated or ungated, and ungated or gated B particles can be fixed at random positions or move randomly. All models are formulated on a d-dimensional regular lattice, and we suppose that the mobile species perform independent, homogeneous, discrete-time lattice random walks. The model involving a single, immobile, ungated target A and a concentration of mobile, gated B particles is solved exactly. For the remaining three models we determine exactly, in the form of rigorous lower and upper bounds showing the same N dependence, the large-N asymptotical behavior of the probability that the A particle survives until the Nth step. We also realize that for all four models studied here the A particle survival probability can be interpreted as the moment generating function of some functionals of random walk trajectories, such as, e. g., the number of self-intersections, the number of sites visited exactly a given number of times, the "residence time" on a random array of lattice sites, etc. Our results thus apply to the asymptotic behavior of corresponding generating functions which are not known as yet.

Diffusion-mediated reactions with a time-dependent absorption rate

Diffusion-mediated reactions models are particularly useful for the characterization of physical, chemical, and biological problems. In this paper we present a theoretical study of the absorption probability density, survival probability, and reaction rate for diffusion-mediated reactions models with a time-dependent finite absorption rate (an extension of a model usually referred to as the "imperfect trap model"). The results are obtained by means of the formalism of continuous time random walk on a lattice and considering a general reaction dynamics upon encounter of the reactives. First jump probability densities are included to take initial conditions into account. Previous results presented by Collins and Kimball [J. Colloid. Sci. 4, 425 (1949)] and Noyes [J. Chem. Phys. 22, 1349 (1954)] are reobtained for the particular case of a time-independent absorptivity. Short and long time behaviors are analyzed resulting, in particular, in that the long time behavior of the absorption probability density exhibits the same time dependence as the first passage time density. The results obtained are illustrated by considering a one-dimensional model with consequent discussion.