No Arabic abstract
Small-world networks describe many important practical systems among which neural networks consisting of excitable nodes are the most typical ones. In this paper we study self-sustained oscillations of target waves in excitable small-world networks. A novel dominant phase-advanced driving (DPAD) method, which is generally applicable for analyzing all oscillatory complex networks consisting of nonoscillatory nodes, is proposed to reveal the self-organized structures supporting this type of oscillations. The DPAD method explicitly explores the oscillation sources and wave propagation paths of the systems, which are otherwise deeply hidden in the complicated patterns of randomly distributed target groups. Based on the understanding of the self-organized structure, the oscillatory patterns can be controlled with extremely high efficiency.
The transition to turbulence via spatiotemporal intermittency is investigated in the context of coupled maps defined on small-world networks. The local dynamics is given by the Chate-Manneville minimal map previously used in studies of spatiotemporal intermittency in ordered lattices. The critical boundary separating laminar and turbulent regimes is calculated on the parameter space of the system, given by the coupling strength and the rewiring probability of the network. Windows of relaminarization are present in some regions of the parameter space. New features arise in small-world networks; for instance, the character of the transition to turbulence changes from second order to a first order phase transition at some critical value of the rewiring probability. A linear relation characterizing the change in the order of the phase transition is found. The global quantity used as order parameter for the transition also exhibits nontrivial collective behavior for some values of the parameters. These models may describe several processes occurring in nonuniform media where the degree of disorder can be continuously varied through a parameter.
We study the strategy to optimally maximize the dynamic range of excitable networks by removing the minimal number of links. A network of excitable elements can distinguish a broad range of stimulus intensities and has its dynamic range maximized at criticality. In this study, we formulate the activation propagation in excitable networks as a message passing process in which the critical state is reached when the largest eigenvalue of the weighted non-backtracking (WNB) matrix is exactly one. By considering the impact of single link removal on the largest eigenvalue, we develop an efficient algorithm that aims to identify the optimal set of links whose removal will drive the system to the critical state. Comparisons with other competing heuristics on both synthetic and real-world networks indicate that the proposed method can maximize the dynamic range by removing the smallest number of links, and at the same time maintain the largest size of the giant connected component.
Population bursts in a large ensemble of coupled elements result from the interplay between the local excitable properties of the nodes and the global network topology. Here collective excitability and self-sustained bursting oscillations are shown to spontaneously emerge in adaptive networks of globally coupled non-excitable units. The ingredients to observe collective excitability are the coexistence of states with different degree of synchronizaton joined to a global feedback acting, on a slow timescale, against the synchronization (desynchronization) of the oscillators. These regimes are illustrated for two paradigmatic classes of coupled rotators: namely, the Kuramoto model with and without inertia. For the bimodal Kuramoto model we analytically show that the macroscopic evolution originates from the existence of a critical manifold organizing the fast collective dynamics on a slow timescale. Our results provide evidence that adaptation can induce excitability by maintaining a network permanently out-of-equilibrium.
Networks of excitable elements are widely used to model real-world biological and social systems. The dynamic range of an excitable network quantifies the range of stimulus intensities that can be robustly distinguished by the network response, and is maximized at the critical state. In this study, we examine the impacts of backtracking activation on system criticality in excitable networks consisting of both excitatory and inhibitory units. We find that, for dynamics with refractory states that prohibit backtracking activation, the critical state occurs when the largest eigenvalue of the weighted non-backtracking (WNB) matrix for excitatory units, $lambda^E_{NB}$, is close to one, regardless of the strength of inhibition. In contrast, for dynamics without refractory state in which backtracking activation is allowed, the strength of inhibition affects the critical condition through suppression of backtracking activation. As inhibitory strength increases, backtracking activation is gradually suppressed. Accordingly, the system shifts continuously along a continuum between two extreme regimes -- from one where the criticality is determined by the largest eigenvalue of the weighted adjacency matrix for excitatory units, $lambda^E_W$, to the other where the critical state is reached when $lambda_{NB}^E$ is close to one. For systems in between, we find that $lambda^E_{NB}<1$ and $lambda^E_W>1$ at the critical state. These findings, confirmed by numerical simulations using both random and synthetic neural networks, indicate that backtracking activation impacts the criticality of excitable networks.
Direct simulations of two-dimensional plane channel flow of a viscoelastic fluid at Reynolds number Re = 3000 reveal the existence of a family of attractors whose structure closely resembles the linear Tollmien-Schlichting (TS) mode, and in particular exhibits strongly localized stress fluctuations at the critical layer position of the TS mode. At the parameter values chosen, this solution branch is not connected to the nonlinear TS solution branch found for Newtonian flow, and thus represents a solution family that is nonlinearly self-sustained by viscoelasticity. The ratio between stress and velocity fluctuations is in quantitative agreement for the attractor and the linear TS mode, and increases strongly with Weissenberg number, Wi. For the latter, there is a transition in the scaling of this ratio as Wi increases, and the Wi at which the nonlinear solution family comes into existence is just above this transition. Finally, evidence indicates that this branch is connected through an unstable solution branch to two-dimensional elastoinertial turbulence (EIT). These results suggest that, in the parameter range considered here, the bypass transition leading to EIT is mediated by nonlinear amplification and self-sustenance of perturbations that excite the Tollmien-Schlichting mode.