Kinetic equations for reaction-subdiffusion systems: Derivation and stability analysis

We derive general kinetic equations for reacting and subdiffusing entities based on a nonlinear continuous time random walk formalism proposed by Vlad and Ross [Phys. Rev. E 66, 061908 (2002)]. Reaction and diffusion processes are separable in a typical reaction-diffusion system, and their combined influence on the evolution of the density of a species is a simple sum. Our derivation shows that this is no longer true for subdiffusive entities undergoing reactions. The strong memory effects in the transport process, i.e., the non-Markovian nature of subdiffusion, results in a nontrivial combination of reactions and spatial dispersal, which we discuss in detail. We carry out a linear stability analysis of the derived reaction-subdiffusion system to understand the effects of memory on pattern formation. We find that the Turing instability persists in the subdiffusive system. However, the memory modifies the Turing threshold and the characteristics of the band of unstable modes close to this threshold.