“…For the time stepping, we use an implicit-explicit scheme (IMEX) [38][39][40], where we implement the diffusion term using the implicit Crank-Nicolson method, and use the explicit Euler method for the nonlinear reaction term, with ∆t = 0.01. The Crank-Nicolson scheme has been shown to be stable even in the presence of noise [41,42]. We have verified that all described behaviour still occurs for a much finer time-step with ∆t = 0.00125 or a finer grid size with ∆x = 0.25, and that the numerical scheme indeed converges as ∆t is decreased.…”