🧠 Modeling Long-Term Disease Dynamics: The SIR Model with Vital Dynamics

The Susceptible–Infectious–Recovered (SIR) model augmented with Vital Dynamics is a foundational epidemiological framework specifically designed to analyze disease spread over temporal scales sufficiently long that demographic events—namely births and natural deaths—cannot be ignored. This inclusion transforms the analysis from acute outbreak prediction (epidemic) to steady-state prevalence assessment (endemic).

---

🧩 Compartmental Structure and Flow Explanation

The model organizes the population, N, into three core health states—Susceptible (S), Infectious (I), and Recovered (R)—while incorporating mechanisms for population renewal and attrition:

• Susceptible (S): Individuals vulnerable to infection. They gain members via births (Λ or B) and lose members through infection (βSI/N) or natural death (μS).
• Infectious (I): Individuals actively spreading the disease. They gain members from successful infection of susceptibles and lose members through recovery (γI) or natural death (μI).
• Recovered (R): Individuals who have gained permanent immunity. They gain members from recovery and lose members through natural death (μR).

In the simplest endemic setting, the birth rate Λ and natural death rate μ are often assumed to satisfy Λ = μN, producing a constant steady-state population. When pathogen-induced mortality is present, this equality does not hold and N varies dynamically.

To summarize the demographic balance explicitly:

N(t) = S(t) + I(t) + R(t),
dN/dt = Λ − μN (in absence of disease-induced death).

---

📊 Mathematical Formulation (ODEs)

Assuming a mass-action incidence rate and constant per-capita natural birth (Λ) and death (μ) rates, the deterministic system describing changes in S(t), I(t), and R(t) is:

dS/dt = Λ − βSI/N − μS
dI/dt = βSI/N − (γ + μ)I
dR/dt = γI − μR

Where:

• N = S + I + R (constant if Λ = μN).
• β = transmission rate.
• γ = recovery rate.
• Λ = recruitment/birth rate.
• μ = natural death rate.

The total removal rate from the infectious class is therefore γ + μ.

---

⚛️ Parameter Definitions and Ranges

Parameter	Definition	Typical Range (Day⁻¹)	Context
β	Transmission rate coefficient.	0.2 – 1.0	Varies with R₀ and contact behavior.
γ	Recovery rate.	0.14 – 0.5	Infectious period ≈ 2–7 days.
μ	Natural mortality rate.	1/27000 – 1/1825	Long-term demographic processes.
Λ	Recruitment/Birth rate.	≈ μN	Entry into susceptible class.
R₀	Basic Reproduction Number.	1.5 – 4.0	Given by β / (γ + μ).

---

🎯 Applicability and Limitations

Applicability

1. Endemic State Prediction:
   Persistence occurs when R₀ > 1. If R₀ ≤ 1, infections decay to zero.
2. Equilibrium Analysis:
   The endemic equilibrium satisfies:
   I* = (Λ/μ)(1 − 1/R₀) − R*,
   S* = N/R₀,
   with R* determined by demographic balance.
3. Recurrent Outbreaks:
   Continuous replenishment of susceptibles (Λ) sustains repeated epidemic waves.
4. Herd Immunity Threshold:
   H = 1 − 1/R₀,
   representing the critical immune fraction necessary for elimination.

Limitations and Assumptions

1. Homogeneous Mixing: Uniform mixing neglects age structure, spatial structure, and heterogeneity.
2. Permanent Immunity: Immunity is assumed lifelong; waning immunity requires R → S transition.
3. Damped Oscillations: Deterministic dynamics produce decaying oscillations unless stochasticity or seasonal forcing is added.

---

📚 Reference Papers

1. Kermack, W. O., & McKendrick, A. G. A contribution to the mathematical theory of epidemics.
2. Hethcote, H. W. The mathematics of infectious diseases.
3. Brauer, F., Castillo-Chavez, C. Mathematical Models in Population Biology and Epidemiology.
4. Anderson, R. M., & May, R. M. Infectious Diseases of Humans: Dynamics and Control.
5. Diekmann, O., Heesterbeek, J. A. P., & Metz, J. A. J. On the definition and computation of the basic reproduction ratio R₀.