Boundary conditions at a thin membrane for the normal diffusion equation which generate subdiffusion

Subdiffusion equation with fractional Caputo time derivative with respect to another function in modeling transition from ordinary subdiffusion to superdiffusion

We use a subdiffusion equation with fractional Caputo time derivative with respect to another function g (g-subdiffusion equation) to describe a smooth transition from ordinary subdiffusion to superdiffusion. Ordinary subdiffusion is described by the equation with the "ordinary" fractional Caputo time derivative, superdiffusion is described by the equation with a fractional Riesz-type spatial derivative. We find the function g for which the solution (Green's function, GF) to the g-subdiffusion equation takes the form of GF for ordinary subdiffusion in the limit of small time and GF for superdiffusion in the limit of long time. To solve the g-subdiffusion equation we use the g-Laplace transform method. It is shown that the scaling properties of the GF for g-subdiffusion and the GF for superdiffusion are the same in the long time limit. We conclude that for a sufficiently long time the g-subdiffusion equation describes superdiffusion well, despite a different stochastic interpretation of the processes. Then, paradoxically, a subdiffusion equation with a fractional time derivative describes superdiffusion. The superdiffusive effect is achieved here not by making anomalously long jumps by a diffusing particle, but by greatly increasing the particle jump frequency which is derived by means of the g-continuous-time random walk model. The g-subdiffusion equation is shown to be quite general, it can be used in modeling of processes in which a kind of diffusion change continuously over time. In addition, some methods used in modeling of ordinary subdiffusion processes, such as the derivation of local boundary conditions at a thin partially permeable membrane, can be used to model g-subdiffusion processes, even if this process is interpreted as superdiffusion.

Renewal equations for single-particle diffusion through a semipermeable interface

Diffusion through semipermeable interfaces has a wide range of applications, ranging from molecular transport through biological membranes to reverse osmosis for water purification using artificial membranes. At the single-particle level, one-dimensional diffusion through a barrier with constant permeability κ_{0} can be modeled in terms of so-called snapping out Brownian motion (BM). The latter sews together successive rounds of partially reflected BMs that are restricted to either the left or right of the barrier. Each round is killed (absorbed) at the barrier when its Brownian local time exceeds an exponential random variable parameterized by κ_{0}. A new round is then immediately started in either direction with equal probability. It has recently been shown that the probability density for snapping out BM satisfies a renewal equation that relates the full density to the probability densities of partially reflected BM on either side of the barrier. Moreover, generalized versions of the renewal equation can be constructed that incorporate non-Markovian, encounter-based models of absorption. In this paper we extend the renewal theory of snapping out BM to single-particle diffusion in bounded domains and higher spatial dimensions. In each case we show how the solution of the renewal equation satisfies the classical diffusion equation with a permeable boundary condition at the interface. That is, the probability flux across the interface is continuous and proportional to the difference in densities on either side of the interface. We also consider an example of an asymmetric interface in which the directional switching after each absorption event is biased. Finally, we show how to incorporate an encounter-based model of absorption for single-particle diffusion through a spherically symmetric interface. We find that, even when the same non-Markovian model of absorption applies on either side of the interface, the resulting permeability is an asymmetric time-dependent function with memory. Moreover, the permeability functions tend to be heavy tailed.

Subdiffusion equation with Caputo fractional derivative with respect to another function

We show an application of a subdiffusion equation with Caputo fractional time derivative with respect to another function g to describe subdiffusion in a medium having a structure evolving over time. In this case a continuous transition from subdiffusion to other type of diffusion may occur. The process can be interpreted as "ordinary" subdiffusion with fixed subdiffusion parameter (subdiffusion exponent) α in which timescale is changed by the function g. As an example, we consider the transition from "ordinary" subdiffusion to ultraslow diffusion. The g-subdiffusion process generates the additional aging process superimposed on the "standard" aging generated by "ordinary" subdiffusion. The aging process is analyzed using coefficient of relative aging of g-subdiffusion with respect to "ordinary" subdiffusion. The method of solving the g-subdiffusion equation is also presented.