Next Article in Journal
Entropy-Based Methods for Motor Fault Detection: A Review
Previous Article in Journal
Revisiting the Surface Brightness Profile of the Stellar Disk with the Statistical Mechanics of the Self-Gravitating System with the Central Body
Previous Article in Special Issue
Decoding the News Media Diet of Disinformation Spreaders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chaos in Opinion-Driven Disease Dynamics

1
Mathematical Institute, University of Koblenz, 56070 Koblenz, Germany
2
Faculty of Information and Communication Technology, Wrocław University of Science and Technology, 50-372 Wrocław, Poland
3
Interdisciplinary Centre for Mathematical and Computational Modelling (ICM), University of Warsaw, 00-927 Warsaw, Poland
4
Faculty of Management, Wroclaw University of Science and Technology, 50-371 Wrocław, Poland
*
Author to whom correspondence should be addressed.
Submission received: 6 February 2024 / Revised: 12 March 2024 / Accepted: 19 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Entropy-Based Applications in Sociophysics)

Abstract

:
During the COVID-19 pandemic, it became evident that the effectiveness of applying intervention measures is significantly influenced by societal acceptance, which, in turn, is affected by the processes of opinion formation. This article explores one among the many possibilities of coupled opinion–epidemic systems. The findings reveal either intricate periodic patterns or chaotic dynamics, leading to substantial fluctuations in opinion distribution and, consequently, significant variations in the total number of infections over time. Interestingly, the model exhibits a protective pattern.

1. Introduction

Understanding the complex dynamics of opinion formation in a society is a challenging cross-disciplinary task where expertise and knowledge from many different fields are required. Mathematical models of opinion dynamics can provide quantitative and qualitative insights into the set of parameters and variables that shape and determine the dynamics and structure of opinions. Thanks to pioneers like Serge Galam in modeling social phenomena with methods borrowed and inspired by statistical physics, whole new sub-fields of research like sociophysics have been established [1] and triggered subsequently work. Meanwhile, there have been a large variety of models designed like the Galam Models [2], Voter model [3], Sznajd model [4], threshold models [5] and bounded confidence models [6], to name just a few.
Mathematical models for the epidemic dynamics of infectious diseases have existed for more than two hundred years, starting with the work of Bernoulli in 1766 on smallpox inoculation [7] and continuing—after a long pause—with Ross and Hudson on Malaria at the beginning of the last century [8].
However, until the COVID-19 pandemics, there was relative little interest among the epidemiological modeling community to incorporate aspects of opinion dynamics into epidemic models. Some notable exceptions are [9,10,11,12,13,14,15,16]. During the COVID-19 pandemics, it became evident that the efficiency of the implementation of so-called Public Health and Social Measures (PHSMs) depend crucially on the level of acceptance in a society and is, therefore, impacted by opinion formation processes. Furthermore, in most countries of the Western hemisphere, we have seen strong polarization dynamics within society and deep digression on what is “the right thing to do” to fight the pandemic. These processes shaped also the attitude toward COVID-19 vaccination and were causal for the large number of deaths in autumn 2021 due to insufficient vaccine uptake in the elderly population in eastern European countries like Poland and Bulgaria (for major works on coupled opinion–epidemic systems since 2008 see [9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]).
In this article, we study one of the many possibilities of coupled opinion–epidemic systems. The opinion dynamics we consider are a slight generalization of the so-called q-voter model (a kind of nonlinear voter process which favors the most prevalent opinion). The number of possible opinions is large and represents a continuum in our setting. For the epidemic dynamics, we consider a simple SIS system. Opinions have a direct impact on the likelihood of becoming infected (one possible interpretation could be to associate opinions with the frequency of mask wearing) and, therefore, cause heterogeneity in the relative share of infections among the different opinion carriers. In the other direction, an infection increases the likelihood of changing an opinion (more precisely, the likelihood of an infected individual to change their opinion is proportional to the the relative proportion of infected of their opinion among all infected). It is easy to show that the coupled system has no stationary solutions in the general case (Section 2.5) despite the fact that in the decoupled system both—the opinion and the epidemic dynamics—converge to stable equilibria.
It is a widespread belief that pure opinion formation systems are unlikely to show chaotic behavior [32], although there are examples of chaos even in simple deterministic opinion models [33]. We observe for the system presented in this article either complex periodic patterns or chaotic dynamics with large fluctuations in the distribution of the opinions causing substantial variations over time in the total number of infected. It has been known that differences in the perceived risk can have dramatic impact on epidemic dynamics [9], and coupled opinion–epidemic systems ‘can exhibit dynamics that do not occur when the two subsystems are isolated from one another’ [17]. For instance, previous works presented the emergence of periodic dynamics in opinion–epidemic models in which behavior dynamics induced an instability of the endemic equilibrium through a supercritical Hopf bifurcation—an effect leading to large oscillations [15,16]. Furthermore, a recent study presented chaotic behavior in a coupled model induced by delayed response of behavior to epidemic variables [31]. Our opinion dynamics and the coupling to the epidemic dynamics differs from the above-mentioned studies and does not involve delayed response or periodic external triggers like seasonality to induce non-stationary dynamics.

2. Materials and Methods—Model Description

The examined dynamical system consists of two coupled dynamical systems: an epidemiological system and an opinion formation system.

2.1. Epidemiological System

In the following, we consider infectious disease dynamics under the impact of opinions that affect the likelihood of becoming infected in a closed population of size N. Opinions are described by a one-dimensional continuous variable x [ 0 , 1 ] . A zero value reflects the point of view in opposition or polarization to the opinion with value 1. Values between 0 and 1 reflect non-extreme opinions. Each individual is required to have one opinion x. Let S ( t , x ) and Z ( t , x ) denote the number of susceptible or infectious individuals with opinion x at time t. A straightforward extension of the classical Kermack–McKendrick SIS model [34] leads to the system
d d t S ( t , x ) = β ( x ) N S ( t , x ) 0 1 Z ( t , y ) d y + γ Z ( t , x ) ,
d d t Z ( t , x ) = β ( x ) N S ( t , x ) 0 1 Z ( t , y ) d y γ Z ( t , x ) .
Here, β ( x ) > 0 denotes an opinion-dependent transmission rate. As a prototypical example, one may consider a linear dependence
β ( x ) = β 0 + ( β 1 β 0 ) x ,
where β 0 and β 1 denote the transmission rates for the two extreme opinions x = 0 and x = 1 . However, it is advised that β 0 β 1 to observe dynamics diverging from standard SIS model. The recovery rate γ > 0 is assumed to be independent of the opinion x.

2.2. Opinion Formation System

