Convergence of a finite volume scheme for a model for ants

We develop and analyse a finite volume scheme for a nonlocal active matter system known to exhibit a rich array of complex behaviours. The model under investigation was derived from a stochastic system of interacting particles describing a foraging ant colony coupled to pheromone dynamics. In this work, we prove that the unique numerical solution converges to the unique weak solution as the mesh size and the time step go to zero. We also show discrete long-time estimates, which prove that certain norms are preserved for all times, uniformly in the mesh size and time step. In particular, we prove higher regularity estimates which provide an analogue of continuum parabolic higher regularity estimates. Finally, we numerically study the rate of convergence of the scheme, and we provide examples of the existence of multiple metastable steady states.