Abstract

We investigate the influence of the surface effects on a diffusive process by considering that the particles may be sorbed or desorbed or undergo a reaction process on the surface with the production of a different substance. Our analysis considers a semi-infinite medium, where the particles may diffuse in contact with a surface with active sites. For the surface effects, we consider integrodifferential boundary conditions coupled with a kinetic equation which takes non-Debye relation process into account, allowing the analysis of a broad class of processes. We also consider the presence of the fractional derivatives in the bulk equations. In this scenario, we obtain solutions for the particles in the bulk and on the surface.

1. Introduction

The description of the dynamic processes present on the surface such as adsorption-desorption or reactions and their influence on the diffusion is very important due to the broadness of applications in several fields such as optimization of industrial processes, electronic devices, the action of pharmaceuticals in the organism, and materials science [1]. Particularly in chemical engineering applications of reactions that occur on solid surfaces are of the great interest and advances have recently been achieved by using the heterogeneous catalysis in the pollutants removal [2], with emphasis on oxidative pathways involving the formation of radicals [35], and in the biodiesel production [68]. In this sense, it is worth not only studying chemical reactions on substrates, but also analyzing different phenomena occurring in a given surface, for example, the adsorption of proteins [912], which can be controlled by nanoroughness on the surface of the material [13, 14] or hybridization of DNA [1517], among others. In these systems, different species may diffuse and react [1821], which implies considering suitable changes in the diffusion equation or in the boundary conditions to account for the processes of interest [18, 19, 22, 23]. These contexts usually are analyzed by using the standard approaches related to Markovian processes. Therefore, extensions of these approaches become very important in order to consider non-Markovian processes where anomalous diffusion or non-Debye relaxations on the kinetic processes are present.

Here, we analyze the surface effects on a diffusion process by considering one-dimensional semi-infinite media in contact with a surface, which may adsorb (and desorb) particles or absorb the particles to perform a reaction process with the formation of another kind of particle. Figure 1 illustrates the scenario analyzed in this manuscript, where two species (particles or substances) and diffuse in the bulk which is in contact with a surface. Species may be adsorbed or sorbed by the surface. For the first process, species is removed from the bulk by the surface and after some time is desorbed with a characteristic time. In the other case, when species is sorbed, a reaction process may occur with the formation of species . In this sense, the reaction only occurs on the surface, for instance, in the presence of an specific catalyst. Thus, species reacts promoting the formation of species , in other words, following a first-order irreversible reaction . This reaction can be described by a kinetic equation [21, 24, 25]. In order to describe these processes, the model considered here is a set of coupled equations that may represent reaction, adsorption, and desorption processes of a substrate on a surface and also the diffusion process in the bulk. These developments are performed in Section 2 by considering generalized diffusion equations [2640] for the species ( and ) in the bulk coupled with kinetic equations as boundary conditions. The discussion and conclusions are presented in Section 3.

2. Diffusion and Surface Effects

Let us start our analysis about the surface effects on the diffusion processes by considering the system governed in the bulk by the following generalized diffusion equationsandwith , where and are the generalized diffusion coefficients related to particles and , respectively. The quantities and represent the density of each particle present in the bulk and , in (1) and (2), is the integrodifferential operator:where defines how the past history of each system influences the time evolution of and depending on the choice of different fractional differential operators may be obtained. For example, the casecorresponds to the well-known Riemann-Liouville fractional operator [41] for , which has been applied in several situations related to anomalous diffusion [4245]. The exponential kernelwith and being a normalization factor, corresponds to the fractional operator of Caputo-Fabrizio [46]. Further possibilities for describing the kernel are discussed by Gómez-Aguilar et al. [47], in particular the kernel , where is the Mittag-Leffler function [41], which leads us to Atangana-Baleanu fractional operator [4850]. It has been considered in the Bloch system [51], by resulting in a different behavior from the one obtained in the usual context where differential operators of integer order are considered. These operators, by taking a variable order into account, have been used to extend the Gray-Scott reaction-diffusion model which describe irreversible reaction between two species. In this scenario, the Atangana-Baleanu-Caputo fractional differential operator has shown a faster stabilization behavior and Liouville-Caputo fractional differential operator presented a stronger memory effect [52]. From these choices for , we observe that the Riemann-Liouville operator has a singularity at the origin (), while the recently proposed Caputo-Fabrizio and Atangana-Baleanu operators are nonsingular [4648, 5355], and the last one may manifest different regimes of diffusion.

