# A General Epidemic Model and Its Application to Mask Design Considering Different Preferences towards Masks

Authors: Chaoqian Wang; Hamdi Kavak
Journal: Complexity
 (2022)

---

## Abstract

While most masks have a limited effect on personal protection, how effective are they for collective protection? How to enlighten the design of masks from the perspective of collective dynamics? In this paper, we assume three preferences in the population: (i) never wearing a mask; (ii) wearing a mask if and only if infected; (iii) always wearing a mask. We study the epidemic transmission in an open system within the Susceptible-Infected-Recovered (SIR) model framework. We use agent-based Monte Carlo simulation and mean-field differential equations to investigate the model, respectively. Ternary heat maps show that wearing masks is always beneficial in curbing the spread of the epidemic. Based on the model, we investigate the potential implications of different mask designs from the perspective of collective dynamics. The results show that strengthening the filterability of the mask from the face to the outside is more effective in most parameter spaces because it acts on individuals with both preferences (ii) and (iii). However, when the fraction of individuals always wearing a mask achieves a critical point, strengthening the filterability from outside to the face becomes more effective because of the emerging hidden reality that the infected individuals become too few to utilize the filterability from their face to outside fully.


---

## Body

## 1. Introduction

As the COVID-19 pandemic is ravaging the world, the protection of masks is a topic of interest. A news article published in Nature indicates that wearing a surgical mask leads to a11% drop in risk, while a 5% drop for cloth [1]. The protective effect of masks on individuals may seem minimal, but it is also necessary to focus on the protective effect on collectives.Since Kermack and McKendrick [2] proposed the Susceptible-Infected-Recovered (SIR) compartment model, various epidemic models have been developed considerably. In the classic SIR model, the population is divided into three compartments: (i) the susceptible S; (ii) the infected I; (iii) the recovered R. Through human-to-human contact or self-healing, individuals flow from one compartment to another. A simple modified version is the SEIR model, which adds an exposed E compartment to the SIR model. Recently, Barlow et al. [3, 4] derived the analytical solutions of the SIR [3] and the SEIR [4] models.Researchers, in recent years, have explored additional factors and mechanisms to the classic epidemic models, such as isolation [5] and vaccination [6–11]. The dynamics of the epidemic transmission can also be applied to the information spreading, creating rumor spreading models [12, 13] or the public opinion dynamics model [14, 15]. From the perspective of verification and validation, the global stability of this class of nonlinear dynamical systems is widely studied [16–21]. In particular, researchers have proved the global stability of endemic equilibria in various epidemic models in multigroup populations [19–21], which are general cases of the model proposed in this work. A common approach to prove global stability is constructing a Lyapunov function (not limited to epidemiology but also widely applied to other complex systems such as evolutionary dynamics [22, 23]), which measures the system’s “energy.” If the energy continues to decay, then the system will stabilize at an equilibrium point.When it comes to the protective effect of masks, several works [24–28] are noticed to have emerged in the COVID-19 period after 2020. Li et al. [24] treated whether people wear masks or not as an evolutionary game. Gondim [25] considered masks in the SEIR model and validated the model by real-world data. Auger and Moussaoui [26] studied the confinement’s release threshold, taking the masks into account. Lasisi and Adeyemo [27] modeled the effect of wearing masks on COVID-19 infection dynamics. Han et al. [28] investigated the effect of three different preferences on wearing a mask.Based on the existing literature, we find the previous works on masks have three shortcomings. First, when classifying the population into three categories with different preferences on wearing masks according to their assumptions, there is no work classifying them into three independent variables. They set only two variables as the fraction of two categories, and the remaining category’s fraction is one minus these two variables. This leads to an inability to ensure constant relative proportions of the other two categories when investigating the effect of the proportion of a certain category. Second, previous work did not carry out a complete analysis of the stability of their models. This makes verification and validation challenging. Third, only focusing on the effect of masks on epidemic spreading, there is no previous work considering providing applications of the epidemic models to the design of masks itself.This paper builds a general epidemic model in an open system considering three different preferences on wearing masks. We start from a set of agent-based rules, and use mean-field analysis to verify and validate the model. In addition to filling in the gaps of previous work by treating three preferences as independent variables and considering global stability analysis, we explore [29, 30] the effect of different preferences towards wearing masks on the epidemic transmission through our model’s eyes. Considering that in the traditional perception, masks are designed at an individual level, we also try to reveal the design strategies of the masks by the collective dynamics based on our model.


## 2. Model

