Advances in Mathematical Physics
Volume 2017 (2017), Article ID 6361598, 7 pages
https://doi.org/10.1155/2017/6361598

Research Article

Anomalous Diffusion with an Irreversible Linear Reaction and Sorption-Desorption Process

Maike A. F. dos Santos,¹ Marcelo K. Lenzi,² and Ervin K. Lenzi¹

¹Departamento de Fisica, Universidade Estadual de Ponta Grossa, Avenida General Carlos Cavalcanti 4748, Ponta Grossa, PR, Brazil
²Departamento de Engenharia Qumica, Universidade Federal do Paraná, Av. Cel. Francisco H. dos Santos, 210 Jardim das Americas, 81531-990 Curitiba, PR, Brazil

Correspondence should be addressed to Ervin K. Lenzi; ervinklenzi@gmail.com

Received 8 May 2017; Accepted 21 June 2017; Published 25 July 2017

Academic Editor: Alkesh Punjabi

Copyright © 2017 Maike A. F. dos Santos et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate the diffusion of two different species in a semi-infinite medium considering the presence of linear reaction terms. The dynamics for these species is governed by fractional diffusion equations. We also consider the presence of an adsorption-desorption boundary condition. The solutions for this system are found in terms of the function of Fox and by analyzing the behavior of the mean square displacement a rich class of diffusion processes is verified. In this sense, we show how the surface effects modify the bulk dynamics and promote an anomalous diffusion of system.

1. Introduction

Recently, a great number experiments on biological, physical, and chemical systems have reported a nonlinear time dependence for the mean square displacement connected to anomalous diffusion. Examples can be found in dynamics of biological cell [1–3], crowding in living cell [4, 5], and diffusion in random fractal geometries [6, 7]. In these situations, ( superdiffusion and subdiffusion), in contrast to the usual diffusion characterized by the linear time growing for the mean square displacement; that is, . A framework usually applied, in these contexts, is the continuous time random walk (CTRW) [8]. It may be characterized by waiting-time distributions with a long-tailed behavior which leads one to the fractional diffusion equations as macroscopic formulation [9, 10] for the diffusion process. These equations extend the Fokker-Planck equation [11] and have been successfully applied in several scenarios such as adsorption phenomena [12, 13], diffusion of ions in liquids [14], and diffusion in porous media [15]. Concerning the reaction process in the context of the anomalous diffusion, several works [16–18] have raised the question about the suitable extension when reaction terms are incorporated. One formulation was proposed in [19], where the authors considered the fractional operator only acting on the diffusive term. However, this formulation may result in an unrealistic solution. A similar attempt was made in [20], without success. Other ways of extending the fractional diffusion equation to reaction contexts can be found in [21, 22]. To overcome possible unrealistic scenarios, a mesoscopic formulation based on CTRW was proposed in [23–25] leading us to some satisfactory results and, consequently, working as a guide to incorporate reaction terms in the fractional diffusion equations. In particular, the presence of linear reaction terms [23] yields

where represents a density function of a species, is the diffusion coefficient, and is related to the reaction process present in the bulk; is Riemann-Liouville fractional derivative [26] defined as follows:

for and for the standard form of the differential operator of first-order is recovered. Equation (1) was also generalized to the case of species in [24] by taking into account linear reaction terms, where solutions for the cases with two species in a reversible and irreversible process were obtained.

Here, we analyze the solutions of these processes governed by fractional diffusion equations for the two species, in semi-infinite space, that is, , with sorption-desorption boundary conditions, where an irreversible process, that is, , is present in the bulk. We consider the following equation:

subjected to the boundary conditions

where may represent an sorption-desorption process. In (3), is related to the numbers of species present in the bulk. and are matrix

respectively. It is worth mentioning that results obtained from (3) were used to investigate biological systems [27, 28]. In particular, they have been applied to analyze experimental results of fluorescence recovery after photobleaching (FRAP) in biological systems [29], which is an experimental method widely used to explore binding interactions in cells in vitro and in vivo [30].

