An effective drift-diffusion model for pandemic propagation and uncertainty prediction

Predicting pandemic evolution involves complex modeling challenges, typically involving detailed discrete mathematics executed on large volumes of epidemiological data. Making them physics based provides added intuition as well as predictive value. Differential equations have the advantage of offering smooth, well-behaved solutions that try to capture overall predictive trends and averages. In this paper, the canonical susceptible-infected-recovered model is simplified, in the process generating quasi-analytical solutions and fitting functions that agree well with the numerics, as well as infection data across multiple countries. The equations provide an elegant way to visualize the evolution of the pandemic spread, by drawing equivalents with the similar dynamics of a particle, whose location over time represents the growing fraction of the population that is infected. This particle slides down a potential whose shape is set by model epidemiological parameters such as reproduction rate. Potential sources of errors and their growth over time are identified, and the uncertainties are mapped into a diffusive jitter that tends to push the particle away from its minimum. The combined physical understanding and analytical expressions offered by such an intuitive drift-diffusion model sets the foundation for their eventual extension to a multi-patch model while offering practical error bounds and could thus be useful in making policy decisions going forward.