Tracer dynamics in single-file system with absorbing boundary
The paper addresses the single-file diffusion in the presence of an absorbing boundary. The emphasis is on an interplay between the hard-core interparticle interaction and the absorption process. The resulting dynamics exhibits several qualitatively new features. First, starting with the exact probability density function for a given particle (a tracer), we study the long-time asymptotics of its moments. Both the mean position and the mean square displacement are controlled by dynamical exponents which depend on the initial order of the particle in the file. Secondly, conditioning on non-absorption, we study the distribution of long-living particles. In the conditioned framework, the dynamical exponents are the same for all particles, however, a given particle possesses an effective diffusion coefficient which depends on its initial order. After performing the thermodynamic limit, the conditioned dynamics of the tracer is subdiffusive, the generalized diffusion coefficient being different from that reported for the system without absorbing boundary.
Consider a single particle diffusing in the semi-infinite one-dimensional channel. The particle escapes from the channel only if it hits the channel boundary situated at the origin. Assuming normal diffusion without any external drift, the mean particle position remains constant in time, its mean escape time is infinite, nevertheless the particle eventually escapes with the probability one. At a given time , starting with an ensemble of all possible particle trajectories, it is interesting to restrict the attention on the subensemble of those paths which have not hit the boundary up to the time . The subensemble is described by the conditional probability density functions (PDFs), the condition being the non-absorption up to the time . The conditioned dynamics exhibits qualitatively different features comparing with the unconditioned one. For instance, the conditioned mean particle position is no longer constant, it growths as . The longer one waits the further from the origin is the typical surviving trajectory. One can say that the conditioning implies an effective force which drags the particle away from the absorbing boundary.
In the present paper we shall investigate the system of hard-core interacting particles diffusing in the semi-infinite one-dimensional channel with the absorbing boundary at the origin. We have three main objectives. First, the hard-core interaction implies an entropic interparticle repulsion and we analyze its effect on the dynamics of the individual particle (a tracer). Second, we are interested in the dynamics of the long-living particles, that is, in the tracer dynamics conditioned on non-absorption. Third, we shall compare the dynamics of the system of particles with that of the corresponding system in thermodynamic limit.
Single particle stochastic dynamics conditioned on non-absorption has been explored extensively in probability theory. A regularly updated bibliography is available in Ref. Pollett (2012). Usually, the conditioning suggests itself in the frame of a biological Hastings (2004), demographic Steinsaltz and Evans (2004) or epidemiological Nåsell (1995) context, where the absorbed diffusion process models the populations undergoing extinction. By the conditioning on non-absorption the focus is shifted on the behavior of the long-surviving paths of the process. It may happen that thus conditioned process converges towards a time-invariant distribution, the so called quasi-stationary distribution. The study of quasi-stationary distributions began with the seminal work of Yaglom on sub-critical Galton-Watson processes Yaglom (1947). For various stochastic processes, the results on existence, uniqueness and other properties of quasi-stationary distributions are reviewed in Cattiaux et al. (2009). Examples of solvable quasi-stationary distributions are the Brownian motion with constant drift absorbed at the origin Mandl (1961); Ferrari et al. (1997), the absorbed logistic Feller diffusion Méléard (2012), and the Wright-Fisher diffusion Huillet (2007). In the demographic context, one specific consequence of the conditioning on non-extinction is the deceleration of the instantaneous mortality rates (mortality rate plateaus) Steinsaltz and Evans (2004).
Many-particle diffusion in one-dimensional channels where the particles are not able to pass each other (the single-file diffusion, SFD) occurs in many systems such as narrow biological pores Hodgkin and Keynes (1955), the channel systems of zeolites Hahn et al. (1996); Kärger and Ruthven (1992), or during the sliding of proteins along the DNA Li et al. (2009); Berg et al. (1981); Bressloff and Newby (2013). In these systems, the diffusion of the tracer is slowed down due to the interparticle interactions. The mean-square displacement of the tracer increases with time as in contrast to its linear increase for the free particle. This was first proved by Harris Harris (1965). Since then, the single-file diffusion was analytically investigated in many different settings including systems in thermodynamic limit Jepsen (1965); Lebowitz and Percus (1967); Levitt (1973), infinite line with a finite number of particles Rödenbeck et al. (1998); Kumar (2008); Aslangul (1998), finite interval Lizana and Ambjörnsson (2008, 2009); Ryabov (2013), particles under the action of external force field Barkai and Silbey (2009); Ryabov and Chvosta (2011). The first-passage problem for a tracer in an infinite system was studied in Sanders and Ambjörnsson (2012). The present paper addresses a different setting: in Ref. Sanders and Ambjörnsson (2012) only the tracer feels the absorbing boundary whereas, in the present paper, each particle can be absorbed. Numerically, the first-passage problem in a driven SFD system was studied in Barma and Ramaswamy (1986).
The present paper continues the study of the model introduced in Ref. Ryabov and Chvosta (2012), where we have investigated the single-file diffusion model including an absorbing boundary. In that paper, we were interested in the time of absorption of the individual particles. In the present paper the central issue is the tracer dynamics. The two prerequisites, i.e., the hard-core interaction and the absorption, are essential for a proper understanding of kinetics of diffusion-limited chemical reactions in crowded environments Liu et al. (2011); Dorsaz et al. (2010); Park et al. (2003); Seki et al. (2011).
The paper is organized as follows. Sec. II contains the definition of the model. In Sec. III, in order to keep the paper self-contained, we present the complete solution of the underlying single-particle case. Secs. IV, V comprise our main results. We first study the unconditioned dynamics of the tracer (Subsecs. IV.1, V.1) and then we condition this dynamics on non-absorption (Subsecs. IV.2, V.2).
Ii Definition of the model
Consider the diffusion of hard-core interacting Brownian particles in a semi-infinite one-dimensional interval with the absorbing boundary at the origin. Initially, particles are distributed along the half-line . During the time evolution, each particle diffuses with the same diffusion constant, . The particles cannot enter the interval from the outside and they are allowed to leave it only through the boundary at the origin. The boundary is perfectly absorbing, i.e., if any particle hits the origin it is absorbed with the probability one. At the initial time , let us label the particles according to ordering of their positions from the left to the right (cf. Fig. 1). We have
where the random variable denotes the position of the -th particle at . The hard-core interaction guarantees that the initial ordering of particles is conserved over time. Notice that the particle No. is the first one that might be absorbed. It is only after this event that the particle No. can approach the origin and be absorbed. Let us denote as the (random) time of the absorption of the -th particle. Then we have
The last inequality, , means that the rightmost particle (and hence any particle) is eventually absorbed with the probability one Ryabov and Chvosta (2012). At the same time, the mean value is infinite Ryabov and Chvosta (2012).
Iii Single diffusing particle
iii.1 Brownian motion absorbed at the origin
Let us take . Suppose that at the initial time the particle is located at the position , . The PDF of the particle’s position at the time conditioned on its initial position is determined by solving the diffusion equation
subject to the absorbing boundary condition at the origin,
and to the initial condition
This problem is readily solved by the method of images Chandrasekhar (1943). The result reads
Having this Green’s function, the time evolution of an arbitrary initial PDF, say , is given by
As for the initial PDF , we only assume that all its moments exist. For the sake of illustrations, we take the particular initial condition
The spatial integral of PDF (7) over the interval equals the survival probability, that is, the probability that the particle has not been absorbed by the time . If we denote by the time of the absorption, we can write
where the last expression stands for the average of the error function Abramowitz and Stegun (1972) with respect to the initial PDF . At the time , the survival probability equals to one. The long-time behavior of may be derived by inserting the power series representation of the the error function into Eq. (9). The expansion is given in Eq. (57) and its first term yields the power-law decay
The prefactor depends on the diffusion constant and on the average initial position of the particle. The sign “” means “is asymptotically equal”.
iii.2 Brownian motion conditioned on non-absorption
According to Eq. (10) the particle will ultimately hit the absorbing boundary at the origin with probability one. Let us now focus on the dynamics of the particle conditioned on non-absorption. By definition
is the PDF for the particle’s position at the time under the condition that the particle has not been absorbed by the time . The power series representation of PDF is given in Eq. (60). It follows that
Notice that this asymptotic representation is non-negative and it is normalized to one on . The distribution with PDF (13) is known as the Rayleigh distribution Krishnamoorthy (2006). Further, the asymptotic PDF (13) is independent of the initial condition and there remains just one length scale associated with the dynamics, the diffusion length . All other length scales which have been introduced by the initial condition are already forgotten.
The first and the second moment of the asymptotic PDF (13) read
The mean position of the surviving trajectories should be compared with the corresponding result for the unconditioned dynamics, that is, with . Provided a trajectory has avoided absorption by the time , it is typically found far from the origin and its typical position growths as . The first-order correction to asymptotic result (14) vanishes as and it depends on the initial condition. Using again Eq. (60), we get
Iv interacting particles
iv.1 Tracer dynamics with absorption
Considering a general number of particles, , the model setting must be completed by the specification of an initial state. We assume the initial joint probability density function for the positions of the particles vanishes outside the domain , and, inside this domain, it is given by
Throughout the paper, all PDFs that have originated in the single-particle problem are denoted by the letter “”. On the other hand, the PDFs in the present -particle problem will be designated by “”. It is a simple consequence of the assumed hard-core interaction that all -particle quantities can be expressed solely through the single-particle PDFs.
The exact PDF for the position of the -th particle, , reads 111See Eq. (35) in Ref. Ryabov and Chvosta (2012)
Apart from the particle labeling, space-time trajectories of hard-core interacting particles are the same as trajectories of noninteracting particles. Hence the probabilistic interpretation behind Eq. (IV.1) can be based on the noninteracting picture. In this picture, the right-hand side (multiplied by ) gives the probability of finding a particle in the interval and, simultaneously, having particles to the right of and particles to the left (including those already absorbed by the boundary). The combinatorial factor accounts for all possible labelings of noninteracting particles.
which can be proved by the direct summation of the expressions (IV.1). In consequence, this equation tells us that the density of particles for the system of interacting particles is the same as that for the system of noninteracting particles. This holds true for all collective properties. However, the dynamics of the individual particles in the two problems is substantially different.
Let us now derive the long-time asymptotics of the PDFs (IV.1). The right-most particle is special. In the long-time limit it behaves in a similar way as the single-diffusing particle Ryabov and Chvosta (2012). In particular, for , the binomial theorem yields
The expression has been written in a way which shows the main asymptotics,
and the first correction, the second term in (IV.1). The correction describes the relaxation towards the main asymptotics and it is negative.
We proceed to the long-time behavior of the -th particle with . We start again with Eq. (IV.1) and rewrite it in the form
If we introduce the renormalized diffusion coefficient,
The only difference is that here we have instead of in (13). The initial order of the particle, , controls, in a decisive way, the main asymptotics. The smaller , the faster is the decay of PDF (for a given ). As a consequence, in the long-time limit, the sum (19) is dominated by the PDF .
Our next goal is analysis of the mean positions of the individual particles. For the rightmost particle, the calculation is based on Eq. (IV.1). We obtain
The main asymptotics is covered by the first term on the right hand side, that is, apart from the multiplication by , the main asymptotics coincides with that for the single particle where we have . The second term describes corrections. As for the remaining particles, the interaction changes the mean-position dynamics. The evaluation of the first moments of the densities (26) yields the main asymptotics
with the prefactor
Thus the initial condition and the total number of particles enters the asymptotics only through the prefactor. Notice that the asymptotics for the -th particle for , differs from that for the rightmost particle (and therefore also from that for the single particle), its mean position asymptotically approaches zero.
In a similar way, we readily obtain the second moments. The results are
with the prefactors
For the rightmost particle, the main asymptotics is proportional to . Interestingly, for , the second moment approaches the nonzero value whereas, for , the second moment decreases towards zero.
valid for any time. In the asymptotic domain, the main asymptotic of the left hand side is covered by the term on the right hand side. Differently speaking, the main asymptotics of the remaining terms in the sum is subdominant with respect to the main asymptotics of the term.
iv.2 Tracer dynamics conditioned on non-absorption
If denotes the time of absorption of the -th particle, its survival probability at the time is defined by
In Ryabov and Chvosta (2012) we have shown that the above integral can be expressed solely through the survival probability of the single-diffusing particle, i.e., through the expression (9). The leading term in the long-time limit is
Being interested in the long-time dynamics of the individual surviving particles, we introduce the conditional PDFs
On the right hand side, the numerator is given in Eq. (IV.1) and the denominator in Eq. (35). In the long-time limit, the fraction greatly simplifies. Dividing the main asymptotics (26) by the leading term in (35) yields
This result is remarkable for its simplicity. The asymptotic dynamics of the -th tracer is the same as the dynamics of a single-diffusing particle with the diffusion coefficients . In other words, the only implication of the hard-core interaction is the renormalization of the diffusion coefficient. The left-most particle diffuses with the smallest effective diffusion coefficients . On the other hand the right-most particle has the same effective diffusion coefficient as a single-diffusing particle, .
The above asymptotic relation means that also the moments of the conditioned dynamics are (except of the value of the diffusion coefficient) simply the moments of the single-diffusing particle. More precisely, using Eq. (37), we get
Finally, notice that in the present conditioned description there exists no simple relationship similar to Eq. (33). The -average of the conditioned -th moments is no more equal to the -th conditioned moment for the single-particle diffusion, i.e.,
V Tracer dynamics in thermodynamic limit
v.1 Tracer dynamics with absorption
We now wish to focus on the dynamics of the tracer in a system of infinite number of particles. First, at the initial time , the particles are distributed randomly on the half-line with a constant density . At the initial instant we activate the absorbing boundary at the origin. Next, the SFD system evolves in time and we are again interested in the dynamics of the -th tagged particle.
Similarly to the previous finite- case, the analysis is based on the exact PDF for the position of the -th particles. The analytical expression which, in the present context, replaces the formula (IV.1) reads
Notice a straightforward probabilistic interpretation of these formulas. The first term on the right-hand side of Eq. (42), , is simply the mean number of the particles absorbed in the time interval . The second term on the right-hand side represents the mean number of particles which are, at the time , diffusing in the space interval . Hence, at the time , stands for the mean number of particles located to the left of the coordinate , including those which were absorbed. In Eq. (41) one recognizes the probability of finding a noninteracting particle in the interval multiplied by the probability that there are particles to the left of (including those already absorbed by the boundary).
The formal derivation of Eq. (41) from Eq. (IV.1) proceeds as follows. At the initial time we assume that particles are homogenously distributed within a finite spatial interval . For a large , the probability of finding a single particle to the right of behaves as
where . We insert this estimation into Eq. (IV.1). The final result (41) follows after performing the thermodynamic limit: , , fixed. Interestingly enough, the passage from Eq. (IV.1) to Eq. (41) is similar in spirit to the well known passage from binomial to the Poisson distribution (the law of rare events).
As for the first two moments of the tracer position we get the asymptotics formulas
v.2 Tracer dynamics conditioned on non-absorption
In the thermodynamic limit, the survival probability of the -th tagged particle still depends on its order and, asymptotically, it assumes the form
which is derived by the spatial integration of PDF (41).
Let us now focus on the dynamics of the tracer which has survived by the time . In the large-time limit, the trace PDF conditioned on non-absorption is given by the ratio of asymptotic expressions (44) and (47):
where we have introduced the generalized diffusion coefficient
Thus the average position of the tracer increases as in contrast to -law as observed for a finite (cf. Eq. (38)). The second moment growths as and hence the tracer dynamics is subdiffusive.
Finally notice that the generalized diffusion coefficient is different as compared to the that obtained in a system without the absorbing boundary Levitt (1973); Leibovich and Barkai (2013). As it was pointed out in Refs. Lizana et al. (2010); Leibovich and Barkai (2013), the coefficient is sensitive to the way the system is prepared. In fact our result (49) indicates that it also depends on boundary conditions used.
Vi Concluding remarks
Returning to the objectives which were outlined in the Introduction, in the long-time limit the following overall picture emerges. First, due to the hard-core repulsion, the particle which possesses a right-hand neighbor feels the (moving) reflecting barrier to the right. The barrier restricts its motion, it reflects the right-moving paths and hence increases the number of left-moving paths. This left-pushing tendency is illustrated by the asymptotic formulae (28), (30), and (31). The mean position and the mean square displacement of the tracer exhibit new dynamical exponents which depend on its initial order. Of course, the rightmost particle has no right-hand neighbor and hence it behaves differently. In the transient regime, its left-hand neighbors still exist and the first particle is pushed to the right. In the asymptotic regime, all other particles have already disappeared and the first one simply undergoes the free diffusion. Second, the conditioning on non-absorption removes a part of the left-moving trajectories from the unconditioned ensemble. Hence it imposes, effectively, the right-oriented drift. Surprisingly enough, the conditioning significantly reduces the effect of the hard-core interaction. The co-operative impact of the both tendencies is behind the asymptotic formulas (38), (39). The conditioned mean position of the tracer growths as regardless its order. The interparticle repulsion manifests itself only through the order-dependent tracer diffusion coefficient. The closer the relative particle position to the boundary the less mobile should that particle be in order to survive for the long times.
The above reasoning holds for the system which initially contains a finite number of particles . In the thermodynamic limit (i.e., assuming initially an infinite number of particles randomly distributed along a half-line with a constant density ), the long-time dynamics of a tracer is rather different. The new features are based on a simple observation that, for any tracer, there are infinite number of particles to the right of it. This implies the -independent exponential damping of the unconditioned moments (45), (46); the initial order appears only in the pre-exponential factor. The conditioned dynamics of a tracer is subdiffusive and independent of (see Eqs. (48), (50), and (51)). This is in parallel to what has been observed in the SFD without an absorbing boundary. Namely, for a finite , Aslangul Aslangul (1998) has shown that, in the long-time limit, the tracer diffusion is normal with the effective diffusion coefficient dependent both on and on . On the other hand, for an infinite , one observes an anomalous diffusion and no -dependence Levitt (1973). The present paper detects the same features in the SFD with an absorbing boundary.
Acknowledgements.This work was supported by the grant SVV-2014-267-305, and by the Charles University Grant Agency (project No. 301311).
Appendix A Asymptotic expansion of the single-particle PDF conditioned on non-absorption
Using the power series representation for the functions inside the averaging brackets, we obtain
The above double sum is treated by the substitution :
where the time-independent coefficients carry all the information concerning the initial condition. Explicitly, they read
into the averaging in (9), we immediately obtain
Interestingly, the numerical factors in the denominators of individual terms of (57) form a sequence (for ), which is A167558 sequence in Sloane’s On-Line Encyclopedia of Integer Sequences Sloane (). This sequence originally emerged in a completely different situation without any obvious connection to the expansion of the error function (see also A167546).
Returning to the main goal of the Appendix, we divide the series (54) by (57). Notice that the prefactor appears in both (54) and (57), therefore it cancels. Representing the fraction by the geometrical series, we finally obtain the sought asymptotic expansion
The asymptotic expansion of the corresponding distribution function, i.e.,
- Pollett (2012) P. Pollett, “Quasi-stationary distributions: a bibliography,” online at http://www.maths.uq.edu.au/ pkp/papers/qsds/qsds.pdf (2012).
- Hastings (2004) A. Hastings, Trends Ecol. Evol. 19, 39 (2004).
- Steinsaltz and Evans (2004) D. Steinsaltz and S. N. Evans, Theor. Pop. Biol. 65, 319 (2004).
- Nåsell (1995) I. Nåsell, Adv. Appl. Probab. 28, 895 (1995).
- Yaglom (1947) A. M. Yaglom, Dokl. Acad. Nauk SSSR (in Russian) 56, 795 (1947).
- Cattiaux et al. (2009) P. Cattiaux, P. Collet, A. Lambert, S. Martínez, S. Méléard, and J. S. Martín, Ann. Probab. 37, 1926 (2009).
- Mandl (1961) P. Mandl, Czechoslovak Mathematical Journal 11, 558 (1961).
- Ferrari et al. (1997) P. A. Ferrari, S. Martínez, and J. S. Martín, J. Stat. Phys. 86, 213 (1997).
- Méléard (2012) S. Méléard, Probability Surveys 9, 340 (2012).
- Huillet (2007) T. Huillet, J. Stat. Mech.-Theory E. 11, 6 (2007).
- Hodgkin and Keynes (1955) A. L. Hodgkin and R. D. Keynes, J. Physiol. 128, 61 (1955).
- Hahn et al. (1996) K. Hahn, J. Kärger, and V. Kukla, Phys. Rev. Lett. 76, 2762 (1996).
- Kärger and Ruthven (1992) J. Kärger and M. Ruthven, Diffusion in Zeolites and in Other Microporous Solids (Wiley, New York, 1992).
- Li et al. (2009) G.-W. Li, O. G. Berg, and J. Elf, Nature Physics 5, 294 (2009).
- Berg et al. (1981) O. G. Berg, R. B. Winter, and P. H. V. Hippel, Biochemistry 20, 6929 (1981).
- Bressloff and Newby (2013) P. C. Bressloff and J. M. Newby, Rev. Mod. Phys. 85, 135 (2013).
- Harris (1965) T. E. Harris, J. Appl. Prob. 2, 323 (1965).
- Jepsen (1965) D. W. Jepsen, J. Math. Phys. 6, 405 (1965).
- Lebowitz and Percus (1967) J. L. Lebowitz and J. K. Percus, Phys. Rev. 155, 122 (1967).
- Levitt (1973) D. G. Levitt, Phys. Rev. A 8, 3050 (1973).
- Rödenbeck et al. (1998) C. Rödenbeck, J. Kärger, and K. Hahn, Phys. Rev. E 57, 4382 (1998).
- Kumar (2008) D. Kumar, Phys. Rev. E 78, 021133 (2008).
- Aslangul (1998) C. Aslangul, Europhys. Lett. 44, 284 (1998).
- Lizana and Ambjörnsson (2008) L. Lizana and T. Ambjörnsson, Phys. Rev. Lett. 100, 200601 (2008).
- Lizana and Ambjörnsson (2009) L. Lizana and T. Ambjörnsson, Phys. Rev. E 80, 051103 (2009).
- Ryabov (2013) A. Ryabov, J. Chem. Phys. 138, 154104 (2013).
- Barkai and Silbey (2009) E. Barkai and R. Silbey, Phys. Rev. Lett. 102, 050602 (2009).
- Ryabov and Chvosta (2011) A. Ryabov and P. Chvosta, Phys. Rev. E. 83, 020106 (2011).
- Sanders and Ambjörnsson (2012) L. P. Sanders and T. Ambjörnsson, J. Chem. Phys. 136, 175103 (2012).
- Barma and Ramaswamy (1986) M. Barma and R. Ramaswamy, J. Stat. Phys. 43, 561 (1986).
- Ryabov and Chvosta (2012) A. Ryabov and P. Chvosta, J. Chem. Phys. 136, 064114 (2012).
- Liu et al. (2011) D.-J. Liu, J. Wang, D. M. Ackerman, I. I. Slowing, M. Pruski, H.-T. Chen, V. S.-Y. Lin, and J. W. Evans, ACS Catal. 1, 751 (2011).
- Dorsaz et al. (2010) N. Dorsaz, C. D. Michele, F. Piazza, P. D. L. Rios, and G. Foffi, Phys. Rev. Lett. 105, 120601 (2010).
- Park et al. (2003) J. Park, H. Kim, and K. J. Shin, J. Chem. Phys. 118, 9697 (2003).
- Seki et al. (2011) K. Seki, M. Wojcik, and M. Tachiya, J. Chem. Phys. 134, 094506 (2011).
- Chandrasekhar (1943) S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943).
- Abramowitz and Stegun (1972) M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th ed. (Dover, New York, 1972).
- Krishnamoorthy (2006) K. Krishnamoorthy, Handbook of Statistical Distributions with Applications, 1st ed. (Chapman & Hall/CRC, New York, 2006).
- (39) See Eq. (35) in Ref. Ryabov and Chvosta (2012).
- Leibovich and Barkai (2013) N. Leibovich and E. Barkai, Phys. Rev. E 88, 032107 (2013).
- Lizana et al. (2010) L. Lizana, T. Ambjörnsson, A. Taloni, E. Barkai, and M. A. Lomholt, Phys. Rev. E 81, 051118 (2010).
- (42) N. J. A. Sloane, “The on-line encyclopedia of integer sequences,” online at https://oeis.org/.