The results obtained from (3) subjected to (4) lead us to comprehend how the processes present on the surface and in the bulk (reaction) influence the dynamic of the species present in system. In particular, these results are presented in the next section, Section 2, and show a rich class of processes which can be related to anomalous diffusion. In Section 3, we present our discussion and conclusions.

2. Subdiffusion and Linear Reaction

Then, we focus our attention on the solutions of (3) for the two species, by considering that initially in the bulk one of the species is present and the other is absent; for example, and . Thus, we may analyze how the reaction process influences the spreading of species and the production of species , as well as the surface effects of both.

After performing some calculations, we can expand (3) in a set of coupled equations as follows:

Notice that these equations are coupled by the reaction term: the fact that and have as a particular case the situation worked out in [24]. This set of equations subjected to the boundary conditions given by (4) can be solved using the Fourier , and Laplace , transforms.

Applying these integral transforms in (6), after changing of variables , , and , we obtain

for the initial conditions and . Performing the inverse of Fourier transform yields

So, performing the Laplace integral transform, we obtain for species

For species , we may use the equation , where

with

In (9) and (10), we have the presence of function of Fox [31],

with the contour defined in [8]. It is worth mentioning that the asymptotic behavior of this function in these equations is essentially characterized by stretched exponentials. In fact, after performing some calculations, it is possible to show that

By using these results, we may obtain the survival probability which is related to the quantity of each species in the bulk. In particular, for species , we have that

and for species

We shall analyze the previous results into a situation characterized by a surface with a variable flux for species . The constant is related to the rate of the particles through the surface to the bulk and the parameter represents a characteristic time. From these equations, we observe that for the survival probability decreases exponentially; that is, with the production of species 2, . The cases and are characterized by the absence of reaction with a flux through the surface in contact with the bulk. The behavior of and is illustrated in Figure 1 for different scenarios.

(a)

(b)

Figure 1: This figure illustrates the behavior of (14) and (15) versus for different scenarios. In (a) the blue (solid), black (dashed-dotted), and green (dashed) lines correspond to the cases with , with , and with . In all cases, we have used, for simplicity, . In (b) the blue dashed line corresponds to the cases , , , and . The black dotted line represents the cases , , , and , respectively. The red solid line is the cases , , , and .

Now, we focus our attention on the time dependence manifested by the mean square displacement . Figure 2 shows the behavior of for species , when (6) are considered. For initial times, it is characterized by a power-law, , and, for long times, it decreases exponentially. In fact, performing some calculations, it is possible to show for this case that

in the asymptotic limit of long times. The behavior manifested by species for long times is different from species . It is characterized in the asymptotic limit of by a power-law, , for (see Figure 3). This behavior is essentially the same as the one obtained for the fractional case in absence of reaction term. In both figures (Figures 2 and 3), we also added straight lines to evidence the behaviors manifested during the spreading of the species. In this sense, it is interesting to note that Figure 3(b) shows for the second species three different regimes. In Figure 4, we have shown that the behavior for intermediated times exhibited by the mean square displacement is depending on parameters values. In fact, the difference between the black and the red lines is the value of which in this case shows a direct influence on the behavior of for intermediate times, before the system reaches the asymptotic behavior.

(a)

(b)

Figure 2: Behavior of the mean square displacement for species versus . In (a) the black dashed-dotted line corresponds to the case with and the red dashed line represents the case . In (b), the black dashed-dotted line is the cases and and the blue dashed line corresponds to the case . In all cases, we used and . We also use straight lines to evidence the asymptotic behavior for small times which for is usual and for is anomalous.

(a)

(b)

Figure 3: This figure illustrates the behavior of the mean square displacement for species versus . In (a) the black dotted line corresponds to the cases , , , and . The green dashed line is the case with and the red line is the cases , , , and with . In (b) the black dotted line is the case and the blue solid line corresponds to , respectively. For these cases, we used , , and . We also use straight lines to evidence the different behavior of the mean square displacement.

Figure 4: This figure illustrates the behavior of the mean square displacement for species versus for . The black solid line corresponds to the case and the red dashed line is the case . We consider, for simplicity, , , and . We also use straight lines (blue) to evidence the different behavior of the mean square displacement.