There is an epidemic disease spreading in the system. To prevent this epidemic, individuals hold different preferences for wearing masks. Concerning the infection state, we divide the population into (i) the susceptiblex; (ii) the infected y; (iii) the recovered z. In terms of different preferences towards wearing masks, we divide the population into (i) those who never wear masks (subscript 0); (ii) those who wear masks if and only if infected (subscript 1); (iii) those who always wear masks (subscript 2). Therefore, we have up to 9 categories according to different combinations of the classification of the two dimensions mentioned above.Before describing evolutionary rules, we list the definition of our mathematical symbols in Table1.Table 1
The definition of mathematical symbols.
SymbolDefinitionx0The number of susceptible individuals never wearing a mask.y0The number of infected individuals never wearing a mask.z0The number of recovered individuals never wearing a mask.x1The number of susceptible individuals wearing a mask if and only if infected.y1The number of infected individuals wearing a mask if and only if infected.z1The number of recovered individuals wearing a mask if and only if infected.x2The number of susceptible individuals always wearing a mask.y2The number of infected individuals always wearing a mask.z2The number of recovered individuals always wearing a mask.nThe number of individuals in the system.ΛThe number of new individuals entering the system within unit time.μThe rate of natural death.rThe rate of recovering.αThe rate of human-to-human infection.ε0The fraction of new individuals never wearing a mask.ε1The fraction of new individuals wearing a mask if and only if infected.ε2The fraction of new individuals always wearing a mask.pIThe protective effect produced when an infected individual wears a mask.pSThe protective effect produced when a susceptible individual wears a mask.


### 2.1. The Agent-Based Rules

Consider an open system containing initiallyn|t=0 agents (i.e., individuals). Within a Monte Carlo step, an agent i is randomly selected, and the following parallel events occur.(1)
If agenti is susceptible, we again select an agent j randomly. If agent j is infected, then agent i is infected with a probability αα>0. If agent j wears a mask, then agent i spares from infection with a probability pI0<pI<1. If agent i wears a mask, then agent i spares from infection with a probability pS0<pS<1.(2)
If agenti is infected, then it recovers with a probability r (the average infection cycle is 1/r). This does not happen at the same Monte Carlo step as the event (1).(3)
Agenti naturally dies with a probability μ (the average lifespan is 1/μ). We do not consider deaths due to the epidemic.To ensure the population remains almost unchanged, we must let new agents enter the system. We set the following very first event in a Monte Carlo step, where the number (0) means it happens before the event (1).(0) A new agent enters the system with a probabilitype. The agent’s personal preference determines it never wears a mask with a probability ε00<ε0<1, wears a mask if and only if infected with a probability ε10<ε1<1, or always wears a mask with a probability ε20<ε2<1, yielding ε0+ε1+ε2=1. The preference of an agent on masks does not change over time.For the populationn remaining almost unchanged with time, we let one time step t contains n|t=0 Monte Carlo steps, such that each agent can be selected once on average. Therefore, the expected number of new agents entering the system within one time step t is Λ=n∗pe. The solution is pe=μ (see Theorem 1).


### 2.2. The Mean-Field Equations

Performing mean-field analysis, we can approximate the agent-based dynamics into a set of differential equations. We do not dwell on the principles of mean-field analysis, but only explain some important points. (i) Within unit time, each agent is selected once on average, such that the number of events descending in a category is the population in the category. (ii) The probability of selecting an infected agent never wearing a mask from the population isy0/n, and y1/n, y2/n for selecting an infected agent holding the other two preferences, respectively. (iii) Thanks to the mask, the probability of sparing from infection is pI or pS, which means the probability of infection is 1−pI or 1−pS. (iv) If there are two layers of protection, they must be both breached for the infection to succeed.We denote the state of the system by a vectorψ, containing the population in nine categories. The mean-field differential equations depicting the agent-based dynamics is(1)Ψ⋅=x˙0y˙0z˙0x˙1y˙1z˙1x˙2y˙2z˙2,where(2)x˙0=ε0Λ−αx0y0+1−pIy1+y2n−μx0,y˙0=αx0y0+1−pIy1+y2n−ry0−μy0,z˙0=ry0−μz0,x˙1=ε1Λ−αx1y0+1−pIy1+y2n−μx1,y˙1=αx1y0+1−pIy1+y2n−ry1−μy1,z˙1=ry1−μz1,x˙2=ε2Λ−α1−pSx2y0+1−pIy1+y2n−μx2,y˙2=α1−pSx2y0+1−pIy1+y2n−ry2−μy2,z˙2=ry2−μz2.The system depicted by (1) is an extended version of the SIR model with a standard incidence rate.


## 2.1. The Agent-Based Rules

Consider an open system containing initiallyn|t=0 agents (i.e., individuals). Within a Monte Carlo step, an agent i is randomly selected, and the following parallel events occur.(1)
If agenti is susceptible, we again select an agent j randomly. If agent j is infected, then agent i is infected with a probability αα>0. If agent j wears a mask, then agent i spares from infection with a probability pI0<pI<1. If agent i wears a mask, then agent i spares from infection with a probability pS0<pS<1.(2)
If agenti is infected, then it recovers with a probability r (the average infection cycle is 1/r). This does not happen at the same Monte Carlo step as the event (1).(3)
Agenti naturally dies with a probability μ (the average lifespan is 1/μ). We do not consider deaths due to the epidemic.To ensure the population remains almost unchanged, we must let new agents enter the system. We set the following very first event in a Monte Carlo step, where the number (0) means it happens before the event (1).(0) A new agent enters the system with a probabilitype. The agent’s personal preference determines it never wears a mask with a probability ε00<ε0<1, wears a mask if and only if infected with a probability ε10<ε1<1, or always wears a mask with a probability ε20<ε2<1, yielding ε0+ε1+ε2=1. The preference of an agent on masks does not change over time.For the populationn remaining almost unchanged with time, we let one time step t contains n|t=0 Monte Carlo steps, such that each agent can be selected once on average. Therefore, the expected number of new agents entering the system within one time step t is Λ=n∗pe. The solution is pe=μ (see Theorem 1).


