Abstract

Resource support between individuals is of particular importance in controlling or mitigating epidemic spreading, especially during pandemics. However, there remains the question of how we can protect ourselves from being infected while helping others by donating resources in fighting against the epidemic. To answer the question, we propose a novel resource allocation model by considering the awareness of self-protection of individuals. In the model, a tuning parameter is introduced to quantify the reaction strength of individuals when they are aware of the disease. And then, a coupled model of resource allocation and disease spreading is proposed to study the impact of self-awareness on resource allocation and its impact on the dynamics of epidemic spreading. Through theoretical analysis and extensive Monte Carlo simulations, we find that in the stationary state, the system converges to two states: the whole healthy or the completely infected, which indicates an abrupt increase in the prevalence when there is a shortage of resources. More importantly, we find that too cautious and too selfless for the people during the outbreak of an epidemic are both not suitable for disease control. Through extensive simulations, we locate the optimal point, at which there is a maximum value of the epidemic threshold, and an outbreak can be delayed to the greatest extent. At last, we study further the effects of the network structure on the coupled dynamics. We find that the degree heterogeneity promotes the outbreak of disease, and the network structure does not alter the optimal phenomenon in behavior response. Based on the results of this study, a constructive suggestion is that in the face of a global pandemic, individuals or countries should strengthen mutual support and cooperation while doing their own prevention to suppress the epidemic optimally.

1. Introduction

Controlling the outbreak of epidemic is one of the most important topics in human history. During the past decades, the onset of several major global health threats such as the 2003 spread of SARS, the H1N1 influenza pandemic in 2009, and the western Africa Ebola outbreaks in 2014 have deprived tens of thousands of lives all around the world [13]. In 2019, a novel coronavirus causing severe acute respiratory disease (COVID-19) emerged. By the end of April 2020, there have been more than two million confirmed COVID-19 infections reported all over the world [4]. The surge in infections has led to a severe shortage of medical resources. Thousands of confirmed and suspected cases await treatment [5]. Facing the rapid outbreak of the disease, the contribution of resources from healthy individuals is necessary, but the self-protection of the susceptible individuals cannot be ignored. Thus, the immediate problem is how can we protect ourselves from being infected while helping others in fighting against the epidemic.

A large number of researchers from various disciplines have made efforts to study the topic of optimal resource allocation in disease suppressing in the past years [610]. For example, Preciado et al. [11] studied the problem of the optimal distribution of vaccination resources to control epidemic spreading based on complex networks. They found the cost-optimal distribution of the vaccination resource when different levels of vaccination are allowed through a convex framework. Further, they studied the optimal allocation of the preventive and corrective resources to achieve the highest level of containment when the budget is given in advance or finding the minimum budget required when the budget is not specified [12]. Chen et al. [13] solved the problem of optimal allocation of a limited medical resource based on mean-field theory. The above works considered the problem from a mathematical perspective, and were solved from the premise that both the number of resources and the spreading state of the epidemic are fixed.

However, the real scenario is more complicated than a static mathematical problem. Multiple dynamical processes always interact and coevolve [14, 15], forming a more realistic starting point. For example, during the outbreak of an epidemic, the quantity of available resources is largely affected by the propagation of the disease, and in turn, the dynamical change of the necessary resources influences the dynamical property of the epidemic spreading. The coevolution of multiple dynamical processes has attracted extensive research in recent years [16]. Böttcher [17] studied the coevolution of resource and epidemics, and they found a critical recovery cost that if the cost is above the critical value, epidemics spiral out of control into “explosive” spread. Chen et al. [18] studied the effects of social support from local connections on the spreading dynamics of the epidemic. They proposed a coevolution spreading model on multiplex networks and found a hybrid phase transition on networks with heterogeneous degree distribution. In this multiplex network framework, Chen et al. [19] further studied the impact of preferential resource allocation on the social subnetwork on the spreading dynamics of the epidemic.

