A spatially adaptive hybrid model in reaction diffusion systems

Many biological reaction-diffusion systems are multiscale: in some regions molecules are abundant, while in others only a few are present. Where numbers are low, intrinsic noise is significant, and a stochastic model such as Gillespie's algorithm is needed to capture the fluctuations and rare events that shape the behaviour. Where numbers are high, this approach is too expensive, and a continuum PDE is sufficient.

Hybrid methods aim to apply each description where it is appropriate, but most require an explicit spatial interface separating the stochastic and deterministic regions. The Spatial Regime Conversion Method (SRCM) avoids this. Each region of space carries both discrete particles and continuous PDE mass, and moves mass between them through conversion events as local concentrations change. The method therefore adapts automatically as the system evolves, resolving stochastic detail wherever intrinsic noise matters and using the cheaper PDE everywhere else, with no fixed interface to track.

In this talk I introduce the method and show how it works, then illustrate it on examples including epidemic spread and a Turing instability driven by noise, where it reproduces the stochastic behaviour that a continuum model alone cannot capture.