Abstract
Understanding the impact of behavior dependent mobility in the spread of epidemics and social disorders is an outstanding problem in computational epidemiology. We present a modelling approach for the study of mobility that adapts dynamically according to individual state, epidemic/socialcontagion state and network topology in accordance with limited data and/or common behavioral models. We demonstrate that even for simple compartmental network processes, our approach leads to complex spatial patterns of infection in the endemic state dependent on individual behavior. Specifically, we characterize the resulting phenomena in terms of phase separation, highlighting phase transitions between distinct spatial states and determining the systems’ phase diagram. The existence of such phases implies that small changes in the populations’ perceptions could lead to drastic changes in the spatial extent and morphology of the epidemic/social phenomena.
Introduction
Understanding and modelling population mobility during epidemics has been highlighted as a key area of importance for epidemiology^{1,2,3}, especially in cases when the epidemics occur over a prolonged time period. Recent Ebola virus disease (EVD) outbreaks in western Africa during 2014–2016, and the even more recent outbreaks in heavily populated regions in the Democratic Republic of the Congo (DRC) emphasized the complexity of the disease spread in volatile and diverse areas. The efforts to contain the virus have been hindered by the high rates of population movement from the affected areas to other areas of the DRC and across borders to neighboring countries, creating complex paths and patterns of the infection spread.
In cases like these, high rates of mobility characterize both infected and uninfected populations due to the general lack of security, longheld mistrust of authority and widespread community distrust. Hence, the nature and dynamics of spatial spread of infectious diseases^{4,5} must be investigated and modelled explicitly, ideally coupled with more established networkbased^{6,7,8,9,10,11,12,13,14,15,16,17} and/or agentbased methods of computational epidemiology^{18,19,20,21}. However, there are several challenges on this path, making the development of refined epidemic models problematic.
One significant obstacle is the lack of highresolution data, both in terms of contact rates and transmission vectors, specified within and across suitable social mixing contexts. Despite successful attempts to utilize proxy data, such as mobile phone usage data, in retrospective modelling of the population mobility in highly volatile regions during the Ebola epidemic in Sierra Leone in 2015^{22}, there are several outstanding questions about spatial epidemic expansion, related to (i) the role of densitydependent transmission, especially in situations where mobility restrictions may increase or decrease close congregation, or in cases with highly local superspreading events (e.g., burials), and (ii) the role of longdistance travel.
Another challenge for spatial modelling is the heterogeneity and sensitivity of spatial epidemic patterns to changes in several underlying parameters (behavioral and disease related)^{22}. In general, the resultant patterns are fairly intricate, demanding more refined mechanisms going beyond simple gravitybased interactions or a wavelike diffusion^{2,20,23,24,25}. when in realworld epidemics human mobility may not agree with these standard predictions^{26}. In many cases these more refined models have considered time variation of the underlying network^{27,28,29}. In related work, a mathematical model of the spatial development of social disorder, treated as a contagion, was shown to simulate emergent spatial patterns, while considering the distances travelled to riot locations and the deterrent effect of policing^{30}. Importantly, this spatial interaction model utilized the concept of the rational offender contemplating the potential utilities (rewards) available at different spatial locations, while minimizing costs and trying to avoid capture. Such examples^{30} highlight that spatial contagion models that address the question of individual mobility would have utility beyond the narrower scope of canonical epidemiological models. Indeed, we may recognise that any evolving social phenomena that shares the underlying contact/interaction basis for spread of a quantity of interest stands to gain from such refinements, with examples of such social dynamics including opinion dynamics, population segregation, social unrest and economic competition, where people may “go mad in herds, while they only recover their senses slowly, one by one”^{31}.
The optimal way to model the behaviour associated with individual mobility in such models, however, remains an open question. In this study we aim to provide a rigorous framework that incorporates such dynamic human decisionmaking into commonly used epidemic/social phenomena models which is capable of explaining and categorizing spatial heterogeneity arising from variable mobility under stress. In doing so usual parameters impacting the evolution of the phenomena (such as the basic reproductive ratio R_{0}^{4,32,33} in epidemic models), are augmented by parameters accounting for varying levels of rationality.
When such behavior is introduced, cases of particular interest arise when the population contains groups with differing preferences and/or misaligned interests. In such cases the resulting behaviour may lead to several qualitatively distinct solutions, due to the possibility of emergent frustrated dynamics. This may be compared with pattern formations studies carried out in the context of vegetation models in ecology^{34,35}, preytaxis and predator avoidance systems^{36,37,38,39,40,41} and other similar diffusive processes. This may lead to the existence of phase transitions between distinct regimes of behavior in a space of control parameters. Population groups with varied compartment dependent^{42,43} and mobility dependent behaviors have appeared in related models^{30,44}, but to our knowledge critical phenomena have not been systematically investigated in such contexts. In identifying such transitions we adopt approaches from statistical physics, such as percolation theory and informationtheoretic measures of criticality.
The study of critical dynamics more broadly is of great importance to epidemiology, complex systems and network processes. For example, slowly varying an underlying control parameter, such as the transmissibility of the disease, may induce a sudden change in an observable order parameter, such as the disease attack rate^{4,45}. Indeed, it has been argued that a maladapted pathogen, initially with \({R}_{0} < 1\), can lead to an epidemic if genetic variations, events such as crossing the species barrier, or changes in the host population cause R_{0} to exceed 1^{46}. Other studies indicate that tracing an approach towards an epidemic threshold can also improve prediction and prevention of epidemics^{47}. Of particular relevance is the study of Balcan and Vespignani^{48}, wherein the consideration of human mobility and activity patterns on the spread of infectious diseases lead to “a phase transition between a regime in which the contagion affects a large fraction of the system and one in which only a small fraction is affected”.
The modelling framework we present is based on the BoltzmannLotkaVolterra (BLV) methodology^{49}. The dynamics underpinning the framework are comprised of a fast and slowscale component reflecting (i) the preference for temporary relocation and minimizing travel cost, and (ii) the underlying dynamics of transmission and recovery (here implemented as a Susceptible–Infectious–Susceptible (SIS) model), respectively. The fastscale dynamics are modeled using the Maximum Entropy (MaxEnt) method which yields the mobility flow with the least bias given incomplete information about the system. Both the travel impedance and the preference for temporary relocation to safer locations are quantified through Lagrange multipliers interpreted as coefficients capturing the bounded rationality of the affected individuals. The mobility flow is then reevaluated dynamically, affecting a nonlinear feedback mechanism in the slowerscale SIS model, defined across a network. Importantly, we identify several distinct phases of spatial patterning in the endemic state of the epidemic, dependent on the rationality parameters of the compartment populations, identifying critical regimes and characterizing the system through phase diagrams.
The proposed approach to modelling population mobility in the presence of a severe disease or social disorder, is highly flexible, straightforwardly adapting to various compartmental models, network topologies and utility functions. It thus serves as an ideal test bed for investigating, within a single setting, various factors affecting individuals’ choice, so that the critical phenomena emerging in response to changes in human behavior can be understood. This addresses questions such as (i) when do small changes in the population’s perception of the epidemic risks, or social disorder threats, trigger abrupt and significant changes in the spatial distribution of people across the affected areas, (ii) what is the extent of mixing within the affected communities, and (iii) under what conditions are there coexisting stable and unstable regions within a larger affected territory.
Methods
Here we present our framework designed to extend common epidemiological/social models. Specifically, we generalise the SISnetwork model^{50,51}, described in SI Appendix Section 1, to incorporate dynamic mobility that responds to the evolving state of the epidemic and which allows for compartmentdependent behavior. This is in contrast to existing models which assume temporal and compartmentally homogeneous population mobility^{50,51}. Although more complex compartmental models such as SusceptibleInfectedRecoveredSusceptible (SIRS) models could be considered under this framework with state dependent mobility occurring only if an individual has been previously infected, the steady state behaviours of these models will be similar, with a well defined endemic equilibrium which is locally asymptotically stable^{52,53}, with the primary difference being the transient behaviour of the model. As our analysis focuses on the potential of this framework and the long term behaviour of the model, there is little additional benefit to considering more complex compartmental models for the analysis presented in this work. However, in the case of analysis to specific data, these extensions would be necessary. The structure of the model is informed by the BoltzmanLotkaVolterra (BLV) framework^{49,54}, consisting of a slow dynamic, the infection dynamics captured by the SISnetwork model^{50,51}, and a fast dynamic, the dynamic, state dependent mobility terms.
The SISnetwork model considers \(M\) locations, each with its own subpopulation of \({N}_{i}\), \(i\in \{1,\ldots ,M\}\), individuals of which \({I}_{i}\) are infectious and \({S}_{i}={N}_{i}{I}_{i}\) are susceptible. The state of the epidemic is then characterized by the set of all \({I}_{i}\), written I. These locations are organized into a network model which is characterized by a matrix C of edge weights \({c}_{ij}\) associated with mobility from location \(i\) to location \(j\). We note that, in the general case, we may have nonzero diagonal elements \({c}_{ii}\) to specify a cost associated with ‘mobility’ to an individual’s home location. These \({c}_{ij}\) may represent physical distance, effective travel cost etc. according to modelling needs.
We then introduce mobility functions for each compartment of the model: infectious and susceptible. These are denoted \({\phi }_{ij}^{I}({\bf{I}},{\bf{C}})\) and \({\phi }_{ij}^{S}({\bf{I}},{\bf{C}})\), explicitly allowing for network, compartment and infection level dependence in the mobilities. Explicitly, \({\phi }_{ij}^{x}({\bf{I}},{\bf{C}})\) is the fraction of the population of compartment \(x\) originating from location \(i\) that will mix in location \(j\), dependent on the topology C and the state of the epidemic I. Since these functions change instantaneously in response to any changes in I they constitute a fast dynamic assumed to equilibrate on a time scale faster than the evolution of the epidemic. These are brought into the standard SISnetwork model such that the progression of the epidemic is given by
where
where \({\hat{N}}_{j}({\bf{I}},{\bf{C}})\) is the total number of individuals mixing at location \(j\) (see Eq. 2), \(\gamma \) is the individual recovery rate and \(\beta \) is the transmission rate between a single infected and susceptible individual within the same location. Consequently, the terms \(\beta {\phi }_{ij}^{S}{\phi }_{kj}^{I}{S}_{i}{I}_{k}/{\hat{N}}_{j}\) reflect the rate of total new infections for the susceptible population of \(i\) caused by the infected population from \(k\) which occur in location \(j\) at time \(t\). We emphasize that the mobilities \({\phi }_{ij}^{x}({\bf{I}},{\bf{C}})\), \(x\in \{S,I\}\), now vary in time as the values of the \({S}_{i}\) and \({I}_{i}\) evolve, introducing a qualitatively distinct mechanism for dynamic population mobility.
The task which remains is to specify the functional form of these \({\phi }_{ij}^{I}({\bf{I}},{\bf{C}})\) and \({\phi }_{ij}^{S}({\bf{I}},{\bf{C}})\) which could arguably take many different forms. The principle and insight behind the BLV methodology^{49,55} is to recognize that not all choices of these functions are equally reasonable given the limited knowledge of the system available to the modeler. Rather, we should be concerned with identifying the minimal constraints that characterize plausible behavior and choose the most likely, or leastbiased, functional form that is consistent with them. This is the reasoning behind the wellknown MaxEnt methodology^{56,57}. Consequently, we postulate two fundamental behaviors: (i) a perceived ‘cost’ associated with the graph structure captured by a function of the \({c}_{ij}\) and (ii) a perceived ‘benefit’ associated with the spatial configuration and extent of the epidemic captured by a function of \({I}_{i}\) and \({S}_{i}\). Specifically, here, we consider the cost of mobility from \(i\) to \(j\) to be given simply by \({c}_{ij}\), whilst the benefit of mobility from \(i\) to \(j\) to be given by \({b}_{j}={N}_{j}^{1}({N}_{j}{I}_{j})\). We note that the cost and benefit functions are not restricted to the functional forms used here and may, in principle, be any reasoned function of the compartment, C and I or even augmented with additional model details. The MaxEnt method then determines the \({\phi }_{ij}^{x}({\bf{I}},{\bf{C}})\), at each point in time, by maximizing the Shannon entropy of the probability distribution of all individual behavior, such that it is consistent with the current state of the epidemic, I, in addition to constraints on the mean benefit due to mobility for infectious and susceptible individuals
and mean cost of mobility for all individuals
Such a solution, detailed in SI Appendix Section 2, is given by
where \({Z}_{x,i}={\sum }_{j}\,\exp \,({\alpha }^{x}{b}_{j}\omega {c}_{ij})\). The functional form of \({\phi }_{ij}^{x}\) is entirely determined by the Lagrange multipliers \({\alpha }^{x}\) and \(\omega \), which correspond to the mean quantities \({B}^{x}\) and \(C\). Methodologically, such a formulation provides two distinct modelling perspectives which we illustrate in (7). In the first, \({B}^{x}\) and \(C\) could be measured from data, from which the \({\alpha }^{x}\) and \(\omega \) would be inferred allowing the least biased fit to data^{58}. In this case the \({\phi }_{ij}^{x}\) would constitute the minimal empirical model consistent with the observables. In contrast, the second, which is considered in this work, considers the relationship in the other direction. In this case \({\alpha }^{x}\) and \(\omega \) are considered to be free model parameters corresponding to different strengths of behavior in response to concepts of cost and benefit.
Different choices thus correspond to distinct, hypothetical, model populations with different behaviours which we can analyze. In thermodynamic contexts, in which MaxEnt techniques were first conceived^{56,57}, the Lagrange multipliers appear as intensive parameters such as temperature. Analogously, the two modelling perspectives would amount to i) inferring the temperature given data on measured energies and ii) generative modelling of hypothetical situations at different temperatures. In practical terms, each of these perspectives amounts to a distinct way to use this modelling framework. Firstly, a theoretical model based on assumptions around the effects of disease presence on mobility could be compared to real data in order to test hypothesis around state dependent mobility or identify likely outcomes of an epidemic. Conversely, if time series for mobility and infection dynamics exist on the same resolution, the mean quantities could be calculated and used in order to infer the behavioural parameters of different countries, indicating the extent to which individuals mobility is impacted by the presence of disease. The combination of these approaches would allow for the inference of parameters to be performed on historical data (second approach) and the parameters generated in this way could then be used for future forecasting of disease within similar communities (first approach).
This formulation is entirely consistent with behavioral models that consider individuals with bounded rationality^{59,60}. In such a formulation an objective payoff function is proposed which is then paired with a rationality parameter. Here the objective payoff functions are the benefit \({b}_{j}\) and cost \({c}_{ij}\) with corresponding rationality parameters \({\alpha }^{x}\) and \(\omega \). Given a fixed \(\omega \) we may consider a two dimensional rationality phase space \(\{{\alpha }^{S},{\alpha }^{I}\}\in {{\mathbb{R}}}^{2}\) such that \({\alpha }^{S}\) and \({\alpha }^{I}\) may be positive or negative capturing a range of possible behaviours. The degree of rationality is captured by \({\alpha }^{x}\) with the benefit \({b}_{j}\) being perceived as negative/positive compartment \(x\) based on the sign of \({\alpha }^{x}\). I.e. when \({\alpha }^{x} > 0\), compartment \(x\) views low infectiousness as a benefit and so higher rationality increases mobility towards locations with low levels of infection. On the other hand when \({\alpha }^{x} < 0\), compartment \(x\) views low susceptibility as a benefit and higher rationality increases mobility towards locations with relative low levels of susceptible individuals, see Fig. 1. When the \({\alpha }^{x}\) is infinite the payoff function entirely determines behavior associated with that payoff, with all individuals at location \(i\) choosing the option(s) which maximize/minimize the payoff \({\sum }_{j}\,{\phi }_{ij}^{x}({\bf{I}},{\bf{C}})\,{b}_{j}\). In contrast, when \({\alpha }^{x}=0\), behavior is indifferent to the payoff such that there is a flat distribution over all options. For the purposes of this study we insist that \(\omega > 0\) such that edge weights on the network are always considered as a (positive) cost. We emphasize that the above formulation insists on homogeneous statistical behavior in the individuals, since the constraints are in terms of averages over the total population, however this assumption can be relaxed as discussed in SI Appendix Section 2.
Finally, we discuss the model parameters utilized in this study. We emphasize that our methodology, as presented, can operate on any graph structure, with topologically inhomogeneous behavior, and incorporated into any existing compartmental metapopulation model on a graph. However, in order to present, in a transparent and straightforward way, the possible richness in behavior such a methodology can produce we restrict ourselves to the SIS model with topologically homogeneous behavior on a simple square lattice of length \(L\), such that we have \(M={L}^{2}\) locations, with periodic boundary conditions. At each location we consider an arbitrary, normalized, population \({N}_{i}=1\). Elements \({c}_{ij}\) of C are taken to be \(1\) for the self interaction term \({c}_{ii}\) in addition to those connecting each location to its four immediate neighbors and \(\infty \) otherwise such that mobility from site \(i\) only occurs to itself and its topological neighbors. Since all finite \({c}_{ij}\) are equal the choice of \(\omega \) is arbitrary so long as \(0 < \omega < \infty \), simplifying the analysis in this case.
Results
Here we investigate the behavior of the model specified in the previous section. Whereas typical epidemic modelling investigates the qualitative extent and reach of epidemics given a phase space of infection and recovery rates (β and γ), we consider a phase space of differing compartmental behavior through varying values of the rationality parameters \({\alpha }^{S}\) and \({\alpha }^{I}\), with a focus on the regimes of emergent spatial structure in the endemic infection levels.
To do so we focus primarily on the postcritical regime (\({R}_{0}=\beta /\gamma > 1\)), comprehensively exploring the rationality phase space defined by \({\alpha }^{S}\) and \({\alpha }^{I}\) with fixed infection and recovery rates \(\beta =10\) and \(\gamma =5\), allowing us to consider a wide range of potential behaviours. Very broadly we can divide these behaviours into four main types which we associate with the four quadrants of the rationality phase space, demonstrated in Fig. 1 where each quadrant is defined by the signs of the \(\alpha \) parameters. The positive quadrant (\({\alpha }^{S},{\alpha }^{I} > 0\)) represents plausible mixing preferences present in epidemics where all individuals predominantly avoid infectious individuals. The second quadrant (\({\alpha }^{I} < 0,{\alpha }^{S} > 0\)) possesses dominant mobility behaviour similar to Schelling style segregation^{61} where individuals prefer to mix with members of their own group (\(S\) or \(I\)). In the third quadrant (\({\alpha }^{I} < 0,{\alpha }^{S} < 0\)) all individuals consider ‘infection’ as a benefit and seek out the infected population. This type of mixing behavior may be of particular use in the modelling of the spread of certain social phenomena, e.g. rumors, social myths etc. Finally, in the fourth quadrant (\({\alpha }^{I} > 0,{\alpha }^{S} < 0\)), individuals seek members of the opposite group. This behaviour is much more dynamic and unstable, suggesting application to more extreme social contexts such as rioting and social unrest.
We integrate Eq. 1 using a forward Euler method with Δ\(t=0.001\). Since the equations are deterministic the ensemble of behavior given specified model parameters derives entirely from the ensemble of initial conditions which we specify to be an i.i.d uniform random initial infectious population, \({I}_{i}\), for each location on \([0,0.05]\). Evolution of the dynamics leads to stable long term solutions with distinct typical spatial structures dependent on the rationality parameters, despite the uniform topology. Figure 2 illustrates these behaviors with data from direct simulation. We observe four qualitatively distinct regions: (i) a flat spatial distribution of infection (coloured green) which we refer to as the ‘unseparated regime’, (ii) large highly connected domains (coloured red) which we refer to as the ‘connected regime’ (iii) many small disconnected domains (coloured blue) which we refer to as the ‘isolated regime’ and (iv) an alternating checkerboard pattern containing longlived defect structures due to initial heterogeneity (coloured orange), which we refer to as the ‘antialigned’ regime. We mention that the regimes observed in the space defined by \({\alpha }^{S} > 0\) and \({\alpha }^{I} < 0\) are reminiscent of the standard sequence of spatial patterns (gaps, labyrinth, spots) observed in vegetation growth models^{34,35}.
It is important to note, however, that these regimes do not map directly onto the four quadrants laid out in Fig. 1. For example, in the first quadrant, expected to be most relevant for conventional epidemic modelling, we find qualitatively different spatial patterns in the endemic state despite the entire population (both \(S\) and \(I\)) sharing a consistent conception of benefit (i.e. \({\alpha }^{S} > 0\) and \({\alpha }^{I} > 0\)), through quantitative differences in their rationality \(({\alpha }^{I}\) vs. \({\alpha }^{I})\).
The existence of these distinct regimes is highly suggestive of a system with several stable phases separated by phase transitions where the system undergoes large reconfigurations in response to very small variations in system parameters. In order to characterize these regimes, we seek to understand the qualitative picture presented in Fig. 2 mathematically. In doing so we offer evidence that such regimes are indeed well characterized as distinct configurational phases separated by discontinuous transitions. We present three order parameters designed to distinguish between the salient features in each of the four identified regimes. We then investigate the effect that changes in the basic reproductive ratio have upon the phase diagram.
First, we consider the distinction between the unseparated regime and the connected/isolated/antialigned regimes. The behavior in the unseparated regime is qualitatively identical to SISnetwork models without state dependent mobility, with the infection levels being topologically homogeneous. However, the connected, isolated and antialigned regimes exhibit behavior not observed in conventional network models, with regions of relative high and low infection levels forming stable domainlike structures. As either rationality parameter increases we observe a rapid reconfiguration from a singular distribution of infection levels over all locations to one with substantial variance with two identifiable peaks. Such a distinction can be compactly characterized by the standard deviation of the distribution, \(p({I}_{i})\), functioning here as an effective order parameter^{62}. Figure 3a illustrates the discontinuous change in standard deviation from zero to positive values at the region in phase space separating the unseparated and connected/isolated and antialigned phases. This discontinuous behavior is extremely sensitive to changes in \({\alpha }^{I}\) and \({\alpha }^{S}\), even at low system sizes. This is indicative of an order/disorder phase transition where the unseparated regime corresponds to the homogeneous ordered phase and the connected/isolated/antialigned regimes correspond to the disordered phase.
Secondly, we consider the distinction between the phase separating behaviour found in the \({\alpha }^{S} < {\alpha }^{I}\) and \({\alpha }^{S} > {\alpha }^{I}\) regimes. We proceed by utilising the staggered magnetisation, commonly used to study antiferromagnetic systems, to identify the underlying checkerboard structure of the antialignment observed in the \({\alpha }^{I} > {\alpha }^{S}\) region. The staggered magnetisation compares a system to an ideal antiferromagnetic system, a Néel antiferromagnet^{63}, which in our case, consists of an \(L\times L\) two dimensional grid of alternating +1 and −1 values. Given a standard raster ordering of the locations, such that the mth and nth coordinates in the \(x\) and \(y\) dimensions map to the \((m+nL)\)th location, the staggered magnetisation \({M}_{N}\) for our epidemic system is thus defined as
Figure 3b demonstrates that the staggered magnetisation \({M}_{N}\) is capable of distinguishing the antialigned phase from the other phase separated regimes. This is due to the strong checkerboard pattern of higher and lower infection levels which most closely resembles the Néel antiferromagnet. We observe this to hold despite the presence of boundary defects between domains of equally stable, but incompatible configurations of the checkerboard pattern. In contrast, in the other phase separated regimes where we observe larger domains forming ‘dotted’ or ‘snaking’ patterns, locations are much more likely to have a similar level of infection to their neighbors, resulting in a significantly lower staggered magnetisation.
The distinction between the connected and isolated regimes is more nuanced since the transition occurs between two phase separated populations with broad distributions of infection levels. Consequently, we instead turn to measures of percolation^{64,65,66} to characterize the regimes using statistics of cluster formations within the emergent structures. For the purposes of this study we define a cluster as a contiguous region of sites with infection levels above the population mean value. In order to avoid artifacts related to the square topology we consider sites to be contiguous to one another if they are laterally or diagonally adjacent such that each location is contiguous to eight others. Once clusters are identified the system is amenable to the methods used to study percolation on binary lattices^{67}. A key quantity in such studies is the probability that a selected site belongs to a cluster of size r, \(p(r)\), which can be normalized by considering \(\hat{r}=r/{L}^{2}\). A related order parameter from the study of percolation on such lattices, \(\langle \hat{r}\rangle \), is illustrated, where relevant, over the rationality phase space in Fig. 3c. Such a measure is undefined for the unseparated regime due to its singular distribution, but the connected and isolated phases show very noticeably different values of this percolation order parameter with a separating boundary running diagonally through the phase space. This transition is much smoother than that bounding the unseparated phase. Consequently, in Fig. 4, we demonstrate increasingly rapid growth of \(\langle \hat{r}\rangle \), and thus emergence of the connected phase, as \({\alpha }^{S}\) is increased along a line of constant \({\alpha }^{I}=\,27.5\) through the rationality phase space, as larger lattice sizes are considered and finite size effects become negligible. Inset are typical configurations associated with the order parameter \(\langle \hat{r}\rangle \) illustrating the rapid onset of system wide connected regions of high infection at higher levels of \({\alpha }^{S}\) from local, isolated, regions at lower \({\alpha }^{S}\).
The existence of a phase transition between these two regimes can also be observed in model free measures of criticality, namely the Fisher information which is proven to diverge at a phase transition^{68,69,70,71}. This divergence is due to its direct relationship to the rate of change of a corresponding order parameter^{68}. The Fisher information is defined as the expectation
where for the present system \(X={\bf{I}}\). However, the full phase probability distribution \(P({\bf{I}};{\alpha }^{x})\) is impractical to work with, due to its high dimensionality, so instead we exploit the inequality \({F}_{T}(\theta )\le {F}_{X}(\theta )\), where \(T\) is a statistic of \(X\) with equality holding if and only if \(T\) is a sufficient statistic of \(X\). Then, under some relatively weak assumptions, a divergence in \({F}_{T}(\theta )\) implies a divergence in \({F}_{X}(\theta )\) and thus the existence of a phase transition. Utilizing \(\hat{r}\) as a statistic of I, we compute the Fisher information of \(\hat{r}\) about \({\alpha }^{S}\), \({F}_{\hat{r}}({\alpha }^{S})\), for the same line of constant \({\alpha }^{I}=\,27.5\) as in Fig. 4 with increasing lattice sizes, illustrated in Fig. 5. As the lattice size increases the Fisher information of the statistic \(\hat{r}\) peaks at larger and larger values indicative of a divergence in the infinite system limit and thus a phase transition. Shown inset are the distributions \(p(\hat{r})\) demonstrating a rapid change in configurations across the transition, with a concentration of probability at small cluster membership at lower \({\alpha }^{S}\) and a concentration of probability at large cluster membership at higher \({\alpha }^{S}\).
Finally, we demonstrate the interaction between the canonical parameters influencing the endemic infection levels of an epidemic under the SIS compartmental model, namely \(\beta \) and \(\gamma \), and the behavioral parameters \({\alpha }^{S}\) and \({\alpha }^{I}\). This is achieved by varying \(\beta \) from precritical to postcritical values given fixed \(\gamma =5\) and determining the resulting phase diagrams in the behavioral space using order parameters \(\sigma \), \(\langle \hat{r}\rangle \) and \({M}_{N}\). The results are shown in Fig. 6. For \(\beta < \gamma \) we observe no spatially inhomogeneous behavior, consistent with the absence of an epidemic at precritical parameter choices. It is known that such systems will have a typical infection phase transition independent of its topology at \(\beta /\gamma =1\)^{50} and accordingly the system sustains endemic infection levels uniformly at \(\beta =5\) for all values of \({\alpha }^{S}\) and \({\alpha }^{I}\).
As \(\beta \) is increased further, the phase separated regimes (isolated/connected/antialigned) grow in size until for very large values of \(\beta \) the phase space is dominated by the unseparated phase due to the high levels of infection which saturate the dynamics. Consequently we observe that at low (but postcritical) and high \(\beta \), a large discrepancy between \({\alpha }^{S}\) and \({\alpha }^{I}\) is required to observe phase separation, but at intermediate values separation appears at much smaller differences in population behavior.
Discussion
In this article we have proposed and investigated a dynamic spatial interaction model applicable to studies of epidemics and social disorders. The model explicitly allows for state dependent mobilities which coevolve with the underlying phenomenon leading, here, to a rich range of behaviors not observed in models without this feature.
This model offers two significant contributions, one methodological and one phenomenological.
Methodologically, we have incorporated the BLV framework, that has independently found success in modelling social phenomena ranging from urban and economic growth^{55,72,73,74,75} to ecological or crime dynamics^{30,76}, into the framework of epidemic modelling. The application to epidemic modelling was first postulated by Wilson^{49}, which has since remained an unexplored avenue of investigation. This methodology has several important benefits to the study of social phenomena where explicit mechanisms remain elusive yet a need for plausible models exist. It achieves this by constructing distributions of behavior that include the minimum number of assumed extraneous features which also match the data and/or underlying phenomena at hand, often characterized as being the ‘leastbiased’. Given a set of constraints from data or well rationalized utility functions, the implication is that any other choice introduces unjustified assumptions into the model. It is also highly flexible due to its close relationship to economic rationality models and the relative simplicity of modifying or augmenting the set of utility functions which determine the values and number of Lagrange multipliers. In turn the framework is trivial to incorporate into any compartment model, not just the SIS dynamics considered here.
Finally, we discuss the phenomenological implications resulting from the results we have presented. We have implemented a model dynamic that allows for modified individual behavior in response to the evolving state of the epidemic or social phenomena. This may be particularly relevant in cases of severe outbreaks, especially those that lead to or coincide with civil unrest where mobility patterns may be maximally influenced by pertinent local factors rather than regular travel patterns. A consequence of such a dynamic is the emergence of distinct phases related to the spatial structure of infectious spread in the endemic steady state. These regimes are separated by sharp transitions in the rationality phase space. As with the sharp onset of disease spread at the epidemic phase transition, such a phenomenon is of great interest and importance as it implies that large scale, qualitatively distinct, behavior can be highly sensitive to very small changes in human decisionmaking. Significantly, we observe such phenomena even in cases where the dominant mixing preferences are identical but only the relative levels of rationality (\({\alpha }^{S}\) and \({\alpha }^{I}\)) differ. For instance our results imply that a population where all individuals exhibit infection avoiding behaviour may still exhibit phase separation should the rationality parameters of the susceptible and infectious individuals differ sufficiently.We emphasize that in the proximity of these thresholds, small changes in the perception of local risks and benefits experienced by human decisionmakers can trigger sudden and significant changes in the global system behavior. These changes can be of importance to real world applications such as the rapid (re)deployment of medical teams however calibration of these models to realworld data sets is required to produce detailed recommendations.
Further work may focus on the extension of such a model to more complicated compartmental models, such as SusceptibleInfectedRecoveredSusceptivle (SEIRS) or SusceptibleExposedInfectedRecoveredSusceptivle (SEIRS) compartmental models, where infection avoidance behaviours could be considered only for individuals who have been infected in the past, allowing for an in depth analysis of the range of transient behaviours present in such a model. Additionally, we may consider more realistic topologies including more complicated spatial structures and/or include long range, migrationlike, transport of individuals across the system.
References
 1.