In addition to the physical resources that can directly mitigate or control the epidemic spreading, the awareness of the epidemic in populations is another type of resource. The public can perceive the threat of the epidemic social network through a platform and mass media and then take measures to protect themselves. Therefore, the interplay between awareness diffusion and epidemic spreading is another topic that have attracted extensive research [2022]. A mass of works addressed the problem from different perspectives considering, for example, the risk perception and behavioral changes [2325]. Granell et al. [14] studied the interplay between the processes of epidemic spreading and awareness diffusion on top of multiplex networks and found a metacritical. Wang et al. [26] studied the coevolution mechanisms using both real online and offline data. Based on the empirical analysis, they proposed a coupled model on a multiplex network to study the coupled dynamics of these two processes.

Although a mass of works on the coevolution of awareness diffusion and disease spreading have been carried out, a question that needs to be addressed is how does awareness affect the behavior of resource donation in populations. To answer this question, a novel resource allocation model that incorporates the influence of the awareness is proposed in this paper. We further consider that the healthy individuals in an outbreak of disease are providers of essential resources, and they can not only produce resources, but also donate them to those who need. The donation behavior in turn leads to less resource for self-protection and a more significant probability to be infected. Thus, when the healthy individuals are aware of the threat of the disease and refuse to donate resources, they would have more resources for self-protection and a lower probability of being infected; however, this is not conducive to the overall prevention and control of the disease. Thus, we are faced with a game-like problem between resource contribution and self-protection [27]. To solve the problem optimally, we should study the complex dynamic mechanisms among self-awareness, resource donation, and epidemic spreading. Thus, a coupled dynamical model on top of the complex networks is proposed, and a dynamic message-passing method is adopted for the theoretical analysis in this paper.

Based on the model, we first investigate the influence of awareness on the coupled dynamics of resource allocation and disease spreading on scale-free networks. To quantify the awareness for self-protection, a parameter is assigned to each individual. Through theoretical analysis and numeric simulations, we find that the system has only two stationary states, namely, the absorb state and the globe outbreak state. With the increase of , the epidemic threshold first increases and then decreases, which indicates that there is an optimal value of . Further, we locate the optimal value through extensive simulations, under which the disease can be suppressed to the greatest extent. Then, we explain qualitatively the optimal phenomenon. And then, we investigate the impact of degree heterogeneity on the coupled dynamics and find that the degree heterogeneity does not alter the optimal phenomenon and the abrupt increase in prevalence with a shortage of resource. At last, we find that the epidemic threshold increases with the decrease in degree heterogeneity, which suggests that network heterogeneity promotes the outbreak of disease.

2. Model Description

2.1. Epidemic Model

A resource-based epidemiological susceptible-infected-susceptible model (r-SIS) is proposed to describe the coupled dynamics of epidemic spreading and resource allocation on a complex network. Individuals are represented by nodes in the network, and an adjacency matrix is introduced to represent the connection between nodes. If there is an edge between nodes and , the element , otherwise . According to this scheme, any individual can be in two different states: susceptible (S) and infected (I). The infection propagates between each pair of I-state and S-state neighbors with an infection rate in one contact, which is assumed to depend on whether the S-state nodes donate resources, see the details in Section 2.2. At any time , each I-state node recovers with a recovery rate . Resources including medical, funds, and food can promote recovery of patients from disease [28, 29]. Thus, we define the recovery rate of each I-state node as a function of the resource quantity received from healthy neighbors in this paper. As each I-state node may get a different amount of resource, the recovery rate varies from node to node. Consequently, the recovery rate of any node at time can be defined aswhere is the resource quantity of node received from healthy neighbors at time and is the basic recovery rate. A parameter is introduced in our model to represent the resource utilization rate [30]. Since in the real scenario, there is the common phenomenon of the waste on resource [31, 32] in medical and other service systems, implying the resource received from healthy neighbors may not be fully utilized on curing and recovery. Without loss of generality, we set throughout this work, i.e., only of the resources received are used.

For the r-SIS model, the value is defined as the probability that any node is in the infected state. The fraction of infected nodes in a network of size at time can be calculated by averaging the overall nodes:

Further, we define the prevalence of the disease in the stationary state as .

2.2. Resource Allocation Model Based on Behavior Response

