Modeling latent infection transmissions through biosocial stochastic dynamics ============================================================================= * Bosiljka Tadić * Roderick Melnik ## Abstract The events of the recent SARS-CoV-02 epidemics have shown the importance of social factors, especially given the large number of asymptomatic cases that effectively spread the virus, which can cause a medical emergency to very susceptible individuals. Besides, the SARS-CoV-02 virus survives for several hours on different surfaces, where a new host can contract it with a delay. These passive modes of infection transmission remain an unexplored area for traditional mean-field epidemic models. Here, we design an agent-based model for simulations of infection transmission in an open system driven by the dynamics of social activity; the model takes into account the personal characteristics of individuals, as well as the survival time of the virus and its potential mutations. A growing bipartite graph embodies this biosocial process, consisting of active carriers (host) nodes that produce viral nodes during their infectious period. With its directed edges passing through viral nodes between two successive hosts, this graph contains complete information about the routes leading to each infected individual. We determine temporal fluctuations of the number of exposed and the number of infected individuals, the number of active carriers and active viruses at hourly resolution. The simulated processes underpin the latent infection transmissions, contributing significantly to the spread of the virus within a large time window. More precisely, being brought by social dynamics and exposed to the currently existing infection, an individual passes through the infectious state until eventually spontaneously recovers or otherwise is moves to a controlled hospital environment. Our results reveal complex feedback mechanisms that shape the dependence of the infection curve on the intensity of social dynamics and other sociobiological factors. In particular, the results show how the lockdown effectively reduces the spread of infection and how it increases again after the lockdown is removed. Furthermore, a reduced level of social activity but prolonged exposure of susceptible individuals have adverse effects. On the other hand, virus mutations that can gradually reduce the transmission rate by hopping to each new host along the infection path can significantly reduce the extent of the infection, but can not stop the spreading without additional social strategies. Our stochastic processes, based on graphs at the interface of biology and social dynamics, provide a new mathematical framework for simulations of various epidemic control strategies with high temporal resolution and virus traceability. ## I. INTRODUCTION Stochastic processes of epidemic spreading in human society comprise a specific type of critical phenomena where the microscopic-scale interactions give raise to collective dynamics. Thus, mathematical modelling approaches are necessary to understand the nature of the process and control parameters that govern the transition from the micro-to global scale behaviours1–3. Recently, the critical phenomena in social systems have been researched using the concepts developed in physics of complex systems and networks4. Some prominent examples are the emotion spreading in online social networks5–8, opinion dynamics9, and constructive engagement for the collective knowledge creation10. A detailed analysis of empirical data of human activity online and related theoretical modelling6,11,12 provided evidence that the prominent dynamical mechanisms enabling these collective phenomena lie in the self-organised criticality13–15. The appropriate agent-based modeling of these social phenomena16 includes the individual emotional5,6 and cognitive properties10 of the interacting agents. In addition to social dynamics, the epidemic spreading processes involve some essential biological factors, such as the biology of pathogens, and certain health factors of individuals and groups. Recent COVID-19 data17 on the outbreak with the new SARS-CoV-2 virus are an excellent example. In the absence of pharmaceutical interventions against the virus, social and sanitary measures remain of primary importance for controlling the epidemics. In this context, there are new traits to deal with, high latency times, rapid transmission, and the potential of the virus to trigger SARS (severe acute respiratory syndrome), which is a medical emergency. For a summary of COVID-19’s unique properties on pathogenic, epidemiological and clinical issues, see18,19 and references therein. Statistically, about 20% of the infected individuals need hospitalisation, some of whom (making about 13% of all infected) have mild and moderate symtoms. However, 7% exhibit severe symptoms and need intensive care; they can further differentiate such that 4.7% have clinical stage of illnesses, ending up with a high fatality rate. Data from around the world confirm this overall picture with some regional variations20–23. On the other hand, about 80 % of those infected show very mild clinical manifestations or remain completely asymptomatic, so that they recover spontaneously after about two weeks; this group of infected individuals often remains undetected and outside formal health management procedures. Therefore, mortality statistics can serve as an indicator of the actual number of infected individuals. The mortality rate ranges from 1.7 % to 9 % of registered cases, depending on the country or region. In addition to organisational health issues, it is hypothesized that such variations in clinical manifestations may be associated with the existence of different strains and potential genetic mutations of the virus20,24. Theoretically, a virus that jumps on a new host can trigger its evolutionary development in the direction of better adaptation to human cells, which can make the pathogen progressively less aggressive towards the host25–27. Such mutations of the virus would be highly desirable to mitigate current epidemics. Currently, various possible mutations in the SARS-CoV-2 virus are being observed, giving rise to an open hot topic of discussion among researchers28–30. In the traditional modelling approaches based on mean-field equations, the standard SIR model1 has been extended to take into account the above-mentioned features of the epidemic manifestation by distinguishing between four SIRU31 or six SEIHR32 infection stages and groups involved. Acronymes indicate the initial letter of “Suspected”, “Infected”, “Recovered”, extended by the groups of “Exposed”, “Undetected”, “Hospitalised”, that have own co-evolutionary dynamics. Recently, due to the characteristics of COVID-19 disease, a model has been extended to take into account eight different groups SIDARTHE33. These models with a large number of phenomenological parameters that are adapted to the actual data, were able to describe the infection curve (increase in the cumulative number of infected individuals) as well assess the effects of social isolation on the flattening of the curve34. A new and promising line of research is opening up through the microscopic agent-based modeling of the epidemic processes35. These type of models are increasingly used for describing several specific issues of COVID-19 like epidemics36–40. Respiratory droplets and contacts are considered to be the primary routes of transmission of SARS-CoV-2 virus41. However, transmission via passive objects (fomits) as well as several other mechanisms (aerosol and fecal-oral transmission) are also reported as highly possible18,42–44. In addition to hospital equipments, various passive objects can get contaminated, e.g., by contacts or respiratory droplets of an infected individual and the infection can be transmitted to a new host. This indirect transmission mechanism is becoming increasingly intriging given the reported long survival time of the SARS CoV-2 virus on different surfaces43,44. Indirect exposure to the virus and the large number of undiagnosed cases in current COVID-19 epidemics underscore the importance of *latent infection transmission* as a new face of epidemic spreading. So far, this problem has remained outside the radar for standard modeling approaches. In this paper, we develop an agent-based model that adequately describes these latent transmissions of the infection at microscopic scale and the emergence of global patterns. The model takes into account survival time of the virus and key personal characteristics of individuals, such as susceptibility to the virus and exposure time, which are crucial for the process. We simulate an open system, where the agents are generated over time through social activity fluctuations. As a proxy for social dynamics we use an empirical time series of fluctuations in activity collected from online social networks, precisely a segment of MySpace45. These time series have some essential features of social dynamics— correlation of fluctuations with a typical daily periodicity (circadian cycle)8,45. Hence, the underlying social structure is implicitly represented through temporal correlations and a variable intensity of time series. Using long time series of the social activity with an hourly resolution, which determines the time unit of the simulation step, our simulations span several weeks of real-time processes. In the time window, which lasts up to fourteen days after its first appearance in the process, an agent changes its state from “Susceptible” to “Ex-posed” and possibly “Infected”, followed by either “Hos-pitalised” or “spontaneously-Recovered”, after which it is removed from the process. During this period, an infected agent generates a number of viruses (contagious spots), which remain infective for other agents within the virus survival time. The most susceptible agents are likely to have severe symptoms; they are hospitalised and thus moved to a controlled environment. Among the remaining agents, most are asymptomatically infected carriers of viruses; thus, they contribute to the latent spread of infection within ongoing social dynamics for a long period before their spontaneous recovery occurs. The process is presented as a growing directed bipartite graph composed of infected agents and viruses that they spread. This graph-based presentation allows us to identify the pathways of infection that lead to each infected individual and the number of hops the virus has carried from its origin to the current host. This mathematical framework enables us to simulate possible scenarios with the virus mutation. The development of the network framework based on the stochastic process is essential also due to multifaceted possible side effects and damages that the virus can have in a long run19. The simulated high-resolution process revealed the features of dynamic feedback in different scenarios that lead to altered course of the infection propagation. We determine the impact of the level of social activity with/out lock-down and exposure of each individual and the mutation of the virus on the shape of the infection curves. ## II. THE MODEL An external input drives the system —a time series of social activity *st*, which introduces new agents in the process and sets the time *t*; the resolution is one hour and the total length is *tmax*. Thus the total number of agents is ![Graphic][1] The process is visualised as an evolving bipartite graph with a growing number of nodes and edges. Two types of nodes are the agent nodes representing humans (*Hnode*), and virus nodes representing contagious spots (*Vnode*). They are connected via directed edges from *Hnode → V node → Hnode* oriented in the direction of infection transmission. Besides its ID and the creation time, each Hnode possesses several other properties that can influence the process: the agent’s state (“Uninfected”, “Infected”, “Hospitalised” and “spontaneously-Recovered”); the individual susceptibility to infections, h-factor, which is fixed by the creation of the node as a random number *hi* ∈ [0, 1], as well as its exposition time *![Graphic][2]*, where *Te* is the maximum exposition time in hours. Besides, we define the virus-host *gvh* variable, which keeps track of the number of hops of the virus from its origin (the first infected node in the system) up to the current *Hnode*. A new *Vnode* can be produced at every time step by an active carrier (infected agent) proportionally to the severity of its infection, which is measured by its susceptibility; the virus node remains active for a fixed number of hours *Tv*. Over time, an agent can experience the transition from the “Uninfected” to “Infected” state, followed by one of the possible scenarios, depending on the agent’s susceptibility level, *hi*. Precisely, the highly susceptible agent (whose *hi >* 0.8) is likely to have a severe illness, and its state will change to “Hospitalised” after a random period between two and seven days. Meanwhile, the less susceptible agents (*hi* ≤ 0.8) represent the asymptomatic cases that will stay mildly infected and unreported until eventually reaching the state “spontaneously-Recovered” after fourteen days. After changing its status to “spontaneously-Recovered”, an agent is removed from the dynamics. Similarly, we remove the agent when its status has changed from “Infected” to “Hospitalised”; even though the hospitalised individuals can transfer infection, e.g., when the so-called nosocomial transmission (via health-care objects) occurs, it takes part in a controlled hospital environment. Mean-while, the free active virus carriers who take part in the social dynamics primarily contribute to the latent infection transmission. See schematic flow in Fig. 1. ![FIG. 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F1.medium.gif) [FIG. 1:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F1) FIG. 1: In the simulation time window between the addition of agents through social dynamics time series and their removal by either spontaneous recovery or hospitalisation, different groups of agents and transition between them are schematically indicated: “Susceptible” (S) → “Exposed” (E) → “spontaneously-Recovered” (sR) or → “Hospitalised” (H). For completeness, further differentiation of the “Hospitalised” group to “Intensive-Care” (IC), “Recovered” (R) and “Deceased” (D) groups are shown in the darker square on the lower right corner; since these latter subgroups do not contribute to the dynamics studied in this work, they are not considered explicitly The simulations start with one infected agent, placed on the active carriers list *Hactive*. At each step, an agent in the *Hactive* list can produce a new *Vnode* with a probability proportional to its susceptibility factor *hi*. Thus, more susceptible agents, who are likely to have severe symptoms, can spread viruses more often than those who are barely ill. The fluctuating time series introduces *st* new uninfected agents at each time step *t*. By creating an agent, we identify its creation time *ti* as current time *t* and fix its individual susceptibility factor *hi* and the exposition time *![Graphic][3]*. The agent is placed on *Hexposed* list where it remains for *![Graphic][4]* hours exposed and can get infected by currently active viruses with the infection rate *λt*. Apart from the agent’s susceptibility factor, the infection transmission rate *λt* depends on several factors and fluctuates in time, as explained below, see Eq. (1). Once infected, the agent is removed from the list of exposed agents and placed to the list of active carriers *Hactive*. We keep the identity of the virus that infected the agent and update the number of hops that the virus passed till that infection event occurs (the virus generation *gv*). The time step ends up with revising the contents of each list. The difference between the current time and the node’s time *ti* is computed for each node on a given list. Then the node is removed according to the criteria described above. The detailed program flow is described in the Appendix 1. The transmission rate *λt*, appart from a constant (empirical value) *λ* is given by ![Formula][5] It combines the impact of the current fluctuation of the number *Va*(*t*) of active viruses, i.e., *f* = *dVa*(*t*)*/dt/Ha*(*t*); it is normalised by the active number of carriers *Ha*(*t*), which represents the upper limit of the new viruses at that instant of time. In addition, the infection rate is pro-protional to the individual susceptibility *hi* of the agent in question. As stated above, our network framework allows us to follow the sequence of the virus transmission along the chain of infection events, see Fig. 2, the virus generation *gv*. Thus, we can consider the impact of hypothetical virus mutation along the chain. For this scenario, we note that by passing through a new agent node, the virus *gv* increases by one starting from *gv* = 1 at the first infected individual. We assume that its sufficient weakening can be described as *g*(*gv*) = 2*/*(1 + *gv*). Alternatively, we simulate the case without the mutation, i.e., by fixing the factor *g* = 1. The developed methodology readily allows the analysis of other scenarios as well. ![FIG. 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F2.medium.gif) [FIG. 2:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F2) FIG. 2: (Colour online) Zoom-in the directed bipartite network for the first 72 hours (*λ*0 = 0.23, with mutations) with *Hnodes* (blue) and created by them *Vnodes* (pale). Infection transmission occurred through virus nodes, which are indicated by red colour, in the direction of the edge. The amount of *Vnodes* emitted by one *Hnode* is proportional to its susceptibility factor and the duration of its stay among active carriers. Among white *Vnodes* are recent ones that still can infect a new exposed agent; otherwise, each *Vnode* is non-infective four hours after its appearance. ## III. RESULTS In the simulations, we fix the parameters *λ* = 0.23, *Tv* = 4 hours, the hospitalisation (2-7 days) and spontaneous recovery time (14 days) as well as the thresh-old susceptibility *h* = 0.8 according to the reported empirical data from SARS-CoV-2 epidemics, described in Introduction. By differentiating between the mutation (“gen”) and non-mutation (“g1”) scenarios, we change the maximum exposition time *Te* and control the intensity of the social dynamics (average number of new uninfected agents) by choosing the corresponding time series. We compute the time-evolving number of created agents and virus nodes, the edges between the infected agents and viruses, the effective transmission rate, the number of active carriers and the number of active virus spots, as well as the number of exposed and the number of infected individuals per time step. We also determine the cumulative number of infected, hospitalised and spontaneously recovered individuals for *tmax* = 1364 hours (56 days). Furthermore, we simulate the lockdown scenario with two types of social dynamics for the total period of 1680 hours, and a lockdown–recovery scenarios for *tmax* = 3044 hours. These results are presented in the following three sections. ### A. Infection transmission network and sampled quantities In Fig. 2, we show a part of the bipartite network that embodies the infection transmissions during the first three days. Agents (blue nodes) are enumerated by order of appearance (addition to the network), starting from the first infected agent. Along with the outgoing links from each agent node, we have *Vnodes* (pale colour) that were emitted by the agent during its active infectious period. Some of the *Vnodes*, shown in red colour, appear to infect another agent along the red edge. Meanwhile, the majority of other *Vnodes* are no longer infectious, excluding recently posted ones, which are still infectious and can change the colour to red by connecting with a new uninfected agent. Given that an agent can become infected only once during the period of interest here, the actual network has a tree structure; thus, the uninfected *Hnodes* remain disconnected from this graph (not shown). Following the creation of new agents from a given time series *st*, we sample several time-dependent quantities, as shown in Fig. 3 and in the following figures. We determine the fluctuations in the number of exposed and infected agents, the number of active carriers and active viruses, as well as the cumulative infectious curve and the number of spontaneously recovered and hospitalised agents. Notably, in Fig. 3, we describe the differences between the cases with and without virus mutations, mean-while, the input time series and the other parameters (*λ* = 0.23,*Te* = 24 and *Tv* = 4 hours) are kept the same. In Fig. 4, we show how the effective transmission rates vary in these two cases. Even though the number of ex-posed individuals is practically the same, the number of infected ones per hour is smaller in the case with virus mutations than without mutations, leading to the two different cumulative infection curves in the panel B. A similar difference then occurs in the fractions of hospitalised and spontaneously recovered cases, as shown in the panel C. Computing the temporal variations in the number of active carriers (infected agents) and the number of active viruses nodes completes the picture. As shown in panel D, these quantities are significantly higher in the case without mutations. It should be stressed that, given eq. (1), the fluctuations in these quantities have dynamical feedback to the effective transmission rate. We note that the obtained infection curves can be fitted by logistic function with different parameters, see more in the following section. ![FIG. 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F3.medium.gif) [FIG. 3:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F3) FIG. 3: Comparative simulations for the case of virus mutations *g*(*gv*) (“gen”) and without the mutations (“g1”) resulting from the same social dynamics–time series depicted in the bottom panel. Time fluctuations of the number of exposed agents and the number of infected agents per hour, panel A, and the corresponding infectious curves (the cumulative number of infected agents), panel B. Fits according to the logistic function (parameters shown in the legend). The cumulative number of spontaneously recovered and hospitalised agents are shown in panel C and the respective fluctuation of the number of active carriers and active viruses, in panel D. In panels C and D, the infection curves from panel A are shown by dotted lines, for comparison. ![FIG. 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F4.medium.gif) [FIG. 4:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F4) FIG. 4: (color online) Top: Effective transmission rate *λt* vs time *t* for the case with the virus mutations (red curve) and without mutations (black curve) for *λ* = 0.23 corresponding As the infection network in Fig. 2 shows, the infection path can be followed forwards along with the directed links via red virus nodes. It appears that the number of agents that get infected from a given previously infected one varies from agent to agent, as it is shown in the bottom panel of Fig. 4. Averaging over the infected agents in a given period, we can get a local measure of the infection spreading, a quantity similar to what is in the epidemiology literature known as “R-factor”46,47. Our network presentation of the process clearly shows that the R-factor is given by the relative ratio between two successive layers of the Hnods. Recently, studies have shown that the reproduction rate is a local measure with a limited predictive value, see the discussion and the empirical data analysed in46–48. ### B. Influence of social dynamics and exposure times on the course of the infection curve To assess the impact of the intensity of social dynamics to the infection curve, here, we simulate the scenarios with the social lockdown, which is modelled by another time series; a representative example is shown in the top panel in Fig. 5. Specifically, starting with a moderately high social activity (we use the same time series as in Fig. 3), the process lasts for six weeks, then the input time series is changed. Here we chose another empirical time series (corresponding to the negative-emotion activity in MySpace data set45), which exhibits about four times smaller intensity but also almost absent correlations (the Hurst exponent is close to 0.5). Some simulation results are shown in the main panel of Fig. 5 for the case with the mutations. These results reveal how the reduced social activity leads to an effective flattening of the infection curve, in qualitative agreement with the overall empirical data. (For easier comparison, the parameters are such that the initial part of the bottom curve corresponds to the lower curve in panel B of Fig. 3). However, the plateau level results from the course of the entire curve from the beginning of the infection. In this way, the impact of social lockdown depends on other factors that are built into the infection process before the intervention. In this figure, we further demonstrate how the increased exposure time of the agents leads to the increased level of the plateau, with other parameters fixed. In the following Fig. 6, we show comparative results with/out mutations and two different lockdown scenarios, while keeping the exposure fixed. ![FIG. 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F5.medium.gif) [FIG. 5:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F5) FIG. 5: Moderate social dynamics with a lockdown: The social activity fluctuations *st* (top panel) and corresponding infection curves *It* vs time *t* for different range of exposition times, indicated in the legend (bottom), for the case with virus mutations. ![FIG. 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F6.medium.gif) [FIG. 6:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F6) FIG. 6: Top: Infection curves for two scenarios with the lockdown: high-intensity dynamics and lockdown after week 4, and low-intensity dynamics with the lockdown after week 6; two lines in each case are the scenarios with/without virus mutations. Fits according to the logistic function (see text). Bottom: Temporal fluctuations of the number of exposed agents, which is the same for both scenarios, and the number of carriers and the number of active virus spots corresponding to the infection curves in the top panel. To predict the extent of the infection33,46,49, the course of the infection curve is standardly fitted by the sigmoid (logistic) function ![Formula][6] The parameter *K* is the so-called “carrying capacity”, and *r* is the rate by which the curve reaches it, while the parameters *m* and *n* relate to the beginning and position of the inflexion point. The rationale behind the occurrence of this functional form (and the related derivative, the bell-like function of the infection rate curve) lies in the very nature of the infectious spreading in a *given population size*, a fraction of susceptible individuals that will be infected. In the beginning, the infection spreads to ever-larger number having practically unlimited resources. When the accelerated growth reaches the maximum infection rate (the inflexion point of the infectious curve), the process starts experiencing the limited space: the number of potential susceptible individuals that are not yet infected is reducing. Consequently, the infection rate starts decreasing while the cumulative number of infected cases asymptotically approaches the final capacity *K*. Moreover, all simulations demonstrate that lowering the intensity of social dynamics, even if for a small factor such as natural day-night fluctuations, will affect the process, but with a delay. A particularly considerable lockdown comes in effect only with a ten days delay, cf. Fig. 5. A closer investigation reveals, see Fig. 6 that, after the lockdown, the number of exposed agents starts decreasing, reaching the corresponding lower level after the period comparable to *Te*. Then the number of active carriers takes about a week to ten days longer to reduce and adjust to the lockdown dynamics. The number of active viruses follows this curve with a small delay (4 hours). ### C. Predicting the course of events after lockdown is lifted In the literature, the impact of the imposed social lockdown during COVID-19 pandemics has been investigated from several different angles. For example, apart from economic issues, diverse social and psychological factors have been reported50, which manifest in altering the social dynamics after a particular lockdown is lifted. In this context, our model allows us to simulate different scenarios. Notably, Fig. 7 presents the results where we simulate the impact of the social activity after the 4-weeks lockdown, meanwhile keeping all other parameters at the same level. In one case, a moderate social activity takes part, which also does not have prominent dynamical correlations (here, we amplified the time series that characterises the lockdown). In the other scenario (“back-to-normal”), we use the same time series as in the period before the lockdown and continue the process for another eight weeks. Note that this time series contains a reasonable level of dynamic correlations, as explained above. For comparison, we also simulate the situation without lockdown keeping the original level of the social activity for the period corresponding to 16 weeks of real-time. The results showing the course of the infection curve are given in the central panel of Fig. 7. The driving social dynamics time series and the corresponding number of exposed agents for the scenario without lockdown are shown in the top panel and for two scenarios with the lockdown, in the bottom panel of Fig. 7. These results reveal that the infection curve increases after the lockdown in all cases, but the increase rate is low for the case of moderate social activity. In the “back-to-normal” case, the increase is much faster even if compared to the corresponding segment of the curve “1”, but slower than at the beginning of infection process before the lockdown. Fitting the segment III of the curve gives the overall capacity which is still lower than the projection from the fit of the curve “1”. It is also essential to notice that the curve “1” has a different course from the possible extension driven from fitting its first part (see black dotted line). We can conclude that even small fluctuations that lower the social dynamics *st* (see top panel) can, through the above-described feedback mechanisms, affect the course of the infection curve, thus reducing the overall projection of the infection. ![FIG. 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F7.medium.gif) [FIG. 7:](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F7) FIG. 7: Central panel: Infection curves *It* vs time for the extended period with the naturally fluctuating social dynamics, curve (1), and for two exiting scenarios with the lockdown intervention followed by moderate-intensity dynamics (2) and “back-tonormal” scenario (3). Periods before, during, and after lockdown are indicated as I, II, III. The fluctuating intensity of the susceptible *st* and exposed *et* agents corresponding to the case (1) are shown in the top panel, and to the curves (2) and (3) in the bottom panel. The last part of the signal *st* in the case “back-to-normal” corresponds to the second half of the signal in the top panel. ## IV. DISCUSSION AND CONCLUSIONS We have designed an agent-based model with the high-resolution dynamics of infection process that explicitly observes the survival time of viruses and personal properties of individuals who produce them. We simulate an open system which is driven by social dynamics of the involved agents, and represented by a growing bipartite graph with the agents and viruses nodes. Considering the actual SARS-CoV-2 epidemics, the model adequately takes into account latent infection transmissions by asymptomatic virus carriers as well as indirect transmissions that can occur, for example, through contaminated surfaces during the virus survival period. These processes comprise a significant part of infection transmissions by each agent occurring in the time window which can last up to two weeks since the agent enters the system through social dynamics until it spontaneously recovers or otherwise moves to a controlled hospital environment. The developed mathematical formalism with the infection network makes it possible to trace a path along which the virus hops over different hosts before infecting a particular individual. It thus allows us to model potential mutations of the virus along its evolution path. The results revealed the key components of this biosocial stochastic process that significantly influence the course of infection spread and the predictions of outcomes. Specifically, we have found that: * *The intensity of social dynamics* in conjunction with the individual susceptibility of each agent is crucial for the latent infection transmission. Hence, a lockdown measure becomes effective, but with a typical delay. The simulated process with a high temporal resolution uncovers the underlying mechanisms at work. Lowering the social activity level gradually reduces the number of exposed individuals until it reaches the level of new social activity; with a delay, it causes a decrease in the number of active viruses carries and the number of active viruses that they produce. The opposite trends occur by lifting the lockdown measures. Depending on the renewed social activity, the outcome can be lower than in the case of the process without lockdown. In particular, in the “back-to-normal” situation, the final projection of the number of infected is still smaller or comparable with the original one. However, it can be reached after a much more extended period. It is interesting to point out that much smaller, natural fluctuations in social dynamics (such as day-night or workdays vs weekends) that appear periodically, as in our time series, can have profound effects on lowering the slope of the infection curves. * *The exposure time* of each individual, is another factor that can considerably increase the course of the infection curve even with a low or moderate social activity level. Thus, modifying the exposure time of individuals or groups is an additional essential characteristic to be considered in conjunction with social dynamics measures. * *Virus mutations scenarios* towards gradually reducing the transmission rate can slow down the growth of the infection curve. Even though these virus mutations are favourable, they have no potential to stop the spread of infection without additional social strategies. A new mathematical framework, developed in this work, provides a robust tool for the analysis of biosocial processes in SARS-CoV-2 like epidemics. Its essential new ingredients are high-resolution dynamics, open system, and the network-assisted process presentation. Meanwhile, the traditional mean-field equations and most agent-based models in the literature have limiting factors due to the fixed size of the system. Given the infection-network representation, however, we do not consider explicit physical distance and the impact of the agent’s mobility (see, for example, a different study in35). Instead, the social dynamics that drive the system in our model can be varied, representing different social structures and situations. Moreover, possible extensions of the model are to include a specific set of connections of each agent, a kind of ego-network, and modified transmission rates inside it. Consequently, a more heterogeneous pattern of the exposure times can emerge. However, such non-random distribution of the exposure time per agent could not significantly influence its impact to the global features of the epidemics found in this work. To reveal the nature of the underlying stochastic process, here we have simulated the infection spreading from a single source. In the meantime, the spreading from different sources or different times, which may increase the slope of the infection curve, can also be considered. Lastly, as we already mentioned, the issue of SARS-CoV-2 virus mutations is an open problem that has received increasing attention of researchers in different fields28–30. The infection network presentation allows us to model different ways of how particular virus strains (virus-host interaction) change over time. These aspects of stochastic biosocial processes, among others highlighted here, could represent a compelling direction for further research. ## Data Availability All data used in this work are publicly available Algorithm 1 Latent Infection Transmissions ![Figure8](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/01/2020.07.30.20164491/F8.medium.gif) [Figure8](http://medrxiv.org/content/early/2020/08/01/2020.07.30.20164491/F8) ## Acknowledgments B.T. work supported by the Slovenian Research Agency (research code funding number P1-0044). R.M. is grateful to the NSERC and the CRC Program for their support and he is also acknowledging the support of the BERC 2018-2021 program and Spanish Ministry of Science, Innovation, and Universities through the Agencia Estatal de Investigacion (AEI) BCAM Severo Ochoa excellence accreditation SEV-2017-0718, and the Basque Government fund “AI in BCAM EXP. 2019/00432” ## Appendix Here, we give a more detailed description of the program flow for the model of infection transmissions driven by the dynamics of social activity, see section II. As explained in the main text, the model appropriately takes into account the latent infection transmissions. * Received July 30, 2020. * Revision received July 30, 2020. * Accepted August 1, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. E; W. O. Kermack, A. G. McKendrick and G. T. Walker. A contribution to the mathematical theory of epidemics. Proc. R. Soc. Lond. A, 115700721, 1927. 2. L.J.S. Allen. An Introduction to Stochastic Epidemic Models. Springer, Berlin, 2008. 3. F. Brauer. Mathematical Epidemiology: past, present and future, Infectious Disease Modelling, 2:113–127, 2017. 4. S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes. Critical phenomena in complex networks. Rev. Mod. Phys., 80(4):1275–1335, 2008. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/RevModPhys.80.1275&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262253500006&link_type=ISI) 5. B. Tadić, V. Gligorijević, M. Mitrović, and M. Šuvakov. Co-evolutionary mechanisms of emotional bursts in online social dynamics and networks. Entropy, 15(12):5084–5120, 2013. 6. M. Šuvakov and B. Tadić. Collective emotion dynamics in chats with agents, moderators and bots. Condensed Matter Physics, 17(3):33801:1–12, 2014. 7. M. Mitrović and B. Tadić. Dynamics of bloggers communities: Bipartite networks from empirical data and agent-based modeling. Physica A: Statistical Mechanics and its Applications, 391(21):5264 – 5278, 2012. 8. 1. J. A. Holyst B. Tadić, Milovan Šuvakov, D. Garcia, and F. Schweitzer. Agent-Based Simulations of Emotional Dialogs in the Online Social Network MySpace, Cyberemotions: Collective Emotions in Cyberspace, J. A. Holyst, editor, pages 207– 229. Springer International Publishing, Cham, 2017. 9. H. Xia, W. Huili, and X. Zhaoguo. Opinion dynamics: A multidisciplinary review and perspective on future research. International Journal of Knowledge and Systems Science (IJKSS), 2, 2011. 10. M. Mitrović Dankulov R. Melnik, and B. Tadić. The dynamics of meaningful social interactions and the emergence of collective knowledge. Scientific Reports, 5:12197, 2015. 11. B. Tadić, M. Mitrović Dankulov, and R. Melnik. Mechanisms of self-organized criticality in social processes of knowledge creation. Phys. Rev. E, 96:032307, 2017. 12. A. Dmitriev, V. Dmitriev, and S. Balybin. Self-organized criticality on twitter: Phenomenological theory and empirical investigations based on data analysis results. Complexity, 2019:8750643, 2019. 13. H. J. Jensen. Self-Organized Criticality: Emergent Complex Behavior in Physical and Biological Systems (Cambridge Lecture Notes in Physics). Cambridge University Press, Cambridge, 1998. 14. D. Dhar. Theoretical studies of self-organized criticality. Physica A, 369:29, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.physa.2006.04.004&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000239159900003&link_type=ISI) 15. M. J. Aschwanden. Self-Organized Criticality Systems. Open Academic Press, Berlin Warsaw, 2013. 16. B. Tadić. Modeling behavior of Web users as agents with reason and sentiment, pages 177–186. Novapublishing, N.Y., 2013. 17. Covid-19 data repository by the center for systems science and engineering (CSSE) at Johns Hopkins University. 18. Y. Wang, Y. Wang, Y. Chen, and Q. Qin. Unique epidemiological and clinical features of the emerging 2019 novel coronavirus pneumonia (covid-19) implicate special control measures. Journal of Medical Virology, 92(6):568–576, 2020. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 19. Y. Zhang, X. Geng, Y. Tan, Q. Li, C. Xu, J. Xu, L. Hao, et al. New understanding of the damage of SARS-CoV-2 infection outside the respiratory system. Biomedicine & Pharmacotherapy, 127:110195, 2020. 20. J. T. Wu, K. Leung, M. Bushman, N. Kishore, R. Niehus, P. M. de Salazar, et al. Estimating clinical severity of covid-19 from the transmission dynamics in Wuhan, China. Nature Medicine, 26:506–510, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0822-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 21. A. J. Kucharski, T. W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, et al. Early dynamics of transmission and control of covid-19: a mathematical modelling study. The Lancet Infectious Diseases, 20(5):553 – 558, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1473-3099(20)30144-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 22. F. Castella. Can the COVID-19 epidemics be controlled on the basis of daily test reports? arxiv:2003.06967v4, 2020. 23. C. Lu, R. Gam, A. P. Pandurangan, and J. Gough. Genetic risk factors for death with SARS-CoV-2 from the UK biobank. medRxiv, 2020.07.01.20144592, 2020. 24. K. A. Lythgoe, M. Hall, L. Ferretti, M. de Cesare, G. MacIntyre-Cockett, A. Trebes, et al. Shared SARS-CoV-2 diversity suggests localised transmission of minority variants. bioRxiv, 2020.05.28.118992, 2020. 25. 1. In Michael G. Katze, 2. Marcus J. Korth, 3. G. Lynn Law, and 4. Neal Nathanson Martin T. Ferris, Mark T. Heise, and Ralph S. Baric. Chapter 13 - host genetics: It is not just the virus, stupid. In Michael G. Katze, Marcus J. Korth, G. Lynn Law, and Neal Nathanson, editors, Viral Pathogenesis (Third Edition), pages 169 – 179. Academic Press, Boston, third edition, 2016. 26. 1. M. G. Katze, 2. M. J. Korth, 3. G. L. Law, and 4. N. Nathanson, editors R. W. Doms. Chapter 3 - basic concepts: A step-by-step guide to viral infection. In M. G. Katze, M. J. Korth, G. L. Law, and N. Nathanson, editors, Viral Pathogenesis (Third Edition), pages 29 – 40. Academic Press, Boston, third edition, 2016. 27. 1. M. G. Katze, 2. M. J. Korth, 3. G. L. Law, and 4. N. Nathanson, editors M. Schneider, J. R. Johnson, N. J. Krogan, and S. K. Chanda. Chapter 12 - the virushost interactome: Knowing the players to understand the game. In M. G. Katze, M. J. Korth, G. L. Law, and N. Nathanson, editors, Viral Pathogenesis (Third Edition), pages 157 – 167. Academic Press, Boston, third edition, 2016. 28. J. Zhao, J. Sun, W-T. He, X. Ji, Q. Gao, X. Zhai, M. A. Suchard, et al. Snapshot of the evolution and mutation patterns of SARS-CoV-2. bioRxiv, 2020.07.04.187435, 2020. 29. N. Vankadari. Overwhelming mutations or snps of SARS-CoV-2: A point of caution. Gene, 752:144792, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.gene.2020.144792&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32445924&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 30. L. van Dorp, M. Acman, D. Richard, L. P. Shaw, C. E. Ford, L. Ormond, C. J. Owen, et al. Emergence of genomic diversity and recurrent mutations in SARS-CoV-2. Infection, Genetics and Evolution, 83:104351, 2020. 31. P. Magal and G. Webb. Predicting the number of reported and unreported cases for the covid-19 epidemic in South Korea, Italy, France and Germany. medRxiv, 2020.03.21.20040154, 2020. 32. Q. Lin, S. Zhao, D. Gao, Y. Lou, S. Yang, S. S. Musa, M. H. Wang, et al. A conceptual model for the coronavirus disease 2019 (covid-19) outbreak in Wuhan, China with individual reaction and governmental action. International Journal of Infectious Diseases, 93:211–216, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.058&link_type=DOI) 33. G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. Di Filippo, A. Di Matteo, and M. Colaneri. Modelling the covid-19 epidemic and implementation of populationwide interventions in italy. Nature Medicine, 26:855–860, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0883-7&link_type=DOI) 34. C. Anastassopoulou, L. Russo, A. Tsakris, and C. Siettos. Data-based analysis, modelling and forecasting of the covid-19 outbreak. PLOS ONE, 15(3):1–21, 03 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0232391&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 35. E. Cuevas. An agent-based model to evaluate the covid-19 transmission risks in facilities. Computers in Biology and Medicine, 121:103827, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.compbiomed.2020.103827&link_type=DOI) 36. M. L Jackson. Low-impact social distancing interventions to mitigate local epidemics of SARS-CoV-2. medRxiv, 2020.06.30.20143735, 2020. 37. Z. Burda. Modelling excess mortality in covid-19-like epidemics. arxiv:2006.15583, 2020. 38. D. Kai, G-P. Goldstein, A. Morgunov, V. Nangalia, and A. Rotkirch. Universal masking is urgent in the covid-19 pandemic: Seir and agent based models, empirical validation, policy recommendations. arxiv:2004.13553v1, 2020. 39. C. Wolfram. An agent-based model of covid-19. Complex Systems, 29:87–105, 2020. 40. S.L. Chang, N. Harding, C. Zacherson, O.M. Cliff, and M. Prokopenko. Modelling transmission and control of the covid-19 pandemic in Australia. arxiv:2003.10218v1, 2020. 41. World Health Organization Team. Modes of transmission of virus causing covid-19: implications for IPC precaution recommendations. Scientific brief, updated on 9 July 2020 42. J. Chen. Pathogenicity and transmissibility of 2019-ncova quick overview and comparison with other emerging viruses. Microbes and Infection, 22(2):69 – 71, 2020. Special issue on the new coronavirus causing the COVID-19 outbreak. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F01%2F2020.07.30.20164491.atom) 43. C.E. Enyoh, A.W. Verla, W. Qingyue, D.K. Yadav, M.A. Hossain Chowdhury, B.O. Isiuku, et al. Indirect exposure to novel coronavirus (SARS-CoV-2): An overview of current knowledge. Preprints 2020, 2020040460, 2020. 44. J. Cai, W. Sun, J. Huang, M. Gamber, J. Wu, and G. He. Indirect virus transmission in cluster of covid-19 cases, Wenzhou, China, 2020. Emerging infectious diseases, 26(6):13431345, 2020. 45. M. Šuvakov, M. Mitrović, V. Gligorijević, and B. Tadić. How the online social networks are used: dialogues-based structure of myspace. Journal of the Royal Society Interface, 10(79):20120819, 2012. 46. D. T Christopoulos. A novel approach for estimating the final outcome of global diseases like covid-19. medRxiv, 2020.07.03.20145672, 2020. 47. A. Grant. The reproduction number r for covid-19 in england: Why hasn’t “lockdown” been more effective? medRxiv, 2020.07.02.20144840, 2020. 48. D. Adam. What reproduction number can’t tell us about managing covid-19. Nature, 583:346, 2020. 49. Narayanan C Viswanath. Prediction of covid-19 active and total cases after a fall and rise of cases. medRxiv, 2020.07.02.20145045, 2020. 50. D. Lahiri, S. Dubey, and A. Ardila. Impact of covid-19 related lockdown on cognition and emotion: A pilot study. medRxiv, 2020.06.30.20138446, 2020. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/graphic-2.gif [6]: /embed/graphic-8.gif