## 2.2. The Mean-Field Equations

Performing mean-field analysis, we can approximate the agent-based dynamics into a set of differential equations. We do not dwell on the principles of mean-field analysis, but only explain some important points. (i) Within unit time, each agent is selected once on average, such that the number of events descending in a category is the population in the category. (ii) The probability of selecting an infected agent never wearing a mask from the population isy0/n, and y1/n, y2/n for selecting an infected agent holding the other two preferences, respectively. (iii) Thanks to the mask, the probability of sparing from infection is pI or pS, which means the probability of infection is 1−pI or 1−pS. (iv) If there are two layers of protection, they must be both breached for the infection to succeed.We denote the state of the system by a vectorψ, containing the population in nine categories. The mean-field differential equations depicting the agent-based dynamics is(1)Ψ⋅=x˙0y˙0z˙0x˙1y˙1z˙1x˙2y˙2z˙2,where(2)x˙0=ε0Λ−αx0y0+1−pIy1+y2n−μx0,y˙0=αx0y0+1−pIy1+y2n−ry0−μy0,z˙0=ry0−μz0,x˙1=ε1Λ−αx1y0+1−pIy1+y2n−μx1,y˙1=αx1y0+1−pIy1+y2n−ry1−μy1,z˙1=ry1−μz1,x˙2=ε2Λ−α1−pSx2y0+1−pIy1+y2n−μx2,y˙2=α1−pSx2y0+1−pIy1+y2n−ry2−μy2,z˙2=ry2−μz2.The system depicted by (1) is an extended version of the SIR model with a standard incidence rate.


## 3. Results and Discussion

In this section, we demonstrate the numerical results of the model from both Monte Carlo simulation and mean-field equations. The algorithm of the Monte Carlo simulation was described in Section2.1. In the numerical simulation of mean-field equations, we use iteration Ψt+Δt=Ψt+Ψ.tΔt, where Δt=0.01.We set three statistical measures: (i) the proportion of susceptible,px=x0+x1+x2/n; (ii) the proportion of infected, py=y0+y1+y2/n; (iii) the proportion of recovered, pz=z0+z1+z2/n.


### 3.1. Impact of Masks on Epidemic Spreading

Figures1 and 2 show the time evolution of px, py, and pz with different parameters and initial conditions.Figure 1
Time evolution ofpx, py, and pz. The results of Monte Carlo simulation and mean-field equations are presented together. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.1, ε1=0.1, ε2=0.8, pI=0.75, pS=0.25. (a) px|t=0=0.9, py|t=0=0.1, pz|t=0=0. (b) px|t=0=0.1, py|t=0=0.9, pz|t=0=0. (a) (b) x0|t=0=x1|t=0=x2|t=0 in x|t=0, as well as y|t=0 and z|t=0.
(a)(b)Figure 2
Time evolution ofpx, py, and pz. The results of Monte Carlo simulation and mean-field equations are presented together. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.3, ε1=0.1, ε2=0.6, pI=0.5, pS=0.05. (a) px|t=0=0.9, py|t=0=0.1, pz|t=0=0. (b) px|t=0=0.1, py|t=0=0.9, pz|t=0=0. (a) (b) x0|t=0=x1|t=0=x2|t=0 in x|t=0, as well as y|t=0 and z|t=0.
(a)(b)From Figures1 and 2, we find the proportions of different individuals always achieve stability after time evolution. The results of the Monte Carlo simulation fluctuate, while the results of mean-field equations are stable. They corroborate each other. In addition, we obverse two phenomena. First, the epidemic may either die out or exist at the end, dependent on different parameters. Second, with the same parameters and different initial conditions, the steady-states are the same.Next, in the heat maps Figures3 and 4, we present the steady-states of py and px as a ternary function of ε0, ε1, ε2, respectively. The results of the Monte Carlo simulation are the average of the last 200 time steps t. The results of mean-field equations are retrieved from the state Ψt+Δt when maxΨtΔt<0.0001.Figure 3
The steady-state ofpy as a ternary function of ε0, ε1, ε2. (a) Monte Carlo simulation. (b) Mean-field equations. α=0.2, μ=0.01, Λ=10, r=0.05, pI=0.5, pS=0.05.
(a)(b)Figure 4
The steady-state ofpx as a ternary function of ε0, ε1, ε2. (a) Monte Carlo simulation. (b) Mean-field equations. α=0.2, μ=0.01, Λ=10, r=0.05, pI=0.5, pS=0.05.
(a)(b)From Figures3 and 4, we find that the results of Monte Carlo simulation and mean-field equations corroborate each other. In Figure 3, we observe that more individuals wearing masks reduce the proportion of infected individuals in the population. In particular, always wearing a mask has a better effect on reducing infected individuals. In Figure 4, we observe that more individuals wearing masks increases the proportion of susceptible individuals in the population. Different from increasing recovered individuals, it means that more people are spared from getting infected once. Also, always wearing a mask has a better effect on increasing susceptible individuals (Figure 4), making more individuals spared from being infected.


