**Advances in Pure Mathematics**

Vol.08 No.01(2018), Article ID:82171,13 pages

10.4236/apm.2018.81007

Imperfect Trapping in a Random Walk with Both Species Mobile

Natalia C. Bustos^{1}, Miguel A. Ré^{1,2}^{ }

^{1}Fa.M.A.F.-U.N.C., Córdoba, Argentina

^{2}CIII-F.R.C.-U.T.N., Córdoba, Argentina

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 21, 2017; Accepted: January 28, 2018; Published: January 31, 2018

ABSTRACT

It is presented here a continuous time random walk model for diffusion mediated reactions with both species mobile. The random walk is carried out over an infinite homogeneouos lattice. They are calculated the probability density for the time of reaction of a pair, the reaction rate and the time evolution of the concentration of the majority species. Analytical results are obtained in the Fourier-Laplace transform representation. Known results for a fixed trap are reobtained with appropriate marginal probabilities. It is thus justified Smoluchowski’s original approximation considering the trap at a fixed position and the majority species diffusing with a coefficient sum of the individual coefficients. The results obtained are illustrated by a one dimensional model with bias.

**Keywords:**

Mobile Trap, Diffusion Mediated Reactions, Continuouos Time Random Walk, Trapping Models

1. Introduction

The dynamics of diffusion mediated reactions processes has been extensively studied [1] for many years. It was the objective of these studies the proposal of theoretical models for describing chemical reactions. The schema so developed has also shown useful for the description of diverse phenomena in physics, biology or echology [2] [3] [4] [5] . The reaction process is modeled as a trapping reaction $A+T\to T$ or $A+T\to C$ and it is said diffusion controlled when the reactives have to diffuse to meet in space and then react.

One hundred years ago, Smoluchowski [6] proposed a simplified model assuming a very dilute species (the minority species denoted here as T, the trap) so that it can be considered as a single isolated spherical particle T sorrounded by a swarm of diffusing pointlike particles A. In Smoluchowski’s model particles A diffusion coeficient is ${D}^{\prime}={D}_{A}+{D}_{T}$ the sum of the individual coeficients of A’s and T’s respectively. In the original model the trapping is inmediate (perfect trapping) upon the collision of an A particle with the sphere T. The reaction rate is computed from the flux of A’s into the sphere.

Since this pioneering work different extensions have been proposed to include diffusion in disordered media [7] [8] [9] or to give a better description for the short time behaviour introducing a finite reaction time [10] (imperfect trapping). In imperfect trapping the particles may even separate after en encounter without reacting [11] .

The discrete formulation of the problem is also of interest in many physical problems [12] [13] and trapping models have also been formulated based on Random Walks on lattices. Among these models we mention imperfect trapping [14] [15] with a finite reaction rate and dynamic or gated trapping [16] [17] [18] [19] [20] when the reaction is also modulated by an another independent reaction $T\leftrightarrow {T}^{\star}$ that switches the trap between an active and an inactive state.

Recently it has also been considered an evanescent Random Walk [21] [22] when the walkers have a finite lifetime.

Another extension of Smoluchowski’s model is to consider mobile traps. Koza and Taitelbaum [13] have proposed a model of a mobile imperfect trap in the continuous space with the same diffusion coeficient for both species. Sánchez et al. [23] and Sánchez [24] have considered a mobile trap that diffuses but with no explicit diffusion model. In these works, the problem is considered in the trap reference system and in the laboratory system respectively.

In this article, we propose a generalized Continuous Time Random Walk for imperfect trapping with both species diffusing with different coeficients. We assume separable Random Walks and a time dependent reaction rate. In Section 2, we present a noninteracting two particles Random Walk to introduce the main magnitudes of the model. In Section 3, we include the reaction process in the model through an statiscally independent process. It is calculated in this section the local and global probability density for the time of reaction or the absorption probability density. In Section 4, an initial uniform distribution of the majority species is considered and we calculate the local and global reaction rate. It is verified that the global reaction rate coincides with the reaction rate for an immobile trap. The time dependent concentration of particles A is also calculated. In Section 5, we present the results from the trap reference system. In this frame we recover some results of the immobile trap model. Finally in Section 6, we illustrate the results of our model with a one dimensional Random Walk. Discussion and conclusios follow.

2. Two Particles Random Walk

Let us consider a CTRW of two non-interacting particles on an homogeneous infinite lattice. The particles will be identified as A and T respectively. The position of each particle is given by an integer vector: $x$ for particle A and $y$ for particle T. Each particle makes a generalized, separable, statistically independent RW, with transitions rates ${\lambda}_{A}{p}_{A}\left(x-{x}^{\prime}\right)$ for a transition ${x}^{\prime}\to x$ of particle A and ${\lambda}_{T}{p}_{T}\left(y-{y}^{\prime}\right)$ for a transition ${y}^{\prime}\to y$ of particle T res- pectively. The sojourn probability of the walkers at a given configuration $\left(x\mathrm{,}y\right)$ is

${\Phi}_{0}\left(t\right)={\text{e}}^{-\lambda t}$ (1)

with $\lambda ={\lambda}_{A}+{\lambda}_{T}$ . In this way the diffusion coeficients for a lattice parameter a are respectively ${D}_{A}={a}^{2}{\lambda}_{A}/2\text{d}$ and ${D}_{T}={a}^{2}{\lambda}_{T}/2\text{d}$ .