To describe the opinion dynamics, we employ a modification of the q-voter model [3]. The classical q-voter model ( q > 1 ) for a population with just two mutually exclusive opinions reads as
v = α ( 1 v ) v q v s . ( 1 v ) q = α v s . ( 1 v ) v q 1 ( 1 v ) q 1 ,
where v denotes the fraction of individuals having one opinion and 1 v denotes the remaining individuals having the alternative opinion. A change of opinion occurs, if an individuals ”meets” q individuals of the other opinion. Then, with rate α , the individuals flips to the other opinion. It is required that α > 0 , since individuals are attracted to the opinion of others. It is obvious that the simple q-model allows for two opinion-polarized equilibria v = 0 and v = 1 and a balanced equilibrium v = 1 2 . The two polarized equilibria are asymptotically stable, whereas the balanced equilibrium is unstable. To see this, we remark that v = 0 , v = 1 and v = 1 / 2 are the roots of the right-hand side f(v) of the above ODE. Computing the derivative at these roots, we observe, that f (0) = f (1) = α < 0 and f (1/2) = ( q 1 ) α ( 1 2 ) q 1 > 0 , showing the local instability or stability of the respective equilibria.
Let U( t , x ) = S ( t , x ) + Z ( t , x ) denote the total number of individuals having opinion x, where x [ 0 , 1 ] denotes a continuous range of opinions, e.g., ranging between complete rejection of a non–pharmaceutical intervention ( x = 0 ) to full agreement with this intervention ( x = 1 ). Then, N = 0 1 U ( t , x ) d x equates to the overall population. When q = 2 and there is a continuous spectrum of opinions x [ 0 , 1 ] , a generalization of the classical q–voter model (2) reads
d d t U ( t , x ) = ( a U ( t , x ) 2 + ϵ ) 0 1 U ( t , y ) k ( x , y ) d y U ( t , x ) 0 1 ( a U ( t , y ) 2 + ϵ ) k ( y , x ) d y .
Here, the non–negative kernel k ( x , y ) : [ 0 , 1 ] 2 R + denotes the confidence (i.e., trust) of individuals of opinion x and y in each others judgment when interacted. With rate a > 0 , this interaction leads to a switch in opinion. Typically, we will assume that k ( x , y ) = ρ ( | x y | ) = ρ ( r ) , depending only on the distance r = | x y | of the two opinions x and y. It seems natural to assume that ρ is decaying with r, i.e., the further apart the two opinions, the less the trust in others judgment. This mechanism is referred to as Bounded Confidence (BC) [6] in this work.
It is important to note that the probability of sampling of an individual with specific opinion in the continuous x equals zero. Therefore, in the system, the interaction of individuals should be interpreted as the attraction of an individual to a particular opinion rather than “meeting” other individuals.
The individual can change their opinion according to q-individuals with rate a. The variable ϵ > 0 denotes a background rate of opinion change independent of encounters with other individuals. The proportion of individuals who take action is defined by parameters a and ϵ with the assumption 0 a + ϵ 1 . The proportion of individuals who do not take any action equals 1 a ϵ .
Some mathematical results for this model can be found in the appendix in Propositions A1 and A2.

2.3. Opinion–Epidemic Model

In this section, we will use z ( t , x ) and u ( t , x ) as a densities of variables Z ( t , x ) and U ( t , x ) , respectively (i.e., z ( t , x ) = Z ( t , x ) 0 1 U ( t , x ) d x and u ( t , x ) = U ( t , x ) 0 1 U ( t , x ) d x ). To couple the opinion dynamics with the infection process, we assume that the rate of changing one’s opinion x scales with the number of infected z ( t , x ) having this opinion. The rational behind this assumption is the following: If an individual of opinion x observes a large number of infected having the very same opinion x, the more likely it is that the individual will change to another opinion y. Furthermore, conversely, if there are only a few infected sharing one’s opinion, then it is less likely for the individual to change that opinion. Scaling all populations with the total population N, incorporating the above idea into Equation (3) and combining it with the the SIS-model (1), we arrive at
d d t z ( t , x ) = β ( x ) u ( t , x ) z ( t , x ) Z ( t ) γ z ( t , x ) ,
Z ( t ) = 0 1 z ( t , y ) d y ,
β ( x ) = β 0 + ( β 1 β 0 ) x ,
d d t u ( t , x ) = 1 Z ( t ) ( a u ( t , x ) 2 + ϵ ) 0 1 z ( t , y ) u ( t , y ) ρ ( | x y | ) d y
u ( t , x ) z ( t , x ) 0 1 ( a u ( t , y ) 2 + ϵ ) ρ ( | x y | ) d y .
Here, Z ( t ) denotes the total number of infected.
A simple model for the interaction kernel ρ ( r ) is the so-called bounded confidence interval
ρ ( | x y | ) = 1 for | x y | τ 0 else . ,
i.e., if the difference between the two opinions x and y exceeds the bounded confidence threshold τ , no interaction occurs. Some analytic results for the asymptotic dynamics of the pure opinion dynamics are provided in Appendix A and Appendix B.

2.4. Discretization of the Opinion Space

To simulate the coupled infection–opinion dynamics (4), we discretize the opinion space [ 0 , 1 ] by n discrete opinions x 1 < < x n . For simplicity, we assume x k = h ( k 1 2 ) for k = 1 , n and h = 1 / n to be equidistantly spaced. Let z i ( t ) = z ( t , x i ) and u i ( t ) = u ( t , x i ) . Then, the opinion-discretized version of (4) reads as
z i = β i u i z i Z γ z i ,
u i = 1 Z ( a u i 2 + ϵ ) k = 1 n z k u k ρ h ( | i k | ) u i z i k = 1 n ( a u k 2 + ϵ ) ρ h ( | i k | ) ,
where
1 1 Z ( t ) = k = 1 n z k ( t ) ,
β i = β 0 + ( β 1 β 0 ) x i ,
z i = d d t z i ( t ) ,
u i = d d t u i ( t )
and ρ h ( d ) = ρ ( h d ) denotes the scaled bounded confidence kernel.
The coupled 2 n -dimensional ODE system (5a,b) can be solved by any standard ODE solver, e.g., a classical Runge–Kutta method with adaptive step sizes. For the simulation implementation details, please see Section 2.6.

2.5. Some Analytic Results for a Simplified Setting

Consider the opinion-discretized system (5) in the simplified setting ρ h 1 and ϵ = 0 . We will show in the following that such a system cannot have a stationary solution except for the case when the likelihood of becoming infected does not depend on the opinion (in that case the system is actually decoupled). Non-trivial equilibria for the infected compartments z i are characterized by the equation
z i = β i u i Z β i Z + γ .
Analogously, the equilibria for the opinion group u i have to satisfy
u i 2 z k u k = u i z i k u k 2 .
Inserting the above relation for z i , yields
k β k u k 2 β k Z + γ = β i β i Z + γ k u k 2 .
Hence, β i / ( β i Z + γ ) has to be a constant independent of i. This is only possible if β i = β is a constant independent of i.
Assuming β i = β to be independent of i, we define A = β Z β Z + γ and hence, z i = A u i . Therefore,
Z = i z i = β Z β Z + γ u i
and thanks to i u i = 1 , we arrive at 1 = β / ( β Z + γ ) , or
Z = 1 γ β .
Feasible solutions for Z < 1 can only exist if β > γ , i.e., if the classical epidemiological threshold condition R 0 = β γ > 1 for the basic reproduction number R 0 is satisfied. The result is likely to hold for the general case ( ρ h 1 and ϵ > 0 ) and values of epsilon that are not too large; however, to date, we have been unable to obtain conclusive analytic results in this direction.

