Kinetics of Infectious Disease Spread: The SIR Model and Chemical Kinetics Analogy

Kinetics of Infectious Disease Spread

Background: The SIR Model

The SIR model is a foundational epidemiological model used to describe the spread of infectious diseases within a population. It divides the population into three compartments:

• Susceptible (S): Individuals who can acquire the disease if exposed to infectious individuals.

• Infectious (I): Individuals who are currently infected and can transmit the disease to susceptibles.

• Removed (R): Individuals who have recovered (or died) and can no longer be infected or transmit the disease.

At any time, the total population is conserved:

• Population conservation: $S + I + R = N$

The SIR model assumes a closed population (no births or deaths unrelated to the disease).

Mathematical Formulation

• Differential equations:

The time evolution of each compartment is described by the following differential equations:

• $\frac{dS}{dt} = -\beta \frac{S I}{N}$

• $\frac{dI}{dt} = \beta \frac{S I}{N} - \gamma I$

• $\frac{dR}{dt} = \gamma I$

Where:

• $\beta$: Transmission rate constant (per-capita rate at which susceptibles become infected upon contact with infectives).

• $\gamma$: Recovery rate constant (rate at which infectives recover or are removed).

Analogy to Chemical Kinetics

The SIR model can be interpreted using chemical kinetics concepts:

• Transmission event: Analogous to a bimolecular reaction ($S + I \rightarrow 2I$), with rate constant $\beta/N$.

• Recovery event: Analogous to a first-order decay ($I \rightarrow R$), with rate constant $\gamma$.

Key Parameters and Their Roles

• Basic reproduction number ($R_0$): The average number of secondary infections produced by a single infection in a fully susceptible population.

It is defined as:

• $R_0 = \frac{\beta}{\gamma}$

If the average recovery time is $t$ days, then $\gamma = 1/t$ (units: days-1).

Example: SIR Model Parameters

• Given:

  • $N = 1000$ (total population)

  • $S = 1000$ (initial susceptibles)

  • $I = 1$ (initial infectives)

  • $R = 0$ (initial removed)

  • $R_0 = 3.0$ (basic reproduction number)

  • $t = 8$ days (average recovery time)

• Calculate:

  • Recovery rate: $\gamma = 1/t = 0.125$ days-1

  • Transmission rate: $\beta = R_0 \times \gamma = 3.0 \times 0.125 = 0.375$ days-1

  • Transmission rate per encounter: $\beta/N = 0.375/1000 = 0.000375$ days-1

Simulating the SIR Model

To simulate the SIR model numerically (e.g., in Excel):

1. Create a table with columns for time ($t$), $S$, $I$, $R$, and their rates of change.

2. Start with initial conditions (e.g., $S_0 = 999$, $I_0 = 1$, $R_0 = 0$).

3. At each time step, update $S$, $I$, and $R$ using Euler's method:

New value = old value + $\Delta t$ × rate

Choose a small time step (e.g., $\Delta t = 0.1$ day).

Analysis and Extensions

• Peak infection: The day when the number of infected individuals is at its maximum.

• Maximum infections: The highest number of infected individuals at any time.

• Parameter roles:

  • $\beta$ controls how quickly the infection spreads (analogous to a rate constant in chemical kinetics).

  • $\gamma$ controls how quickly individuals recover (analogous to a decay constant).

• Model extension for vaccination: To include vaccination, add a new compartment or modify the equations so that a fraction of susceptibles are moved directly to the removed class (immune).

Modified equations with vaccination rate $v$:

• $\frac{dS}{dt} = -\beta \frac{S I}{N} - v S$

• $\frac{dI}{dt} = \beta \frac{S I}{N} - \gamma I$

• $\frac{dR}{dt} = \gamma I + v S$

Summary Table: SIR Model Parameters and Equations

Parameter/Symbol	Meaning	Typical Units
$S$	Number of susceptible individuals	persons
$I$	Number of infectious individuals	persons
$R$	Number of removed (recovered/immune) individuals	persons
$N$	Total population ($S + I + R$)	persons
$\beta$	Transmission rate constant	days-1
$\gamma$	Recovery rate constant	days-1
$R_0$	Basic reproduction number ($\beta/\gamma$)	dimensionless

Key Takeaways

• The SIR model uses differential equations to describe the kinetics of disease spread, analogous to chemical reaction kinetics.

• Key parameters ($\beta$, $\gamma$, $R_0$) determine the speed and extent of an epidemic.

• Numerical simulation (e.g., in Excel) allows for visualization and analysis of epidemic dynamics.

• Extensions such as vaccination can be incorporated by modifying the model equations.