Global stability of epidemic models with uniform susceptibility

Transmission dynamics of infectious diseases are often studied using compartmental mathematical models, which are commonly represented as systems of autonomous ordinary differential equations. A key step in the analysis of such models is to identify equilibria and find conditions for their stability. Local stability analysis reduces to a problem in linear algebra, but there is no general algorithm for establishing global stability properties. Substantial progress on global stability of epidemic models has been made in the last 20 years, primarily by successfully applying Lyapunov’s method to specific systems. Here, we show that any compartmental epidemic model in which susceptible individuals cannot be distinguished and can be infected only once, has a globally asymptotically stable (GAS) equilibrium. If the basic reproduction number ${R}_0$ satisfies ${R}_0 > 1$, then the GAS fixed point is an endemic equilibrium (i.e., constant, positive disease prevalence). Alternatively, if ${R}_0 \le 1$, then the GAS equilibrium is disease-free. This theorem subsumes a large number of results published over the last century, strengthens most of them by establishing global rather than local stability, avoids the need for any stability analyses of these systems in the future, and settles the question of whether co-existing stable solutions or non-equilibrium attractors are possible in such models: they are not.

💡 Research Summary

This groundbreaking paper provides a definitive answer to a long-standing question in mathematical epidemiology: what is the long-term fate of disease trajectories in a broad class of compartmental models? The authors prove that any epidemic model characterized by “uniform susceptibility” – meaning there is a single, homogeneous susceptible class (S) and infected individuals never return to it – possesses a Globally Asymptotically Stable (GAS) equilibrium. The dynamics are entirely dictated by the basic reproduction number R0: if R0 ≤ 1, all solutions converge to the disease-free equilibrium; if R0 > 1, all solutions (except the trivial case with no initial infection) converge to a unique endemic equilibrium.

The research formalizes a general S x1…xn model class, where S is the susceptible compartment and x1,…,xn are any number of non-susceptible compartments (e.g., exposed, infectious, recovered, hospitalized). The model allows for arbitrary flows among these non-susceptible compartments, described by a Post-Infection Transfer (PIT) matrix M. A key insight is that the structural assumptions (no flow back to S, positive mortality) guarantee M is a non-singular M-matrix, a property central to the analysis.

The major technical achievement is the construction of a Lyapunov function that works for this entire model class, not just specific instances. The authors demonstrate that building such a function reduces to solving a system of linear equations and verifying a set of algebraic inequalities, which are always satisfied due to the M-matrix properties of the system. This elegantly bypasses the need for ad-hoc, model-specific stability proofs that have populated the literature for decades.

The implications are profound. The theorem unifies and strengthens stability results for thousands of existing models (SIR, SEIR, multi-stage models, models with arbitrarily distributed stage durations) by establishing global, not just local, stability. It conclusively rules out the possibility of complex long-term dynamics like sustained oscillations, chaos, or multiple stable equilibria within this vast class. The discussion carefully delineates the boundaries of the result: extensions like maternally-derived immunity or vaccination compartments (with no flow to infected classes) preserve global stability, but models where infection-derived immunity wanes (allowing flow back to S) fall outside the theorem’s scope and can indeed exhibit periodic solutions.

In essence, this work provides a powerful “global stability guarantee” for a core family of epidemic models, significantly simplifying future analyses and solidifying our theoretical understanding of their predictable asymptotic behavior.