In turn the probability for a transition from configuration $\left({x}^{\prime}\mathrm{,}{y}^{\prime}\right)$ to con- figuration $\left(x\mathrm{,}y\right)$ between t and $t+\text{d}t$ after a sojourn time t is given by ${\psi}_{0}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)\text{d}t$ with the probability density

${\psi}_{0}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)=\left[{\delta}_{y,{y}^{\prime}}{\lambda}_{A}{p}_{A}\left(x-{x}^{\prime}\right)+{\delta}_{x,{x}^{\prime}}{\lambda}_{T}{p}_{T}\left(y-{y}^{\prime}\right)\right]{\text{e}}^{-\lambda t}$ (2)

Here we are neglecting simultaneous transitions of both particles for being of second order in $dt$ . This is the main magnitude in our model since all the results will be expressed in terms of this density.

Let ${G}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\text{d}t$ be the conditional probability of reaching con- figuration $\left(x\mathrm{,}y\right)$ between t and $t+\text{d}t$ given that the walkers started at the configuration $\left({x}_{0}\mathrm{,}{y}_{0}\right)$ . The probability density ${G}_{0}$ satisfies the recursive rela- tion

$\begin{array}{l}{G}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\\ ={\displaystyle \underset{{x}^{\prime}\mathrm{,}{y}^{\prime}}{\sum}}{\psi}_{0}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)\star {G}_{0}\left({x}^{\prime}\mathrm{,}{y}^{\prime}\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)+{\delta}_{x,{x}^{\prime}}{\delta}_{y,{y}^{\prime}}\delta \left(t-{0}^{+}\right)\end{array}$ (3)

since the walkers may reach the configuration $\left(x\mathrm{,}y\right)$ by a transition from a previous configuration $\left({x}^{\prime}\mathrm{,}{y}^{\prime}\right)$ or $\left(x\mathrm{,}y\right)$ is the initial configuration, these being mutually exclusive events.

We are representing here by the symbol $\star $ the convolution product

$f\left(t\right)\star g\left(t\right)={\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{d}{t}^{\prime}\text{\hspace{0.05em}}f\left(t-{t}^{\prime}\right)g\left({t}^{\prime}\right)$ (4)

${G}_{0}$ is the Green’s function for the diffusion problem of two particles.

An analytic solution for the Green’s function is obtained by taking Fourier transform in the spatial coordinates and Laplace transform in the time variable

${\stackrel{^}{G}}_{0}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)=\frac{{\text{e}}^{i{k}_{A}\cdot {x}_{0}}{\text{e}}^{i{k}_{T}\cdot {y}_{0}}}{1-{\stackrel{^}{\psi}}_{0}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u\right)}$ (5)

with

${\stackrel{^}{\psi}}_{0}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u\right)=\frac{1}{u+\lambda}\left[{\lambda}_{A}{p}_{A}\left({k}_{A}\right)+{\lambda}_{T}{p}_{T}\left({k}_{T}\right)\right]$ (6)

Hereafter the superindex L represents the Laplace transform in the temporal variable

${f}^{L}\left(u\right)={\displaystyle {\int}_{0}^{\infty}}\text{\hspace{0.05em}}\text{d}t\text{\hspace{0.05em}}{\text{e}}^{-ut}f\left(t\right)$ (7)

and the symbol ${}^{\wedge}$ denotes the Fourier transform in spatial variables

$\stackrel{^}{f}\left(k\right)={\displaystyle \underset{x}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\text{e}}^{ik\cdot x}f\left(x\right)$ (8)

In this way the Green’s function Fourier-Laplace ( $\mathcal{F}\mathcal{L}$ ) transform in (5) is given in terms of the $\mathcal{F}\mathcal{L}$ transform of the transition probability density (2) only.

Jointly the conditional probability of a given configuraion at time t of the particles is

${P}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)={\Phi}_{0}\left(t\right)\star {G}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)$ (9)

after reaching the confiugration $\left(x\mathrm{,}y\right)$ at time ${t}^{\prime}$ the walkers must stay in these positions at least for a time $t-{t}^{\prime}$ .

The marginal probability for the A particle position

${P}_{A}\left(x\mathrm{;}t|{x}_{0}\right)={\displaystyle \underset{y}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)$

is in the $\mathcal{F}\mathcal{L}$ transform representation

${\stackrel{^}{P}}_{A}^{L}\left({k}_{A}\mathrm{;}u|{x}_{0}\right)={\Phi}_{A}^{L}\left(u\right)\frac{{\text{e}}^{i{k}_{A}\cdot {x}_{0}}}{1-{\stackrel{^}{\psi}}_{A}^{L}\left({k}_{A}\mathrm{;}u\right)}$ (10)

with

${\psi}_{A}\left(x-{x}_{0}\mathrm{;}t\right)={\lambda}_{A}{p}_{A}\left(x-{x}_{0}\right){\text{e}}^{-{\lambda}_{A}t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Phi}_{A}\left(t\right)={\text{e}}^{-{\lambda}_{A}t}$ (11)

the transition probability density and sojourn probability for the A walker. This result is consistent with the assumption of no interaction between the particles: particle A diffuses with diffusion coeficient ${D}_{A}={a}^{2}{\lambda}_{A}/2\text{d}$ on an infinite lattice of parameter a.

A similar result is obtained for the position marginal probability of particle T position