In the real scenario, healthy individuals can produce resources. For simplicity, we consider that each individual (node) in the network can generate one unit resource at a time step. During an outbreak of a disease, the susceptible individuals can perceive the threat of the disease intuitively by acquiring the information from neighbors. Generally, the more infected neighbors of an individual, the deeper he/she is aware of the disease [33, 34]. People aware of the disease may have different reactions [33]. Also, to quantify the reaction strength of an individual to the local information of disease, a tuning parameter is introduced. Based on the description above, we can define the probability that a healthy individual with infected neighbors donates resource aswhere is a basic donation probability. When , all healthy nodes have the same donation probability . Besides, we consider that a healthy node can donate one unit resource equally to its I-state neighbors at a time. Based on the resource allocation scheme, the amount of resource that node with infected neighbors donate to one of its I-state neighbor, , can be expressed as

When disease breaks out in the human population, people aware of the disease in their proximity can take measures to reduce their susceptibility, leading to a reduction in the effective rate of infection [14, 35]. We consider that if an individual is aware and refuses to donate resource for self-protection, it reduces its infectivity by a factor . The basic infection rate is denoted as , and the annealing infection rate is denoted as , which can thus be expressed as

If a healthy individual donates resource to infected neighbors, the individual has a larger probability to be infected; on the contrary, if the individual does not donate resource, there is a relatively smaller probability to be infected. The annealing infection rate of any node can also be expressed as a function of :

3. Dynamic Message-Passing Method

In order to theoretically analyze the dynamic processes, we develop a generated dynamic message-passing method (GDMP) [36, 37]. In this method, the message is defined on the directed edges of a network to carry causal information of the flow of contagion, which can only transfer one way along directed links. represents the probability that node is infectious because it was infected by one of its neighbors other than node . In computing , we only take into account the contributions to that come from the neighbors other than . The higher order process of being infected by and then passes the infection back to is neglected. Combining and equation (3) for resource allocations, the resources that an infected node receive from its healthy neighbors can be expressed aswhere is the expected number of I-state neighbors of node at time , which is expressed aswhere the plus one takes into account that node is infected at this moment. The factor in equation (7) stands for the probability that node is susceptible at time . With the definition above, the discrete-time version of the evolution of in a time interval reads [38]where is the probability that the node is not infected by any neighbor with the product being over the set of the neighbors of node . The expression of is as follows:

Note that the first term on the right-hand side of equation (9) stands for the probability that node is in S-state and infected by at least one of its neighbors. The second term is the probability that node is in I-state and does not recover. Similarly, we can get the time evolution of aswhere is the probability that node is not infected by any of its neighbors excluding node , which can be expressed as

The product in equation (12) is over the set of the neighbors of excluding . Further, by setting and considering situation in a stationary state, equations (9) and (11) become

Through numerical iteration, we can compute the infection probability of any node at any time and prevalence in the stationary state for different values of and . However, the equations can only be solved numerically, except for the trivial solutions of and , for all , which leads to an overall phase of an all-healthy population.

Due to nonlinearities in equations (7)–(12), they do not have a closed analytic form, and this disallows obtaining the epidemic threshold for fixed values of , such that if and when . The calculation of can be performed by considering that when , and , and the number of infected neighbors of any healthy node is approximately zero in the thermodynamic limit. Then, prior to reaching , the expression is valid. We can get a physical picture that the isolated infected nodes are well separated and surrounded by healthy nodes, and any infected node can receive all the resource from each of its neighbors. By adding these assumptions to equation (7), resource becomes . By linearizing equation (1) and neglecting the second-order terms for small , we obtain

Equation (15) suggests that the recovery rate is proportional to the node degree and inversely proportional to when . For the sake of clarity, the basic recovery rate is set at in this paper. Further equations (10) and (12) can also be linearized using aswhere is the set of directed edges and is the nonbacktracking matrix [39] of the network with the elements labelled by the edges:with being the Dirac delta function. Substituting equation (17) into equation (14) and ignoring the higher order terms of give

Finally, considering that when , equation (6) becomes

To estimate the epidemic threshold, we calculate the average recovery rate as

By inserting equations (20) and (21) into equation (19), we get

The system of equations in equation (22) has a nontrivial solution if and only if is an eigenvalue of the matrix [38]. The lowest value is then given bywhere is the largest eigenvalue of [15, 37, 40].

4. Numerical Verification and Simulation Results