Halloran, M. E. et al. Ebola: Mobility data. Sci. 346, 433–433 (2014).
 2.
Wesolowski, A. et al. Commentary: Containing the ebola outbreak — the potential and challenge of mobile network data. PLoS Curr. 6 (2014).
 3.
Berge, T., Lubuma, J.S., Moremedi, G., Morris, N. & KonderaShava, R. A simple mathematical model for ebola in africa. J. Biol. Dyn. 11, 42–74 (2017).
 4.
Kermack, W. O. & McKendrick, A. G. A contribution to the mathematical theory of epidemics. Proc. Royal Soc. Lond. Ser. A, Containing papers a mathematical physical character 115, 700–721 (1927).
 5.
Anderson, R. M. & May, R. M. Infectious diseases of humans: Dynamics and control (Oxford university press, 1992).
 6.
Newman, M. E. J. & Watts, D. J. Scaling and percolation in the smallworld network model. Phys. Rev. E 60, 7332–7342 (1999).
 7.
Moore, C. & Newman, M. E. J. Epidemics and percolation in smallworld networks. Phys. Rev. E 61, 5678–5682 (2000).
 8.
Newman, M. E. J. Spread of epidemic disease on networks. Phys. Rev. E 66, 016128 (2002).
 9.
Keeling, M. J. & Eames, K. T. Networks and epidemic models. J. Royal Soc. Interface 2, 295–307 (2005).
 10.