${\stackrel{^}{P}}_{T}^{L}\left({k}_{A}\mathrm{;}u|{y}_{0}\right)={\Phi}_{T}^{L}\left(u\right)\frac{{\text{e}}^{i{k}_{T}\cdot {y}_{0}}}{1-{\stackrel{^}{\psi}}_{T}^{L}\left({k}_{T}\mathrm{;}u\right)}$ (12)

with

${\psi}_{T}\left(x-{x}_{0}\mathrm{;}t\right)={\lambda}_{T}{p}_{T}\left(x-{x}_{0}\right){\text{e}}^{-{\lambda}_{T}t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Phi}_{T}\left(t\right)={\text{e}}^{-{\lambda}_{T}t}$ (13)

If we have an uniform distribution of A particles, ${c}_{0}$ in each lattice site, the particle concentration at site $x$ at time t is given by

${c}_{A}\left(x;t\right)={\displaystyle \underset{{x}_{0}}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{0}\left(x;t|{x}_{0}\right){c}_{0}$ (14)

and, replacing (5) in (9) and then in (14), we get in the $\mathcal{F}\mathcal{L}$ transform representation

${\stackrel{^}{c}}_{A}^{L}\left(k;u\right)={\left(2\text{\pi}\right)}^{d}\delta \left(k\right)\frac{{c}_{0}}{u}$ (15)

with d the dimension of space.

So it is verified that the uniform distribution is the equilibrium distribution for the noninteracting particles system.

3. The Reaction Process

In this section we include in the model the possibility of the reaction $A+T\to T$ when both particles are at the same position in the lattice. In this case the particles may separate by one of them changing position or they may react with a time dependent reaction rate $\gamma \left(t\right)$ . Assuming diffusion and reaction as statiscally independent processes, the sojourn probability of the walkers at the same position is [15]

${\Phi}_{1}\left(t\right)=\mathrm{exp}\left[-\lambda t-{\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\gamma \left(u\right)\text{d}u\right]$ (16)

Concurrently the probability density for the time of separation of the particles is

${\psi}_{1}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)={\psi}_{0}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right){\Xi}_{r}\left(t\right)$ (17)

and the probability density for the time of reaction is

${\xi}_{A}\left(t\right)={\xi}_{r}\left(t\right){\text{e}}^{-\lambda t}$ (18)

with

${\Xi}_{r}\left(t\right)=exp\left[-{\displaystyle {\int}_{0}^{t}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\gamma \left(s\right)\text{d}s\right]$ (19)

and

${\xi}_{r}\left(t\right)=\gamma \left(t\right){\Xi}_{r}\left(t\right)$ (20)

When particles A and T are not at the same position the state of the system can only change by one of the particles jumping to another position with probability density as given by (2).

We denote by ${G}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\text{d}t$ the Green’s function for the diffusion of the particles when the reaction is allowed. It must satisfy the recursive relation

$\begin{array}{l}{G}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\\ ={\displaystyle \underset{{x}^{\prime}\mathrm{,}{y}^{\prime}}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\psi \left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)\star {G}_{r}\left({x}^{\prime}\mathrm{,}{y}^{\prime}\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)+{\delta}_{x,{x}^{\prime}}{\delta}_{y,{y}^{\prime}}\delta \left(t-{0}^{+}\right)\end{array}$ (21)

This equation is similar to (3) but here

$\psi \left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)=\{\begin{array}{l}{\psi}_{1}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}^{\prime}={y}^{\prime}\\ {\psi}_{0}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}t\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}^{\prime}\ne {y}^{\prime}\end{array}$ (22)

We obtain an analytic solution of Equation (21) in the $\mathcal{F}\mathcal{L}$ transform representation by an extension of the local inhomogeneity method [15] [17] and making use of the identity [12]

${\delta}_{x\mathrm{,}y}={\displaystyle {\int}_{0}^{\text{2\pi}}}\dots {\displaystyle {\int}_{0}^{\text{2\pi}}}\frac{{\text{d}}^{d}\varphi}{{\left(2\text{\pi}\right)}^{d}}exp\left[i\left(x-y\right)\cdot \left(\varphi -k\right)\right]$ (23)

with d the dimension of the space.

The probability density for the time of both particles arriving at the same position $E\left({x}_{e};t|{x}_{0},{y}_{0}\right)={G}_{r}\left({x}_{e},{x}_{e};t|{x}_{0},{y}_{0}\right)$ is in the $\mathcal{F}\mathcal{L}$ transform repre- sentation

${\stackrel{^}{E}}_{r}^{L}\left(k\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)=\frac{I\left(k\mathrm{;}u|{x}_{0}-{y}_{0}\right){\text{e}}^{ik\cdot {y}_{0}}}{{\Xi}_{r}^{L}\left(u+\lambda \right)+{\xi}_{r}^{L}\left(u+\lambda \right)I\left(k\mathrm{;}u|0\right)}$ (24)

with the definition

$I\left(k\mathrm{;}u|{x}_{0}-{y}_{0}\right)=\frac{1}{\left(u+\lambda \right)}{\displaystyle {\int}_{0}^{\text{2\pi}}}\dots {\displaystyle {\int}_{0}^{\text{2\pi}}}\frac{{\text{d}}^{d}\varphi}{{\left(2\text{\pi}\right)}^{d}}\frac{exp\left[i\varphi \cdot \left({x}_{0}-{y}_{0}\right)\right]}{1-{\stackrel{^}{\psi}}_{0}^{L}\left(\varphi \mathrm{,}k-\varphi \mathrm{;}u\right)}$ (25)

