American Institute of Mathematical Sciences

March  2017, 10(1): 1-32. doi: 10.3934/krm.2017001

Opinion dynamics over complex networks: Kinetic modelling and numerical methods

 1 TU München, Faculty of Mathematics, Boltzmannstra$\mathcal{B}$ e 3, D-85748, Garching (München), Germany 2 University of Ferrara, Department of Mathematics and Computer Science, Via N. Machiavelli 35,44121, Ferrara, Italy

Received  March 2016 Revised  May 2016 Published  November 2016

In this paper we consider the modeling of opinion dynamics over time dependent large scale networks. A kinetic description of the agents' distribution over the evolving network is considered which combines an opinion update based on binary interactions between agents with a dynamic creation and removal process of new connections. The number of connections of each agent influences the spreading of opinions in the network but also the way connections are created is influenced by the agents' opinion. The evolution of the network of connections is studied by showing that its asymptotic behavior is consistent both with Poisson distributions and truncated power-laws. In order to study the large time behavior of the opinion dynamics a mean field description is derived which allows to compute exact stationary solutions in some simplified situations. Numerical methods which are capable to describe correctly the large time behavior of the system are also introduced and discussed. Finally, several numerical examples showing the influence of the agents' number of connections in the opinion dynamics are reported.

Citation: Giacomo Albi, Lorenzo Pareschi, Mattia Zanella. Opinion dynamics over complex networks: Kinetic modelling and numerical methods. Kinetic & Related Models, 2017, 10 (1) : 1-32. doi: 10.3934/krm.2017001
References:

show all references