Meyers, L. A. Contact network epidemiology: Bond percolation applied to infectious disease prediction and control. Bull. The Am. Math. Soc. 44, 63–87 (2006).
 11.
Piraveenan, M., Prokopenko, M. & Hossein, L. Percolation centrality: Quantifying graphtheoretic impact of nodes during percolation in networks. PLoS One 8 (2013).
 12.
Miller, J. C. & Kiss, I. Z. Epidemic spread in networks: Existing methods and current challenges. Math. Model. Nat. Phenom. 9, 4–42 (2014).
 13.
PastorSatorras, R., Castellano, C., Van Mieghem, P. & Vespignani, A. Epidemic processes in complex networks. Rev. Mod. Phys. 87, 925 (2015).
 14.
Keeling, M. The implications of network structure for epidemic dynamics. Theor. Popul. Biol. 67, 1–8 (2005).
 15.
PastorSatorras, R. & Vespignani, A. Epidemic dynamics and endemic states in complex networks. Phys. Rev. E 63, 066117 (2001).
 16.
PastorSatorras, R. & Vespignani, A. Epidemic spreading in scalefree networks. Phys. Rev. Lett. 86, 3200 (2001).
 17.
Erten, E., Lizier, J., Piraveenan, M. & Prokopenko, M. Criticality and information dynamics in epidemiological models. Entropy 19 (2017).
 18.
