We report finite-size numerical investigations and mean-field analysis of a Kuramoto model with inertia for fully coupled and diluted systems. In particular, we examine, for a gaussian distribution of the frequencies, the transition from incoherence to coherence for increasingly large system size and inertia. For sufficiently large inertia the transition is hysteretic, and within the hysteretic region clusters of locked oscillators of various sizes and different levels of synchronization coexist. A modification of the mean-field theory developed by Tanaka, Lichtenberg, and Oishi [Physica D 100, 279 (1997)] allows us to derive the synchronization curve associated to each of these clusters. We have also investigated numerically the limits of existence of the coherent and of the incoherent solutions. The minimal coupling required to observe the coherent state is largely independent of the system size, and it saturates to a constant value already for moderately large inertia values. The incoherent state is observable up to a critical coupling whose value saturates for large inertia and for finite system sizes, while in the thermodinamic limit this critical value diverges proportionally to the mass. By increasing the inertia the transition becomes more complex, and the synchronization occurs via the emergence of clusters of whirling oscillators. The presence of these groups of coherently drifting oscillators induces oscillations in the order parameter. We have shown that the transition remains hysteretic even for randomly diluted networks up to a level of connectivity corresponding to a few links per oscillator. Finally, an application to the Italian high-voltage power grid is reported, which reveals the emergence of quasiperiodic oscillations in the order parameter due to the simultaneous presence of many competing whirling clusters.

Abstract. -We study the dynamics of two symmetrically coupled populations of identical leaky integrate-and-fire neurons characterized by an excitatory coupling. Upon varying the coupling strength, we find symmetry-breaking transitions that lead to the onset of various chimera states as well as to a new regime, where the two populations are characterized by a different degree of synchronization. Symmetric collective states of increasing dynamical complexity are also observed. The computation of the the finite-amplitude Lyapunov exponent allows us to establish the chaoticity of the (collective) dynamics in a finite region of the phase plane. The further numerical study of the standard Lyapunov spectrum reveals the presence of several positive exponents, indicating that the microscopic dynamics is high-dimensional.

Information transmission in the human brain is a fundamentally dynamic network process. In partial epilepsy, this process is perturbed and highly synchronous seizures originate in a local network, the so-called epileptogenic zone (EZ), before recruiting other close or distant brain regions. We studied patient-specific brain network models of 15 drug-resistant epilepsy patients with implanted stereotactic electroencephalography (SEEG) electrodes. Each personalized brain model was derived from structural data of magnetic resonance imaging (MRI) and diffusion tensor weighted imaging (DTI), comprising 88 nodes equipped with region specific neural mass models capable of demonstrating a range of epileptiform discharges. Each patient’s virtual brain was further personalized through the integration of the clinically hypothesized EZ. Subsequent simulations and connectivity modulations were performed and uncovered a finite repertoire of seizure propagation patterns. Across patients, we found that (i) patient-specific network connectivity is predictive for the subsequent seizure propagation pattern; (ii) seizure propagation is characterized by a systematic sequence of brain states; (iii) propagation can be controlled by an optimal intervention on the connectivity matrix; (iv) the degree of invasiveness can be significantly reduced via the proposed seizure control as compared to traditional resective surgery. To stop seizures, neurosurgeons typically resect the EZ completely. We showed that stability analysis of the network dynamics, employing structural and dynamical information, estimates reliably the spatiotemporal properties of seizure propagation. This suggests novel less invasive paradigms of surgical interventions to treat and manage partial epilepsy.

The emergence of dynamical abrupt transitions in the macroscopic state of a system is currently a subject of the utmost interest. Given a set of phase oscillators networking with a generic wiring of connections and displaying a generic frequency distribution, we show how combining dynamical local information on frequency mismatches and global information on the graph topology suggests a judicious and yet practical weighting procedure which is able to induce and enhance explosive, irreversible, transitions to synchronization. We report extensive numerical and analytical evidence of the validity and scalability of such a procedure for different initial frequency distributions, for both homogeneous and heterogeneous networks, as well as for both linear and nonlinear weighting functions. We furthermore report on the possibility of parametrically controlling the width and extent of the hysteretic region of coexistence of the unsynchronized and synchronized states.

We investigate the onset of collective oscillations in a excitatory pulse-coupled network of leaky integrateand-fire neurons in the presence of quenched and annealed disorder. We find that the disorder induces a weak form of chaos that is analogous to that arising in the Kuramoto model for a finite number N of oscillators ͓O. V. Popovych et al., Phys. Rev. E 71 065201͑R͒ ͑2005͔͒. In fact, the maximum Lyapunov exponent turns out to scale to zero for N → ϱ, with an exponent that is different for the two types of disorder. In the thermodynamic limit, the random-network dynamics reduces to that of a fully homogeneous system with a suitably scaled coupling strength. Moreover, we show that the Lyapunov spectrum of the periodically collective state scales to zero as 1 / N 2 , analogously to the scaling found for the "splay state."

Two symmetrically coupled populations of N oscillators with inertia m display chaotic solutions with broken symmetry similar to experimental observations with mechanical pendulums. In particular, we report evidence of intermittent chaotic chimeras, where one population is synchronized and the other jumps erratically between laminar and turbulent phases. These states have finite lifetimes diverging as a power law with N and m. Lyapunov analyses reveal chaotic properties in quantitative agreement with theoretical predictions for globally coupled dissipative systems.

The microscopic and macroscopic dynamics of random networks is investigated in the strong-dilution limit (i.e., for sparse networks). By simulating chaotic maps, Stuart-Landau oscillators, and leaky integrate-and-fire neurons, we show that a finite connectivity (of the order of a few tens) is able to sustain a nontrivial collective dynamics even in the thermodynamic limit. Although the network structure implies a nonadditive dynamics, the microscopic evolution is extensive (i.e., the number of active degrees of freedom is proportional to the number of network elements).

A synaptic theory of Working Memory (WM) has been developed in the last decade as a possible alternative to the persistent spiking paradigm. In this context, we have developed a neural mass model able to reproduce exactly the dynamics of heterogeneous spiking neural networks encompassing realistic cellular mechanisms for short-term synaptic plasticity. This population model reproduces the macroscopic dynamics of the network in terms of the firing rate and the mean membrane potential. The latter quantity allows us to gain insight of the Local Field Potential and electroencephalographic signals measured during WM tasks to characterize the brain activity. More specifically synaptic facilitation and depression integrate each other to efficiently mimic WM operations via either synaptic reactivation or persistent activity. Memory access and loading are related to stimulus-locked transient oscillations followed by a steady-state activity in the β-γ band, thus resembling what is observed in the cortex during vibrotactile stimuli in humans and object recognition in monkeys. Memory juggling and competition emerge already by loading only two items. However more items can be stored in WM by considering neural architectures composed of multiple excitatory populations and a common inhibitory pool. Memory capacity depends strongly on the presentation rate of the items and it maximizes for an optimal frequency range. In particular we provide an analytic expression for the maximal memory capacity. Furthermore, the mean membrane potential turns out to be a suitable proxy to measure the memory load, analogously to event driven potentials in experiments on humans. Finally we show that the γ power increases with the number of loaded items, as reported in many experiments, while θ and β power reveal non monotonic behaviours. In particular, β and γ rhythms are crucially sustained by the inhibitory activity, while the θ rhythm is controlled by excitatory synapses.