Collective behaviors of animal groups may stem from visual lateralization-Tending to obtain information through one eye.

Animal groups exhibit various captivating movement patterns, which manifest as intricate interactions among group members. Several models have been proposed to elucidate collective behaviors in animal groups. These models achieve a certain degree of efficacy; however, inconsistent experimental findings suggest insufficient accuracy. Experiments have shown that some organisms employ a single information channel and visual lateralization to glean knowledge from other individuals in collective movements. In this study, we consider individuals' visual lateralization and a single information channel and develop a self-propelled particle model to describe the collective behavior of large groups. The results suggest that homogeneous visual lateralization gives the group a strong sense of cohesiveness, thereby enabling diverse collective behaviors. As the overlapping field grows, the cohesiveness gradually dissipates. Inconsistent visual lateralization among group members can reduce the cohesiveness of the group, and when there is a high degree of heterogeneity in visual lateralization, the group loses their cohesiveness. This study also examines the influence of visual lateralization heterogeneity on specific formations, and the results indicate that the directional migration formation is responsive to such heterogeneity. We propose an information network to portray the transmission of information within groups, which explains the cohesiveness of groups and the sensitivity of the directional migration formation.

Composite spiral waves in discrete-time systems.

Spiral waves are a type of typical pattern in open reaction-diffusion systems far from thermodynamic equilibrium. The study of spiral waves has attracted great interest because of its nonlinear characteristics and extensive applications. However, the study of spiral waves has been confined to continuous-time systems, while spiral waves in discrete-time systems have been rarely reported. In recent years, discrete-time models have been widely studied in ecology because of their appropriateness to systems with nonoverlapping generations and other factors. Therefore, spiral waves in discrete-time systems need to be studied. Here, we investigated a novel type of spiral wave called a composite spiral wave in a discrete-time predator-pest model, and we revealed the formation mechanism. To explain the observed phenomena, we defined and quantified a move state effect of multiperiod states caused by the coupling of adjacent stable multiperiod orbits, which is strictly consistent with the numerical results. The other move state effect is caused by an unstable focus, which is the state of the local points at the spiral center. The combined effect of these two influences can lead to rich dynamical behaviors of spiral waves, and the specific structure of the composite spiral waves is the result of the competition of the two effects in opposite directions. Our findings shed light on the dynamics of spiral waves in discrete-time systems, and they may guide the prediction and control of pests in deciduous forests.

Irregular spots on body surfaces of vertebrates induced by supercritical pitchfork bifurcations.

The classical Turing mechanism containing a long-range inhibition and a short-range self-enhancement provides a type of explanation for the formation of patterns on body surfaces of some vertebrates, e.g., zebras, giraffes, and cheetahs. For other type of patterns (irregular spots) on body surfaces of some vertebrates, e.g., loaches, finless eels, and dalmatian dogs, the classical Turing mechanism no longer applies. Here, we propose a mechanism, i.e., the supercritical pitchfork bifurcation, which may explain the formation of this type of irregular spots, and present a method to quantify the similarity of such patterns. We assume that, under certain conditions, the only stable state of "morphogen" loses its stability and transitions to two newly generated stable states with the influence of external noise, thus producing such ruleless piebald patterns in space. The difference between the competitiveness of these two states may affect the resulting pattern. Moreover, we propose a mathematical model based on this conjecture and obtain this type of irregular patterns by numerical simulation. Furthermore, we also study the influence of parameters in the model on pattern structures and obtain the corresponding pattern structures of some vertebrates in nature, which verifies our conjecture.

Non-Markovian majority-vote model.

Non-Markovian dynamics pervades human activity and social networks and it induces memory effects and burstiness in a wide range of processes including interevent time distributions, duration of interactions in temporal networks, and human mobility. Here, we propose a non-Markovian majority-vote model (NMMV) that introduces non-Markovian effects in the standard (Markovian) majority-vote model (SMV). The SMV model is one of the simplest two-state stochastic models for studying opinion dynamics, and displays a continuous order-disorder phase transition at a critical noise. In the NMMV model we assume that the probability that an agent changes state is not only dependent on the majority state of his neighbors but it also depends on his age, i.e., how long the agent has been in his current state. The NMMV model has two regimes: the aging regime implies that the probability that an agent changes state is decreasing with his age, while in the antiaging regime the probability that an agent changes state is increasing with his age. Interestingly, we find that the critical noise at which we observe the order-disorder phase transition is a nonmonotonic function of the rate ß of the aging (antiaging) process. In particular the critical noise in the aging regime displays a maximum as a function of ß while in the antiaging regime displays a minimum. This implies that the aging/antiaging dynamics can retard/anticipate the transition and that there is an optimal rate ß for maximally perturbing the value of the critical noise. The analytical results obtained in the framework of the heterogeneous mean-field approach are validated by extensive numerical simulations on a large variety of network topologies.