In this section, we study systematically the effects of self-awareness and the network structure on the coupled dynamics of resource allocation and disease spreading, respectively, through numerical verification and Monte Carlo simulations. In the simulation, the synchronous updating method [41, 42] is applied to the disease infection and resource allocation processes. Within each time increment , where in this paper, infection propagates from any I-state node to S-state node with probability and any I-state node recovers to S-state with a probability . With the spreading of disease, the resource allocation process co-occurs. The dynamics terminate once it enters a steady state in which the number of infected nodes only fluctuates within a small range. Note that, we fix the factor at a constant value throughout the paper, such that if any healthy individual chooses to reserve their resource, the probability that they are infected in one contact with an infected neighbor reduces to .

4.1. Effects of Self-Awareness on the Spreading Dynamics

In this section, we investigate the effects of awareness for self-protection on the spreading dynamics. We consider that the coupled processes of resource allocation and disease spreading takes place on a scale-free network, as many real-world networks have skewed degree distributions [4346]. To build the network, we adopt the uncorrelated configuration model (UCM) [47, 48], according to a given degree distribution with maximum degree [49] and minimum degree , which assures no degree correlation of the network when is sufficiently large. To avoid the influence of the network structure on the result, the degree exponent is set at , the network size is set at , and the average degree is set at in the simulations. In addition, we leverage the susceptibility measure to determine the epidemic threshold through simulations [50], which is expressed aswhere represents the ensemble average over all realizations. The epidemic threshold can then be determined when the value of exhibits diverging peaks at the certain infection rate [50, 51].

We first investigate the effects of self-awareness on the spreading dynamics using Monte Carlo simulations. Initially, a fraction of nodes are selected randomly as seeds, and the remaining nodes are in the susceptible state. To present different reaction strength of individuals when they are aware of a certain disease from local information, we select eight typical values of from to in the simulation. In Figures 1(a) and 1(c), we plot the prevalence in the stationary state as a function of basic infection rate for different . Symbols in Figures 1(a) and 1(c) represent the results obtained by Monte Carlo simulations, and lines are the theoretical results obtained from numeric iterations, respectively. From the curves in Figures 1(a) and 1(c), we observe that the system converges to two possible stationary states: either the whole population is healthy, or it becomes completely infected for any , which tells us that when there is a shortage in resource, the disease breaks out abruptly.

Besides, we can observe from Figures 1(a) and 1(b) that with the increase of from to , the epidemic threshold increases gradually, see the peaks of for the corresponding . It reveals that the stronger the individual’s sense of self-protection, the more delayed the outbreak of the disease within this parameter interval (see the right arrow). On the contrary, we observe from Figures 1(c) and 1(d) that when increases from to , the threshold decreases gradually, which reveals that the disease breaks out more easily with a stronger sense of self-protection within this parameter interval (see the left arrow). The phenomenon suggests that too cautious or too selfless for the people during the outbreak of an epidemic are both not suitable for disease control, and there is an optimal value of the reaction strength, at which an epidemic outbreak is postponed to the greatest extent.

We further study systematically the effects of behavior response and the basic infection rate on the spreading dynamics. In Figure 2, we exhibit the full phase diagram of the coupled dynamics of resource allocation and disease spreading. Colors in Figure 2(a) encode the fraction of infected nodes in the stationary state. The epidemic threshold , marked by red circles, rises monotonically until it reaches the maximum at (indicated by the blue dotted line) and then falls gradually with the increase of . Besides, we observe that there are only two possible stationary states: the whole healthy (marked by blue color) and the whole infected of the population (marked by yellow color).

Figure 3(a) plots the time evolution of for six typical values of when the basic infection rate is fixed at . We find that when the value of is small, the system converges to a stationary state rapidly, such as for . With the increase of , it takes a longer time for the system to reach a stationary state. Further, to exhibit the effects of on the dynamics more intuitively, we plot the fraction of infected nodes at a fixed time as a function of in Figure 3(b), which is denoted as for the sake of clarity. We observe that the value of decreases continuously with until reaching the minimum value at (marked by red circle in Figure 3(b)) and then increases gradually with .

Next, we qualitatively explain the optimal phenomena by studying the time evolution of the critical quantities.