while the Green’s function is in the $\mathcal{F}\mathcal{L}$ transform representation

$\begin{array}{l}{\stackrel{^}{G}}_{r}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)\\ ={\stackrel{^}{G}}_{0}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)-\left[\frac{1}{1-{\stackrel{^}{\psi}}_{0}^{L}\left({k}_{A}\mathrm{,}{k}_{T}\mathrm{;}u\right)}-1\right]{\psi}_{r}^{L}\left(u+\lambda \right){\stackrel{^}{E}}_{r}^{L}\left(k\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)\end{array}$ (26)

More details of this calculation is given in Appendix A.

For the reaction to take place the walkers have to be at the same position. So the probability density for the time of reaction at a particular position ${x}_{e}$ or local absorption probability density is given by the convolution product

${A}_{l}\left({x}_{e}\mathrm{;}t|{x}_{0}-{y}_{0}\right)={\psi}_{A}\left(t\right)\star E\left({x}_{e}\mathrm{;}t|{x}_{0}-{y}_{0}\right)$ (27)

since for the reaction to take place both particles have to meet at the same lattice site and react before they move away.

In the $\mathcal{F}\mathcal{L}$ representation we get from (18) and (24)

${\stackrel{^}{A}}_{l}^{L}\left(k\mathrm{;}u|{x}_{0}-{y}_{0}\right)=\frac{I\left(k\mathrm{;}u|{x}_{0}-{y}_{0}\right){\text{e}}^{ik\cdot {y}_{0}}}{\frac{{\Xi}_{r}\left(u+\lambda \right)}{{\xi}_{r}\left(u+\lambda \right)}+I\left(k\mathrm{;}u|0\right)}$ (28)

With this result we may acknowledge Equation (26) as the $\mathcal{F}\mathcal{L}$ transform of

$\begin{array}{l}{G}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\\ ={G}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)-{\displaystyle \underset{{x}^{\prime}}{\sum}}{G}_{1}\left(x\mathrm{,}y\mathrm{;}t|{x}^{\prime}\mathrm{,}{x}^{\prime}\right)\star {A}_{l}\left({x}^{\prime}\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)\end{array}$ (29)

i.e. the Green’s function for the reaction process is the Green’s function of the two particles RW (3) minus the contribution of those realizations with a reaction at a previous time. Here

${G}_{1}\left(x\mathrm{,}y\mathrm{;}t|{x}^{\prime}\mathrm{,}{x}^{\prime}\right)={G}_{0}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)-{\delta}_{x\mathrm{,}{x}^{\prime}}{\delta}_{y,{y}^{\prime}}\delta \left(t-{0}^{+}\right)$ (30)

is the probability density for the time of reaching the configuration $\left(x\mathrm{,}y\right)$ after making at least one transition. For $x=y$ and ${x}_{0}={y}_{0}$ this density is the equivalent of the probability density of the time of return to the origin for one particle random walk.

It is interesting to note that in the case of inmediate reaction upon encounter $\left(\gamma \left(t\right)\to \delta \left(t-{0}^{+}\right)\right)$

${\Xi}_{r}\left(t\right)=\{\begin{array}{l}1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}t=0\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0\end{array},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\xi}_{r}\left(t\right)=\delta \left(t-{0}^{+}\right)$ (31)

The probability density of the meeting time (in this case the first meeting due to inmediate reaction) independent of the lattice position is

${\stackrel{\u02dc}{E}}_{r}\left(t|{x}_{0},{y}_{0}\right)={\displaystyle \underset{{x}_{e}}{\sum}}E\left({\stackrel{\to}{x}}_{e};t|{x}_{0},{y}_{0}\right)={\stackrel{^}{E}}_{r}\left(k=0;t|{x}_{0},{y}_{0}\right)$ (32)

Taking Laplace transform, evaluating Equation (24) at $k=0$ and using (31) we get

${\stackrel{\u02dc}{E}}_{r}^{L}\left(u|{x}_{0},{y}_{0}\right)=\frac{I\left(0;u|{x}_{0}-{y}_{0}\right)}{I\left(0;u|0\right)}$ (33)

${\stackrel{\u02dc}{E}}_{r}^{L}\left(u|{x}_{0},{y}_{0}\right)$ may be identified with Siegert’s formula for the First Passage Time for the lattice position ${y}_{0}$ of a walker that started a RW at ${x}_{0}$ with transition probability density

${\psi}_{E}\left(x-{x}^{\prime}\mathrm{;}t\right)=\left[{\lambda}_{A}{p}_{A}\left(x-{x}^{\prime}\right)+{\lambda}_{T}{p}_{T}\left({x}^{\prime}-x\right)\right]{\text{e}}^{-\lambda t}$ (34)

by recognizing the identity

$I\left(\mathrm{0;}u|{x}_{0}-{y}_{0}\right)={G}^{L}\left({x}_{0}-{y}_{0}\mathrm{;}u\right)$

the Green’s function of one particle Random Walk with transition probability density (34) on an infinite lattice [15] . Thus it is verified Smoluchowski’s assumption that for the global reaction process we may consider one of the particles fixed on the lattice while the other diffuses with coeficient $D={D}_{A}+{D}_{T}$ the sum of the individual diffusion coeficients.

