Abstract:We consider an epidemiological SIR model with an infection rate depending on the recovered population. We establish sufficient conditions for existence, uniqueness, and stability (local and global) of endemic equilibria and consider also the stability of the disease-free equilibrium. We show that, in contrast with classical SIR models, a system with a recovery-dependent infection rate can have multiple endemic stable equilibria (multistability) and multiple stable and unstable saddle points of equilibria. We establish conditions for the occurrence of these phenomena and illustrate the results with some examples.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to introduce an infection rate that depends on the recovered population in the epidemiological SIR model and study the new characteristics of this model in terms of equilibrium point stability. Specifically, the paper focuses on the following aspects:
1. **Existence and Uniqueness**: Determine the conditions for the existence of equilibrium points in the model, especially the existence and uniqueness of the endemic equilibrium points. These equilibrium points describe the long - term state of the disease in the population.
2. **Stability Analysis**: Study the local and global stabilities of these equilibrium points. In particular, the paper explores whether this new model can have multiple stable equilibrium points (multistability) and multiple stable and unstable saddle - point equilibrium points compared to the classical SIR model.
3. **Multistability Phenomenon**: Explore whether the system can exhibit the multistability phenomenon, that is, the co - existence of multiple stable equilibrium points, under specific conditions. This is of great significance in epidemic control strategies because different initial conditions may lead to completely different long - term behaviors.
### Main Contributions
- **Theoretical Results**: By defining the auxiliary function \(g(R)\), the paper establishes sufficient conditions for the existence, uniqueness, and stability of equilibrium points. These conditions are expressed based on the relationship between the infection rate function \(f(R)\) and the auxiliary function \(g(R)\).
- **Numerical Examples**: Through specific numerical examples, the paper demonstrates the existence of the multistability phenomenon and illustrates the possible practical impacts of these phenomena.
### Model Background
- **SIR Model**: The classical SIR model divides the population into three groups: susceptibles \(S\), infectives \(I\), and recovered \(R\). The model usually assumes that the infection rate is a constant.
- **Improved Model**: In the model proposed in this paper, the infection rate \(\beta\) is a function \(f(R)\) of the proportion of recovered \(R\), which reflects the impact of human behavior on disease transmission.
### Mathematical Expressions
- **Basic Equations**:
\[
\frac{dS}{d\tau}=1 - f(R)SI - S
\]
\[
\frac{dI}{d\tau}=f(R)SI - kI
\]
\[
\frac{dR}{d\tau}=(k - 1)I - R
\]
where \(f(R)\) is a positive function representing the change of the infection rate with the proportion of recovered.
- **Auxiliary Function**:
\[
g(R)=\frac{k - 1}{k - R}
\]
### Key Conclusions
- **Existence of Equilibrium Points**: When \(f(R)>g(R)\), the model has an endemic equilibrium point.
- **Stability**: If \(\frac{df}{dR}(R^*)<\frac{dg}{dR}(R^*)\), then the endemic equilibrium point \((I^*, R^*)\) is locally stable; conversely, if \(\frac{df}{dR}(R^*)>\frac{dg}{dR}(R^*)\), then \((I^*, R^*)\) is a local saddle point.
- **Multistability**: In some cases, the model can have multiple stable equilibrium points, which depends on the relationship between \(f(R)\) and \(g(R)\).
Through these studies, the paper provides a new perspective for understanding the impact of complex human behavior on disease transmission and theoretical support for formulating more effective epidemic control strategies.