Hybrid multiscale coarse-graining for dynamics on complex networks.

We propose a hybrid multiscale coarse-grained (HMCG) method which combines a fine Monte Carlo (MC) simulation on the part of nodes of interest with a more coarse Langevin dynamics on the rest part. We demonstrate the validity of our method by analyzing the equilibrium Ising model and the nonequilibrium susceptible-infected-susceptible model. It is found that HMCG not only works very well in reproducing the phase transitions and critical phenomena of the microscopic models, but also accelerates the evaluation of dynamics with significant computational savings compared to microscopic MC simulations directly for the whole networks. The proposed method is general and can be applied to a wide variety of networked systems just adopting appropriate microscopic simulation methods and coarse graining approaches.

Large deviation induced phase switch in an inertial majority-vote model.

We theoretically study noise-induced phase switch phenomena in an inertial majority-vote (IMV) model introduced in a recent paper [Chen et al., Phys. Rev. E 95, 042304 (2017)]. The IMV model generates a strong hysteresis behavior as the noise intensity f goes forward and backward, a main characteristic of a first-order phase transition, in contrast to a second-order phase transition in the original MV model. Using the Wentzel-Kramers-Brillouin approximation for the master equation, we reduce the problem to finding the zero-energy trajectories in an effective Hamiltonian system, and the mean switching time depends exponentially on the associated action and the number of particles N. Within the hysteresis region, we find that the actions, along the optimal forward switching path from the ordered phase (OP) to disordered phase (DP) and its backward path show distinct variation trends with f, and intersect at f = fc that determines the coexisting line of the OP and DP. This results in a nonmonotonic dependence of the mean switching time between two symmetric OPs on f, with a minimum at fc for sufficiently large N. Finally, the theoretical results are validated by Monte Carlo simulations.

First-order phase transition in a majority-vote model with inertia.

We generalize the original majority-vote model by incorporating inertia into the microscopic dynamics of the spin flipping, where the spin-flip probability of any individual depends not only on the states of its neighbors, but also on its own state. Surprisingly, the order-disorder phase transition is changed from a usual continuous or second-order type to a discontinuous or first-order one when the inertia is above an appropriate level. A central feature of such an explosive transition is a strong hysteresis behavior as noise intensity goes forward and backward. Within the hysteresis region, a disordered phase and two symmetric ordered phases are coexisting and transition rates between these phases are numerically calculated by a rare-event sampling method. A mean-field theory is developed to analytically reveal the property of this phase transition.

Critical noise of majority-vote model on complex networks.

The majority-vote model with noise is one of the simplest nonequilibrium statistical model that has been extensively studied in the context of complex networks. However, the relationship between the critical noise where the order-disorder phase transition takes place and the topology of the underlying networks is still lacking. In this paper, we use the heterogeneous mean-field theory to derive the rate equation for governing the model's dynamics that can analytically determine the critical noise f(c) in the limit of infinite network size N→∞. The result shows that f(c) depends on the ratio of 〈k〉 to 〈k(3/2)〉, where 〈k〉 and 〈k(3/2)〉 are the average degree and the 3/2 order moment of degree distribution, respectively. Furthermore, we consider the finite-size effect where the stochastic fluctuation should be involved. To the end, we derive the Langevin equation and obtain the potential of the corresponding Fokker-Planck equation. This allows us to calculate the effective critical noise f(c)(N) at which the susceptibility is maximal in finite-size networks. We find that the f(c)-f(c)(N) decays with N in a power-law way and vanishes for N→∞. All the theoretical results are confirmed by performing the extensive Monte Carlo simulations in random k-regular networks, Erdös-Rényi random networks, and scale-free networks.

Mobility and density induced amplitude death in metapopulation networks of coupled oscillators.

We investigate the effects of mobility and density on the amplitude death of coupled Landau-Stuart oscillators and Brusselators in metapopulation networks, wherein each node represents a subpopulation occupied any number of mobile individuals. By numerical simulations in scale-free topology, we find that the systems undergo phase transitions from incoherent state to amplitude death, and then to frequency synchronization with increasing the mobility rate or density of oscillators. Especially, there exists an extent of intermediate mobility rate and density that can lead to global oscillator death. Furthermore, we show that such nontrivial phenomena are robust to diverse network topologies. Our findings may invoke further efforts and attentions to explore the underlying mechanism of collective behaviors in coupled metapopulation systems.

Explosive synchronization transitions in complex neural networks.

