Abstract: We proposed an improved next subvolume method in this paper. Numerical results show that the proposed method can be applied to biochemical reactions systems with complex spatial structure, indicating significant improvement in efficiency compared to existing approaches under some conditions.

Key words: gene regulatory network, reaction-diffusion master equation, Gillespie algorithm, the next subvolume method