## Abstract

With more than 60 million confirmed cases and more than 1.5 million deaths, SARS-CoV-2 has paralyzed our societies, leading to self isolation and quarantine for several months. A COVID-19 vaccine remains a critical element in the eventual solution to this public health crisis. Two vaccines are ready to be mass produced and eventually supplied to the population. Here, we develop an epidemiological network model able to represent COVID-19 pandemic dynamics of different countries such as in Italy. Stochastic computational simulations identify the necessary number of vaccines and vaccine efficacy thresholds capable of preventing an epidemic whilst adhering to lockdown guidelines. Assuming a vaccine efficacy of 100% in a mass vaccination program, at least 60% of a given population should be vaccinated to obtain herd immunity. Otherwise, new waves of COVID-19 would be reported. However, reaching 60% of the population will be an unprecedented mass immunisation campaign that will take several months and millions of dollars. If the vaccine efficacy reports lower levels of efficacy in practice, then the coverage of vaccination would be needed to be even higher. Simulations suggest that the “Ring of Vaccination” strategy, vaccinating susceptible contact and contact of contacts, would be a more tractable vaccine strategy to prevent the new waves of COVID-19 meanwhile a high percent of the population is vaccinated.

## 1 Introduction

The novel Coronavirus SARS-CoV-2 epidemic which first emerged from Wuhan, China in December 2019 has spread globally, causing high levels of mortality and morbidity worldwide. To curb the spread, governments across the world have implemented measures ranging from quarantining, social distancing, wearing of face masks among others. Amidst this crisis, national health care systems such as in Italy and the United States of America have been overwhelmed by the ever-increasing number of infection cases [36].

COVID-19 epidemiological models have been formulated to understand and curb the spread of the disease. Many of these models follow an SIR framework either in the deterministic or stochastic form or both [2, 59, 53, 55, 9, 10]. Other variations and modifications to this general model have been considered including SEIR [47, 30, 50] and SIRD [13] compartmental models. Some models also include parameters such as age-heterogeneity [56], guiding the flow of users in supermarkets [19], and governmental policies [50, 4, 32]. In addition, a few studies incorporate the dynamics of the disease within an individual host [21, 1]. However, only a few of these models, currently predicting the SARS-CoV-2 pandemic, consider the structure of the population and the underlying interactions between individuals [15, 61, 22].

The assumption of random homogeneous mixing in epidemiological models has been doc-umented to be unrealistic in nature as populations have underlying structural properties and individuals tend to interact with each other [14]. Increasingly, network theory is being used in epidemiology [35, 26, 51, 16]. In particular, social networks have gained popularity in conceptualising the effects of social interaction during epidemics in a given population [38, 25, 54]. Contacts between individuals in a given population can be captured in a network where nodes represent individuals and the edges represent the connections between them [20]. Social networks are thus important determinants of infectious disease transmission as for example, infections transmitted by close contact can easily spread along the paths of a network.

With the continuation of SARS-CoV-2 in many parts of the world, the push for a vaccine has become highly necessary. Pharmaceutical companies are in a race to develop suitable vaccines as there is a lack of other alternatives. As of October 2020, there were 17 candidate vaccines undergoing trial at various stages. Owing to the fact that it is a novel viral disease, it is still unclear what levels of vaccine efficacies will be sufficient to curb the spread of the virus. Identifying such efficacies earlier can direct vaccine development and administration in the population [23]. Previous studies suggest that vaccination would be effective for protecting the host against SARS-CoV-2 [28], however, there have not been studies to evaluate the potential effects of different vaccination programs over a network model.

A key question to be answered is, how much vaccine is required to create herd immunity to block SARS-CoV-2 transmission? [3]. In other words, how many people need to be vaccinated in order to reach herd immunity? Here, we employ a network-based approach to explore the potentials of two vaccination schemes, classical mass vaccination and ring vaccination, in minimizing the spread of SARS-CoV-2, see Figure 1. Ring vaccination is a vaccination strategy in which infected cases and contacts of cases are identified and vaccinated [18, 39]. Our analysis uses a stochastic network simulation model of SARS-CoV-2 transmission to examine its control by mass and ring vaccination strategies with varying efficacies in the presence of non-pharmaceutical disease control interventions. Given that vaccinating millions of people will require a lot of time, this study implements a lockdown period as a further control measure during vaccination strategies in our simulations. The main objective is to identify the necessary vaccine efficacy thresholds capable of preventing an epidemic whilst adhering to lockdown instructions.