It has been recently reported that explosive synchronization transitions can take place in networks of phase oscillators [Gómez-Gardeñes et al. Phys. Rev. Lett. 106, 128701 (2011)] and chaotic oscillators [Leyva et al. Phys. Rev. Lett. 108, 168702 (2012)]. Here, we investigate the effect of a microscopic correlation between the dynamics and the interacting topology of coupled FitzHugh-Nagumo oscillators on phase synchronization transition in Barabási-Albert (BA) scale-free networks and Erdös-Rényi (ER) random networks. We show that, if natural frequencies of the oscillations are positively correlated with node degrees and the width of the frequency distribution is larger than a threshold value, a strong hysteresis loop arises in the synchronization diagram of BA networks, indicating the evidence of an explosive transition towards synchronization of relaxation oscillators system. In contrast to the results in BA networks, in more homogeneous ER networks, the synchronization transition is always of continuous type regardless of the width of the frequency distribution. Moreover, we consider the effect of degree-mixing patterns on the nature of the synchronization transition, and find that the degree assortativity is unfavorable for the occurrence of such an explosive transition.

Nucleation pathways on complex networks.

Identifying nucleation pathway is important for understanding the kinetics of first-order phase transitions in natural systems. In the present work, we study nucleation pathway of the Ising model in homogeneous and heterogeneous networks using the forward flux sampling method, and find that the nucleation processes represent distinct features along pathways for different network topologies. For homogeneous networks, there always exists a dominant nucleating cluster to which relatively small clusters are attached gradually to form the critical nucleus. For heterogeneous ones, many small isolated nucleating clusters emerge at the early stage of the nucleation process, until suddenly they form the critical nucleus through a sharp merging process. Moreover, we also compare the nucleation pathways for different degree-mixing networks. By analyzing the properties of the nucleating clusters along the pathway, we show that the main reason behind the different routes is the heterogeneous character of the underlying networks.

Strategy to suppress epidemic explosion in heterogeneous metapopulation networks.

We propose an efficient strategy to suppress epidemic explosion in heterogeneous metapopulation networks, wherein each node represents a subpopulation with any number of individuals and is assigned a curing rate that is proportional to kα with the node degree k and an adjustable parameter α. We perform stochastic simulations of the dynamical reaction-diffusion processes associated with the susceptible-infected-susceptible model in scale-free networks. We find that the epidemic threshold reaches a maximum when α is tuned at αopt≃1.3. This nontrivial phenomenon is robust to the change of the network size and the average degree. In addition, we carry out a mean field analysis to further validate our scheme, which also demonstrates that epidemic explosion follows different routes for α larger or less than αopt. Our work suggests that in order to efficiently suppress epidemic spreading on heterogeneous complex networks, subpopulations with higher degrees should be allocated more resources than just being linearly dependent on k.

Coarse-grained Monte Carlo simulations of the phase transition of the Potts model on weighted networks.

Developing an effective coarse-grained (CG) approach is a promising way for studying dynamics on large size networks. In the present work, we have proposed a strength-based CG (s-CG) method to study critical phenomena of the Potts model on weighted complex networks. By merging nodes with close strengths together, the original network is reduced to a CG network with much smaller size, on which the CG Hamiltonian can be well defined. In particular, we make an error analysis and show that our s-CG approach satisfies the condition of statistical consistency, which demands that the equilibrium probability distribution of the CG model matches that of the microscopic counterpart. Extensive numerical simulations are performed on scale-free networks and random networks, without or with strength correlation, showing that this s-CG approach works very well in reproducing the phase diagrams, fluctuations, and finite-size effects of the microscopic model, while the d-CG approach proposed in our recent work [Phys. Rev. E 82, 011107 (2010)] does not.

Nucleation in scale-free networks.

We have studied nucleation dynamics of the Ising model in scale-free networks whose degree distribution follows a power law with the exponent γ by using the forward flux sampling method and focusing on how the network topology would influence the nucleation rate and pathway. For homogeneous nucleation, the new phase clusters grow from those nodes with smaller degree, while the cluster sizes follow a power-law distribution. Interestingly, we find that the nucleation rate R{Hom} decays exponentially with network size and, accordingly, the critical nucleus size increases linearly with network size, implying that homogeneous nucleation is not relevant in the thermodynamic limit. These observations are robust to the change of γ and are also present in random networks. In addition, we have also studied the dynamics of heterogeneous nucleation, wherein w impurities are initially added either to randomly selected nodes or to targeted ones with the largest degrees. We find that targeted impurities can enhance the nucleation rate R{Het} much more sharply than random ones. Moreover, ln(R{Het}/R{Hom}) scales as w{(γ-2)/(γ-1)} and w for targeted and random impurities, respectively. A simple mean-field analysis is also present to qualitatively illustrate the above simulation results.