2.6. Simulation Setup

2.6.1. Software

We run simulations in Julia [35] version 1.9.4. To resolve systems numerically, we use DynamicalSystems [36,37] and OrdinaryDiffEq [38]. For the postprocessing and analysis, ChaosTools [37], StatsBase, and FFTW [39] packages are used.

2.6.2. Hardware

Computations are performed on a Cray XC40 (Okeanos), which is part of the ICM computing infrastructure. The system is composed of 1084 computing nodes. Each node has 24 Intel Xeon E5-2690 v3 CPU cores with a two-way Hyper Threading (HT) with 2.6 GHz clock frequency.

2.6.3. Simulations

Simulations were conducted in a discretized opinion space x with space points n. Simulations run for 30,000 time steps with resolution 1. To allow for the system to stabilize its dynamics, the initial 20,000 steps are discarded and only the last 10,000 steps are used for further investigation.
The coupled 2 n -dimensional ODE system (5a) is solved using Verner’s “Most Efficient” Runge–Kutta method with order 6(5) [40].
To analyze the system and find parameter sets with chaotic dynamics, we run simulations with parameter sampling in a grid search manner. We decided to alter values of parameters n, ϵ , and τ in Equation (4). The parameters a, γ , β 0 , and β 1 we kept constant. This setup allowed us to vary and evaluate the influence of the number of discretization space points n, the proportion of individuals who change their opinion in unstructured manner ϵ , and the bounded confidence interval τ that controls the limits of mixing of a particular opinion in the population. We found these parameters to be the most important for the process. We set a to 0.6 to account for the fact that most people are conformist and follow the majority opinion. To make the system over-critical, we set the limits of the function β ( x ) (i.e., β 0 and β 1 ) to be above the γ value. This way, the classical epidemiological threshold condition R 0 = β γ > 1 is satisfied. Values of β 0 and β 1 equal to 0.11 and 0.225 , respectively, were chosen to ensure an increase in transmission with increasing x. The initial opinion distribution was uniform u ( t = 0 , x ) = 1.0 . The distribution of initially infected was uniform and equal to 0.01 (i.e., z ( t = 0 , x ) = 0.01 ).
Fixed parameters and their values used in the grid search are compiled in Table 1. Varied parameters and their range characteristics are compiled in Table 2.

2.6.4. Analysis Methods

Since we observed either periodic and chaotic behavior, we prepared a set of analysis methods that would permit an evaluation regardless of whether the system is periodic or chaotic and provide an insight into dynamics heterogeneity. We used the following methods:
  • autocorrelation;
  • maximum Lyapunov exponent (MLE) [41,42];
  • spectral Shannon; entropy [43,44,45,46];
  • standard Shannon entropy [46];
  • Poincaré maps;
  • Fourier Transform.
Autocorrelation is a standard tool for detecting periodic signals and their frequencies. In the setup, we compute the maximal value of autocorrelation of arbitrarily chosen lags ranging from 150 to 300. Autocorrelation values should be close to 1.0 for the periodic dynamics. Lyapunov exponents measure rates of separation of nearby trajectories in the flow of a dynamical system. In the analysis, we use the maximum Lyapunov exponent (MLE) to account for the maximal exponential separation. A positive MLE indicates chaotic dynamics. Spectral Shannon entropy is a measure of the heterogeneity in a signal and is computed on the basis of the square of the amplitudes of its Fourier transform. The values are than normalized to 1 and the standard Shannon entropy is computed. The closer the signal is to white noise, i.e., chaotic, the higher the entropy. Standard Shannon entropy is used as a measure of the heterogeneity in opinion dynamics. It is computed based on the densities of opinions in each time step. Then, by obtaining the maximum, minimum, mean, and range length of all entropies (i.e., difference of maximum and minimum), we gain an insight into the diversity of opinions. A Poincaré map is a standard descriptive method used to determine the periodicity in a system. In our work, we plot u i ( t ) vs. u i ( t + 1 ) and z i ( t ) vs. z i ( t + 1 ) to visualize orbits. When a system is periodic, it has evident closed loops. The Fast Fourier Transform (FFT) is a standard mathematical method for spectral analysis of discrete signals. It is a tool for describing time series by the frequencies and amplitudes of components.

3. Results and Discussion