### 3.2. Potential Implication of Different Mask Designs

Based on our model, we can reveal the potential implication of different mask designs. When an infected individual wears a mask, it is the filterability from the face to the outside that provides protection (to the population), and when a susceptible individual wears a mask, the filterability from the outside to the face matters.We show in Figures5(a) and 5(b) the steady infected fraction py as a binary function of the protective effect produced when an infected individual wears a mask pI and when a susceptible individual wears a mask pS.Figure 5
(a) Monte Carlo simulation. The steady-state ofpy as a binary function of pI, pS. (b) Mean-field equations. The steady-state of py as a binary function of pI, pS. (c) Mean-field equations. The steady-state of ∂py/∂pI and ∂py/∂pS as a binary function of pI, pS. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.3, ε1=0.1, ε2=0.6, pI=0.5, pS=0.05.
(a)(b)(c)We can observe that in both the Monte Carlo simulation (Figure5(a)) and mean-field equations (Figure 5(b)), an increase in protective effect pI or pS leads to a decrease in the steady infected fraction py. On this basis, we further ask which one in increasing pI or pS is more effective in reducing the infected fraction?We show in Figure5(c) the steady-state of ∂py/∂pI and ∂py/∂pS as a binary function of pI and pS. If ∂py/∂pI<∂py/∂pS, increasing the unit protective effect from the infected side is more conducive to reducing the infected fraction, and vice versa. Intuitively, increasing pI should have always been more conducive than increasing pS because the former acts on individuals with two preferences, (ii) those who wear masks if and only if infected and (iii) those who always wear masks. In contrast, the latter only acts on individuals with one preference, (iii) those who always wear masks. Increasing pI obviously has a broader scope of action than the increasing pS and covers the latter’s population. However, Figure 5(c) presents a different phenomenon. This indicates that we can provide the designs of the masks with different insights from the group dynamics. We will give this further analysis in Section 4.5.


## 3.1. Impact of Masks on Epidemic Spreading

Figures1 and 2 show the time evolution of px, py, and pz with different parameters and initial conditions.Figure 1
Time evolution ofpx, py, and pz. The results of Monte Carlo simulation and mean-field equations are presented together. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.1, ε1=0.1, ε2=0.8, pI=0.75, pS=0.25. (a) px|t=0=0.9, py|t=0=0.1, pz|t=0=0. (b) px|t=0=0.1, py|t=0=0.9, pz|t=0=0. (a) (b) x0|t=0=x1|t=0=x2|t=0 in x|t=0, as well as y|t=0 and z|t=0.
(a)(b)Figure 2
Time evolution ofpx, py, and pz. The results of Monte Carlo simulation and mean-field equations are presented together. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.3, ε1=0.1, ε2=0.6, pI=0.5, pS=0.05. (a) px|t=0=0.9, py|t=0=0.1, pz|t=0=0. (b) px|t=0=0.1, py|t=0=0.9, pz|t=0=0. (a) (b) x0|t=0=x1|t=0=x2|t=0 in x|t=0, as well as y|t=0 and z|t=0.
(a)(b)From Figures1 and 2, we find the proportions of different individuals always achieve stability after time evolution. The results of the Monte Carlo simulation fluctuate, while the results of mean-field equations are stable. They corroborate each other. In addition, we obverse two phenomena. First, the epidemic may either die out or exist at the end, dependent on different parameters. Second, with the same parameters and different initial conditions, the steady-states are the same.Next, in the heat maps Figures3 and 4, we present the steady-states of py and px as a ternary function of ε0, ε1, ε2, respectively. The results of the Monte Carlo simulation are the average of the last 200 time steps t. The results of mean-field equations are retrieved from the state Ψt+Δt when maxΨtΔt<0.0001.Figure 3
The steady-state ofpy as a ternary function of ε0, ε1, ε2. (a) Monte Carlo simulation. (b) Mean-field equations. α=0.2, μ=0.01, Λ=10, r=0.05, pI=0.5, pS=0.05.
(a)(b)Figure 4
The steady-state ofpx as a ternary function of ε0, ε1, ε2. (a) Monte Carlo simulation. (b) Mean-field equations. α=0.2, μ=0.01, Λ=10, r=0.05, pI=0.5, pS=0.05.
(a)(b)From Figures3 and 4, we find that the results of Monte Carlo simulation and mean-field equations corroborate each other. In Figure 3, we observe that more individuals wearing masks reduce the proportion of infected individuals in the population. In particular, always wearing a mask has a better effect on reducing infected individuals. In Figure 4, we observe that more individuals wearing masks increases the proportion of susceptible individuals in the population. Different from increasing recovered individuals, it means that more people are spared from getting infected once. Also, always wearing a mask has a better effect on increasing susceptible individuals (Figure 4), making more individuals spared from being infected.


