  Spatial reaction-diffusion models have been employed to describe many
emergent phenomena in biological systems. The modelling technique most commonly
adopted in the literature implements systems of partial differential equations
(PDEs), which assumes there are sufficient densities of particles that a
continuum approximation is valid. However, due to recent advances in
computational power, the simulation, and therefore postulation, of
computationally intensive individual-based models has become a popular way to
investigate the effects of noise in reaction-diffusion systems in which regions
of low copy numbers exist.
  The stochastic models with which we shall be concerned in this manuscript are
referred to as `compartment-based'. These models are characterised by a
discretisation of the computational domain into a grid/lattice of
`compartments'. Within each compartment particles are assumed to be well-mixed
and are permitted to react with other particles within their compartment or to
transfer between neighbouring compartments.
  We develop two hybrid algorithms in which a PDE is coupled to a
compartment-based model. Rather than attempting to balance average fluxes, our
algorithms answer a more fundamental question: `how are individual particles
transported between the vastly different model descriptions?' First, we present
an algorithm derived by carefully re-defining the continuous PDE concentration
as a probability distribution. Whilst this first algorithm shows strong
convergence to analytic solutions of test problems, it can be cumbersome to
simulate. Our second algorithm is a simplified and more efficient
implementation of the first, it is derived in the continuum limit over the PDE
region alone. We test our hybrid methods for functionality and accuracy in a
variety of different scenarios by comparing the averaged simulations to
analytic solutions of PDEs for mean concentrations.