Another interesting behavior present in these systems is illustrated in Figure 5, where plateaus may be exhibited in some time intervals due to the choice of the boundary conditions considered for species . The plateau for this species indicates that during this time interval the flux by the surface related to species is essentially in equilibrium with the reaction process in the bulk, leading us to a stationary state for the distribution related to this species while the species is produced.

(a)

(b)

Figure 5: Behavior of the mean square displacement for species versus , when different choices for the boundary conditions are performed. In (a) the black dashed-dotted line corresponds to the case with and the blue solid line considers . For these cases, for simplicity, we used and . In (b) the black dashed-dotted line is the case with and the blue solid line considers . We also use straight lines to evidence the different behavior of the mean square displacement.

3. Conclusions

We have investigated the behavior of a system composed of two species ( and ) which are governed by fractional diffusion equations in a semi-infinity region. In presence of reaction, we have considered a term characterizing an irreversible process and a boundary condition which can be related to a sorption-desorption process. These equations are coupled by the reaction term which has a direct consequence on the spreading of these species. This point is shown by the behavior of the mean square displacement obtained from (6) for species and which lead us to an anomalous diffusion for both. In particular, for species , we obtained different diffusive behaviors. We have also shown that the intermediated behavior exhibited by the mean square displacement before reaching the asymptotic limit depends on the choice of the parameters present in (3). We have shown the influence of the boundary condition on the system and, in this sense, we also analyzed the behavior of the survival probability for these situations and showed that it depends on the choice of the boundary conditions, that is, , which is directly proportional to the flux of particles through the surface.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors thank the CNPq (Brazilian agency) for partial financial support.

References

F. Höfling and T. Franosch, “Anomalous transport in the crowded world of biological cells,” Reports on Progress in Physics, vol. 76, no. 4, Article ID 046602, 2013.

M. J. Saxton, “A biological interpretation of transient anomalous subdiffusion. II. Reaction kinetics,” Biophysical Journal, vol. 94, no. 3, pp. 760–771, 2008.

A. J. Ellery, R. E. Baker, S. W. McCue, and M. J. Simpson, “Modeling transport through an environment crowded by a mixture of obstacles of different shapes and sizes,” Physica A. Statistical Mechanics and its Applications, vol. 449, pp. 74–84, 2016.

R. Metzler, J.-H. Jeon, and A. G. Cherstvy, “Non-Brownian diffusion in lipid membranes: experiments and simulations,” Biochimica et Biophysica Acta - Biomembranes, vol. 1858, no. 10, pp. 2451–2467, 2016.

S. K. Ghosh, A. G. Cherstvy, D. S. Grebenkov, and R. Metzler, “Anomalous, non-Gaussian tracer diffusion in crowded two-dimensional environments,” New Journal of Physics, vol. 18, no. 1, Article ID 013027, 2016.

Y. Mardoukhi, J.-H. Jeon, and R. Metzler, “Geometry controlled anomalous diffusion in random fractal geometries: looking beyond the infinite cluster,” Physical Chemistry Chemical Physics, vol. 17, no. 44, pp. 30134–30147, 2015.

T. Miyaguchi and T. Akimoto, “Anomalous diffusion in a quenched-trap model on fractal lattices,” Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, vol. 91, no. 1, Article ID 010102, 2015.

R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,” Physics Reports. A Review Section of Physics Letters, vol. 339, no. 1, 77 pages, 2000.

I. M. Sokolov and J. Klafter, “Field-induced dispersion in subdiffusion,” Physical Review Letters, vol. 97, no. 14, Article ID 140602, 2006.

R. Metzler and J. Klafter, “The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics,” Journal of Physics. A. Mathematical and General, vol. 37, no. 31, pp. R161–R208, 2004.

H. Risken and T. Frank, The FokkerPlanck Equation: Method of Solution and Applications, Springer, 2010.

E. K. Lenzi, L. R. Evangelista, G. Barbero, and F. Mantegazza, “Anomalous diffusion and the adsorption-desorption process in anisotropic media,” EPL, vol. 85, no. 2, Article ID 28004, 2009.