## 2 Modelling Epidemic Spread on Networks

### 2.1 Model Setup

#### Network Generation

To study the impact of vaccination on the epidemic dynamics of SARS-CoV-2, the underlying structure of human interactions is represented by a typical complex network. As heterogeneous networks are often used to explore epidemic spread, we consider the infection dynamics on an Erdos-Renyi (ER) network due to its tractability and practicability [27, Chapter 1]. The ER network is based on the *G*(*n, M*) random graph model characterized by two parameters; the network size *n* and the number of edges *M* which assigns exactly *M* number of edges to each graph. In the ER networks used for this study, *n* = *N, M* = 5*N*, that is *G*(*N*, 5*N*), which yields an average degree distribution ⟨ *k*⟩ = 10.

#### Epidemic Spread

SARS-CoV-2 is a disease which spreads primarily through close contact with an infected person. Following the exposure to the SARS-CoV-2 virus and before symptom onset, individuals go through an incubation period of about 14 days with the average being about 5 − 6 days [33, 52, 7]. After this incubation period, infected individuals become symptomatic and are able to transmit to others through respiratory droplets or by direct contact [33].

#### Epidemic States

We simulated a Susceptible - Exposed - Infected - Recovered - Vaccinated (SEIRV) process on an Erdos-Renyi network with *N* = 10^{6} nodes. At any instant during the infection process, the status of a node can be in any one of five possible states: Susceptible (*S*, not infected but can be infected), Exposed (*E*, infectious, may not show symptoms and can transmit to others), Infected (*I*, infectious and symptomatic, capable of transmission to others), Recovered (*R*, recovered and immune to the disease) and Vaccinated (*V*, vaccinated). Therefore, at any time, *t*, in the infection process, *N* (*t*) = *S*(*t*) + *E*(*t*) + *I*(*t*) + *R*(*t*) + *V* (*t*).

#### Epidemic Process