## 3.2. Potential Implication of Different Mask Designs

Based on our model, we can reveal the potential implication of different mask designs. When an infected individual wears a mask, it is the filterability from the face to the outside that provides protection (to the population), and when a susceptible individual wears a mask, the filterability from the outside to the face matters.We show in Figures5(a) and 5(b) the steady infected fraction py as a binary function of the protective effect produced when an infected individual wears a mask pI and when a susceptible individual wears a mask pS.Figure 5
(a) Monte Carlo simulation. The steady-state ofpy as a binary function of pI, pS. (b) Mean-field equations. The steady-state of py as a binary function of pI, pS. (c) Mean-field equations. The steady-state of ∂py/∂pI and ∂py/∂pS as a binary function of pI, pS. α=0.2, μ=0.01, Λ=10, r=0.05, ε0=0.3, ε1=0.1, ε2=0.6, pI=0.5, pS=0.05.
(a)(b)(c)We can observe that in both the Monte Carlo simulation (Figure5(a)) and mean-field equations (Figure 5(b)), an increase in protective effect pI or pS leads to a decrease in the steady infected fraction py. On this basis, we further ask which one in increasing pI or pS is more effective in reducing the infected fraction?We show in Figure5(c) the steady-state of ∂py/∂pI and ∂py/∂pS as a binary function of pI and pS. If ∂py/∂pI<∂py/∂pS, increasing the unit protective effect from the infected side is more conducive to reducing the infected fraction, and vice versa. Intuitively, increasing pI should have always been more conducive than increasing pS because the former acts on individuals with two preferences, (ii) those who wear masks if and only if infected and (iii) those who always wear masks. In contrast, the latter only acts on individuals with one preference, (iii) those who always wear masks. Increasing pI obviously has a broader scope of action than the increasing pS and covers the latter’s population. However, Figure 5(c) presents a different phenomenon. This indicates that we can provide the designs of the masks with different insights from the group dynamics. We will give this further analysis in Section 4.5.


## 4. Verification and Validation

This section verifies and validates the properties that we concluded in numerical results by analyzing them at a mathematical level.


### 4.1. The Total Population Dynamics

The total population dynamics follows Theorem1.Theorem 1.
Fort⟶∞, we have n⟶Λ/μ.Proof.
According to (1),(3)n˙=x0.+y0.+z0.+x1.+y1.+z1.+x2.+y2.+z2.=Λ−μn.
Solving (3), we get(4)n=n|t=0−Λμe−μt+Λμ.
From (4), we complete the proof that n⟶Λ/μ for t⟶∞. We give this significant value a symbol n∗,(5)n∗=limt⟶∞n=Λμ.
In the same way, we can also prove thatx0+y0+z0⟶ε0Λ/μ, x1+y1+z1⟶ε1Λ/μ, x2+y2+z2⟶ε2Λ/μ for t⟶∞.
Theorem1 gives us another important insight: when discussing the steady-state, we can substitute n for n∗=Λ/μ in the system of (1).


### 4.2. The Basic Reproduction Number

The basic reproduction numberℛ0 is one of the most important measures in epidemiology. It can assist in analyzing both the stability of the system and the effect of parameters on the epidemic transmission.We letΨ.=0 and solve for the epidemic-free equilibrium, denoted by Ψ∗,(6)Ψ∗=Λμε0,0,0,ε1,0,0,ε2,0,0T.Then, we can follow the meWthod in [31] to find the basic reproduction number (see Appendix A):(7)ℛ0=αr+με0+1−pIε1+1−pS1−pIε2.The basic reproduction number reveals the following theorem.Theorem 2.
The epidemic-free equilibriumΨ∗ is locally asymptotically stable if ℛ0<1, and the epidemic-free equilibrium Ψ∗ is not stable if ℛ0>1.
(See [31] for proof).
Substituting the parameters in Figure1 into Equation (A.5), we can calculate ℛ0=0.9167<1, which means the epidemic-free equilibrium is stable, consistent with that shown in Figure 1. Similarly, substituting parameters in Figure 2 produces ℛ0=2.1167>1, such that the epidemic-free equilibrium is not stable, which is also consistent with that shown in Figure 2.


### 4.3. Epidemic-Free and Endemic Equilibria

