Monte Carlo simulation and linear stability analysis of Turing pattern formation in reaction-subdiffusion systems

Subdiffusion is an important physical phenomenon observed in many systems. However, numerical techniques to study it, especially when coupled to reactions, are lacking. In this paper, we develop an efficient Monte Carlo algorithm based on the Gillespie algorithm and the continuous-time random walk to simulate reaction-subdiffusion systems. Using this algorithm, we investigate Turing pattern formation in the Schnakenberg model with subdiffusion. First, we show that, as the system becomes more subdiffusive, the homogeneous state becomes more difficult to destablize and Turing patterns form less easily. Second, we show that, as the number of particles in the system decreases, the magnitude of fluctuations increases and again the Turing patterns form less easily. Third, we show that, as the system becomes more subdiffusive, the ratio between the two diffusive constants must be higher in order to observe Turing patterns. Finally, we also carry out linear stability analysis to validate the results obtained from our algorithm.