Eubank, S. et al. Modelling disease outbreaks in realistic urban social networks. Nat. 429, 180 (2004).
 19.
Longini, I. M. et al. Containing pandemic influenza at the source. Sci. 309, 1083–1087 (2005).
 20.
Cliff, O. M. et al. Investigating spatiotemporal dynamics and synchrony of influenza epidemics in Australia: An agentbased modelling approach. Simul. Model. Pract. Theory 87, 412–431 (2018).
 21.
Zachreson, C. et al. Urbanization affects peak timing, prevalence, and bimodality of influenza pandemics in Australia: Results of a censuscalibrated model. Sci. Adv. 4 (2018).
 22.
Peak, C. M. et al. Population mobility reductions associated with travel restrictions during the ebola epidemic in Sierra Leone: Use of mobile phone data. Int. J. Epidemiol. 47, 1562–1570 (2018).
 23.
Viboud, C. et al. Synchrony, waves, and spatial hierarchies in the spread of influenza. Sci. 312, 447–451 (2006).
 24.
Belik, V., Geisel, T. & Brockmann, D. Natural human mobility patterns and spatial spread of infectious diseases. Phys. Rev. X 1, 011001 (2011).
 25.
D’Silva, J. P. & Eisenberg, M. C. Modeling spatial invasion of ebola in west Africa. J. Theor. Biol. 428, 65–75 (2017).
 26.