References:
Stationary states of (21) with relaxation coefficients ${V_r}={V_a}=1$, mean density of connectivity $\gamma=30$, ${{c}_{\max }}=1500$ and several values of the attraction parameters $\alpha$, and having fixed $\beta = 0$. Left: convergence toward the Poisson distribution for big values of $\alpha$. Right: convergence toward a power-law distribution in the limit $\alpha\rightarrow 0$, we indicated with $p_{\infty}^{(-k)}, k=1,2,3$ the $\alpha-$dependent stationary solutions for $\alpha=10^{-1},10^{-2},10^{-3}$, respectively.
Stationary solutions of type $f_\infty(w,c) = g_\infty(w)p_\infty(c)$, where $g_\infty(w)$ is given by (51) with $\kappa =1$, $m_w = 0$, $\sigma^2 = 0.05$ and $p_\infty(c)$ defined by (25), with ${V_r}={V_a}=1$, $\gamma=30$ and $\alpha = 10$ on the left and $\alpha = 0.1$ on the right.
Stationary solutions captured via Monte Carlo simulations, with $N_s=2\times10^4$ samples. Parameters of the model are chosen as follows $\sigma^2 = 0.05$, ${V_r}={V_a}=1$, $\beta = 0$, $\alpha = 10$ on the right hand side and $\alpha =0.1$ on the left hand side.
Test $\#1$. One-dimensional setting: on the left, convergence of (58) to the stationary solution (51), of the Fokker-Planck equation, for decreasing values of the parameter $\varepsilon$, $g^N_0$ represents the initial distribution. On the right, convergence of the Monte-Carlo (55) to the reference solution (24) for increasing values of the the number of samples $N_s$.
Test $\#1.$ One-dimensional setting: on the left, the solution of the Chang-Cooper type scheme with the flux (70) is indicated with $g^N_T$ and compared with the stationary solution (51), also the initial data $g^N_0$ (72) is reported. On the right we report the decay of the $L^1$ relative error (73) for different choices of the quadrature rule, mid-point rule (67) and Milne's rule (respectively of $2^{nd}$ and $4^{th}$ order).
Test $\#2.$ From left to right and from the top to the bottom: evolution of the density $f(w,c,t)$ at different time steps. The plot $(a)$ represents the initial data $f_0(w,c)$ (74) and plot $(d)$ the stationary solution. On the plane $(z,c)$ we depict with a blu line the marginal distribution $p(c,t)$ of the solution at time $t$, with red line we represent the reference marginal distribution of the stationary solution.
Test $\#2.$ Decay of the $L^1$ relative error with respect to the stationary solution (51). On the left, fixed characteristic rates $V=\{10^3,10^4,10^5\}$, on the right, variable characteristic rates defined as in (76) with $U =\{10^3,10^4,10^5\}$. In both cases for increasing values of the characteristic rate $V$ and $U$ the stationary state is reached faster.
Test $\#2.$ Evolution at time $t = 1$ of the initial data $f_0(w, c)$ (74) as isoline plot. On the left in the case of constant characteristic rate on the right variable characteristic rates defined as in (76). The right plot shows that for lower opinion's density the evolution along the connection is faster and slower where the opinions are more concentrated.
Test $\#3.$ From left to right and from the first row to the second row, evolution of the initial data (78) in time frame $[0,T]$ with $T = 2$. The evolution shows how a small portion of density with high connectivity can bias the majority of the population towards their position. (Note: The density is scaled according to the marginal distribution $\rho(c,t)$ in order to better show its evolution, the actual marginal density $\rho(c,t)$ is depicted in the background, scaled by a factor 10).
Test $\#3.$ On the left-hand side final and initial state of the marginal distribution $g(w,t)$ of the opinion, the green line represents the evolution of the average opinion $\bar{m}(t)$, the red and blue dashed lines represent respectively the opinions $\bar{w}_L=0.75$ and $\bar{w}_F=-0.5$, which are the two leading opinions of the initial data (78).
Test $\#4.$ Evolution of the Fokker-Planck model (41) where the interaction are described by (80) with $\Delta=0.25$, in the time frame $[0,T]$ with $T = 100$. The evolution shows the emergence of three main opinion clusters, which are not affected by the connectivity variable. (Note: In order to better show its evolution, we represent the solution as $\log(f(w,c,t)+\epsilon)$, with $\epsilon = 0.001$.)
Test $\#4.$ Evolution of the solution of the Fokker-Planck model (41), where the interaction are described by (79) with $\Delta(c)=d_0c/{c_{\max }}$ and $d_0=1.01$, in the time frame $[0,T]$ with $T = 100$. The choice of $\Delta(c)$ reflects in the heterogeneous emergence of clusters with respect to the connectivity level: for higher level of connectivity consensus is reached, instead for lower levels of connectivity multiple opinion clusters are present. (Note: In order to better show its evolution, we represent the solution as $\log(f(w,c,t)+\epsilon)$, with $\epsilon = 0.001$.
Parameters in the various test cases
 Test $\sigma^2$ $\sigma^2_F$ $\sigma^2_L$ ${{c}_{\max }}$ ${V_r}$ ${V_a}$ $\gamma_0$ $\alpha$ $\beta$ #1 $5\times10^{-2}$ $6\times10^{-2}$ $-$ $250$ $1$ 1 $30$ $1\times10^{-1}$ $0$ #2 $5\times10^{-2}$ $6\times10^{-2}$ $-$ $250$ $-$ $-$ $30$ $1\times10^{-1}$ $0$ #3 $5\times10^{-3}$ $4\times10^{-2}$ $2.5\times10^{-2}$ 250 $1$ $1$ 30 $1\times10^{-4}$ $0$ #4 $1\times10^{-3}$ $-$ $-$ 250 $1$ $1$ 30 $1\times10^{-1}$ $0$
 Test $\sigma^2$ $\sigma^2_F$ $\sigma^2_L$ ${{c}_{\max }}$ ${V_r}$ ${V_a}$ $\gamma_0$ $\alpha$ $\beta$ #1 $5\times10^{-2}$ $6\times10^{-2}$ $-$ $250$ $1$ 1 $30$ $1\times10^{-1}$ $0$ #2 $5\times10^{-2}$ $6\times10^{-2}$ $-$ $250$ $-$ $-$ $30$ $1\times10^{-1}$ $0$ #3 $5\times10^{-3}$ $4\times10^{-2}$ $2.5\times10^{-2}$ 250 $1$ $1$ 30 $1\times10^{-4}$ $0$ #4 $1\times10^{-3}$ $-$ $-$ 250 $1$ $1$ 30 $1\times10^{-1}$ $0$
 [1] Xuefei He, Kun Wang, Liwei Xu. Efficient finite difference methods for the nonlinear Helmholtz equation in Kerr medium. Electronic Research Archive, 2020, 28 (4) : 1503-1528. doi: 10.3934/era.2020079 [2] Fabio Camilli, Giulia Cavagnari, Raul De Maio, Benedetto Piccoli. Superposition principle and schemes for measure differential equations. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2020050 [3] Anton A. Kutsenko. Isomorphism between one-Dimensional and multidimensional finite difference operators. Communications on Pure & Applied Analysis, 2021, 20 (1) : 359-368. doi: 10.3934/cpaa.2020270 [4] Xin Guo, Lei Shi. Preface of the special issue on analysis in data science: Methods and applications. Mathematical Foundations of Computing, 2020, 3 (4) : i-ii. doi: 10.3934/mfc.2020026 [5] Juan Pablo Pinasco, Mauro Rodriguez Cartabia, Nicolas Saintier. Evolutionary game theory in mixed strategies: From microscopic interactions to kinetic equations. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2020051 [6] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [7] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [8] Min Xi, Wenyu Sun, Jun Chen. Survey of derivative-free optimization. Numerical Algebra, Control & Optimization, 2020, 10 (4) : 537-555. doi: 10.3934/naco.2020050 [9] Hoang The Tuan. On the asymptotic behavior of solutions to time-fractional elliptic equations driven by a multiplicative white noise. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020318 [10] Lin Shi, Xuemin Wang, Dingshi Li. Limiting behavior of non-autonomous stochastic reaction-diffusion equations with colored noise on unbounded thin domains. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5367-5386. doi: 10.3934/cpaa.2020242 [11] D. R. Michiel Renger, Johannes Zimmer. Orthogonality of fluxes in general nonlinear reaction networks. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 205-217. doi: 10.3934/dcdss.2020346 [12] Bernold Fiedler. Global Hopf bifurcation in networks with fast feedback cycles. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 177-203. doi: 10.3934/dcdss.2020344 [13] Predrag S. Stanimirović, Branislav Ivanov, Haifeng Ma, Dijana Mosić. A survey of gradient methods for solving nonlinear optimization. Electronic Research Archive, 2020, 28 (4) : 1573-1624. doi: 10.3934/era.2020115 [14] Neng Zhu, Zhengrong Liu, Fang Wang, Kun Zhao. Asymptotic dynamics of a system of conservation laws from chemotaxis. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 813-847. doi: 10.3934/dcds.2020301 [15] Emre Esentürk, Juan Velazquez. Large time behavior of exchange-driven growth. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 747-775. doi: 10.3934/dcds.2020299 [16] Yue Feng, Yujie Liu, Ruishu Wang, Shangyou Zhang. A conforming discontinuous Galerkin finite element method on rectangular partitions. Electronic Research Archive, , () : -. doi: 10.3934/era.2020120 [17] Shun Zhang, Jianlin Jiang, Su Zhang, Yibing Lv, Yuzhen Guo. ADMM-type methods for generalized multi-facility Weber problem. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020171 [18] Dan Zhu, Rosemary A. Renaut, Hongwei Li, Tianyou Liu. Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020076 [19] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [20] Wei Feng, Michael Freeze, Xin Lu. On competition models under allee effect: Asymptotic behavior and traveling waves. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5609-5626. doi: 10.3934/cpaa.2020256

2019 Impact Factor: 1.311