The conditional probability of finding the walkers in configuration $\left(x\mathrm{,}y\right)$ at time t assuming that they started at configuration $\left({x}_{0},{y}_{0}\right)$ is given by the convolution product

${P}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)=\Phi \left(t\right)\star {G}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)$ (35)

with

$\Phi \left(t\right)=\{\begin{array}{ll}{\Phi}_{1}\left(t\right)\hfill & x=y\hfill \\ {\Phi}_{0}\left(t\right)\hfill & x\ne y\hfill \end{array}$

i.e. as in the case of (9), after reaching the configuration $\left(x\mathrm{,}y\right)$ at time ${t}^{\prime}$ the walkers must stay in this configuration at least for a time $t-{t}^{\prime}$ .

4. Reaction Rate

We consider in this section a distribution of particles A with initial uniform concentration ${c}_{0}$ on the lattice. This is the equilibrium concentration in the absence of reaction for the diffusion model under consideration. We assume that the A particles do not interact among them. Then the local reaction rate at ${x}_{e}$ at time t when T starts at ${y}_{0}$ ; i.e. the number of particles that react at ${x}_{e}$ with T between t and $t+\text{d}t$ is ${\mathcal{R}}_{l}\left({x}_{e};t|{y}_{0}\right)\text{d}t$ with the probability density

${\mathcal{R}}_{l}\left({x}_{e};t|{y}_{0}\right)={\displaystyle \underset{{x}_{0}}{\sum}}A\left({x}_{e};t|{y}_{0}\right){c}_{0}$ (36)

In the $\mathcal{F}\mathcal{L}$ transform representation we get from (28)

${\stackrel{^}{\mathcal{R}}}_{l}^{L}\left(k\mathrm{;}u|{y}_{0}\right)={c}_{0}\frac{{\text{e}}^{ik\cdot {y}_{0}}}{u+{\lambda}_{T}\left(1-{\stackrel{^}{p}}_{T}\left(k\right)\right)}{\left[\frac{{\Xi}_{r}^{L}\left(u+\lambda \right)}{{\xi}_{r}^{L}\left(u+\lambda \right)}+I\left(k\mathrm{;}u|0\right)\right]}^{-1}$ (37)

Here

${\stackrel{^}{G}}_{1}^{L}\left(k\mathrm{;}u|{y}_{0}\right)=\frac{{\text{e}}^{ik\cdot {y}_{0}}}{u+{\lambda}_{T}\left(1-{\stackrel{^}{p}}_{T}\left(k\right)\right)}$ (38)

is the $\mathcal{F}\mathcal{L}$ transform of Green’s function of the isolated particle T random walk.

The global reaction rate, defined as ${\mathcal{R}}_{g}\left(t\right)={\displaystyle {\sum}_{x}}{\mathcal{R}}_{l}\left(x;t|{y}_{0}\right)$ is in the Laplace transform representation

${\mathcal{R}}_{g}^{L}\left(u\right)=\frac{{c}_{0}}{u}{\left[\frac{{\Xi}_{r}^{L}\left(u+\lambda \right)}{{\xi}_{r}^{L}\left(u+\lambda \right)}+\frac{1}{u}\right]}^{-1}$ (39)

Note that the global reaction rate does not depend on the initial trap position as should be in a homogeneous lattice.

This reaction rate is coinciding with that of a set of walkers in the presence of a fixed trap when the walkers diffuses with the equivalent transition density (34).

In the presence of reaction the marginal probability of an A particle position is

${P}_{r,A}\left(x;t|{x}_{0},{y}_{0}\right)={\displaystyle \underset{y}{\sum}}\text{\hspace{0.05em}}{P}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)$ (40)

with ${P}_{r}\left(x\mathrm{,}y\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right)$ in (35). From this marginal probability we calculate the concentration of A’s at a particular site

${c}_{r}\left(x\mathrm{;}t|{y}_{0}\right)={\displaystyle \underset{{x}_{0}}{\sum}}\text{\hspace{0.05em}}{P}_{r\mathrm{,}A}\left(x\mathrm{;}t|{x}_{0}\mathrm{,}{y}_{0}\right){c}_{0}$ (41)

$\mathcal{F}\mathcal{L}$ transforming this equation, using (15) and (36) and noting that

${\Phi}_{1}^{L}\left(u\right)-{\Phi}_{0}^{L}\left(u\right)=-{\xi}_{r}^{L}\left(u+\lambda \right){\Phi}_{0}^{L}(\; u\; )$

we get

${\stackrel{^}{c}}_{r}{}^{L}\left(k;u|{y}_{0}\right)={\left(2\text{\pi}\right)}^{d}\delta \left(k\right)\frac{{c}_{0}}{u}-{\Phi}_{0}\left(u\right){\stackrel{^}{\mathcal{R}}}_{l}{}^{L}\left(k;u|{y}_{0}\right)$ (42)

In the space-time representation this result corresponds to

$c\left(x\mathrm{;}t|{y}_{0}\right)={c}_{0}-{\Phi}_{0}\left(t\right)\star {\mathcal{R}}_{l}\left(x\mathrm{;}t|{y}_{0}\right)$ (43)

the local concentration is the uniform initial concentration deducted the contri- bution of those A particles that have reacted at a previous time ${t}^{\prime}<t$ .