We separateΨ into(8)Φ1=x0,y0,x1,y1,x2,y2T,Φ2=z0,z1,z2T.From (1), we can assert that the steady-state of Φ1 can determine the steady-state of Φ2, and Φ2 does not affect the evolution of Φ1. Therefore, the stability of Ψ is equivalent to (i) the stability of Φ1 and (ii) the stability of Φ2 when Φ1 achieves stability.We first prove the global stability of the epidemic-free equilibrium by constructing the Lyapunov function.Theorem 3.
The epidemic-free equilibriumΨ∗ is global asymptotically stable in ℝ≥09 if ℛ0<1.Proof.
Consider the Lyapunov functionℒΦ1 in ℝ≥06,(9)ℒΦ1=x0−x0∗22x0∗+y0+1−pIx1−x1∗22x1∗+y1+1−pIx2−x2∗22x2∗+y2.
We can conclude that, (i)ℒΦ1=0 when Φ1=Φ1∗, (ii) ℒΦ1>0 when Φ1≠Φ1∗. Therefore, Φ1 is positive definite in the neighborhood of Φ1∗. Secondly, we have (see Appendix B),(10)ℒ.Φ1≤0,when Φ1≠Φ1∗. Note that ℒ.Φ1=0 can be confirmed by Equation (B.1) when Φ1=Φ1∗. Therefore, Φ1 is negative semidefinite in the neighborhood of Φ1∗. Hence, according to Lasalle’s Invariance Principle [32], Φ1∗ is globally asymptotically stable in ℝ≥06. Given Φ1∗ stable, it is easy to prove the global asymptotic stability of Φ2∗ in ℝ≥03 by constructing Lyapunov function ℒzΦ2=z0+z1+z2. Therefore, Ψ∗ is global asymptotically stable in ℝ≥09.
The equationΨ.=0 has two solutions. We denote the second solution by Ψ∗∗. The infected population is nonzero; thus, we call it the endemic equilibrium. This equilibrium Ψ∗∗ corresponds to the results shown in Figure 2. It is not easy to express it analytically, but we can still show its existence condition.Theorem 4.
The endemic equilibriumΨ∗∗ exists and is unique in ℝ≥09 if ℛ0>1.Proof.
First, we show the relationship between the existence and uniqueness of positiveyi∗∗, i=0,1,2 (y0∗∗>0, y1∗∗>0, y2∗∗>0) and ℛ0>1 (see Appendix C). Then, the existence and uniqueness of xi∗∗, zi∗∗, i=0,1,2 can be naturally confirmed, hence the existence and uniqueness of Ψ∗∗.
Theorem3 validates that in Figure 1, the steady-state with the same parameters is independent of the initial conditions. Theorem 4 is consistent with Figure 2.


### 4.4. Robustness Analysis for the Effect of Wearing Masks

The basic reproduction number measures the average number of individuals that an infected individual can transmit the epidemic. The higher the basic reproduction number, the more severe the epidemic.Analyzing Equation (A.5), we can see that the coefficient beforeε0 is 1, and 1−pI for ε1, and 1−pS1−pI for ε2. Since pS>0, pI>0, we have 1−pS1−pI<1−pI<1. Considering the constraints: 0≤ε0≤1, 0≤ε1≤1, 0≤ε2≤1, ε0+ε1+ε2=1, we know the following facts. (i) ℛ0 takes the minimum when ε0=0, ε1=0, ε2=1. (ii) ℛ0 takes the maximum when ε0=1, ε1=0, ε2=0. Therefore, everyone always wearing a mask minimizes the epidemic severity, while no one wearing masks maximizes the epidemic severity.


### 4.5. Application to Mask Design

The basic reproduction numberℛ0 can play the same role as the infected fraction py in measuring the outbreak severity. From Equation (A.5), we see that the coefficient 1−pI acts on both ε1 and ε2 while 1−pS only acts on ε2, which means the mask protection from the infected side acts on two population categories and that from the susceptible side acts on only one. This brings us a misleading intuition that increasing the protective effect from the infected side is always more conducive. However, we can reveal a hidden different reality by group dynamics.We write the partial derivatives ofℛ0 with respect to pI and pS in (11) and (12).(11)∂ℛ0∂pI=−αr+με1+1−pSε2,(12)∂ℛ0∂pS=−αr+μ1−pIε2.Increasing the protective effect from the infected side is better than that of the susceptible one means(13)∂ℛ0∂pI<∂ℛ0∂pS,or(14)ε1ε2>pS−pI,and increasing the protective effect from the susceptible side is better than that of the infected one means(15)∂ℛ0∂pS<∂ℛ0∂pI,or(16)ε1ε2<pS−pI.We can discuss the parameter space in two cases. First, ifpS<pI, then (13) always holds. Second, if pS>pI, then (15) does not always hold. This suggests that the “intuitive” phenomena (increasing pI is more effective) occupy more parameter space, which can be verified by Figure 5(c).For the latter case,pS>pI, we can transform (16) into ε2>ε1/pS−pI. This indicates that if the fraction of individuals always wearing a mask ε2 exceeds a critical point, ε1/pS−pI, then, increasing pS is more effective than increasing pI, even if pI can act on both category ε1 and ε2. We can interpret this result in daily language by taking into account the fraction of existing infected individuals. As we showed in Section 4.4, ℛ0 decreases (i.e., infected individuals increasing) with an increase in ε2. In this way, the critical point of ε2 is rational to exist, over which the infected individuals are too few to exert the protective effect that the mask produces on their side.