In this section, we evaluate the coupled model dynamics according to their chaotic or periodic nature. We performed 2870 simulations using grid search to investigate the consistency, heterogeneity, variety, and quality of the results. We present the numerical results along with their qualitative and quantitative analyses.
In Figure 1, we can see a visualization of the simulation in which the MLE and autocorrelation indicate a chaotic regime. In the top two subplots, we can see visualizations of the raw u and z data. The absence of clear recurrent patterns is a visual indicator for chaotic behavior (in Figure 2, we see also that for the chosen parameter values, we have positive Lyapunov exponents). This is also apparent for the three bottom subplots, which show aggregate Z, entropy of z, and entropy of u. We can see that Z exhibits a deviation in values within ∼ 20 % . Entropy of u exhibit a greater range of fluctuations than entropy of z. This might be due to the fact that u does not translate to the values of z as the nonuniform β ( x ) function is used. Nevertheless, there is little noticeable order in the timelines.
In Figure 3, we can see a visualization of a simulation leading to a periodic regime. In the top two subplots, we can see visualizations of the raw u and z data. A recurrent pattern is readily observable in this example—specifically, u and z are oscillating in a descending order. Even though the pattern seems simple, it exhibits complex dynamics, observable in the plots of Z and both entropies. Similar to the first panel Figure 1, entropy of u exhibits a greater range of values than entropy of z. However, contrary to the previous panel, we can easily distinguish periodicity in each subplot.
In Figure 4, we present extended disordered Z dynamics from Figure 1. Notwithstanding, it is evident that the dynamics of the system exhibit pseudo-periodic intervals that are intermittently disrupted by chaotic regimes. While the root cause of this chaotic behavior remains unexplored in this particular study, it serves as a potential hint for further investigations into the underlying pathways leading to chaos.
Up until this point, particular examples of numerical patterns of chaotic and periodic behavior have been presented. Simulations have had easily observable highly nontrivial chaotic behavior and nontrivial periodic behavior. To give a comprehensive perspective on the complexity of dynamics involved in the system, the assessment of disorder in simulations encompassing a diverse array of parameters are illustrated in subsequent figures.
The heatmaps in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 present the variety of metrics described in Section 2.6.4.
The panel in Figure 2 consists of four heatmaps showing the values for Autocorellation, MLE, Spectral Entropy, and Wavelet Entropy. From Figure 1 and Figure 3, we know that u displays greater entropy variation than z, therefore, we use it to compute autocorrelation and entropy values. The missing tiles in the Autocorrelation heatmap represent the NaN values that mark stationary solutions. The chaotic scenarios are easy to distinguish in the Autocorrelation (low values), MLE (high values), and Spectral Entropy (high values) heatmaps. We can see three nonuniform regions: two in the top and one in the cenere. Overlap of all three is evident with some divergence on the perimeter of chaotic regimes. Convergence of these three measures is a strong computational indicator of chaos in the system. Maximum Wavelet Entropy is an outlier displaying a level of consistency less aligned with the other measures. The chaotic regions in this instance are subtly and delicately outlined. To better visualize the fluctuations and correlation of Autocorrelation and MLE, we include cross-section plots along the x and y axes in the Figure A1.
In Figure 5, there is a panel picture with four heatmaps with statistical properties of entropy of u: minimum entropy, maximum entropy, range length of entropy (difference of maximum and minimum value), and mean entropy. We can see that the larger the BC threshold, the larger the oscillations in the opinion range. For high epsilon and low BC, we have a very consistent and narrow opinion distribution. The highest oscillations in entropy match with the highest MLE in Figure 2.
In Figure 6, there are two heatmaps: mean Z and mean opinion x. The large BC threshold correlates with the low mean of infected and the low mean opinion. Inversely, the lower the ϵ , the lower the mean of the number of infected. This suggests a rational response to epidemics, wherein individuals prioritize self-protection, aligning their opinions with those in favor of a protective regime (lower values of opinion x). Notably, when both mechanisms are involved with the strongest influence (largest ϵ smallest τ ), we can see the highest mean infected, indicating the worst epidemic outcomes. Surprisingly, there is no visible influence of the chaoticity mode on the mean epidemics results. It remains uncertain whether the mode is irrelevant in general or if the chaotic fluctuations are too subtle to have a noticeable impact on the epidemics.
Figure 7 Descriptive statistical analysis of three simulations: one chaotic and two periodic. The chaotic mode simulation is presented in the top row. The middle and bottom rows illustrate simulations with only the τ parameter changed (0.55 → 0.85) and with only the ϵ parameter changed (0.25 → 0.31), respectively. In the two columns on the left-hand side, we have Poincaré maps of z ( x = 5 ) and u ( x = 5 ) . In the two columns on the right-hand side, there are the Fast Fourier Transform results of z ( x = 5 ) and u ( x = 5 ) . For the sake of clear visualization, we present only absolute amplitudes and half of the frequencies (to avoid a symmetric picture). The difference is evident in the panel plots. Poincaré maps have closed smooth loops in the middle and bottom rows, contrary to the behavior in the top row. The spectral analysis in the top row is highly noisy in comparison to those below it. Both periodic and chaotic simulations show highly complex dynamics. Periodic regimes demonstrate nontrivial dynamics even when there is a lack of any noise. We can clearly see that using descriptive statistical methods, we can easily distinguish chaotic from periodic modes and detect the evident noisiness of the chaotic regime.
All the simulations presented up to now had the number of discretization points n set to 10. In Figure 8, we present heatmaps for the autocorrelation and MLE measures according to various discretization resolution ranging from 4 to 20 points. The minimal resolution required for the chaotic regime to occur is 5. The missing tiles in autocorrelation heatmap correspond to NaN values that mark stationary solutions. We can see that the number of confined parameter spaces with chaotic regimes varies between n values and the size of the parameter space revealing chaotic dynamics increases with increasing n. Most of the chaotic parameter regions have unstructured shapes and sizes. The pattern of formation of these areas is noisy and unclear. The positive MLEs marking chaotic space take up to ∼ 90 % of the total space for n = 20 .

Discussion

The results show that a minimum of five spatial points is required to observe chaotic behavior. The manifestation of chaotic dynamics is conditioned upon a sufficient number of states and bounded confidence that facilitates the mixing of individuals with at least three different opinions. This outcome aligns with expectations, as, when mixing only two states, there is just one possible opinion transition route (to the opposite opinion). As a result, the opinion formation dynamics are strictly limited and periodicity is enforced.
Subsequently, when the bounded confidence mechanism is deactivated (i.e., τ = 1.05 ), chaotic behavior cases are still observed. This suggests that while bounded confidence may contribute to chaotic dynamics, it is not an obligatory factor for their occurrence. On the contrary, all instances of chaotic dynamics emerge upon the presence of the ϵ > 0 . This underscores the essential role of opinion change independent of encounters as the mandatory mechanism for chaotic dynamics to develop.
An intriguing aspect of the model is the protective behavior embedded within it. We observe that the strong influence of independent opinion change and BC correlates with worse epidemic outcomes. Interestingly, in scenarios where communication across the population is widespread and uniform, and individuals follow the majority, the average number of infected cases is the lowest. Although we provided numerical evidence for chaotic dynamics only for systems with a finite number n > 4 of opinions (up to n = 20 ), we see no reason why chaotic patterns should become absent for large n, which is to say, a continuum of opinions—Figure 8 indicates that the opposite is true. A systematic numerical investigation of the continuous opinion system is challenging, especially when search through the parameter space is required (we observed runtimes of about 2 days for n = 20 ).

4. Conclusions and Outlook

This study introduces a new coupled model for opinion and epidemic dynamics. Despite the absence of factors like seasonal effects, delayed response, or contrarians, which are known to trigger chaotic or complex dynamical pattern, our simulations and analyses reveal compelling evidence for both chaotic and complex periodic behavior. This observation is unexpected given that the constituent of the decoupled system is deterministic and exhibits only stationary dynamics and the coupling mechanism is quite simple (infected individuals are just more likely to change their opinion than non-infected).
A closer examination of the chaotic timeline unveils pseudo-periodic intervals disrupted by a noisy signal. This phenomenon might be a hint for a possible route to chaos via intermittency. We did not find any evidence for period-doubling bifurcations.
One implication of chaotic behavior in opinion–epidemics coupled systems is the obvious difficulty in forecasting dynamics beyond a time horizon larger than the inverse of the largest Lyapunov exponent [47]. There are several natural extensions of our system. On the epidemic side, it would be interesting to look at SIRS dynamics. Obviously, many more ways of coupling the opinion and the epidemic dynamics are possible, e.g., the opinion could not only impact the likelihood of becoming infected but also the infectivity of an individual. Finally, there are a large number of opinion models studying the impact of heterogeneity of types of individuals in a society on opinion dynamics, like the contrarians and conformists, and it would be interesting to see how such systems coupled with epidemic dynamics behave.

Author Contributions

Conceptualization, T.G., T.K., K.N., R.P., M.S. and J.S.; methodology, T.G., T.K., K.N. and J.S.; software, K.N., R.P. and M.S.; validation, K.N.; formal analysis, T.G., T.K. and M.S.; investigation, K.N. and M.S.; resources, K.N.; data curation, K.N.; writing—original draft preparation, K.N.; writing—review and editing, T.G., T.K., K.N., R.P. and M.S.; visualization, K.N.; supervision, T.G. and T.K.; project administration, T.G. and T.K.; funding acquisition, T.G. and T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This collaboration between the groups in Koblenz, Germany, and in Poland was funded by the DAAD–NAWA joint project “MultiScale Modelling and Simulation for Epidemics”—MSS4E, DAAD project number: 57602790, NAWA grant number: PPN/B-DE/2021/1/00019/U/DRAFT/00001.