Wesolowski, A. et al. Impact of human mobility on the emergence of dengue epidemics in Pakistan. Proc. Natl. Acad. Sci. 112, 11887–11892 (2015).
 27.
Valdano, E., Ferreri, L., Poletto, C. & Colizza, V. Analytical computation of the epidemic threshold on temporal networks. Phys. Rev. X 5, 021005 (2015).
 28.
Ren, G. & Wang, X. Epidemic spreading in timevarying community networks. Chaos: An Interdiscip. J. Nonlinear Sci. 24, 023116 (2014).
 29.
Nadini, M. et al. Epidemic spreading in modular timevarying networks. Sci. Reports 8, 1–11 (2018).
 30.
Davies, T. P., Fry, H. M., Wilson, A. G. & Bishop, S. R. A mathematical model of the London riots and their policing. Sci. Reports 3, 1303 (2013).
 31.
Mackay, C. Extraordinary popular delusions (Templeton Foundation Press, 2015).
 32.
Diekmann, O., Heesterbeek, J. A. P. & Metz, J. A. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 28, 365–382 (1990).
 33.
LloydSmith, J. O., Schreiber, S. J., Kopp, P. E. & Getz, W. M. Superspreading and the effect of individual variation on disease emergence. Nat. 438, 355 (2005).
 34.