On the surface, we consider that the processes are governed by the following equations:andHere is related to the rate of particle sorption by the surface which by a reaction process produces substance . In addition to (6) and (7), we also have the boundary conditions and to solve (1) and (2). In (6), represents the density of particles which are adsorbed by the surface. For the adsorption and desorption processes on the surface, we assume that they may be modeled by the following kinetic equation [56]:In (8), is the bulk density just in front of the surface, which may be adsorbed by the surface. The parameter is connected to the adsorption phenomenon, being related to a characteristic adsorption time , and is a kernel that governs the desorption phenomenon. Thus, the surface density of adsorbed particles depends on the bulk density of particles just in front of the membrane, and on the surface density of particles already sorbed [56]. Equation (8) also extends the usual kinetic equations of first order to situations characterized by unusual relaxations, i.e., non-Debye relaxations for which a nonexponential behavior of the densities can be obtained, depending on the choice of the kernels [56, 57].

From (1), (2), (6), and (7), it is possible to show thatandIn (9) the term implies the removal of the particles from the bulk by the surface to promote the production of species by a reaction process. Thus, (11) shows that the mass (number of particles) variation on species is connected to the variations on species . In particular, the negative sign shows that the variation of particles of one species produces an opposite variation on the other. The particles produced by a reaction process on the surface are being released from the surface to the bulk. Equations (9) and (10) imply which consequently yieldsi.e., a direct consequence of the conservation of the total number of particles present in the system.

These systems of coupled equations can be solved by using the Laplace transform and the Green function approach. In fact, by applying the Laplace transform ( and ) it is possible to simplify the previous equations and obtain thatandwhich are subjected, in the Laplace domain, to the boundary conditions:with and , where . By using the Green function approach, the solutions for and can found by considering the previous equations and they are given byandwhere and correspond to the Green function for each species and represents the quantity of particles which may initially be present on the surface. These Green functions are obtained by solving the following equation:by taking into account the suitable boundary conditions for each species related to (17) and (18). Thus, the solutions given by (17) and (18) are obtained from the combination of (15), (16), and (19) by taking into account their boundary conditions, following the procedure employed in [58]. It is also interesting to note that the boundary conditions chosen for the Green function were performed in order to evidence the adsorption-desorption process manifested by the surface on the distribution . For species , by performing some calculations, it possible to show that the Green function is given bywhen the boundary conditions, which are consistent with (17),and are considered to solve (21). Note that the boundary condition used to solve the equation for the Green function incorporates the reaction on the surface. The term related to the adsorption-desorption process was not incorporated in (21) to be evident on the solution as an additional term. In (20), the first term corresponds to the spreading of the initial condition and the second term shows the influence of the surface on the spreading of the system, which in this case corresponds to a reaction sorption process where species is removed from the system. It is worth mentioning that depending on the choice of (or ) the previous Green function may present a stationary solution. In particular, for case with , for (i.e., for long times), we have, for the Green function, the following stationary behavior:which in this case is independent of , where . This feature is interesting and shows that the surface effects with this characteristic play an important role at intermediate times. For the case (or ), we also have an stationary solution and it is given bywhich is different from (22) due to the boundary conditions used to obtain the Green function. The scenario is typical, for example, of the Caputo-Fabrizio operator, which is characterized by an exponential kernel as mentioned before.

For species , we have thatwhen the boundary conditions, which are consistent with Eq. (18),are employed to solve (21). Equation (24) may also present a stationary solution depending on the choice of , similar to (20).

In order to perform the inverse of Laplace transform and obtain the time dependent behavior of (20) and (24), we first consider the Riemann-Liouville fractional time derivative and after that the Caputo-Fabrizio fractional time derivative. We also consider that ; i.e., the sorption of species by the surface occurs in a constant rate. For the Riemann-Liouville fractional time derivative, we have that and, consequently,with , where is the Fox function [59] and is the generalized Mittag-Leffler function [59] (see Figure 2). For the Caputo-Fabrizio fractional time operator, we have thatwith . By substituting the previous equation in (20) and performing some calculations with , we obtain thatfor the Caputo-Fabrizio fractional time operator with :where is the complementary error function and . Figure 3 shows the behavior of (29) for different values of . In Figure 4, we present the behavior of (29) for different times in order to show that for long times a stationary state is reached. For (24), we havefor the Riemann-Liouville operator (see Figure 5) andfor the Caputo-Fabrizio fractional time operator (see Figure 6).