Data Availability Statement

The code used in the research is available under https://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.10624921 (accessed on 11 March 2024).

Acknowledgments

We would like to thank the University of Koblenz and the Wroclaw University of Science and Technology for providing the necessary infrastructure and scientific environment for several meetings of the authors in Koblenz and Wroclaw within the DAAD–NAWA joint project. Karol Niedzielewski and Jan Schneider would like to thank the Infodemics Pandemics Summer School 2023 Lübeck for stimulating their interest in coupled opinion and epidemics models.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BCBounded Confidence
FFTFast Fourier Transform
MLEMaximum Lyapunov Exponent
MDPIMultidisciplinary Digital Publishing Institute
DOAJDirectory of Open Access Journals
TLAThree-Letter Acronym
LDLinear Dichroism

Appendix A

Proposition A1.
Consider the q-voter model
d d t u ( t , x ) = u q ( t , x ) u ( t , y ) d y u ( t , x ) u q ( t , y ) d y .
Let x ˜ = argmax u ( t = 0 , x ) . Then x ˜ = argmax u ( t , x ) for all times t 0 .
Proof. 
Let the solution u in the point x ˜ be defined as u ˜ ( t ) : = u ( t , x ˜ ) , and the solution at an arbitrary point x * be defined as u * ( t ) : = u ( t , x * ) . We assume u * ( 0 ) < u ˜ ( 0 ) , so that Δ u ( t ) : = u ˜ ( t ) u * ( t ) > 0 , as well as ( u * ( t ) + Δ u ( t ) ) q = u ˜ q ( t ) , resulting in the equation u * q ( t ) u ˜ q ( t ) = k = 1 q u * q k ( t ) Δ u k ( t ) . It thus holds that
u * ( t ) u ˜ ( t ) = ( u * q ( t ) u ˜ q ( t ) ) u ( t , y ) d y ( u * ( t ) u ˜ ( t ) ) u q ( t , y ) d y = k = 1 q u * q k ( t ) Δ u k ( t ) u ( t , y ) d y + Δ u ( t ) u q ( t , y ) d y
It holds true that u * ( t ) u ˜ ( t ) < 0 for q = 2 , since after division by Δ u ( t ) 0 , it holds that
k = 1 2 u * 2 k ( t ) Δ u k 1 ( t ) = k = 1 2 u * 2 k ( t ) ( u * ( t ) u ˜ ( t ) ) k 1 = u ˜ ( t ) u ( t , y ) d y > u 2 ( t , y ) d y .
So, it also holds that u * ( t ) < u ˜ ( t ) for all q 2 , and thus, u * ( t ) < u ˜ ( t ) for all t 0 . For q 2 , it can be shown by induction that it holds that
k = 1 q u * q k ( t ) ( u * ( t ) u ˜ ( t ) ) k 1 > u q ( t ) d y ,
so that the statement holds for a general q 2 -voter model as well.  □
Proposition A2
(Stability of the equilibrium in the pure 2-voter model). Regarding the system (3), if ϵ is sufficiently large, i.e., ϵ > a , the uniform equilibrium u ( t , x ) = 1 is stable.
Proof. 
We consider a perturbation ansatz u ( t , x ) = 1 + δ e λ t v ( x ) , where 0 1 v ( x ) d x = 0 for δ 1 . If λ < 0 , then, the uniform equilibrium u 1 is stable, otherwise not. Inserting the ansatz into system (3), we obtain
d d t u ( t , x ) = a ( u 2 + ϵ a ) 0 1 u d y u 0 1 ( u 2 + ϵ a ) d y = a 1 + ϵ a + 2 δ e λ t v ( x ) + δ 2 e 2 λ t v 2 ( x ) 1 + δ e λ t 0 1 v ( y ) d y 1 + δ e λ t v ( x ) 0 1 1 + ϵ a + 2 δ e λ t v ( y ) + δ 2 e 2 λ t v 2 ( y ) d y .
Linearizing with respect to δ and using v d y = 0 , we get
d d t u ( t , x ) = δ a ϵ e λ t v s . + O ( δ 2 ) ,
and thanks to d d t u ( t , x ) = δ λ e λ t v , we finally arrive at
1 1 λ = a ϵ + O ( δ ) .
For ϵ sufficiently large, i.e., ϵ > a , we get λ < 0 and, therefore, the stability of the uniform equilibrium u ( t , x ) = 1 .
As one can see from the above derivation, this result also holds true for the generalization to a q 2 -voter model.  □
The above result shows the local stability of the equilibrium in the 2-voter model. Proving global stability is a different issue, particularly for the integro–differential model. The lack of knowledge of a suitable Lyapunov function is the main obstacle when trying to analyze global stability. Therefore, this issue is still the subject of current research and beyond the scope of this paper.

Appendix B

Figure A1. Cross sections along the Epsilon and τ axes of the Autocorrelation and MLE heatmaps in Figure 2. Panel of four scatter plots: Maximum Autocorrelation along the τ axis (top left), MLE along the τ axis (bottom left), Maximum Autocorrelation along the ϵ axis (top right), and Maximum Autocorrelation along the ϵ axis (bottom right). We can see large fluctuations in each of the four plots. The low Autocorrelation values correlate with high MLE values, as a comparison of the top and bottom rows shows. This is expected behavior and increases our confidence in the existence of chaoticity in these regions.
Figure A1. Cross sections along the Epsilon and τ axes of the Autocorrelation and MLE heatmaps in Figure 2. Panel of four scatter plots: Maximum Autocorrelation along the τ axis (top left), MLE along the τ axis (bottom left), Maximum Autocorrelation along the ϵ axis (top right), and Maximum Autocorrelation along the ϵ axis (bottom right). We can see large fluctuations in each of the four plots. The low Autocorrelation values correlate with high MLE values, as a comparison of the top and bottom rows shows. This is expected behavior and increases our confidence in the existence of chaoticity in these regions.
Entropy 26 00298 g0a1