The model (depicted in Algorithm 1 proceeds in discrete one-day time steps for a period of 360 days to determine disease dynamics. Simulation codes can be found in Supplemental Materials. At the initial state of the epidemic process, all individuals in the network are susceptible except one (patient zero) which is in the exposed state. On each day during the epidemic process, there is interaction between individuals and infected persons can potentially transmit to their susceptible contacts. If a susceptible individual comes into contact with someone who has the virus (that is, a person in the *E* or *I* state), a Binomial trial is used to determine if the contact results in an infection. If yes, the newly infected susceptible individual moves from the *S* state to the *E* state. Any individual exposed to the disease remains in the *E* state for the duration of the incubation period after which it moves to the *I* state. Infectious individuals can transmit to their neighbours, when they come into contact, in Binomial trial with a given probability, *β*(*t*) and then move into the *R* state after the infectious period.

In the course of an epidemic, the rate of infection is never constant. As interventions are being executed, the rate of infection decreases, otherwise it increases. To model this response, we utilise a double logistic function to model the various phases of the infection dynamics leading to a decline in cases when interventions are initiated and a resurgence in cases when interventions cease. The double logistic function used to define a time dependent probability of infection, *β*(*t*) is as follows:
where *b*_{1} is the first boundary (i.e. function value at time zero), *b*_{2} is the second boundary, *r*_{1} is the rate of change of first period, *r*_{2} is rate of change of second period, *m*_{1} is the midpoint of the first period (start of interventions), *m*_{2} is the midpoint of the second period (end of interventions) and *t* is time. Note that if *b*_{1} *> b*_{2} the function increases first and then decreases, and vice versa. To get a good set of parameters for *β*(*t*), we applied Eqn (1) to emulate infection data (from February 22 to September 1 2020) from Italy, one of the worst hit countries during the SARS-CoV-2 pandemic.

The main goal of vaccination in this model is to block transmission. Every day, there is contact between infected and vaccinated individuals. Vaccinated individuals that become exposed due to contact with an exposed individual move into the *E* state with probability (1 − *η*)*β*(*t*), where *η* is the efficacy of the vaccine. Individuals for whom the vaccine is effective remain in the *V* state whilst those for whom the vaccine is not effective move into the *E* state and later into the *I* state after which they proceed to the *R* state. Furthermore, depending on the vaccination strategy being modelled, the time of vaccination as well as population coverage varies.

For the first strategy which we refer to as the **classical mass vaccination** strategy, a fraction of individuals start in the *V* state (denoted %*vac*) with one index case (patient zero) and the rest in the *S* state. The population coverage in this case ranges from 10%, 20%, …, 100%. The modelling process follows as in Algorithm 1.

For the second strategy which we refer to as **ring vaccination**, vaccination occurs after a percentage (1% or 3%) of the population has been exposed to the virus during the epidemic. In this case, we simulate the epidemics as described above (in Algorithm 1) with no vaccination and only begin vaccine administration after a proportion of the population (%*exposed*) has been exposed (See Algorithm 2 for more details). We assume that once a case is diagnosed, all the contacts or neighbours are traceable and vaccinated. To initiate the ring vaccination, susceptible contacts of exposed individuals can be vaccinated in a binomial trial with probability (1 − *η*)*β*(*t*). This vaccination is further extended to susceptible contacts of these first contacts with the assumption that contacts are identified through contact tracing. Contacts of contacts are also vaccinated with the same probability. In our model, we also assume that traced and identified susceptible contacts and contacts of contacts are vaccinated. This process is described in Algorithm 2 and detailed implementation can be found at Code and examples. In both strategies, we assume that the vaccine does not have an effect on infected or exposed individuals. Table 1 summarizes all parameters and key terms used in this study.

### 2.2 Simulation Scenarios

Several scenarios are explored here. First, we consider an epidemic implemented with no vaccination. This scenario is evident of the epidemic outcome without the discovery of a vaccine. Subsequent scenarios consist of the overall performance of the vaccine by considering its efficacy in preventing transmission and symptomatic infection using either the classical mass vaccination or the ring vaccination strategy. In our simulations, vaccine efficacy varied between 40% − 100% whereas population coverage for classical vaccination varied between 10% − 100%. For ring vaccination, the proportion of exposed individuals (%*exposed*) varied as 1% and 3% (see Table 2). For all simulations: *N* = 10^{6}, *b*_{1} = 0.028, *b*_{2} = 0.001, *r*_{1} = 0.09, *r*_{2} = 0.04, *m*_{1} = day 50, *m*_{2} = day 126 and *T* = 360 days. Each scenario is repeated 50 times and the mean of infection cases taken for analysis.

### 2.3 High-Performance Computing

Implementations of this network model was computationally demanding and challenging using conventional resources. For instance, one simulation can take up to several hours or days to complete in a modern desktop computer. Thus, due to these limitations, we employed a High-Performance Computing (HPC) cluster for our simulations. The cluster used is FUCHS-CSC from the Center for Scientific Computing (CSC, Frankfurt, Germany). It is based on 72 dual-socket AMD Magny-Cours CPU compute nodes with 64 GB of RAM, 250 dual-socket AMD Istanbul compute nodes with 32 GB of RAM and 36 quad-socket AMD Magny-Cours compute nodes with 128 GB of RAM each. A simulation in the HPC takes about 5h to complete for the classical vaccination scenarios and about 5 − 8 days for the ring vaccination scenarios. This project is coded in Python.

## 3 Results

Our parameter fitting for *β*(*t*) shows infection peaks at comparable time points with varying population percentages but essentially a qualitative fit of the data (see Figure 2). The parameter values used to fit *β*(*t*) are: *b*_{1} = 0.028, *b*_{2} = 0.001, *r*_{1} = 0.09, *r*_{2} = 0.04, *m*_{1} = 50, *m*_{2} = 126.

### 3.1 COVID-19 dynamics - control scenario

In a completely susceptible population, the introduction of one exposed individual leads to the spread of the infection with more than one peak of cases of infection after some months. In Figure 3, we see that the mean number of infected cases first peaks in the second month after infection onset with about 12.8% of the population infected. By the fifth month during the infection course, there arises a second peak of infected cases which closely matches the first.

### 3.2 Mass Vaccination Scenarios

These scenarios determine the outcome of vaccination at the initial stages of the epidemic before infected cases peak. When vaccine efficacy is 40%, a population coverage of 50% or more is needed to achieve infection peak with less than 10% infection cases. Table 3 reveals that with vaccine efficacy of 60%, a coverage more than 70% and above keeps mean infection cases on the low with elimination of the peak occurring when coverage is more than 80%. Furthermore, low cases of infection are observed when 70% or more of the population is vaccinated and vaccine efficacy is 80%. On the other hand, a vaccine which is 100% efficacious requires 60% or more of the population to be vaccinated to ensure there is no peak of infections. As seen in Figure 4 (see also Figure S1 - S4 in Supplementary Material), peak(s) of infection cases are observed when a small percentage of the population is vaccinated. For instance, in all scenarios, two infection peaks are observed when %*vac* is 40%. Similarly, there are peaks of infection for vaccine coverage between 20% and 60% even when *η* varies. Therefore, vaccinating a small proportion of the population is not useful in these instances as cases can still peak even with a vaccine with 100% efficacy. Thus, the key to eliminating infection peaks is to administer very efficacious vaccines to many individuals.

### 3.3 Ring vaccination scenario

We carried out simulations using Algorithm 2 to determine the infection outcome when ring vaccination is used. In Figure 5 (see also Figure S5 and Table 4), we show these outcomes with varying scenarios of vaccine efficacy and when 1% or 3% of the population is already exposed to the disease before the onset of vaccination. We see that there are less cases realised with a 1% exposed population before the start of vaccination as compared to a 3% exposed population. In addition, the number of infected cases in these scenarios are considerably lower than that of the classical vaccination method. Also, even with a 100% efficacious vaccine, total eradication of the peak is not achieved irrespective of the exposed population prior to vaccination (See also Figure S5a and Figure S5b). Furthermore, only one peak of infection is seen using this strategy compared to the sometimes two peaks of the mass vaccination strategy.

In comparison to mass vaccination, a lower percentage of the population has to be vaccinated when using a ring vaccination protocol in order to attain low infection cases (see Table 4). This is especially so as with effective contract tracing, more individuals can be vaccinated and thus decreasing the number of infections. It is worth noting also that even with the above results, the percentage of vaccinated individuals is lesser when only 1% of the population is exposed prior to vaccination than when the prior exposed population is 3%. Also, from Table 4 we see that the vaccinated populations are very similar in both scenarios respectively irrespective of the efficacy of the vaccine.

## 4 Discussion

With the ongoing spread of SARS-CoV-2 worldwide, pharmaceutical companies have also been in a race to produce safe and highly effective vaccines to counter transmission of the disease. As of December 2020, there were 52 vaccine candidates in clinical trials in humans, 13 of which are in Phase 3 trials (see Table 5) [60]. In November 2020, some pharmaceutical companies and institutes, including Pfizer Inc and BioNTech, Moderna, the University of Oxford (in collaboration with AstraZeneca), announced positive results from the first interim analyses of their Phase 3 vaccine trials [46, 48, 37]. Initial data released from these trials report that the vaccines manufactured by Pfizer Inc/BioNTech and Moderna both yielded 95% efficacy whereas that by University of Oxford (in collaboration with AstraZeneca) yielded 70% efficacy. On 2 December 2020, the UK medicines regulator MHRA granted a temporary regulatory approval for the Pfizer-BioNTech vaccine [43] which was under evaluation for emergency use authorization (EUA) status by the United States Food and Drug Administration [57] and approved for use on 11 December 2020 [43] in the United States.

Although the Pfizer, Moderna and University of Oxford/AstraZeneca vaccines are expected to be used worldwide in coming months, in other countries such as China and Russia, other vaccines (from Sinovac in China and Gameleya Research Institute in Russia) have been developed and are currently in limited use in these countries. Table 5 summarizes the candidate vaccines under phase 3 trial and their current usage status as well as their manufacturers and vaccine platform (that is, a system which uses certain basic components as the backbone but is relatively flexible and can be adapted quickly to be used against different pathogens [24]).

In this study, we modelled the spread of SARS-CoV-2 infection using an SEIRV model structure on a social network. We run stochastic simulations to determine vaccine efficacy and population coverage limits capable of extinguishing the disease. We considered two vaccination strategies, each with varying scenarios regarding vaccine efficacy and population coverage. We found that, the introduction of a single infected person into a completely susceptible population leads to the spread of infection giving rise to more infected cases and subsequently more than one infection peak. This is an indication that in the absence of a good enough vaccine, there is a much higher chance of more than one infection peak which is already being observed in some countries [8, 6].

In addition, the introduction of vaccination lowered the number of infected cases irrespective of the type of vaccination. Our simulations and analysis show that when using a mass vaccination strategy, the lower the efficacy of the vaccine, the more people needed to be vaccinated against the disease in order to eliminate infection peak(s). For instance, a vaccine with an efficacy of 40% will require more than 80% of the entire population to be vaccinated in order to reach low infection peak values whereas a 100% effective vaccine will require more than 60% of the population to be vaccinated. Given the experimental nature and limited initial supply of vaccines, a classical mass vaccination campaign might not be feasible in many countries. However, ring vaccination of likely case contacts and contacts of cases could provide an effective alternative in distributing the vaccine to ensure low levels of infection and subsequently preventing infection peaks. This is with the assumption of effective contract tracing of infected people.

Researchers, policy makers and the general public are of the opinion that the discovery of a vaccine will allow the return to normality of life before the SARS-CoV-2 pandemic. It is worth noting that, the discovery of a proven-to-be effective vaccine however, might not reduce transmission completely. This is due to the fact that a vaccine which effectively reduces the severity of transmission does not necessarily reduce virus transmission to a comparable degree.

It is also important to consider the potential impact of voluntary mass vaccination in the efforts of clearing the epidemic [45, 49, 31, 11]. In years past, the roll-out of vaccination has been faced with declining vaccine confidence in the general public which could possibly lead to hesitancy in getting vaccinated [12]. Such instances could easily lead to a disruption of people receiving a vaccine voluntarily thus, having a detrimental effect on efforts to eradicate the disease and hence such situations should not be underestimated.

## Limitations and future work

Though exploring the effects of SARS-CoV-2 transmission on social networks, in this study we have limited our analysis to one type of social network which is the Erdos-Renyi network. To proceed towards increasing practicality, the analysis performed in this study can be extended to consider disease spread and vaccination on other social networks, such as the Barabási–Albert network [5] and small-world network [58] models, using similar scenarios and protocols. Such network models can be compared with each other and the outcomes analysed. Data on known social networks such as contact matrices and mixing patterns can also be used for further analysis and evaluation.

Simulations in this study were performed with the assumption of an optimistic condition in which there was vaccination of all contacts and contacts of contacts. In reality, contact tracing is especially difficult in the course of an ongoing epidemic and thus reduces the impact of vaccination efforts [34, 29]. Also, the model used here assumes there is equal mixing of individuals and their neighbours in the network whereas in reality this is not the case and consequently infection cases could be lower and population coverage may be reduced.

The model also assumes that all infected people recover from the disease and are immune. The effect of mortality on the dynamics of the disease is not considered as this model aims to study the general transmission dynamics and the effects of varying vaccination strategies. Within-host dynamics such as the viral dynamics and immune responses would be important to be included in a framework with higher resolution [44].

Epidemiological models of disease spread and transmission generally take into account the reproduction number (*R*_{0}) of the disease to gain knowledge on the transmission process. Estimates of *R*_{0} vary widely as data continue to emerge in an ongoing epidemic. This study however, does not examine this parameter since it does not directly impact the transmission dynamics.

In summary, this study describes the transmission dynamics of SARS-CoV-2 on a social contact network and the vaccine efficacy thresholds needed to prevent new waves of the disease. For future work, we will focus on other network models and their analysis with respect to SARS-CoV-2. We will also consider the transmission potential based on *R*_{0} and the dynamics in this regard.

## Data Availability

The code and the data-sets supporting this article have been uploaded as electronic supplementary material.

https://github.com/Josephine-Tetteh/COVID-19-Network-Model.git

## Supplementary materials

## References

- [1].↵
- [2].↵
- [3].↵
- [4].↵
- [5].↵
- [6].↵
- [7].↵
- [8].↵
- [9].↵
- [10].↵
- [11].↵
- [12].↵
- [13].↵
- [14].↵
- [15].↵
- [16].↵
- [17].
- [18].↵
- [19].↵
- [20].↵
- [21].↵
- [22].↵
- [23].↵
- [24].↵
- [25].↵
- [26].↵
- [27].↵
- [28].↵
- [29].↵
- [30].↵
- [31].↵
- [32].↵
- [33].↵
- [34].↵
- [35].↵
- [36].↵
- [37].↵
- [38].↵
- [39].↵
- [40].
- [41].
- [42].
- [43].↵
- [44].↵
- [45].↵
- [46].↵
- [47].↵
- [48].↵
- [49].↵
- [50].↵
- [51].↵
- [52].↵
- [53].↵
- [54].↵
- [55].↵
- [56].↵
- [57].↵
- [58].↵
- [59].↵
- [60].↵
- [61].↵