## 4.1. The Total Population Dynamics

The total population dynamics follows Theorem1.Theorem 1.
Fort⟶∞, we have n⟶Λ/μ.Proof.
According to (1),(3)n˙=x0.+y0.+z0.+x1.+y1.+z1.+x2.+y2.+z2.=Λ−μn.
Solving (3), we get(4)n=n|t=0−Λμe−μt+Λμ.
From (4), we complete the proof that n⟶Λ/μ for t⟶∞. We give this significant value a symbol n∗,(5)n∗=limt⟶∞n=Λμ.
In the same way, we can also prove thatx0+y0+z0⟶ε0Λ/μ, x1+y1+z1⟶ε1Λ/μ, x2+y2+z2⟶ε2Λ/μ for t⟶∞.
Theorem1 gives us another important insight: when discussing the steady-state, we can substitute n for n∗=Λ/μ in the system of (1).


## 4.2. The Basic Reproduction Number

The basic reproduction numberℛ0 is one of the most important measures in epidemiology. It can assist in analyzing both the stability of the system and the effect of parameters on the epidemic transmission.We letΨ.=0 and solve for the epidemic-free equilibrium, denoted by Ψ∗,(6)Ψ∗=Λμε0,0,0,ε1,0,0,ε2,0,0T.Then, we can follow the meWthod in [31] to find the basic reproduction number (see Appendix A):(7)ℛ0=αr+με0+1−pIε1+1−pS1−pIε2.The basic reproduction number reveals the following theorem.Theorem 2.
The epidemic-free equilibriumΨ∗ is locally asymptotically stable if ℛ0<1, and the epidemic-free equilibrium Ψ∗ is not stable if ℛ0>1.
(See [31] for proof).
Substituting the parameters in Figure1 into Equation (A.5), we can calculate ℛ0=0.9167<1, which means the epidemic-free equilibrium is stable, consistent with that shown in Figure 1. Similarly, substituting parameters in Figure 2 produces ℛ0=2.1167>1, such that the epidemic-free equilibrium is not stable, which is also consistent with that shown in Figure 2.


## 4.3. Epidemic-Free and Endemic Equilibria

We separateΨ into(8)Φ1=x0,y0,x1,y1,x2,y2T,Φ2=z0,z1,z2T.From (1), we can assert that the steady-state of Φ1 can determine the steady-state of Φ2, and Φ2 does not affect the evolution of Φ1. Therefore, the stability of Ψ is equivalent to (i) the stability of Φ1 and (ii) the stability of Φ2 when Φ1 achieves stability.We first prove the global stability of the epidemic-free equilibrium by constructing the Lyapunov function.Theorem 3.
The epidemic-free equilibriumΨ∗ is global asymptotically stable in ℝ≥09 if ℛ0<1.Proof.
Consider the Lyapunov functionℒΦ1 in ℝ≥06,(9)ℒΦ1=x0−x0∗22x0∗+y0+1−pIx1−x1∗22x1∗+y1+1−pIx2−x2∗22x2∗+y2.
We can conclude that, (i)ℒΦ1=0 when Φ1=Φ1∗, (ii) ℒΦ1>0 when Φ1≠Φ1∗. Therefore, Φ1 is positive definite in the neighborhood of Φ1∗. Secondly, we have (see Appendix B),(10)ℒ.Φ1≤0,when Φ1≠Φ1∗. Note that ℒ.Φ1=0 can be confirmed by Equation (B.1) when Φ1=Φ1∗. Therefore, Φ1 is negative semidefinite in the neighborhood of Φ1∗. Hence, according to Lasalle’s Invariance Principle [32], Φ1∗ is globally asymptotically stable in ℝ≥06. Given Φ1∗ stable, it is easy to prove the global asymptotic stability of Φ2∗ in ℝ≥03 by constructing Lyapunov function ℒzΦ2=z0+z1+z2. Therefore, Ψ∗ is global asymptotically stable in ℝ≥09.
The equationΨ.=0 has two solutions. We denote the second solution by Ψ∗∗. The infected population is nonzero; thus, we call it the endemic equilibrium. This equilibrium Ψ∗∗ corresponds to the results shown in Figure 2. It is not easy to express it analytically, but we can still show its existence condition.Theorem 4.
The endemic equilibriumΨ∗∗ exists and is unique in ℝ≥09 if ℛ0>1.Proof.
First, we show the relationship between the existence and uniqueness of positiveyi∗∗, i=0,1,2 (y0∗∗>0, y1∗∗>0, y2∗∗>0) and ℛ0>1 (see Appendix C). Then, the existence and uniqueness of xi∗∗, zi∗∗, i=0,1,2 can be naturally confirmed, hence the existence and uniqueness of Ψ∗∗.
Theorem3 validates that in Figure 1, the steady-state with the same parameters is independent of the initial conditions. Theorem 4 is consistent with Figure 2.


## 4.4. Robustness Analysis for the Effect of Wearing Masks