G. Barbero and L. R. Evangelista, Adsorption Phenomena and Anchoring Energy in Nematic Liquid Crystals, Taylor and Francis, 2005.

P. A. Santoro, J. L. De Paula, E. K. Lenzi, and L. R. Evangelista, “Anomalous diffusion governed by a fractional diffusion equation and the electrical response of an electrolytic cell,” Journal of Chemical Physics, vol. 135, no. 11, Article ID 114704, 2011.

P. Brault, C. Josserand, J.-M. Bauchire, A. Caillard, C. Charles, and R. W. Boswell, “Anomalous diffusion mediated by atom deposition into a porous substrate,” Physical Review Letters, vol. 102, no. 4, Article ID 045901, 2009.

K. Seki, M. Wojcik, and M. Tachiya, “Fractional reaction-diffusion equation,” Journal of Chemical Physics, vol. 119, no. 4, pp. 2165–2170, 2003.

S. B. Yuste, L. Acedo, and K. Lindenberg, “Reaction front in an reaction-subdiffusion process,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 69, no. 3, Article ID 036126, 2004.

T. A. M. Langlands, B. I. Henry, and S. L. Wearne, “Turing pattern formation with fractional diffusion and fractional reactions,” Journal of Physics Condensed Matter, vol. 19, no. 6, Article ID 065115, 2007.

B. I. Henry and S. L. Wearne, “Fractional reaction-diffusion,” Physica A. Statistical Mechanics and its Applications, vol. 276, no. 3-4, pp. 448–455, 2000.

V. V. Gafiychuk and B. Y. Datsko, “Stability analysis and oscillatory structures in time-fractional reaction-diffusion systems,” Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, vol. 75, no. 5, Article ID 055201, 2007.

E. K. Lenzi, M. A. F. Dos Santos, D. S. Vieira, R. S. Zola, and H. V. Ribeiro, “Solutions for a sorption process governed by a fractional diffusion equation,” Physica A: Statistical Mechanics and its Applications, vol. 443, pp. 32–41, 2016.

X. Li and M. Xu, “A model for reversible reaction in a subdiffusive regime,” Journal of Mathematical Physics, vol. 50, no. 10, Article ID 102708, 102708, 9 pages, 2009.

B. I. Henry, T. A. Langlands, and S. L. Wearne, “Anomalous diffusion with linear reaction dynamics: from continuous time random walks to fractional reaction-diffusion equations,” Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, vol. 74, no. 3, Article ID 031116, 2006.

T. A. Langlands, B. I. Henry, and S. L. Wearne, “Anomalous subdiffusion with multispecies linear reaction dynamics,” Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, vol. 77, no. 2, Article ID 021111, 2008.

T. A. Langlands and B. I. Henry, “Fractional chemotaxis diffusion equations,” Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, vol. 81, no. 5, Article ID 051102, 2010.

I. Podlubny, Fractional Differential Equations, Academic Press, 1999.

E. P. Ipiña and S. P. Dawson, “The effect of reactions on the formation and readout of the gradient of Bicoid,” Physical Biology, vol. 14, no. 1, p. 016002, 2017.

M. Kang, E. Di Benedetto, and A. K. Kenworthy, “Proposed correction to Feder's anomalous diffusion frap equations,” Biophysical Journal, vol. 100, no. 3, pp. 791-792, 2011.

S. B. Yuste, E. Abad, and K. Lindenberg, “A reaction-subdiffusion model of fluorescence recovery after photobleaching (FRAP),” Journal of Statistical Mechanics: Theory and Experiment, vol. 2014, no. 11, Article ID 11014, 2014.

J. Wu, Y. Fang, V. I. Zarnitsyna, T. P. Tolentino, M. L. Dustin, and C. Zhu, “A coupled diffusion-kinetics model for analysis of contact-area FRAP experiment,” Biophysical Journal, vol. 95, no. 2, pp. 910–919, 2008.

A. A. Kilbas, H-Transforms: Theory and Applications, CRC Press, 2004.