Modeling partial lockdowns in multiplex networks using partition strategies.

National stay-at-home orders, or lockdowns, were imposed in several countries to drastically reduce the social interactions mainly responsible for the transmission of the SARS-CoV-2 virus. Despite being essential to slow down the COVID-19 pandemic, these containment measures are associated with an economic burden. In this work, we propose a network approach to model the implementation of a partial lockdown, breaking the society into disconnected components, or partitions. Our model is composed by two main ingredients: a multiplex network representing human contacts within different contexts, formed by a Household layer, a Work layer, and a Social layer including generic social interactions, and a Susceptible-Infected-Recovered process that mimics the epidemic spreading. We compare different partition strategies, with a twofold aim: reducing the epidemic outbreak and minimizing the economic cost associated to the partial lockdown. We also show that the inclusion of unconstrained social interactions dramatically increases the epidemic spreading, while different kinds of restrictions on social interactions help in keeping the benefices of the network partition. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s41109-021-00366-7.

Interplay between [Formula: see text]-core and community structure in complex networks.

The organisation of a network in a maximal set of nodes having at least k neighbours within the set, known as [Formula: see text]-core decomposition, has been used for studying various phenomena. It has been shown that nodes in the innermost [Formula: see text]-shells play a crucial role in contagion processes, emergence of consensus, and resilience of the system. It is known that the [Formula: see text]-core decomposition of many empirical networks cannot be explained by the degree of each node alone, or equivalently, random graph models that preserve the degree of each node (i.e., configuration model). Here we study the [Formula: see text]-core decomposition of some empirical networks as well as that of some randomised counterparts, and examine the extent to which the [Formula: see text]-shell structure of the networks can be accounted for by the community structure. We find that preserving the community structure in the randomisation process is crucial for generating networks whose [Formula: see text]-core decomposition is close to the empirical one. We also highlight the existence, in some networks, of a concentration of the nodes in the innermost [Formula: see text]-shells into a small number of communities.

Inferring directed networks using a rank-based connectivity measure.

Inferring the topology of a network using the knowledge of the signals of each of the interacting units is key to understanding real-world systems. One way to address this problem is using data-driven methods like cross-correlation or mutual information. However, these measures lack the ability to distinguish the direction of coupling. Here, we use a rank-based nonlinear interdependence measure originally developed for pairs of signals. This measure not only allows one to measure the strength but also the direction of the coupling. Our results for a system of coupled Lorenz dynamics show that we are able to consistently infer the underlying network for a subrange of the coupling strength and link density. Furthermore, we report that the addition of dynamical noise can benefit the reconstruction. Finally, we show an application to multichannel electroencephalographic recordings from an epilepsy patient.

Mean field phase synchronization between chimera states.

We study two-layer networks of identical phase oscillators. Each individual layer is a ring network for which a non-local intra-layer coupling leads to the formation of a chimera state. The number of oscillators and their natural frequencies is in general different across the layers. We couple the phases of individual oscillators in one layer to the phase of the mean field of the other layer. This coupling from the mean field to individual oscillators is done in both directions. For a sufficient strength of this inter-layer coupling, the phases of the mean fields lock across the two layers. In contrast, both layers continue to exhibit chimera states with no locking between the phases of individual oscillators across layers, and the two mean field amplitudes remain uncorrelated. Hence, the networks' mean fields show phase synchronization which is analogous to the one between low-dimensional chaotic oscillators. The required coupling strength to achieve this mean field phase synchronization increases with the mismatches in the network sizes and the oscillators' natural frequencies.

Robustness and versatility of a nonlinear interdependence method for directional coupling detection from spike trains.

The detection of directional couplings between dynamics based on measured spike trains is a crucial problem in the understanding of many different systems. In particular, in neuroscience it is important to assess the connectivity between neurons. One of the approaches that can estimate directional coupling from the analysis of point processes is the nonlinear interdependence measure L. Although its efficacy has already been demonstrated, it still needs to be tested under more challenging and realistic conditions prior to an application to real data. Thus, in this paper we use the Hindmarsh-Rose model system to test the method in the presence of noise and for different spiking regimes. We also examine the influence of different parameters and spike train distances. Our results show that the measure L is versatile and robust to various types of noise, and thus suitable for application to experimental data.

Measures of spike train synchrony for data with multiple time scales.

BACKGROUND: Measures of spike train synchrony are widely used in both experimental and computational neuroscience. Time-scale independent and parameter-free measures, such as the ISI-distance, the SPIKE-distance and SPIKE-synchronization, are preferable to time scale parametric measures, since by adapting to the local firing rate they take into account all the time scales of a given dataset. NEW METHOD: In data containing multiple time scales (e.g. regular spiking and bursts) one is typically less interested in the smallest time scales and a more adaptive approach is needed. Here we propose the A-ISI-distance, the A-SPIKE-distance and A-SPIKE-synchronization, which generalize the original measures by considering the local relative to the global time scales. For the A-SPIKE-distance we also introduce a rate-independent extension called the RIA-SPIKE-distance, which focuses specifically on spike timing. RESULTS: The adaptive generalizations A-ISI-distance and A-SPIKE-distance allow to disregard spike time differences that are not relevant on a more global scale. A-SPIKE-synchronization does not any longer demand an unreasonably high accuracy for spike doublets and coinciding bursts. Finally, the RIA-SPIKE-distance proves to be independent of rate ratios between spike trains. COMPARISON WITH EXISTING METHODS: We find that compared to the original versions the A-ISI-distance and the A-SPIKE-distance yield improvements for spike trains containing different time scales without exhibiting any unwanted side effects in other examples. A-SPIKE-synchronization matches spikes more efficiently than SPIKE-synchronization. CONCLUSIONS: With these proposals we have completed the picture, since we now provide adaptive generalized measures that are sensitive to firing rate only (A-ISI-distance), to timing only (ARI-SPIKE-distance), and to both at the same time (A-SPIKE-distance).

Generalized synchronization between chimera states.

Networks of coupled oscillators in chimera states are characterized by an intriguing interplay of synchronous and asynchronous motion. While chimera states were initially discovered in mathematical model systems, there is growing experimental and conceptual evidence that they manifest themselves also in natural and man-made networks. In real-world systems, however, synchronization and desynchronization are not only important within individual networks but also across different interacting networks. It is therefore essential to investigate if chimera states can be synchronized across networks. To address this open problem, we use the classical setting of ring networks of non-locally coupled identical phase oscillators. We apply diffusive drive-response couplings between pairs of such networks that individually show chimera states when there is no coupling between them. The drive and response networks are either identical or they differ by a variable mismatch in their phase lag parameters. In both cases, already for weak couplings, the coherent domain of the response network aligns its position to the one of the driver networks. For identical networks, a sufficiently strong coupling leads to identical synchronization between the drive and response. For non-identical networks, we use the auxiliary system approach to demonstrate that generalized synchronization is established instead. In this case, the response network continues to show a chimera dynamics which however remains distinct from the one of the driver. Hence, segregated synchronized and desynchronized domains in individual networks congregate in generalized synchronization across networks.