In (17) and (18), the first term promotes the spreading of the initial condition and the other terms represent the influence of the surface on the diffusive process in the bulk. Thus, we need to quantify these terms, in order to obtain the solution for each species. From (17), after some calculations, it is also possible to show thatwhere . By using the previous result, it is possible to show that the quantity of adsorbed particles by the surface is given by

By using (33) and (34) and by applying the inverse Laplace transform in (13) and (14), it is possible to obtain and . We consider, as performed for the Green functions, the fractional operators characterized by (4) and (5) in order to obtain the inverse Laplace transform of (33) and (34). For the first case, i.e., the Riemann-Liouville fractional time derivative, we obtain for (33) and (34) thatFor the case characterized by (4), i.e., the Caputo-Fabrizio fractional time operator, we obtain thatThe quantities and presented in previous equations are defined as follows:with

By using the previous equations, the inverse of Laplace transform of (17) and (18)andEquations (41) and (42) are the solutions for the scenario analyzed here, which is characterized by a diffusion with processes occurring on a surface. By using these results, it is possible to obtain the survival probability, which is defined as and it is related to the quantity of particles each species present in the bulk. In Figure 7, we illustrate the time behavior of the survival probability for species and by considering different fractional time operators. We observe that for the fractional operator related to the exponential kernel a stationary state is reached. Another interesting quantity to be analyzed is the mean square displacement, i.e., , related to the spreading of the distributions of each species. For example, in absence of adsorption-desorption process ( and ) with , species manifests two different regimes, when the Riemann-Liouville fractional time derivative is considered. In particular, for small times and for long times . In this context, for the Caputo-Fabrizio fractional time operator, we have for small times and for long times. For species , the behavior of the mean square displacement is illustrated in Figure 8 and, similar to species , the Riemann-Liouville fractional time derivative presents two different regimes. For the Caputo-Fabrizio fractional time operator, we observe a superdiffusion (see Figure 8(a)) or a transient (see Figure 8(b)) before reaching a stationary state. The initial behavior observed for species is a consequence of the initial condition considered in this application.

3. Discussions and Conclusions

We have investigated a diffusion process in semi-infinite media in contact with a surface, where the particles may be adsorbed (and desorbed) or sorbed by the surface in order to promote though a reaction process the formation of another species of particle. For particles in the bulk, we have considered generalized diffusion equations, which may be related to a different types of fractional differential operator. In particular, we have considered the Riemann-Liouville and Caputo-Fabrizio fractional time derivative. For each one of these operators, we have found the solutions and obtained the quantity of particles adsorbed by the surface as a function of time. The reaction process considered here occurs only on the surface and it is irreversible, i.e., . The particles obtained from the reaction process are released from the surface to the bulk. In the bulk, we have considered that reaction between the particles is absent; i.e., only the diffusion governs the dynamics of the particles in the bulk. The solutions for these cases have been obtained in terms of the Fox H-functions and generalized Mittag-Leffler functions which are usually related to anomalous diffusion. They have exhibited different diffusive behaviors depending on the conditions employed to analyze the system. In particular, for the situation discussed in Figure 8(a), we have verified that species manifests a superdiffusive behavior for small times when the Riemann-Liouville fractional time operator is considered. For long times, the behavior is essentially governed by the diffusion equations and, consequently, subdiffusive behaviors were obtained. In Figure 8(b), we have analyzed the behavior of the mean square displacement for the Caputo-Fabrizio fractional time operator. For small times, we have a transient until the system reaches a stationary solution, i.e., a time independent solution, for long times. The stationary solution manifested by the system, when the Caputo-Fabrizio is employed, may be related to the resetting process connected to this operator as pointed out in [60]. Finally, we hope that the results presented here may be useful to discussion of situations related to anomalous diffusion.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

E. K. Lenzi acknowledges CNPq for partial financial support.