5. Description from the Trap Reference System

We present the previous results as described from the trap position. Let us start calculating the conditional probability of the relative position to the trap T of an A particle given by the coordinate $z=x-y$ . From (35) by the variable trans- formation method [25] [26]

${P}_{Z}\left(z;t|{z}_{0}\right)={\displaystyle \underset{y}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{P}_{r}\left(z+y,y;t|{z}_{0}+{y}_{0},{y}_{0}\right)$ (44)

and in the $\mathcal{F}\mathcal{L}$ representation

${\stackrel{^}{P}}_{Z}^{L}\left(k\mathrm{;}u|{z}_{0}\right)={\stackrel{^}{P}}_{0}^{L}\left(k\mathrm{,}-k\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)-\frac{{\Phi}_{0}\left(u\right)}{1-{\stackrel{^}{\psi}}_{0}^{L}\left(k\mathrm{,}-k\mathrm{;}u\right)}{A}_{g}^{L}\left(u\right)$ (45)

From this probablity we calculate the concentration of A’s, that in the $\mathcal{F}\mathcal{L}$ representation is

${\stackrel{^}{c}}_{Z}^{L}\left(k\mathrm{;}u\right)={\left(2\text{\pi}\right)}^{d}\delta \left(k\right)\frac{{c}_{0}}{u}-{\stackrel{^}{P}}_{Z}{}^{L}\left(k\mathrm{;}u|0\right){\mathcal{R}}_{g}^{L}\left(u\right)$ (46)

Going back to the space time representation we may understand this result

${c}_{Z}\left(z\mathrm{;}t\right)={c}_{0}-{P}_{Z}\left(z\mathrm{;}t|0\right)\star {\mathcal{R}}_{g}\left(t\right)$ (47)

in the following way: the concentration of A’s at a given distance from the trap is given by the initial concentration subtracting the contribution of those realizations with a reaction at a previous time.

Note that the concentration as seen from the trap position is coinciding with the concentration for a fixed trap position as could be expected.

From this result we may conclude that the description of the reaction process from the trap position is comparable to the model of the trap at a fixed position on the lattice.

The main difference among the results at the reference systems appears in the A’s concentration as given in (43) in the laboratory frame when compared with (47) in the trap frame.

6. One Dimensional Random Walk

We present here the results of our model when we consider an one dimensional Random Walk as schematically shown in Figure 1. The probability density for the change of configuration is

$\begin{array}{c}{\psi}_{0}\left(x-{x}^{\prime},y-{y}^{\prime};t\right)=[{\delta}_{y,{y}^{\prime}}{\lambda}_{A}\left({p}_{d}{\delta}_{x,{x}^{\prime}+1}+{p}_{i}{\delta}_{x,{x}^{\prime}-1}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{\delta}_{x,{x}^{\prime}}{\lambda}_{T}{p}_{T}\left(\frac{1}{2}{\delta}_{y,{y}^{\prime}+1}+\frac{1}{2}{\delta}_{y,{y}^{\prime}-1}\right)]{\text{e}}^{-\lambda t}\end{array}$ (48)

so that particle A makes a biased Random Walk. With this probability density the integral in (25) is

$I\left(k;u|{x}_{0}-{y}_{0}\right)=\frac{{R}_{c}^{{x}_{0}-{y}_{0}}{\zeta}^{\left|{x}_{0}-{y}_{0}\right|}}{{R}_{p}{R}_{1}}$ (49)

where we have defined

Figure 1. One dimensional Random Walk scheme. Transition rate for the trap T is ${\lambda}_{T}$ and for particle A is ${\lambda}_{A}$ . Particle T makes a symmetric Random Walk while particle A makes a biased Random Walk: ${p}_{i}$ is the probability of jumping to the left and ${p}_{d}$ is the probability of jumping to the right.

$\begin{array}{l}{\Lambda}_{d}=2{\lambda}_{A}{p}_{d}+{\lambda}_{T}{\text{e}}^{-ik}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Lambda}_{i}=2{\lambda}_{A}{p}_{i}+{\lambda}_{T}{\text{e}}^{ik}\\ {R}_{p}=\sqrt{{\Lambda}_{i}{\Lambda}_{d}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{R}_{c}=\sqrt{\frac{{\Lambda}_{i}}{{\Lambda}_{d}}}\\ {R}_{1}=\sqrt{{\left(\frac{u+\lambda}{{R}_{p}}\right)}^{2}-1}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\zeta =\frac{u+\lambda}{{R}_{p}}-{R}_{1}\end{array}$ (50)

For the reaction process we assume an exponential dynamics with a constant reaction rate $\gamma $

${\Xi}_{r}\left(t\right)={\text{e}}^{-\gamma t},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\xi}_{r}\left(t\right)=\gamma {\text{e}}^{-\gamma t}$ (51)

For this one dimensional model the local reaction rate in the $\mathcal{F}\mathcal{L}$ repre- sentation is

${\stackrel{^}{\mathcal{R}}}_{l}^{L}\left(k;u|{y}_{0}\right)=\frac{{\text{e}}^{ik{y}_{0}}}{u+{\lambda}_{T}\left(1-\mathrm{cos}k\right)}\frac{{R}_{p}{R}_{1}}{1+\frac{{R}_{p}{R}_{1}}{\gamma}}{c}_{0}$ (52)

and the A’s concentration in the laboratory frame is

${\stackrel{^}{c}}_{r}^{L}\left(k;u|{y}_{0}\right)=2\text{\pi}\delta \left(k\right)\frac{{c}_{0}}{u}-\frac{{c}_{0}}{u+\lambda}\frac{{\text{e}}^{ik{y}_{0}}}{u+{\lambda}_{T}\left(1-\mathrm{cos}k\right)}\frac{{R}_{p}{R}_{1}}{1+\frac{{R}_{p}{R}_{1}}{\gamma}}$ (53)

In turn the A’s concentration in the trap frame is

${\stackrel{^}{c}}_{Z}^{L}\left(k\mathrm{;}u\right)=2\text{\pi}\delta \left(k\right)\frac{{c}_{0}}{u}-\frac{{c}_{0}}{\left(u+\lambda \right)\left[1-{\stackrel{^}{\psi}}_{A}^{L}\left(k,k;u\right)\right]}\frac{\gamma}{u+\gamma}$ (54)

In Figure 2 we plot the normalized concentration $c/{c}_{0}$ in both reference frames: the laboratory frame and the trap frame as a function of the distance to the initial trap position, $x-{y}_{0}$ , and as a function of the relative position to the trap, $z=x-y$ , respectively. We include several curves with an adimensional time $\lambda t$ as a paramenter.

(a) (b)

Figure 2. Majority species concentration in presence of a mobile imperfect trap vs. relative position to the trap, trap reference frame, (filled symbols in (a), thick lines in (b)) or vs. relative position to the initial position of the trap, laboratory reference frame, (hollowed symbols in (a), thin lines in (b)). It can be appreciated a deeper and narrower depletion zone around the trap in the trap reference frame.

As can be seen there is a greater scavenging at the trap position when the concentration is calculated in the reference frame of the trap. A more extended scavenging region is observed around the initial trap position in the laboratory frame. This apparent disagreement can be understood if we take into account that in the laboratory frame we are considering an average realization of the trap Random Walk while in the trap frame we are just on the trap.

7. Discussion and Conclusions

We have presented a theoretical study of diffusion mediated reaction processes to extend previous treatments of imperfect trapping. Here both species, the minority and majority species, diffuse with arbitrary constants. We have calculated global and local absorption probability densities, reaction rates and concentration of the majority species in the laboratory and trap reference frames. All the results are expressed in terms of the relevant magnitudes of the model, the transition probability density between configurations of the system and the reaction dynamics. We have got analytical results in the Fourier Laplace trans- form representation.

Fixed trap absorption probability density and reaction rate are reobtained when we consider the global behaviour in the laboratory frame or in the trap frame of reference. The main difference between mobile and fixed trap models emerges when we consider the majority species concentration as can be appreciated in the one dimensional case of Section 6. In the trap frame, there is a deeper and spatially more reduced scavenging zone when compared with the concentration in the laboratory reference frame where the scavenging is not so deep but more extended spatially. This difference can be understood when we take into account that in the laboratory frame we are considering an average behaviour of the trap diffusion process, while in the trap frame we are con- sidering a particular realization.

Acknowledgements

The authors thank SeCyT-UTN for partial support of this project.

Cite this paper

Bustos, N.C. and Ré, M.A. (2018) Imperfect Trapping in a Random Walk with Both Species Mobile. Advances in Pure Mathematics, 8, 141-153. https://doi.org/10.4236/apm.2018.81007

References

- 1. Rice, S.A. (1985) Diffusion-Limited Reactions. Elsevier, Amsterdam.
- 2. Blumen, A., Klafter, J. and Zumofen, G. (1986) Optical Spectroscopy of Glasses. Reidel, Dordrecht.
- 3. Goel, N.S. and Richter-Dyn, N. (1974) Stochastic Models in Biology. Academic, New York.
- 4. Alberts, B., Bray, D., Lewis, J., Raff, M., Roberts, K. and Watson, J. (1992) Biologa Molecular de la Célula. Omega, Barcelona.
- 5. Viswanathan, G.N., Buldyrev, S.V., Havlin, S., da Luz, M.G.E., Raposo, E.P. and Stanley, H.E. (1999) Optimizing the Success of Random Searches. Nature, 401, 911-914. https://doi.org/10.1038/44831
- 6. von Smoluchowski, R. (1917) Versucheiner Mathematischen Theorie der Koagulations Kinetic Kolloider Lousungen. Zeitschrift für Physikalische Chemie, 92.
- 7. Scher, H. and Lax, M. (1973) Stochastic Transport in a Disordered Solid. I. Theory. Phys. Rev. B, 7. https://doi.org/10.1103/PhysRevB.7.4491
- 8. Scher, H. and Lax, M. (1973) Stochastic Transport in a Disordered Solid. II. Impurity Conduction, Physical Review B, 7, 4491.
- 9. Bendler, J.T. and Shlesinger, M.F. (1985) The Wonderful World of Stochastics. Elsevier, Amsterdam.
- 10. Collins, F.C. and Kimball, G.E. (1949) Diffusion-Controlled Reaction Rates. Journal of Colloid Science, 4, 425-437. https://doi.org/10.1016/0095-8522(49)90023-9
- 11. Noyes, R.M. (1954) A Treatment of Chemical Kinetics with Special Applicability to Diffusion Controlled Reactions. The Journal of Chemical Physics, 22, 1349-1359. https://doi.org/10.1063/1.1740394
- 12. Aslangul, C. (1999) Mutual Annihilation of Two Diffusing Particles in One- and Two-Dimensional Lattices. Journal of Statistical Physics, 94, 219-240.
- 13. Koza, Z. and Taitelbaum (1998) Spatiotemporal Properties of Diffusive Systems with a Mobile Imperfect Trap. Physical Review E, 57, 237-243. https://doi.org/10.1103/PhysRevE.57.237
- 14. Condat, C.A. (1989) Solution to the Glarum Model with a Finite Relaxation Rate. Zeitschrift für Physik B Condensed Matter, 77, 313-320.
- 15. Ré, M. and Budde, C.E. (2000) Diffusion-Mediated Reactions with a Time Dependent Absorption Rate. Physical Review E, 61, 1110-1120.
- 16. Budde, C.E., Cáceres, M.O. and Ré, M.A. (1995) Transient Behaviour in the Absorption in the Presence of a Non-Markovian Probability Distribution Dynamic Trap. Europhysics Letters, 32, 205.
- 17. Cáceres, M.O., Budde, C.E. and Ré, M.A. (1995) Theory of the Absorption Probability Density of Diffusing Particles in the Presence of a Dynamic Trap. Physical Review E, 52, 3462.
- 18. Ré, M.A., Budde, C.E. and Cáceres, M.O. (1996) Survival Probability in the Presence of a Dynamic Trap. Physical Review E, 54, 4427.
- 19. Spouge, J.L., Szabo, A. and Weiss, G.H. (1996) Single-Particle Survival in Gated Trapping. Physical Review E, 54, 2248. https://doi.org/10.1103/PhysRevE.54.2248
- 20. Bressloff, P.C., Karamched, B.R., Lawley, S.D. and Levien, E. (2017) Diffusive Transport in the Presence of Stochastically Gated Absorption. Physical Review E, 96, Article ID: 022102. https://doi.org/10.1103/PhysRevE.96.022102
- 21. Yuste, S.B., Abad, E. and Linenberg, K. (2013) Exploration and Trapping of Mortal Random Walkers. Physical Review Letters, 110, Article ID: 220603. https://doi.org/10.1103/PhysRevLett.110.220603
- 22. Ré, M. and Bustos, N.C. (2015) Reaction Rate in an Evanescent Random Walkers System. Papers in Physics, 7, Article ID: 070003. https://doi.org/10.4279/pip.070003
- 23. Sánchez, A.D., Rodrguez, M.A. and Wio, H.S. (1998) Results in Trapping Reactions for Mobile Particles and a Single Trap. Physical Review E, 57, 6390.
- 24. Sánchez, A.D. (1999) Trapping Reactions for Mobile Particles and a Trap in the Laboratory Frame. Physical Review E, 59, 5021. https://doi.org/10.1103/PhysRevE.59.5021
- 25. Papoulis, A. (1991) Probability, Random Variables and Stochastic Processes. McGraw-Hill, Singapore.
- 26. Gillespie, D.T. (1982) A Theorem for Physicist in the Theory of Random Variables. American Journal of Physics, 51, 520-533.

Appendix A: Local Inhomogeneity Method

Starting with the recursive relation (21) in the Laplace representation we work out the sum ${\sum}_{{x}^{\prime},{y}^{\prime}}}{G}_{0}^{L}\left(x\mathrm{,}y\mathrm{;}u|{x}^{\prime},{y}^{\prime}\right){G}_{r}^{L}\left({x}^{\prime},{y}^{\prime}\mathrm{;}u|{x}_{0}\mathrm{,}{y}_{0}\right)$ . Noting that in the Laplace representation

