Dynamic Opinion Model and Invasion Percolation

Editors' Suggestion

Dynamic Opinion Model and Invasion Percolation

Jia Shao, Shlomo Havlin, and H. Eugene Stanley

Phys. Rev. Lett. 103, 018701 – Published 1 July 2009

Abstract

We propose a “nonconsensus” opinion model that allows for stable coexistence of two opinions by forming clusters of agents holding the same opinion. We study this nonconsensus model on lattices, several model complex networks, and a real-life social network. We find that the model displays a phase transition behavior characterized by a large spanning cluster of nodes holding the same opinion appearing when the concentration of nodes holding the same opinion (even minority) is above a certain threshold. Because of the clustering (community support) of agents holding the same opinion, these clusters cannot be invaded by the other opinion (similar to incompressible fluids). Our extensive simulations show that the nonconsensus opinion model appears to belong to the same universality class as invasion percolation.

Received 25 March 2009

DOI:https://doi.org/10.1103/PhysRevLett.103.018701

Authors & Affiliations

Jia Shao¹, Shlomo Havlin², and H. Eugene Stanley¹

¹Center for Polymer Studies and Department of Physics, Boston University, Boston, Massachusetts 02215, USA
²Minerva Center and Department of Physics, Bar-Ilan University, 52900 Ramat-Gan, Israel

Comments & Replies

Comment on “Dynamic Opinion Model and Invasion Percolation”

Arsalan Sattari, Maya Paczuski, and Peter Grassberger

Phys. Rev. Lett. 109, 079801 (2012)

Shao, Havlin, and Stanley Reply:

Jia Shao, Shlomo Havlin, and H. Eugene Stanley

Phys. Rev. Lett. 109, 079802 (2012)

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 103, Iss. 1 — 3 July 2009

Reuse & Permissions

Access Options

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
Dynamics of the NCO model showing the approach to a stable state on a network with $N = 9$ nodes. For simplicity, we assume $p = 0$ . (a) At $t = 0$ , five nodes are randomly assigned to be $σ_{+}$ (filled circle). The remaining four nodes are assigned $σ_{-}$ (open circle). In the set comprising of node $A$ and its 4 neighboring nodes (dashed box), node $A$ is in a local minority opinion (2 $σ_{+}$ nodes and 3 $σ_{-}$ nodes), while the remaining nodes are not. At the end of simulation step $t = 0$ , node $A$ is converted into $σ_{-}$ . (b) At $t = 1$ , in the set of nodes comprising node $B$ and its 6 neighboring nodes (dotted box), node $B$ becomes in a local minority opinion (3 $σ_{+}$ nodes and 4 $σ_{-}$ nodes), while the remaining nodes are not. Node $B$ is converted into $σ_{-}$ at the end of simulation step $t = 1$ . (c) At $t = 2$ , the nine-nodes system reaches a stable state.Reuse & Permissions
Figure 2
Plot of the normalized size of the largest cluster $s_{1}$ (dotted line), the second largest cluster $s_{2}$ (full line), and the fraction of $σ_{-}$ nodes $F$ (dashed line) in the stable state as a function of $f$ for (a) a SF network with $λ = 2.5$ and $N = 10^{5}$ , (b) an ER network with $⟨ k ⟩ = 4$ and $N = 10^{5}$ , (c) a HEP network with $λ \approx 2.9$ , and (d) a HX lattice of size $1000 \times 1000$ [20]. Each curve represents an average over 100 realizations. The sharp increase of $s_{1}$ and the peak of $s_{2}$ at $f^{*}$ indicate a second-order phase transition. The insets in (a) and (d) show $g (S_{1}, N)$ [Eq. (1)] as a function of $N$ at the critical threshold $f^{*}$ (●), $f^{*} + Δ$ (◇), and $f^{*} - Δ$ (□), where $Δ = 0.01$ for SF and $Δ = 0.005$ for HX. At $f^{*}$ , $g (S_{1}, N)$ approaches a constant, indicating the size of the largest cluster $S_{1}$ proportional to $N^{ζ}$ at $f^{*}$ [Eq. (1)], which is another characteristic of a second-order phase transition.Reuse & Permissions
Figure 3
(a) The cumulative distribution function of cluster sizes $P (S^{'} > S)$ at $f^{*}$ for a NCO model on a SQ lattice with $N = 9 \times 10^{6}$ ( $3000 \times 3000$ ) and a SF network with $λ = 2.5$ and $N = 10^{5}$ . $P (S^{'} > S) \sim S^{1 - τ}$ for both SQ and SF, where $τ \approx 1.89$ (SQ) and $τ \approx 2.5$ (SF). $P (S^{'} > S)$ of the sizes of pores of TIP on SQ with $N = 9 \times 10^{6}$ is also shown, which also takes the form $P (S^{'} > S) \sim S^{1 - τ}$ with $τ \approx 1.90$ . Averages over 100 realizations are shown for all curves. (b) For the NCO model, at criticality, $R_{g} / S^{1 / d_{f}}$ as functions of $S$ for SQ and TR with $N = 9 \times 10^{6}$ . For the trial value $d_{f} = 1.84$ , $R_{g} / S^{1 / d_{f}}$ approaches a constant. For $d_{f} = 1.896$ (regular percolation fractal dimension), $R_{g} / S^{1 / d_{f}}$ is an increasing function of $S$ . (c) For the NCO model, at criticality, $ℓ / S^{1 / d_{ℓ}}$ as a function of $S$ for SQ and TR with $N = 9 \times 10^{6}$ . For the trial value $d_{ℓ} = 1.54$ , $ℓ / S^{1 / d_{ℓ}}$ approaches a constant. For $d_{ℓ} = 1.678$ (regular percolation fractal dimension), $ℓ / S^{1 / d_{ℓ}}$ is an increasing function of $S$ . In (b) and (c), each curve is averaged over 1000 realizations.Reuse & Permissions
Figure 4
(a) The average degree $⟨ k (f) ⟩$ of $σ_{-}$ nodes in the stable state as a function of $f$ for a SF network with $λ = 2.5$ and $N = 10^{5}$ and the HEP network. It is seen that, for $f < 0.5$ , $⟨ k (f) ⟩$ is significantly smaller than that for $f > 0.5$ , since the high degree nodes join the majority opinion [18]. (b) The cumulative degree distribution $P_{f} (k^{'} > k)$ of $σ_{-}$ nodes in the stable state for a SF network with $λ = 2.5$ and $N = 10^{5}$ . $P_{f} (k^{'} > k)$ for $f = 0.450$ and 0.5 are plotted. Notice that $f^{*} \approx 0.450$ . The absence of high degree $σ_{-}$ nodes in a stable state for $f < 0.5$ is again confirmed.Reuse & Permissions

Physical Review Letters