References

  1. Galam, S. Sociophysics: A Physicist’s Modeling of Psycho-Political Phenomena; Understanding Complex Systems; Springer: Boston, MA, USA, 2012. [Google Scholar] [CrossRef]
  2. Galam, S. Sociophysics: A Review of Galam models. Int. J. Mod. Phys. C 2008, 19, 409–440. [Google Scholar] [CrossRef]
  3. Castellano, C.; Muñoz, M.A.; Pastor-Satorras, R. Nonlinear q-voter model. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2009, 80, 041129. [Google Scholar] [CrossRef]
  4. Sznajd-Weron, K.; Sznajd, J. Opinion evolution in closed community. Int. J. Mod. Phys. C 2000, 11, 1157–1165. [Google Scholar] [CrossRef]
  5. Granovetter, M. Threshold Models of Collective Behavior. Am. J. Sociol. 1978, 83, 1420–1443. [Google Scholar] [CrossRef]
  6. Hegselmann, R.; Krause, U. Opinion dynamics and bounded confidence: Models, analysis and simulation. J. Artif. Soc. Soc. Simul. 2002, 5, 1–33. [Google Scholar]
  7. Dietz, K.; Heesterbeek, J.A.P. Bernoulli was ahead of modern epidemiology. Nature 2000, 408, 513–514. [Google Scholar] [CrossRef]
  8. Ross, R.; Hudson, H.P. An application of the theory of probabilities to the study of a priori pathometry—Part III. Proc. R. Soc. Lond. Ser. A Contain. Pap. A Math. Phys. Charact. 1917, 93, 225–240. [Google Scholar] [CrossRef]
  9. Epstein, J.M.; Parker, J.; Cummings, D.; Hammond, R.A. Coupled Contagion Dynamics of Fear and Disease: Mathematical and Computational Explorations. PLoS ONE 2008, 3, e3955. [Google Scholar] [CrossRef]
  10. Wu, Y.; Li, M.; Xiao, J.; Orgun, M.A.; Xue, L. The Impact of Discrimination on the Spread of Infectious Diseases in Complex Networks. New Gener. Comput. 2014, 32, 193–211. [Google Scholar] [CrossRef]
  11. Wu, J.; Ni, S.; Shen, S. Dynamics of public opinion under the influence of epidemic spreading. Int. J. Mod. Phys. C 2016, 27, 1650079. [Google Scholar] [CrossRef]
  12. Pires, M.A.; Crokidakis, N. Dynamics of epidemic spreading with vaccination: Impact of social pressure and engagement. Phys. A Stat. Mech. Appl. 2017, 467, 167–179. [Google Scholar] [CrossRef]
  13. Ni, S.; Weng, W.; Zhang, H. Modeling the effects of social impact on epidemic spreading in complex networks. Phys. A Stat. Mech. Appl. 2011, 390, 4528–4534. [Google Scholar] [CrossRef]
  14. Durham, D.P.; Casman, E.A.; Albert, S.M. Deriving Behavior Model Parameters from Survey Data: Self-Protective Behavior Adoption During the 2009–2010 Influenza A(H1N1) Pandemic. Risk Anal. 2012, 32, 2020–2031. [Google Scholar] [CrossRef] [PubMed]
  15. d’Onofrio, A.; Manfredi, P. Information-related changes in contact patterns may trigger oscillations in the endemic prevalence of infectious diseases. J. Theor. Biol. 2009, 256, 473–478. [Google Scholar] [CrossRef]
  16. Greenhalgh, D.; Rana, S.; Samanta, S.; Sardar, T.; Bhattacharya, S.; Chattopadhyay, J. Awareness programs control infectious disease—Multiple delay induced mathematical model. Appl. Math. Comput. 2015, 251, 539–563. [Google Scholar] [CrossRef]
  17. Sooknanan, J.; Comissiong, D.M.G. Trending on Social Media: Integrating Social Media into Infectious Disease Dynamics. Bull. Math. Biol. 2020, 82, 86. [Google Scholar] [CrossRef] [PubMed]
  18. Sooknanan, J.; Seemungal, T.A.R. FOMO (fate of online media only) in infectious disease modeling: A review of compartmental models. Int. J. Dyn. Control 2023, 11, 892–899. [Google Scholar] [CrossRef]
  19. Fang, F.; Ma, J.; Li, Y. The coevolution of the spread of a disease and competing opinions in multiplex networks. Chaos Solitons Fractals 2023, 170, 113376. [Google Scholar] [CrossRef]
  20. Ali, R.N.; Sarkar, S. Impact of opinion dynamics on the public health damage inflicted by COVID-19 in the presence of societal heterogeneities. Front. Digit. Health 2023, 5, 1146178. [Google Scholar] [CrossRef]
  21. Zanella, M. Kinetic Models for Epidemic Dynamics in the Presence of Opinion Polarization. Bull. Math. Biol. 2023, 85, 36. [Google Scholar] [CrossRef]
  22. Agusto, F.B.; Numfor, E.; Srinivasan, K.; Iboi, E.A.; Fulk, A.; Saint Onge, J.M.; Peterson, A.T. Impact of public sentiments on the transmission of COVID-19 across a geographical gradient. PeerJ 2023, 11, e14736. [Google Scholar] [CrossRef]
  23. Kastalskiy, I.A.; Pankratova, E.V.; Mirkes, E.M.; Kazantsev, V.B.; Gorban, A.N. Social stress drives the multi-wave dynamics of COVID-19 outbreaks. Sci. Rep. 2021, 11, 22497. [Google Scholar] [CrossRef] [PubMed]
  24. Bernardes, A.T.; Ribeiro, L.C. Information, opinion and pandemic. Phys. A Stat. Mech. Appl. 2021, 565, 125586. [Google Scholar] [CrossRef] [PubMed]
  25. Sooknanan, J.; Mays, N. Harnessing Social Media in the Modelling of Pandemics—Challenges and Opportunities. Bull. Math. Biol. 2021, 83, 57. [Google Scholar] [CrossRef] [PubMed]
  26. Jankowski, R.; Chmiel, A. Role of Time Scales in the Coupled Epidemic-Opinion Dynamics on Multiplex Networks. Entropy 2022, 24, 105. [Google Scholar] [CrossRef]
  27. Du, E.; Chen, E.; Liu, J.; Zheng, C. How do social media and individual behaviors affect epidemic transmission and control? Sci. Total Environ. 2021, 761, 144114. [Google Scholar] [CrossRef]
  28. Peng, K.; Lu, Z.; Lin, V.; Lindstrom, M.R.; Parkinson, C.; Wang, C.; Bertozzi, A.L.; Porter, M.A. A multilayer network model of the coevolution of the spread of a disease and competing opinions. Math. Model. Methods Appl. Sci. 2021, 31, 2455–2494. [Google Scholar] [CrossRef]
  29. Carballosa, A.; Mussa-Juane, M.; Muñuzuri, A.P. Incorporating social opinion in the evolution of an epidemic spread. Sci. Rep. 2021, 11, 1772. [Google Scholar] [CrossRef]
  30. Epstein, J.M.; Hatna, E.; Crodelle, J. Triple contagion: A two-fears epidemic model. J. R. Soc. Interface 2021, 18, 20210186. [Google Scholar] [CrossRef]
  31. Wagner, J.; Bauer, S.; Contreras, S.; Fleddermann, L.; Parlitz, U.; Priesemann, V. Societal feedback induces complex and chaotic dynamics in endemic infectious diseases. arXiv 2023, arXiv:2305.15427. [Google Scholar] [CrossRef]
  32. Lim, C.; Zhang, W. Social opinion dynamics is not chaotic. Int. J. Mod. Phys. B 2016, 30, 1541006. [Google Scholar] [CrossRef]
  33. Borghesi, C.; Galam, S. Chaotic, staggered, and polarized dynamics in opinion forming: The contrarian effect. Phys. Rev. E 2006, 73, 066118. [Google Scholar] [CrossRef]
  34. Kermack, W.O.; McKendrick, A.G. Contributions to the mathematical theory of epidemics—I. Bull. Math. Biol. 1991, 53, 33–55. [Google Scholar] [CrossRef]
  35. Bezanson, J.; Edelman, A.; Karpinski, S.; Shah, V.B. Julia: A Fresh Approach to Numerical Computing. SIAM Rev. 2017, 59, 65–98. [Google Scholar] [CrossRef]
  36. Datseris, G. DynamicalSystems.jl: A Julia software library for chaos and nonlinear dynamics. J. Open Source Softw. 2018, 3, 598. [Google Scholar] [CrossRef]
  37. Datseris, G.; Parlitz, U. Nonlinear Dynamics: A Concise Introduction Interlaced with Code; Undergraduate Lecture Notes in Physics; Springer International Publishing: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  38. Rackauckas, C.; Nie, Q. DifferentialEquations.jl—A Performant and Feature-Rich Ecosystem for Solving Differential Equations in Julia. J. Open Res. Softw. 2017, 5, 15. [Google Scholar] [CrossRef]
  39. Frigo, M.; Johnson, S. The Design and Implementation of FFTW3. Proc. IEEE 2005, 93, 216–231. [Google Scholar] [CrossRef]
  40. Verner, J.H. Numerically optimal Runge—Kutta pairs with interpolants. Numer. Algorithms 2010, 53, 383–396. [Google Scholar] [CrossRef]
  41. Lyapunov, A.M. The general problem of the stability of motion. Int. J. Control 1992, 55, 531–534. [Google Scholar] [CrossRef]
  42. Benettin, G.; Galgani, L.; Strelcyn, J.M. Kolmogorov entropy and numerical experiments. Phys. Rev. A 1976, 14, 2338–2345. [Google Scholar] [CrossRef]
  43. Haaga, K.A.; Datseris, G.; Kottlarz, I.; White, A.; Martinuzzi, F.; HeineRugland; Johnson, S.G. JuliaDynamics/ComplexityMeasures.jl: V2.8.0. Zenodo. 2023. Available online: https://zenodo.org/records/8186595 (accessed on 5 February 2024).
  44. Llanos, F.; Alexander, J.M.; Stilp, C.E.; Kluender, K.R. Power spectral entropy as an information-theoretic correlate of manner of articulation in American English. J. Acoust. Soc. Am. 2017, 141, EL127–EL133. [Google Scholar] [CrossRef] [PubMed]
  45. Tian, Y.; Zhang, H.; Xu, W.; Zhang, H.; Yang, L.; Zheng, S.; Shi, Y. Spectral Entropy Can Predict Changes of Working Memory Performance Reduced by Short-Time Training in the Delayed-Match-to-Sample Task. Front. Hum. Neurosci. 2017, 11, 437. [Google Scholar] [CrossRef] [PubMed]
  46. Shannon, C.E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  47. Alfaro, M.; Fuertes, G.; Vargas, M.; Sepúlveda, J.; Veloso-Poblete, M. Forecast of Chaotic Series in a Horizon Superior to the Inverse of the Maximum Lyapunov Exponent. Complexity 2018, 2018, 1452683. [Google Scholar] [CrossRef]
