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.

In reaction-transport systems with non-standard diffusion, the memory of the transport causes a coupling of reactions and transport. We investigate the effect of this coupling for systems with Fisher-type kinetics and obtain a general analytical expression for the front speed. We apply our results to the specific case of subdiffusion.

We determine the conditions for the occurrence of Turing instabilities in activator-inhibitor systems, where one component undergoes subdiffusion and the other normal diffusion. If the subdiffusing species has a nonlinear death rate, then coupling between the nonlinear kinetics and the memory effects of the non-Markovian transport process advances the Turing instability if the inhibitor subdiffuses and delays the Turing instability if the activator subdiffuses. We apply the results of our analysis to the Schnakenberg model, the Gray-Scott model, the Oregonator model of the Belousov-Zhabotinsky reaction, and the Lengyel-Epstein model of the chlorine dioxide-iodine-malonic acid reaction.

We study the Ising-Bloch bifurcation in two systems, the Complex Ginzburg Landau equation (CGLE) and a FitzHugh Nagumo (FN) model in the presence of spatial inhomogeneity introduced by Dirichlet boundary conditions. It is seen that the interaction of fronts with boundaries is similar in both systems, establishing the generality of the Ising-Bloch bifurcation. We derive reduced dynamical equations for the FN model that explain front dynamics close to the boundary. We find that front dynamics in a highly non-adiabatic (slow front) limit is controlled by fixed points of the reduced dynamical equations, that occur close to the boundary.

We present a general method of analyzing the influence of finite size and boundary effects on the dynamics of localized solutions of non-linear spatially extended systems. The dynamics of localized structures in infinite systems involve solvability conditions that require projection onto a Goldstone mode. Our method works by extending the solvability conditions to finite sized systems, by incorporating the finite sized modifications of the Goldstone mode and associated nonzero eigenvalue. We apply this method to the special case of non-equilibrium domain walls under the influence of Dirichlet boundary conditions in a parametrically forced complex Ginzburg Landau equation, where we examine exotic nonuniform domain wall motion due to the influence of boundary conditions.