Gowda, K., Riecke, H. & Silber, M. Transitions between patterned states in vegetation models for semiarid ecosystems. Phys. Rev. E 89, 022701 (2014).
 35.
Gowda, K., Chen, Y., Iams, S. & Silber, M. Assessing the robustness of spatial pattern sequences in a dryland vegetation model. Proc. Royal Soc. A: Math. Phys. Eng. Sci. 472, 20150893 (2016).
 36.
Yong, K., Herrera, E. D. & CastilloChavez, C. From bee species aggregation to models of disease avoidance: The benhur effect. In Mathematical and Statistical Modeling for Emerging and Reemerging Infectious Diseases, 169–185 (Springer, 2016).
 37.
Lee, J., Hillen, T. & Lewis, M. Pattern formation in preytaxis systems. J. Biol. Dyn. 3, 551–573 (2009).
 38.
Bate, A. M. & Hilker, F. M. Preytaxis and travelling waves in an ecoepidemiological model. Bull. Math. Biol. 81, 995–1030 (2019).
 39.
Tulumello, E., Lombardo, M. C. & Sammartino, M. Crossdiffusion driven instability in a predatorprey system with crossdiffusion. Acta Appl. Math. 132, 621–633 (2014).
 40.
Wang, J., Shi, J. & Wei, J. Dynamics and pattern formation in a diffusive predator–prey system with strong allee effect in prey. J. Differ. Equations 251, 1276–1304 (2011).
 41.