Figure 1. Example of a chaotic timeline for the evolution of the dynamics with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 , where MLE and autocorrelation values indicate chaotic behavior. The timeline begins with 20,000+ time steps and finishes with 20,000 + 2000. We show heatmaps of opinion u ( x , t ) , infected z ( x , t ) , sum of infected Z ( t ) , and entropies of z ( x , t ) and u ( x , t ) , respectively. Entropies of z and u are computed with base 10. In each plot, there are fluctuations and irregularities in the data, especially in u and entropy of u.
Figure 1. Example of a chaotic timeline for the evolution of the dynamics with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 , where MLE and autocorrelation values indicate chaotic behavior. The timeline begins with 20,000+ time steps and finishes with 20,000 + 2000. We show heatmaps of opinion u ( x , t ) , infected z ( x , t ) , sum of infected Z ( t ) , and entropies of z ( x , t ) and u ( x , t ) , respectively. Entropies of z and u are computed with base 10. In each plot, there are fluctuations and irregularities in the data, especially in u and entropy of u.
Entropy 26 00298 g001
Figure 2. Panel of heatmaps with measures of disorder of dynamics of system with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Maximum Autocorrelation (top left), MLE (bottom left), Maximum Spectral Entropy (top right), and Maximum Wavelet Entropy (bottom right) are displayed. In each plot, there are large fluctuations in values. The missing tiles in the Autocorrelation heatmap represent NaN values that mark stationary solutions. The parameter spaces with chaotic simulations are easy to distinguish in the Autocorrelation (low values), MLE (high values) and Spectral Entropy (high values) heatmaps. The low Autocorrelation values, high MLE values, and high Spectral Entropy overlap. This is expected behavior that increases confidence in the existence of chaoticity in these regions. Maximum Wavelet Entropy is an outlier, less consistent with other measures.
Figure 2. Panel of heatmaps with measures of disorder of dynamics of system with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Maximum Autocorrelation (top left), MLE (bottom left), Maximum Spectral Entropy (top right), and Maximum Wavelet Entropy (bottom right) are displayed. In each plot, there are large fluctuations in values. The missing tiles in the Autocorrelation heatmap represent NaN values that mark stationary solutions. The parameter spaces with chaotic simulations are easy to distinguish in the Autocorrelation (low values), MLE (high values) and Spectral Entropy (high values) heatmaps. The low Autocorrelation values, high MLE values, and high Spectral Entropy overlap. This is expected behavior that increases confidence in the existence of chaoticity in these regions. Maximum Wavelet Entropy is an outlier, less consistent with other measures.
Entropy 26 00298 g002
Figure 3. Example of a periodic timeline for the evolution of the dynamics with parameters n = 10 , ϵ = 0.25 , and τ = 0.85 where MLE and autocorrelation values indicate periodic behavior. The timeline begins with 20,000+ time steps and finishes with 20,000 + 200. We show heatmaps of opinion u ( x , t ) , infected z ( x , t ) , sum of infected Z ( t ) , and entropies of z ( x , t ) and u ( x , t ) , respectively. Entropies of z and u are computed with base 10. In each plot, there are fluctuations in the data, especially in u and entropy of u. One can clearly distinguish periodic behavior in each plot.
Figure 3. Example of a periodic timeline for the evolution of the dynamics with parameters n = 10 , ϵ = 0.25 , and τ = 0.85 where MLE and autocorrelation values indicate periodic behavior. The timeline begins with 20,000+ time steps and finishes with 20,000 + 200. We show heatmaps of opinion u ( x , t ) , infected z ( x , t ) , sum of infected Z ( t ) , and entropies of z ( x , t ) and u ( x , t ) , respectively. Entropies of z and u are computed with base 10. In each plot, there are fluctuations in the data, especially in u and entropy of u. One can clearly distinguish periodic behavior in each plot.
Entropy 26 00298 g003
Figure 4. Example of a chaotic timeline for the sum of infected Z ( t ) from Figure 1. A simulation with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 with extended time up to 20,000 + 5000 steps is displayed.
Figure 4. Example of a chaotic timeline for the sum of infected Z ( t ) from Figure 1. A simulation with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 with extended time up to 20,000 + 5000 steps is displayed.
Entropy 26 00298 g004
Figure 5. Panel of heatmaps with measures of dispersion in the dynamics of opinion u ( x , t ) with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Minimum Entropy (top left), Maximum Entropy (bottom left), difference in Maximum and Minimum Entropy (top right), and Mean Entropy (bottom right) are displayed. In each plot, there are large fluctuations in values. We can see that the larger the oscillations in the opinion range, the larger the BC threshold. For high epsilon and low BC, we have a very consistent and narrow opinion distribution. The highest oscillations in entropy match with the highest MLE values in Figure 2.
Figure 5. Panel of heatmaps with measures of dispersion in the dynamics of opinion u ( x , t ) with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Minimum Entropy (top left), Maximum Entropy (bottom left), difference in Maximum and Minimum Entropy (top right), and Mean Entropy (bottom right) are displayed. In each plot, there are large fluctuations in values. We can see that the larger the oscillations in the opinion range, the larger the BC threshold. For high epsilon and low BC, we have a very consistent and narrow opinion distribution. The highest oscillations in entropy match with the highest MLE values in Figure 2.
Entropy 26 00298 g005
Figure 6. Panel of heatmaps with measures of expected values of the sum of infected Z ( t ) and opinion u ( x , t ) with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Mean Z (left) and average opinion (right) are displayed. In each plot, there are large fluctuations in values. The higher the τ , the lower the mean number of infected and mean opinion. The higher the Epsilon, the higher the mean number of infected and mean opinion. When Epsilon is at its largest value and τ is at it smallest value (bottom right corner), we can see the highest mean number of infected.
Figure 6. Panel of heatmaps with measures of expected values of the sum of infected Z ( t ) and opinion u ( x , t ) with n = 10 . The results are from the grid search simulations for ϵ and τ parameters that are on the x and y axis, respectively. Mean Z (left) and average opinion (right) are displayed. In each plot, there are large fluctuations in values. The higher the τ , the lower the mean number of infected and mean opinion. The higher the Epsilon, the higher the mean number of infected and mean opinion. When Epsilon is at its largest value and τ is at it smallest value (bottom right corner), we can see the highest mean number of infected.
Entropy 26 00298 g006
Figure 7. Panel of descriptive statistical analysis plots of three simulations: one chaotic (top row) and two periodic (middle and bottom row). The top row presents results from a simulation with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 . The middle illustrates a simulation with only the τ parameter changed (0.55 → 0.85). The bottom row of pictures illustrates a simulation with only the ϵ 1 parameter changed (0.25 → 0.31). In the first column, there is a Poincaré map of z ( x = 5 ) . In the second column, there is a Poincaré map of u ( x = 5 ) . In the third column, there is a Fast Fourier Transform of z ( x = 5 ) . In the fourth column, there is aFast Fourier Transform of u ( x = 5 ) . For the sake of clear visualization, we present only absolute amplitudes and half of the frequencies (to avoid a symmetric picture). The Poincaré maps have closed loops in the middle and bottom rows, which differs to the behavior in the top row. The spectral analysis in the top row is highly noisy compares to the analyses below it.
Figure 7. Panel of descriptive statistical analysis plots of three simulations: one chaotic (top row) and two periodic (middle and bottom row). The top row presents results from a simulation with parameters n = 10 , ϵ = 0.25 , and τ = 0.55 . The middle illustrates a simulation with only the τ parameter changed (0.55 → 0.85). The bottom row of pictures illustrates a simulation with only the ϵ 1 parameter changed (0.25 → 0.31). In the first column, there is a Poincaré map of z ( x = 5 ) . In the second column, there is a Poincaré map of u ( x = 5 ) . In the third column, there is a Fast Fourier Transform of z ( x = 5 ) . In the fourth column, there is aFast Fourier Transform of u ( x = 5 ) . For the sake of clear visualization, we present only absolute amplitudes and half of the frequencies (to avoid a symmetric picture). The Poincaré maps have closed loops in the middle and bottom rows, which differs to the behavior in the top row. The spectral analysis in the top row is highly noisy compares to the analyses below it.
Entropy 26 00298 g007
Figure 8. Panel of heatmaps of Autocorrelation and MLE measures according to various discretization resolutions. The first row consists of Autocorrelation heatmaps. The second row consists of MLE. In the columns, we see results for resolutions ranging from n = 4 to n = 20 . The missing tiles in the Autocorrelation heatmaps correspond to NaN values that mark stationary solutions. We can see that the number of confined parameter spaces with chaotic regimes varies between n values, and the size of area of unordered dynamics increases with increasing n. Chaotic space (positive MLEs) takes up to ∼ 90 % of the total space for n = 20 .
Figure 8. Panel of heatmaps of Autocorrelation and MLE measures according to various discretization resolutions. The first row consists of Autocorrelation heatmaps. The second row consists of MLE. In the columns, we see results for resolutions ranging from n = 4 to n = 20 . The missing tiles in the Autocorrelation heatmaps correspond to NaN values that mark stationary solutions. We can see that the number of confined parameter spaces with chaotic regimes varies between n values, and the size of area of unordered dynamics increases with increasing n. Chaotic space (positive MLEs) takes up to ∼ 90 % of the total space for n = 20 .
Entropy 26 00298 g008
Table 1. Fixed parameters and their values used in a grid search.
Table 1. Fixed parameters and their values used in a grid search.
Parametera γ β 0 β 1
Value0.60.10.110.225
Table 2. Varied parameters and their ranges used in a grid search.
Table 2. Varied parameters and their ranges used in a grid search.
Parametern ϵ τ
Initial values40.00.15
Final value10 10.41.05
Step size10.010.1
1 20 is appended to the sequence in the end.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Götz, T.; Krüger, T.; Niedzielewski, K.; Pestow, R.; Schäfer, M.; Schneider, J. Chaos in Opinion-Driven Disease Dynamics. Entropy 2024, 26, 298. https://0-doi-org.brum.beds.ac.uk/10.3390/e26040298

AMA Style

Götz T, Krüger T, Niedzielewski K, Pestow R, Schäfer M, Schneider J. Chaos in Opinion-Driven Disease Dynamics. Entropy. 2024; 26(4):298. https://0-doi-org.brum.beds.ac.uk/10.3390/e26040298

Chicago/Turabian Style

Götz, Thomas, Tyll Krüger, Karol Niedzielewski, Radomir Pestow, Moritz Schäfer, and Jan Schneider. 2024. "Chaos in Opinion-Driven Disease Dynamics" Entropy 26, no. 4: 298. https://0-doi-org.brum.beds.ac.uk/10.3390/e26040298

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop