Super-critical and sub-critical bifurcations in a reaction-diffusion Schnakenberg model with linear cross-diffusion

In this paper the Turing pattern formation mechanism of a two components reaction-diffusion system modeling the Schnakenberg chemical reaction is considered. In Ref. (Madzavamuse et al., J Math Biol 70(4):709–743, 2015) it was shown how the presence of linear cross-diffusion terms favors the destabilization of the constant steady state. We perform the weakly nonlinear multiple scales analysis to derive the equations for the amplitude of the Turing patterns and to show how the cross-diffusion coefficients influence the occurrence of super-critical or sub-critical bifurcations. We present a numerical exploration of far from equilibrium regimes and prove the existence of multistable stationary solutions.