Encounter-based reaction-subdiffusion model II: partially absorbing traps and the occupation time propagator

Paul C Bressloff
2023-03-19
Abstract:In this paper we develop an encounter-based model of reaction-subdiffusion in a domain $\Omega$ with a partially absorbing interior trap $\calU\subset \Omega$. We assume that the particle can freely enter and exit $\calU$, but is only absorbed within $\calU$. We take the probability of absorption to depend on the amount of time a particle spends within the trap, which is specified by a Brownian functional known as the occupation time $A(t)$. The first passage time (FPT) for absorption is identified with the point at which the occupation time crosses a random threshold $\widehat{A}$ with probability density $\psi(a)$. Non-Markovian models of absorption can then be incorporated by taking $\psi(a)$ to be non-exponential. The marginal probability density for particle position $\X(t)$ prior to absorption depends on $\psi$ and the joint probability density for the pair $(\X(t),A(t))$, also known as the occupation time propagator. In the case of normal diffusion one can use a Feynman-Kac formula to derive an evolution equation for the propagator. However, care must be taken when combining fractional diffusion with chemical reactions in the same medium. Therefore, we derive the occupation time propagator equation from first principles by taking the continuum limit of a heavy-tailed CTRW. We then use the solution of the propagator equation to investigate conditions under which the mean FPT (MFPT) for absorption within a trap is finite. We show that this depends on the choice of threshold density $\psi(a)$ and the subdiffusivity. Hence, as previously found for evanescent reaction-subdiffusion models, the processes of subdiffusion and absorption are intermingled.
Statistical Mechanics
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to develop a reaction - sub - diffusion model based on encounters and study the behavior of particles in partially absorbing internal traps \( U \subset \Omega \). Specifically, the paper focuses on the following key issues: 1. **Particle behavior in partially absorbing traps**: - Particles can freely enter and exit the trap \( U \), but are only absorbed within \( U \). - The probability of absorption depends on the time the particle stays in the trap, which is described by a Brownian functional (called the occupation time \( A(t) \)). 2. **Definition of the first - passage time (FPT)**: - The first - passage time (FPT) of absorption is defined as the moment when the occupation time \( A(t) \) exceeds a random threshold \( \hat{A} \), where the probability density function of \( \hat{A} \) is \( \psi(a) \). 3. **Non - Markovian absorption models**: - By choosing different \( \psi(a) \), non - Markovian absorption models can be incorporated, which are more common in many practical applications. 4. **Derivation of the occupation - time propagator equation**: - For normal diffusion, the evolution equation of the occupation - time propagator can be derived through the Feynman - Kac formula. - For sub - diffusion, the occupation - time propagator equation needs to be derived from first principles through the continuous limit of continuous - time random walks (CTRW) with heavy tails. 5. **Finiteness conditions of the mean first - passage time (MFPT)**: - Explore under what conditions the mean first - passage time (MFPT) within the trap is finite, which depends on the threshold density \( \psi(a) \) and sub - diffusivity. ### Main contributions - **Theoretical framework**: Proposed a reaction - sub - diffusion model based on encounters, which is applicable to the behavior of particles in partially absorbing traps. - **Mathematical derivation**: Derived in detail the occupation - time propagator equation and explored its solutions under different conditions. - **Application example**: Taking spherical traps as an example, studied the finiteness conditions of the mean first - passage time under specific geometric configurations. ### Conclusions Through theoretical analysis and mathematical derivation, the paper reveals the interaction between sub - diffusion and absorption processes, providing a new perspective for understanding reaction - diffusion processes in complex media. This is of great significance for practical problems such as the transport of neurotransmitter receptors in neurons.