Potts, J. R. & Lewis, M. A. How memory of direct animal interactions can lead to territorial pattern formation. J. Royal Soc. Interface 13, 20160059 (2016).
 42.
Rizzo, A., Frasca, M. & Porfiri, M. Effect of individual behavior on epidemic spreading in activitydriven networks. Phys. Rev. E 90, 042801 (2014).
 43.
Zino, L., Rizzo, A. & Porfiri, M. An analytical framework for the study of epidemic models on activity driven networks. J. Complex Networks 5, 924–952 (2017).
 44.
Sander, E. & Topaz, C. M. The zombie swarm: Epidemics in the presence of social attraction and repulsion’. Math. Model. Zombies, Univ. Ott. Press. Ottawa, Can. 265–300 (2014).
 45.
Porta, M. A Dictionary of Epidemiology (Oxford University Press, 2014).
 46.
Antia, R., Regoes, R. R., Koella, J. C. & Bergstrom, C. T. The role of evolution in the emergence of infectious diseases. Nat. 426, 658 (2003).
 47.
O’Regan, S. M. & Drake, J. M. Theory of early warning signals of disease emergence and leading indicators of elimination. Theor. Ecol. 6, 333–357 (2013).
 48.
Balcan, D. & Vespignani, A. Phase transitions in contagion processes mediated by recurrent mobility patterns. Nat. Phys. 7, 581 (2011).
 49.