$\begin{array}{l}{\psi}_{1}^{L}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}u\right)-{\psi}_{0}^{L}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}u\right)\\ =-{\xi}_{r}^{L}\left(u+\lambda \right){\psi}_{0}^{L}\left(x-{x}^{\prime}\mathrm{,}y-{y}^{\prime}\mathrm{;}u\right)\end{array}$ (55)

and making use of the recursive relation (3) we get

$\begin{array}{l}{G}_{r}^{L}\left(x,y;u|{x}_{0},{y}_{0}\right)\\ ={G}_{0}^{L}\left(x,y;u|{x}_{0},{y}_{0}\right)-{\xi}_{r}\left(u+\lambda \right){\displaystyle \underset{{x}^{\prime},{y}^{\prime}}{\sum}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\delta}_{{x}^{\prime},{y}^{\prime}}\left[{G}_{0}^{L}\left(x,y;u|{x}^{\prime},{y}^{\prime}\right)-{\delta}_{x,{x}^{\prime}}{\delta}_{y,{y}^{\prime}}\right]\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\times {G}_{r}^{L}\left({x}^{\prime},{y}^{\prime};u|{x}_{0},{y}_{0}\right)\end{array}$ (56)

It should be noted that in the right hand side it appears ${G}_{r}^{L}\left({x}^{\prime},{x}^{\prime};u|{x}_{0},{y}_{0}\right)={E}_{r}^{L}\left({x}^{\prime};u|{x}_{0},{y}_{0}\right)$ as defined in (24). Then we evaluate (56) at $x=y$ and take Fourier transform obtaining (24). For this calculation use have been made of identity (23) and definition (25).