We begin by studying the case when is small, for example, . We observe in Figure 4 that in the initial stage, the donation probability for is the highest (see the blue line in the top panel of Figure 4(a)) since a smaller value of means a higher willingness of healthy individuals to allocate resources. Although the resource of healthy individuals can improve the recovery probability of infected neighbors to a certain extent, it also makes themselves more likely to be infected. We can observe in Figures 4(b) and 4(c) that the average recovery rate and the infection rate is the highest for ; meanwhile, there is a lowest value of the effective infection rate , as shown in Figure 4(d). However, with the high probability of being infected for the healthy nodes, the number of infected individuals increases at a high rate (see the blue line in the bottom pane of Figure 4(a)). When people are aware of the increment of the infected neighbors, they reduce their donation willingness, which leads to a reduction in infection rate , as shown in Figures 4(a) and 4(b). Consequently, with less resource received from healthy neighbors, the recovery rate of infected nodes reduces accordingly, see Figure 4(c), which leads to an increase of the effective infection rate [15]. The increase in has led to a further increase in the number of infected nodes. Then, people become more aware of the threat of disease and thus reduce the probability of resource donation further, which leads to a further decrease in the infection rate and the recovery rate and finally the increase of the effective infection rate .

Specifically, we observe from Figure 4(d) that when it surpasses a critical time , indicated by the dotted line in the figure, the value of proliferates, which suggests that in this stage, the infection of healthy individuals is much faster than the recovery of infected individuals. With more newly infected nodes, the donation probability and the infection rate decrease further, which results in less resources donated to support the recovery of infected nodes. Thus, the recovery rate of infected nodes drops abruptly, which in turn promotes the increases of the effective infection rate further, and then more and number of infected nodes appear. Consequently, the cascading failure of the entire system occurs.

Based on the above analysis for a small value of , i.e., , we can reasonably explain why people are more willing to contribute a resource while the disease is more likely to break out.

Secondly, we study the case when is significant, for example, . As a larger value of means more sensitive of the individuals to the disease and a lower willingness to allocate resources. Thus, we observe from Figure 4(a) that, initially, there is a smallest value of (see the green stars in top pane of Figure 4(a)) and the infection rate . Consequently, the infected nodes receive the lowest value of the resource to recover, which leads to the smallest value of the recovery rate , as shown by the green stars in Figure 4(c). Then, the recovery of infected nodes is delayed leading to a high effective infection rate. We can observe in Figure 4(d) that, when , there is a highest value of . The high effective infection rate leads to a rapid increase in the number of infected nodes. We can observe in the bottom pane of Figure 4(a) that, in the early stage, there is a second largest value of for , as denoted by the green stars. The large value of can further reduce the willingness of resource donation for the healthy individuals; thus, we can observe a continuous decline in and . The worse thing is that the recovery rate of infected nodes keeps declining with less and less resource (see the curve in Figure 4(c)), which leads to a rapid growth of (see the curve in Figure 4(d)).

Thus, we can explain the reason why a higher sense of self-protection of the population cannot suppress the disease effectively.

At last, we observe in Figure 4 that, when the value of is around the optimal value , there is a relatively lower value of comparing to the case of in the initial stage, which results in a lower value of (see the yellow squares in Figures 4(a) and 4(b)). The lower willingness of resource donation induces to a relatively smaller value of the recovery rate , as shown in Figure 4(c). However, we can observe from Figure 4(d) that the effective infection rate keeps the lowest value in the early stage, which suggests that the disease propagates slowly in the population, and the number of infected nodes increases slowly, which is verified by the curve in the bottom pane of Figure 4(a). Further, the small value of promotes the increase of (see the curve in the top pane of Figure 4(a)), which results in the increase of the recovery rate . And finally, the effective infection rate decreases further, as shown in Figure 4(d). Thus, the disease can be suppressed to the greatest extend.

Through the three steps, we explain the optimal phenomena in the coupled dynamics of resource allocation and disease spreading.

Finally, we further verify our explanation by studying the critical quantities as the function of parameter at a fixed time and basic infection rate . Figures 5(a) to 5(d) plot the values of , , , , and as a function of when and . For the sake of clarity, we denote the local minimum and maximum value as and and the global minimum and maximum value as and , respectively, where . We observe that although when is around , there is a local maximum of and . The recovery rate reaches maximum , and the effective infection rate reaches the lowest , which indicates that the disease can be optimally suppressed at this point.

