Integrodifferential formulations of the continuous-time random walk for solute transport subject to bimolecular A+B→0 reactions: From micro- to mesoscopic

Subdiffusion-limited fractional reaction-subdiffusion equations with affine reactions: Solution, stochastic paths, and applications

In contrast to normal diffusion, there is no canonical model for reactions between chemical species which move by anomalous subdiffusion. Indeed, the type of mesoscopic equation describing reaction-subdiffusion systems depends on subtle assumptions about the microscopic behavior of individual molecules. Furthermore, the correspondence between mesoscopic and microscopic models is not well understood. In this paper, we study the subdiffusion-limited model, which is defined by mesoscopic equations with fractional derivatives applied to both the movement and the reaction terms. Assuming that the reaction terms are affine functions, we show that the solution to the fractional system is the expectation of a random time change of the solution to the corresponding integer order system. This result yields a simple and explicit algebraic relationship between the fractional and integer order solutions in Laplace space. We then find the microscopic Langevin description of individual molecules that corresponds to such mesoscopic equations and give a computer simulation method to generate their stochastic trajectories. This analysis identifies some precise microscopic conditions that dictate when this type of mesoscopic model is or is not appropriate. We apply our results to several scenarios in cell biology which, despite the ubiquity of subdiffusion in cellular environments, have been modeled almost exclusively by normal diffusion. Specifically, we consider subdiffusive models of morphogen gradient formation, fluctuating mobility, and fluorescence recovery after photobleaching (FRAP) experiments. We also apply our results to fractional ordinary differential equations.

Anomalous reaction-diffusion equations for linear reactions

Deriving evolution equations accounting for both anomalous diffusion and reactions is notoriously difficult, even in the simplest cases. In contrast to normal diffusion, reaction kinetics cannot be incorporated into evolution equations modeling subdiffusion by merely adding reaction terms to the equations describing spatial movement. A series of previous works derived fractional reaction-diffusion equations for the spatiotemporal evolution of particles undergoing subdiffusion in one space dimension with linear reactions between a finite number of discrete states. In this paper, we first give a short and elementary proof of these previous results. We then show how this argument gives the evolution equations for more general cases, including subdiffusion following any fractional Fokker-Planck equation in an arbitrary d-dimensional spatial domain with time-dependent reactions between infinitely many discrete states. In contrast to previous works which employed a variety of technical mathematical methods, our analysis reveals that the evolution equations follow from (1) the probabilistic independence of the stochastic spatial and discrete processes describing a single particle and (2) the linearity of the integro-differential operators describing spatial movement. We also apply our results to systems combining reactions with superdiffusion.

Continuous time random walk with local particle-particle interaction

The continuous time random walk (CTRW) is often applied to the study of particle motion in disordered media. Yet most such applications do not allow for particle-particle (walker-walker) interaction. In this paper, we consider a CTRW with particle-particle interaction; however, for simplicity, we restrain the interaction to be local. The generalized Chapman-Kolmogorov equation is modified by introducing a perturbation function that fluctuates around 1, which models the effect of interaction. Subsequently, a time-fractional nonlinear advection-diffusion equation is derived from this walking system. Under the initial condition of condensed particles at the origin and the free-boundary condition, we numerically solve this equation with both attractive and repulsive particle-particle interactions. Moreover, a Monte Carlo simulation is devised to verify the results of the above numerical work. The equation and the simulation unanimously predict that this walking system converges to the conventional one in the long-time limit. However, for systems where the free-boundary condition and long-time limit are not simultaneously satisfied, this convergence does not hold.

Time-fractional characterization of brine reaction and precipitation in porous media

Brine reaction and precipitation in porous media sometimes occur in the presence of a strong fluid flowing field, which induces the mobilization of the precipitated salts and distorts their spatial distribution. It is interesting to investigate how the distribution responds to such mobilization. We view these precipitates as random walkers in the complex inner space of the porous media, where they make stochastic jumps among locations and possibly wait between successive transitions. In consideration of related experimental results, the waiting time of the precipitates at a particular position is allowed to range widely from short sojourn to permanent residence. Through the model of a continuous-time random walk, a class of time-fractional equations for the precipitate's concentration profile is derived, including that in the Riemann-Liouville formalism and the Prabhakar formalism. The solutions to these equations show the general pattern of the precipitate's spatiotemporal evolution: a coupling of mass accumulation and mass transport. And the degree to which the mass is mobilized turns out to be monotonically correlated to the fractional exponent α. Moreover, to keep the completeness of the model, we further discuss how the interaction among the precipitates influences the precipitation process. In doing so, a time-fractional non-linear Fokker-Planck equation with source term is introduced and solved. It is shown that the interaction among the precipitates slightly perturbs their spatial distribution. This distribution is largely dominated by the brine reaction itself and the interaction between the precipitates and the porous media.

Reaction–diffusion with stochastic decay rates

Microscopic physical and chemical fluctuations in a reaction–diffusion system lead to anomalous chemical kinetics and transport on the mesoscopic scale. Emergent non-Markovian effects lead to power-law reaction times and localization of reacting species.