Wilson, A. Boltzmann, lotka and volterra and spatial structural evolution: an integrated methodology for some dynamical systems. J. The Royal Soc. Interface 5, 865–871 (2008).
 50.
Arino, J. & Van den Driessche, P. A multicity epidemic model. Math. Popul. Stud. 10, 175–193 (2003).
 51.
Stolerman, L. M., Coombs, D. & Boatto, S. SIRnetwork model and its application to dengue fever. SIAM J. on Appl. Math. 75, 2581–2609 (2015).
 52.
Chen, J. An SIRS epidemic model. Appl. Math. J. Chin. Univ. 19, 101–108 (2004).
 53.
Liu, Q., Li, B. & Sun, M. Global dynamics of an SIRS epidemic model with distributed delay on heterogeneous network. Math. Probl. Eng. 2017 (2017).
 54.
Wilson, A. Entropy in Urban and Regional Modelling, vol. 1 (Routledge, 2011).
 55.
Wilson, A. The Science of Cities and Regions: Lectures on Mathematical Model Design (Springer Science & Business Media, 2012).
 56.
Jaynes, E. T. Information theory and statistical mechanics. Phys. Rev. 106, 620 (1957).
 57.
Jaynes, E. T. Information theory and statistical mechanics. II. Phys. Rev. 108, 171 (1957).
 58.
Bialek, W. et al. Statistical mechanics for natural flocks of birds. Proc. Natl. Acad. Sci. 109, 4786–4791 (2012).
 59.
Wolpert, D. H., Harré, M., Olbrich, E., Bertschinger, N. & Jost, J. Hysteresis effects of changing the parameters of noncooperative games. Phys. Rev. E 85, 036102 (2012).
 60.
McKelvey, R. D. & Palfrey, T. R. Quantal response equilibria for normal form games. Games Econ. Behav. 10, 6–38 (1995).
 61.
Schelling, T. C. Dynamic models of segregation. J. Math. Sociol. 1, 143–186 (1971).
 62.
Luque, B., Ballesteros, F. J. & Fernández, M. Variances as order parameter and complexity measure for random boolean networks. J. Phys. A: Math. Gen 38, 1031 (2005).
 63.
Néel, M. L. Propriétés magnétiques des ferrites; ferrimagnétisme et antiferromagnétisme. Annales de Physique 12, 137–198 (1948).
 64.
Bollobás, B. & Riordan, O. Percolation (Cambridge University Press, 2006).
 65.
Stauffer, D. & Aharony, A. Introduction to Percolation Theory (Taylor and Francis, 1991).
 66.
Sander, L. M. & Warren, C. P. Percolation on heterogeneous networks as a model for epidemics. Math. Biosci. 180, 293–305 (2002).
 67.
Newman, M. E. Power laws, pareto distributions and zipf’s law. Contemp. physics 46, 323–351 (2005).
 68.
Prokopenko, M., Lizier, J. T., Obst, O. & Wang, X. R. Relating fisher information to order parameters. Phys. Rev. E 84, 041116 (2011).
 69.
Wang, X. R., Lizier, J. T. & Prokopenko, M. Fisher information at the edge of chaos in random boolean networks. Artif. Life 17, 315–329 (2011).
 70.
Harding, N., Nigmatullin, R. & Prokopenko, M. Thermodynamic efficiency of contagions: A statistical mechanical analysis of the sis epidemic model. Interface focus 8, 20180036 (2018).
 71.
Harding, N., Spinney, R. E. & Prokopenko, M. Phase transitions in spatial connectivity during influenza pandemics. Entropy 22, 133 (2020).
 72.
Wilson, A. & Dearden, J. Phase transitions and path dependence in urban evolution. J. Geogr. Syst. 13, 1–16 (2011).
 73.
Ellam, L., Girolami, M., Pavliotis, G. & Wilson, A. Stochastic modelling of urban structure. Proc. Royal Soc. A: Math. Phys. Eng. Sci. 474, 20170700 (2018).
 74.
Crosato, E., Nigmatullin, R. & Prokopenko, M. On critical dynamics and thermodynamic efficiency of urban transformations. Royal Soc. Open Sci. 5, 180863 (2018).
 75.
Slavko, B., Glavatskiy, K. & Prokopenko, M. Dynamic resettlement as a mechanism of phase transitions in urban configurations. Phys. Rev. E 99, 042143 (2019).
 76.
Wilson, A. G. Ecological and urban systems models: Some explorations of similarities in the context of complexity theory. Environ. Plan. A 38, 633–646 (2006).
Author information
Affiliations
Contributions
N.H., R.E.S. and M.P. designed research; N.H. and R.E.S. performed research; N.H., R.E.S. and M.P. analyzed data; N.H., R.E.S. and M.P. wrote the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Harding, N., Spinney, R.E. & Prokopenko, M. Population mobility induced phase separation in SIS epidemic and social dynamics. Sci Rep 10, 7646 (2020). https://doi.org/10.1038/s41598020641831
Received:
Accepted:
Published:
Further reading

Modelling transmission and control of the COVID19 pandemic in Australia
Nature Communications (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.