4.2. Effects of Network Structure on Spreading Dynamics

In this section, we investigate the effects of the network structure on the coupled dynamics of resource allocation and disease spreading. To avoid the impact of reaction strength on the result, the parameter is fixed at . In addition, we adopt the UCM model to generate scale-free networks with different degree distributions . As the degree heterogeneity decreases with the increase of the power exponent [52, 53], it approaches to random regular networks () when [18].

Figure 6 plots the prevalence in the stationary state as a function of the basic infection rate for networks with four typical values of : (blue circles), (upper triangles), (purple squares), and (red rhombus). We observe that there are only two stationary states of the system: all healthy or completely infected for all networks, which implies that the network structure does not alter the first-order transition of . Besides, we find that, with an increase of , the outbreak of disease is delayed gradually. It suggests that the degree heterogeneity enhances the disease spreading, which is consistent with the existing research conclusions [54].

In the end, we study the effects of behavior response on the spreading dynamics systematically. Figure 7 is the phase diagram in the parameter plane on RRNs. Colors encode the prevalence in the stationary state . We find that there is also an optimal value , at which the epidemic threshold reaches the maximum, indicated by the blue dotted line in Figure 7. The results suggest that the network structure does not alert the optimal phenomenon in behavior response.

5. Discussion

In this paper, we have focused on the problem of how can we protect ourselves from being infected while helping others by donating resources during an outbreak of an epidemic. To answer this question, we have proposed a novel resource allocation model in controlling the epidemic spreading by considering the following two facts, namely, the healthy individuals are the providers of essential resources and there is a kind of game between individual’s self-protection and resource contribution. To quantify the awareness for self-protection, a parameter has been assigned to each individuals in the model. Besides, to study the coupled dynamics of resource allocation and disease spreading, a resource-based SIS model has been proposed. First of all, we have theoretically analyzed the model by using a generated dynamic message-passing method and then carried out extensive Monte Carlo simulations on both scale-free and random regular networks. Through theoretical analysis and simulations, we have found that the coupled dynamics converges to two stationary states: the whole infected or all healthy, which indicates that a shortage of resource can induce an abrupt outbreak of the epidemic. More importantly, we have found that too cautious or too selfless for the people during the outbreak of an epidemic are both not suitable for epidemic containment. There is an optimal (balance) point where the epidemic spreading can be controlled to the greatest extent. It also suggests that one can donate resource appropriately to support the people in need, but at the same time, they should reserve the right amount of resources for self-protection. Further, we have located the optimal point. At last, we have investigated the effects of the network structure on the coupled dynamics and found that the degree heterogeneity promotes the outbreak of disease, and the network structure does not alter the optimal phenomenon in behavior response.

Our research is of practical significance in the context of the global outbreak of COVID-19. It will guide us to make the most reasonable choice between resource contribution and self-protection when perceiving the threat of disease and also have a direct application in the development of strategies to suppress the outbreaks of epidemics. Moreover, our suggestions that in the face of a global pandemic, individuals or countries should strengthen mutual support and cooperation while doing their own prevention are consistent with the current measures taken by most individuals and countries in combating the epidemic. At present, not only the individuals but also the nations are donating resources to support each other while ensuring its own prevention and control needs. For example, when the outbreak in China is effectively contained, it announces assistance to many COVID-19 countries by donating medical resources, such as respirator, mask, nucleic acid testing reagent, and sending the medical staffs [55].

There is still much more work need to be done. For example, the SIS model adopted in this work has its own limitations, and it cannot fully describe the characteristics of most real epidemics. As we all know that there is an incubation period in COVID-19, so an SEIR model may be more suitable. Moreover, in some other epidemics, the recovered individuals would obtain a short acquired immune and then turn into the susceptible state again, so the SIRS model is more suitable. Therefore, the research of the dynamical properties when the present mechanisms are applied in various epidemic models would be the future directions.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (nos. JBK190972, JBK171113, and JBK170505), National Natural Science Foundation of China (nos. 61903266, 71671141, 71873108, and 61703292), the Financial Intelligence & Financial Engineering Key Lab of Sichuan Province, China Postdoctoral Science Foundation (no. 2018M631073), and China Postdoctoral Science Special Foundation (no. 2019T120829).