Advancing Mathematical Epidemiology and Chemical Reaction Network Theory via Synergies Between Them

Abstract

Our paper reviews some key concepts in chemical reaction network theory and mathematical epidemiology, and examines their intersection, with three goals. The first is to make the case that mathematical epidemiology (ME), and also related sciences like population dynamics, virology, ecology, etc., could benefit by adopting the universal language of essentially non-negative kinetic systems as developed by chemical reaction network (CRN) researchers. In this direction, our investigation of the relations between CRN and ME lead us to propose for the first time a definition of ME models, stated in Open Problem 1. Our second goal is to inform researchers outside ME of the convenient next generation matrix (NGM) approach for studying the stability of boundary points, which do not seem sufficiently well known. Last but not least, we want to help students and researchers who know nothing about either ME or CRN to learn them quickly, by offering them a Mathematica package “bootcamp”, including illustrating notebooks (and certain sections below will contain associated suggested notebooks; however, readers with experience may safely skip the bootcamp). We hope that the files indicated in the titles of various sections will be helpful, though of course improvement is always possible, and we ask the help of the readers for that.

Keywords

• mathematical epidemiology
• essentially non-negative ODE systems
• chemical reaction networks
• symbolic computation
• algebraic biology