The basic reproduction number measures the average number of individuals that an infected individual can transmit the epidemic. The higher the basic reproduction number, the more severe the epidemic.Analyzing Equation (A.5), we can see that the coefficient beforeε0 is 1, and 1−pI for ε1, and 1−pS1−pI for ε2. Since pS>0, pI>0, we have 1−pS1−pI<1−pI<1. Considering the constraints: 0≤ε0≤1, 0≤ε1≤1, 0≤ε2≤1, ε0+ε1+ε2=1, we know the following facts. (i) ℛ0 takes the minimum when ε0=0, ε1=0, ε2=1. (ii) ℛ0 takes the maximum when ε0=1, ε1=0, ε2=0. Therefore, everyone always wearing a mask minimizes the epidemic severity, while no one wearing masks maximizes the epidemic severity.


## 4.5. Application to Mask Design

The basic reproduction numberℛ0 can play the same role as the infected fraction py in measuring the outbreak severity. From Equation (A.5), we see that the coefficient 1−pI acts on both ε1 and ε2 while 1−pS only acts on ε2, which means the mask protection from the infected side acts on two population categories and that from the susceptible side acts on only one. This brings us a misleading intuition that increasing the protective effect from the infected side is always more conducive. However, we can reveal a hidden different reality by group dynamics.We write the partial derivatives ofℛ0 with respect to pI and pS in (11) and (12).(11)∂ℛ0∂pI=−αr+με1+1−pSε2,(12)∂ℛ0∂pS=−αr+μ1−pIε2.Increasing the protective effect from the infected side is better than that of the susceptible one means(13)∂ℛ0∂pI<∂ℛ0∂pS,or(14)ε1ε2>pS−pI,and increasing the protective effect from the susceptible side is better than that of the infected one means(15)∂ℛ0∂pS<∂ℛ0∂pI,or(16)ε1ε2<pS−pI.We can discuss the parameter space in two cases. First, ifpS<pI, then (13) always holds. Second, if pS>pI, then (15) does not always hold. This suggests that the “intuitive” phenomena (increasing pI is more effective) occupy more parameter space, which can be verified by Figure 5(c).For the latter case,pS>pI, we can transform (16) into ε2>ε1/pS−pI. This indicates that if the fraction of individuals always wearing a mask ε2 exceeds a critical point, ε1/pS−pI, then, increasing pS is more effective than increasing pI, even if pI can act on both category ε1 and ε2. We can interpret this result in daily language by taking into account the fraction of existing infected individuals. As we showed in Section 4.4, ℛ0 decreases (i.e., infected individuals increasing) with an increase in ε2. In this way, the critical point of ε2 is rational to exist, over which the infected individuals are too few to exert the protective effect that the mask produces on their side.


## 5. Conclusion

Although most masks have little to no effect on personal protection [1], we are still interested in the protective effects of masks on a population. We proposed a general epidemic model in the classic SIR framework considering three different preferences towards wearing masks. Some individuals never wear masks; others wear masks if and only if infected, and some always wear masks. We started from agent-based rules and used a set of mean-field differential equations to approximate the model. The results of the two corroborate each other. In this work, the three preferences are independent of each other.The first aspect is the effect of masks on epidemics. The ternary heat maps revealed that wearing masks can reduce the number of infected individuals and increase the number of susceptible individuals. We provided the global stability analysis of the results and showed the robustness of the effectiveness of masks by analyzing the basic reproduction number of the epidemic. We concluded that wearing masks are beneficial to the control of epidemics.The second aspect is the application of the epidemic model to mask design. The protective effect from the infected sidepI can be understood as the filterability of the mask from the face to the outside against viruses, while the protective effect to the susceptible side pS can be interpreted as the filterability from the outside to the face. This can be influenced by the material and design of the mask [33], and we analyzed which side strengthening would provide better results. We showed that strengthening the infected side is more effective in most parameter spaces. This is intuitive since strengthening the infected side acts on two categories of individuals (those wearing masks only if infected and those always wearing masks), while strengthening the susceptible side acts on only one category (those always wearing masks). However, there is a hidden reality from the perspective of group dynamics. We found that once the fraction of individuals always wearing masks exceeds a critical point, ε2>ε1/pS−pI, then, strengthening the susceptible side becomes more effective. This is because the preference of always wearing masks reduces the infected fraction in the population, so that the infected individuals are too few to exert the protective effect of masks produced on their side. In the daily language, both the cases above seem to make sense. However, noticing the latter case from the group perspective and further giving the mask design strategies according to parameter spaces are not straightforward without the help of system dynamics.Real-world situations may have more complexity and different insights. For instance, the underlying assumptions—people’s preferences do not change with time, ignores human subjectivity, which has the potential to reveal more insights. In fact, people can change their preference on whether to wear masks by either estimating the epidemic severity (evolutionary games) or being affected by the propaganda of the effectiveness of masks (opinion dynamics). In this way, future work may consider time-dependent preferences, and apply any modified model to mask design.



---

*Source: 1626008-2022-10-30.xml*