Search

Reservoir computing has been shown to be a useful framework for predicting critical transitions of a dynamical system if the bifurcation parameter is also provided as an input. Its utility is significant because in real-world scenarios, the exact model equations are unknown. This Letter shows how the theory of dynamical system provides the underlying mechanism behind the prediction. Using numerical methods, by considering dynamical systems which show Hopf bifurcation, we demonstrate that the map produced by the reservoir after a successful training undergoes a Neimark-Sacker bifurcation such that the critical point of the map is in immediate proximity to that of the original dynamical system. In addition, we have compared and analyzed different structures in the phase space. Our findings provide insight into the functioning of machine learning algorithms for predicting critical transitions.

Impact of noise in coupled oscillators with pairwise interactions has been extensively explored. Here, we study stochastic second-order coupled Kuramoto oscillators with higher-order interactions, and show that as noise strength increases the critical points associated with synchronization transitions shift toward higher coupling values. By employing the perturbation analysis, we obtain an expression for the forward critical point as a function of inertia and noise strength. Further, for overdamped systems we show that as noise strength increases, the first-order transition switches to second-order even for higher-order couplings. We include a discussion on nature of critical points obtained through Ott-Antonsen ansatz.

Finite-size systems of Kuramoto model display intricate dynamics, especially in the presence of multi-stability where both coherent and incoherent states coexist. We investigate such scenario in globally coupled populations of Kuramoto phase oscillators with higher-order interactions, and observe that fluctuations inherent to finite-size systems drives the transition to the synchronized state occurring before the critical point in the thermodynamic limit. Using numerical methods, we plot the first exit time distribution of the magnitude of complex order parameter and obtain numerical transition probabilities across various system sizes. Further, we extend this study to a two-population oscillator system, and using velocity field of the associated order parameters, show the emergence of a new fixed point corresponding to a partially synchronized state arising due to the finite-size effect which is absent in the thermodynamics limit.

Explosive synchronization refers to an abrupt (first order) transition to non-zero phase order parameter in oscillatory networks, underpinned by the bistability of synchronous and asynchronous states. Growing evidence suggests that this phenomenon might be no less general then the celebrated Kuramoto scenario that belongs to the second order universality class. Importantly, the recent examples demonstrate that explosive synchronization can occur for certain network topologies and coupling types, like the global higher-order coupling, without specific requirements on the individial oscillator dynamics or dynamics-network correlations. Here we demonstrate a rich picture of explosive synchronization and desynchronization transitions in multiplex networks, where it is sufficient to have a single random sparsly connected layer with higher-order coupling terms (and not necessarily in the synchronization regime on its own), the other layer being a regular lattice without own phase transitions at all. Moreover, explosive synchronization emerges even when the random layer has only low-order pairwise coupling, althoug the hysteresis interval becomes narrow and explosive desynchronization is no longer observed. The relevance to the normal and pathological dynamics of neural-glial networks is pointed out., Comment: 8 pages, 6 figures

Coupled limit cycle oscillators with pairwise interactions depict phase transitions to amplitude or oscillation death. This Letter introduces a scheme for higher-order interactions, which can not be decomposed into pairwise interactions. We investigate Stuart Landau oscillators' dynamical evolution under the impression of such a coupling scheme and discover a particular type of oscillator death where a coupling-dependent stable death state, away from the origin, arises in isolation without being accompanied by any other stable state. We call such a state a Solitary death state. Moreover, the explosive transition to the death state is preceded by a surge in amplitude, followed by the revival of the oscillations. Such versatile dynamical states are further enriched with sensitivity to initial conditions. Finally, we point out the resemblance of the results with different dynamical states associated with epileptic seizures.

The Ethereum blockchain and its ERC20 token standard have revolutionized the landscape of digital assets and decentralized applications. ERC20 tokens developed on the Ethereum blockchain have gained significant attention since their introduction. They are programmable and interoperable tokens, enabling various applications and token economies. Transaction graphs, representing the flow of the value between wallets within the Ethereum network, have played a crucial role in understanding the system's dynamics, such as token transfers and the behavior of traders. Here, we explore the evolution of daily transaction graphs of ERC20 token transactions, which sheds light on the trader's behavior during the Black Swan Events -- 2018 crypto crash and the COVID-19 pandemic. By using the tools from network science and differential geometry, we analyze 0.98 billion of ERC20 token transaction data from November 2015 to January 2023. Our analysis reveals that ERC20 financial ecosystem has evolved from a localized wealth formation period to a more mature financial ecosystem where wealth has dispersed among the traders in the network after the crypto crash and during the pandemic period. Before the crash, most sellers only sell the tokens, and buyers only buy the tokens. However, after the crash and during the pandemic period, sellers and buyers both performed buying and selling activities. In addition, we observe no significant negative impact of the COVID-19 pandemic on user behavior in the financial ecosystem., Comment: 12 pages, 9 figures

The effect of phase-lag parameter in pairwise interactions has been a topic of great interest for long. However, real-world systems often have interactions that are beyond pairwise and can be modeled using simplicial complexes. We investigate the effect of the inclusion of phase-lag in coupled Kuramoto oscillators with simplicial interactions and find that it shifts the critical points at which first-order transition from cluster synchronized state to incoherent state occurs. In the thermodynamic limit, using the Ott-Antonsen approach we derive a reduced equation for order parameter measuring cluster synchronization. Further, we progress through the self-consistency method to achieve a closed form of the order parameter measuring global synchronization which was lacking in Ott-Antonsen approach. Moreover, considering polar coordinates framework we obtain rotation frequency of the clusters which comes out to be a function of the phase-lag parameter further indicating that phase-lag can be used as a control parameter to achieve a desired cluster frequency.

We investigate the dynamical evolution of globally connected Stuart-Landau oscillators coupled through conjugate or dis-similar variables on simplicial complexes. We report a first-order explosive phase transition from oscillatory state to death state, with 2-simplex (triadic) interactions, as opposed to the second-order transition with only 1-simplex (dyadic) interactions. Moreover, the system displays four distinct homogeneous steady states in the presence of triadic interactions, in contrast to the two homogeneous steady states observed with dyadic interactions. We calculate the backward transition point analytically, confirming the numerical results and providing the origin of the dynamical states in the transition region. The study will be useful in understanding complex systems, such as ecological and epidemiological, having higher-order interactions and coupling through conjugate variables.

We investigate the impact of contrarians (via negative coupling) in a multilayer network of phase oscillators having higher-order interactions. We show that the multilayer framework facilitates synchronization onset in the negative pairwise coupling regime. The multilayering strength governs the onset of synchronization and the nature of the phase transition, whereas the backward critical couplings depend on higher-order interaction strength. The system does not synchronize below a critical value of multilayering strength. The numerical results agree with the analytical predictions using the Ott-Antonsen approach. The results presented here may be useful for understanding emergent behaviors in real-world complex systems with contrarians and higher-order interactions, such as the brain and society.

We show how multiplexing influences propagating fronts in multilayer networks of coupled bistable oscillators. Using numerical simulation, we investigate both deterministic and noise-sustained propagation. In particular, we demonstrate that the multiplexing allows to reduce the intra-layer dynamics to a common regime where the front propagation speed in all the interacting layers attains the same fixed value. In the presence of noise the dynamics is more complicated and is characterized by the ability of the system to adjust to the common propagation speed for varying the multiplexing strength. In addition, we find that the noise-induced stabilization of wavefront propagation in multilayer networks allows to obtain less pronounced deviations of the wavefront compared to the stabilization achieved in the isolated layer. Finally, we demonstrate that the reduction of the wavefront deviations can be enhanced by increasing the number of interacting layers., Comment: 7 pages, 4 figures

Most real-world networks are endowed with the small-world property, by means of which the maximal distance between any two of their nodes scales logarithmically rather than linearly with their size. The evidence sparkled a wealth of studies trying to reveal possible mechanisms through which the pairwise interactions amongst the units of a network are structured in a way to determine such observed regularity. Here we show that smallworldness occurs also when interactions are of higher order. Namely, by considering Q-uniform hypergraphs and a process through which connections can be randomly rewired with given probability p, we find that such systems may exhibit prominent clustering properties in connection with small average path lengths for a wide range of p values, in analogy to the case of dyadic interactions. The nature of small-world transition remains the same at different orders Q of the interactions, however, the increase in the hyperedge order reduces the range of rewiring probability for which smallworldness emerge., Comment: 4 pages, 5 figures

The inclusion of inertia in the Kuramoto model has been long reported to change the nature of phase transition, providing a fertile ground to model the dynamical behaviors of interacting units. More recently, higher-order interactions have been realized as essential for the functioning of real-world complex systems ranging from the brain to disease spreading. Yet, analytical insights to decipher the role of inertia with higher-order interactions remain challenging. Here, we study the Kuramoto model with inertia on simplicial complexes, merging two research domains. We develop an analytical framework in a mean-field setting using self-consistent equations to describe the steady-state behavior, which reveals a prolonged hysteresis in the synchronization profile. Inertia and triadic interaction strength exhibit isolated influence on system dynamics by predominantly governing, respectively, the forward and backward transition points. This work sets a paradigm to deepen our understanding of real-world complex systems such as power grids modeled as the Kuramoto model with inertia., Comment: 5+ pages, 2 figures. Comments welcome

An incorporation of higher-order interactions is known to lead an abrupt first-order transition to synchronization in otherwise smooth second-order one for pair-wise coupled systems. Here, we show that adaptation in higher-order coupling strength may yield completely different phenomena, notably, second-order transition to synchronization and tiered synchronization. Using the Ott-Antonsen approach, we perform rigorous theoretical calculations illustrating the origin of these emerging phenomena along with a complete description of all (un)stable states. Numerical simulations for limited-size networks are in agreement with the analytical predictions. These results would be important to comprehend dynamical behaviors of real-world complex systems with inherent higher-order interactions and adaptation through feedback coupling.

Eigenvalues statistics of various many-body systems have been widely studied using the nearest neighbor spacing distribution under the random matrix theory framework. Here, we numerically analyze eigenvalue ratio statistics of multiplex networks consisting of directed Erdos-Renyi random networks layers represented as, first, weighted non-Hermitian random matrices and then weighted Hermitian random matrices. We report that the multiplexing strength rules the behavior of average spacing ratio statistics for multiplexing networks represented by the non-Hermitian and Hermitian matrices, respectively. Additionally, for both these representations of the directed multiplex networks, the multiplexing strength appears as a guiding parameter for the eigenvector delocalization of the entire system. These results could be important for driving dynamical processes in several real-world multilayer systems, particularly, understanding the significance of multiplexing in comprehending network properties.

Quenching of oscillations, namely amplitude and oscillations death, is an emerging phenomenon exhibited by many real-world complex systems. Here, we introduce a scheme that combines dissimilar couplings and repulsive feedback links for the interactions of Stuart Landau oscillators and analytically derives the conditions required for the amplitude death. Importantly, this analysis is independent of the network size, presents a generalized approach to calculate the stability conditions for various different coupling schemes, and befits for non-identical oscillators as well. Last, we discuss the similarities of the quenching of oscillations phenomenon with the postictal generalized EEG suppression in convulsive seizures.

Correlation matrices contain a wide variety of spatio-temporal information about a dynamical system. Predicting correlation matrices from partial time series information of a few nodes characterizes the spatio-temporal dynamics of the entire underlying system. This information can help to predict the underlying network structure, e.g., inferring neuronal connections from spiking data, deducing causal dependencies between genes from expression data, and discovering long spatial range influences in climate variations. Traditional methods of predicting correlation matrices utilize time series data of all the nodes of the underlying networks. Here, we use a supervised machine learning technique to predict the correlation matrix of entire systems from finite time series information of a few randomly selected nodes. The accuracy of the prediction validates that only a limited time series of a subset of the entire system is enough to make good correlation matrix predictions. Furthermore, using an unsupervised learning algorithm, we furnish insights into the success of the predictions from our model. Finally, we employ the machine learning model developed here to real-world data sets., Comment: 17 pages, 7 figures

Localization behaviours of Laplacian eigenvectors of complex networks provide understanding to various dynamical phenomena on the corresponding complex systems. We numerically investigate role of hyperedges in driving eigenvector localization of hypergraphs Laplacians. By defining a single parameter \gamma which measures the relative strengths of pair-wise and higher-order interactions, we analyze the impact of interactions on localization properties. For, \gamma < 1 there exists no impact of pairwise links on eigenvector localization while the higher-order interactions instigate localization in the larger eigenvalues. For \gamma > 1, pair-wise interactions cause localization of eigenvector corresponding to small eigenvalues, where as higherorder interactions, despite being much lesser than the pair-wise links, keep driving localization of the eigenvectors corresponding to larger eigenvalues. The results will be useful to understand dynamical phenomena such as diffusion, and random walks on a range of real-world complex systems having higher-order interactions., Comment: 8 pages, 7 figures

The presence of higher-order interactions (simplicial complexes) in networks and certain types of multilayer networks has shown to lead to the abrupt first-order transition to synchronization. We discover that simplicial complexes on multilayer networks can yield multiple basins of attraction, leading to multiple routes to the abrupt first-order transition to synchronization. Using the Ott-Antonsen approach, we develop an analytical framework for simplicial complexes on multilayer networks, which thoroughly explains the origin and stability of all possible dynamical states, including multiple synchronization transitions, of the associated coupled dynamics. The study illustrating rich dynamical behaviours could be pivotal to comprehending the impacts of higher-order interactions on dynamics of complex real-world networks, such as brain, social and technological, which have inherent multilayer networks architecture., Comment: 7 pages, 5 figures

The distribution of the ratios of consecutive eigenvalue spacings of random matrices has emerged as an important tool to study spectral properties of many-body systems. This article numerically investigates the eigenvalue ratios distribution of various model networks, namely, small-world, Erd\H{o}s-R\'enyi random, and (dis)assortative random having a diagonal disorder in the corresponding adjacency matrices. Without any diagonal disorder, the eigenvalues ratio distribution of these model networks depict Gaussian orthogonal ensemble (GOE) statistics. Upon adding diagonal disorder, there exists a gradual transition from the GOE to Poisson statistics depending upon the strength of the disorder. The critical disorder (wc) required to procure the Poisson statistics increases with the randomness in the network architecture. We relate wc with the time taken by the maximum entropy random walker to reach the steady-state. These analyses will be helpful to understand the role of eigenvalues other than the principal one for various network dynamics such as transient behaviour., Comment: 10 pages, 11 Figures

This Letter investigates the transition to synchronization of oscillator ensembles encoded by simplicial complexes in which pairwise and higher-order coupling weights alter with time through a rate-based adaptive mechanism inspired by the Hebbian learning rule. These simultaneously evolving disparate adaptive coupling weights lead to a phenomenon in which the in-phase synchronization is completely obliterated; instead, the anti-phase synchronization is originated. In addition, the onsets of antiphase synchronization and desynchronization are manageable through both dyadic and triadic learning rates. The theoretical validation of these numerical assessments is delineated thoroughly by employing Ott-Antonsen dimensionality reduction. The framework and results of the Letter could help in the understanding of the underlying synchronization behavior of a range of real-world systems, such as brain functions and social systems where interactions evolve with time., Comment: 5 figures, 5 pages + supplemental material

This Letter investigates the upshots of adaptive development of pure 2- and 3- simplicial complexes (triad and tetrad) on the nature of the transition to desynchrony of the oscillator ensembles. The adaptation exercised in the pure simplicial coupling takes a cue from the Hebbian learning rule, i.e., the coupling weight of a triad (tetrad) is prone to increase if the oscillators forming it are in phase and decrease if they are out of phase. The coupling weights in these pure simplicial complexes experiencing such adaptation give rise to first-order routes to desynchronization, whose onsets are entirely characterized by respective Hebbian learning parameters. Mean-field analyses presented for the order parameters for the adaptive 2- and 3- simplicial complexes strongly corroborate with the respective numerical assessments., Comment: 5+ pages, 9 figures

We investigate the spectra of adjacency matrices of multiplex networks under random matrix theory (RMT) framework. Through extensive numerical experiments, we demonstrate that upon multiplexing two random networks, the spectra of the combined multiplex network exhibit superposition of two Gaussian orthogonal ensemble (GOE)s for very small multiplexing strength followed by a smooth transition to the GOE statistics with an increase in the multiplexing strength. Interestingly, randomness in the connection architecture, introduced by random rewiring to 1D lattice, of at least one layer may govern nearest neighbor spacing distribution (NNSD) of the entire multiplex network, and in fact, can drive to a transition from the Poisson to the GOE statistics or vice versa. Notably, this transition transpires for a very small number of the random rewiring corresponding to the small-world transition. Ergo, only one layer being represented by the small-world network is enough to yield GOE statistics for the entire multiplex network. Spectra of adjacency matrices of underlying interaction networks have been contemplated to be related with dynamical behaviour of the corresponding complex systems, the investigations presented here have implications in achieving better structural and dynamical control to the systems represented by multiplex networks against structural perturbation in only one of the layers., Comment: 20 pages, 14 figures

This Letter investigates the nature of synchronization in multilayered and multiplexed populations in which the interlayer interactions are randomly pinned. First, we show that a multilayer network constructed by setting up all-to-all interlayer connections between the two populations leads to explosive synchronization in the two populations successively, leading to the coexistence of coherent and incoherent populations forming chimera states. Second, a multiplex formation of the two populations in which only the mirror nodes are interconnected espouses explosive transitions in the two populations concurrently. The occurrence of both explosive synchronization and chimera are substantiated with rigorous theoretical mean-field analysis. The random pinning in the interlayer interactions concerns the practical problems where the impact of dynamics of one network on that of other interconnected networks remains elusive, as is the case for many real-world systems., Comment: 6 pages, 6 figures

We show that an introduction of a phase parameter ($\alpha$), with $0 \le \alpha \le \pi/2$, in the interlayer coupling terms of multiplex networks of Kuramoto oscillators can induce explosive synchronization (ES) in the multiplexed layers. Along with the {\alpha} values, the hysteresis width is determined by the interlayer coupling strength and the frequency mismatch between the mirror (inter-connected) nodes. A mean-field analysis is performed to support the numerical results. Similar to the earlier works, we find that the suppression of synchronization is accountable for the origin of ES. The robustness of ES against changes in the network topology and frequency distribution is tested. Finally, taking a suggestion from the synchronized state of the multiplex networks, we extend the results to the classical concept of the single-layer networks in which some specific links are assigned a phase-shifted coupling. Different methods have been introduced in the past years to incite ES in coupled oscillators; our results indicate that a phase-shifted coupling can also be one such method to achieve ES., Comment: 9 pages, 8 figures

Chimera and Solitary states have captivated scientists and engineers due to their peculiar dynamical states corresponding to the co-existence of coherent and incoherent dynamical evolution in coupled units in various natural and artificial systems. It has been further demonstrated that such states can be engineered in systems of coupled oscillators by the suitable implementation of communication delays. Here, using supervised machine learning, we predict (a) the precise value of delay which is sufficient for engineering chimera and solitary states for a given set of system parameters, as well as (b) the intensity of incoherence for such engineered states. The results are demonstrated for two different examples consisting of single layer and multi layer networks. First, the chimera states (solitary states) are engineered by establishing delays in the neighboring links of a node (the interlayer links) in a 2-D lattice (multiplex network) of oscillators. Then, different machine learning classifiers, KNN, SVM and MLP-Neural Network are employed by feeding the data obtained from the network models. Once a machine learning model is trained using a limited amount of data, it makes predictions for a given unknown systems parameter values. Testing accuracy, sensitivity, and specificity analysis reveal that MLP-NN classifier is better suited than Knn or SVM classifier for the predictions of parameters values for engineered chimera and solitary states. The technique provides an easy methodology to predict critical delay values as well as the intensity of incoherence for designing an experimental setup to create solitary and chimera states., Comment: 11 Pages, 9 Figures, Contains revised abstract and publication details

Adaptation plays a pivotal role in the evolution of natural and artificial complex systems, and in the determination of their functionality. Here, we investigate the impact of adaptive inter-layer processes on intra-layer synchronization in multiplex networks. The considered adaptation mechanism is governed by a Hebbian learning rule, i.e., the link weight between a pair of interconnected nodes is enhanced if the two nodes are in phase. Such adaptive coupling induces an irreversible first-order transition route to synchronization accompanied with a hysteresis. We provide rigorous analytic predictions of the critical coupling strengths for the onset of synchronization and de-synchronization, and verify all our theoretical predictions by means of extensive numerical simulations., Comment: 9 figures, 15 pages

Many real-world complex systems have small-world topology characterized by the high clustering of nodes and short path lengths.It is well-known that higher clustering drives localization while shorter path length supports delocalization of the eigenvectors of networks. Using multifractals technique, we investigate localization properties of the eigenvectors of the adjacency matrices of small-world networks constructed using Watts-Strogatz algorithm. We find that the central part of the eigenvalue spectrum is characterized by strong multifractality whereas the tail part of the spectrum have Dq->1. Before the onset of the small-world transition, an increase in the random connections leads to an enhancement in the eigenvectors localization, whereas just after the onset, the eigenvectors show a gradual decrease in the localization. We have verified an existence of sharp change in the correlation dimension at the localization-delocalization transition, Comment: 8 pages, 7 figures

The phenomenon of explosive synchronization, which originates from hypersensitivity to small perturbation caused by some form of frustration prevailed in various physical and biological systems, has been shown to lead events of cascading failure of the power grid to chronic pain or epileptic seizure in the brain. Furthermore, networks provide a powerful model to understand and predict the properties of a diverse range of real-world complex systems. Recently, a multilayer network has been realized as a better suited framework for the representation of complex systems having multiple types of interactions among the same set of constituents. This article shows that by tuning the properties of one layer (network) of a multilayer network, one can regulate the dynamical behavior of another layer (network). By taking an example of a multiplex network comprising two different types of networked Kuramoto oscillators representing two different layers, this article attempts to provide a glimpse of opportunities and emerging phenomena multiplexing can induce which is otherwise not possible for a network in isolation. Here we consider explosive synchronization to demonstrate the potential of multilayer networks framework. To the end, we discuss several possible extensions of the model considered here by incorporating real-world properties., Comment: 11 pages, 8 figures

Investigation of eigenvector localization properties of complex networks is not only important for gaining insight into fundamental network problems such as network centrality measure, spectral partitioning, development of approximation algorithms, but also is crucial for understanding many real-world phenomena such as disease spreading, criticality in brain network dynamics. For a network, an eigenvector is said to be localized when most of its components take value near to zero, with a few components taking very high values. In this article, we devise a methodology to construct a principal eigenvector (PEV) localized network from a given input network. The methodology relies on adding a small component having a wheel graph to the given input network. By extensive numerical simulation and an analytical formulation based on the largest eigenvalue of the input network, we compute the size of the wheel graph required to localize the PEV of the combined network. Using the susceptible-infected-susceptible model, we demonstrate the success of this method for various models and real-world networks consider as input networks. We show that on such PEV localized networks, the disease gets localized within a small region of the network structure before the outbreaks. The study is relevant in controlling spreading processes on complex systems represented by networks., Comment: 6 pages, 5 figures

Chimera state refers to coexistence of coherent and non-coherent phases in identically coupled dynamical units found in various complex dynamical systems. Identification of Chimera, on one hand is essential due to its applicability in various areas including neuroscience, and on other hand is challenging due to its widely varied appearance in different systems and the peculiar nature of its profile. Therefore, a simple yet universal method for its identification remains an open problem. Here, we present a very distinctive approach using machine learning techniques to characterize different dynamical phases and identify the chimera state from given spatial profiles generated using various different models. The experimental results show that the performance of the classification algorithms varies for different dynamical models. The machine learning algorithms, namely random forest, oblique random forest based on tikhonov, parallel-axis split and null space regularization achieved more than $96\% $ accuracy for the Kuramoto model. For the logistic-maps, random forest and tikhonov regularization based oblique random forest showed more than $90\%$ accuracy, and for the H\'enon-Map model, random forest, null-space and axis-parallel split regularization based oblique random forest achieved more than $80\%$ accuracy. The oblique random forest with null space regularization achieved consistent performance (more than $83\%$ accuracy) across different dynamical models while the auto-encoder based random vector functional link neural network showed relatively lower performance. This work provides a direction for employing machine learning techniques to identify dynamical patterns arising in coupled non-linear units on large-scale, and for characterizing complex spatio-temporal patterns in real-world systems for various applications., Comment: 20 Pages, 4 Figures; Comments welcome

It is known that intra-layer adaptive coupling among connected oscillators instigates explosive synchronization (ES) in multilayer networks. Taking an altogether different cue in the present work, we consider inter-layer adaptive coupling in a multiplex network of phase oscillators and show that the scheme gives rise to ES with an associated hysteresis irrespective of the network architecture of individual layers. The hysteresis is shaped by the inter-layer coupling strength and the frequency mismatch between the mirror nodes. We provide rigorous mean-field analytical treatment for the measure of global coherence and manifest they are in a good match with respective numerical assessments. Moreover, the analytical predictions provide a complete insight into how adaptive multiplexing suppresses the formation of a giant cluster, eventually giving birth to ES. The study will help in spotlighting the role of multiplexing in the emergence of ES in real-world systems represented by multilayer architecture. Particularly, it is relevant to those systems which have limitations towards change in intra-layer coupling strength., Comment: 7 Pages, 6 Figures

One of the most challenging problems in biomedicine and genomics is the identification of disease biomarkers. In this study, proteomics data from seven major cancers were used to construct two weighted protein-protein interaction (PPI) networks i.e., one for the normal and another for the cancer conditions. We developed rigorous, yet mathematically simple, methodology based on the degeneracy at -1 eigenvalues to identify structural symmetry or motif structures in network. Utilising eigenvectors corresponding to degenerate eigenvalues in the weighted adjacency matrix, we identified structural symmetry in underlying weighted PPI networks constructed using seven cancer data. Functional assessment of proteins forming these structural symmetry exhibited the property of cancer hallmarks. Survival analysis refined further this protein list proposing BMI, MAPK11, DDIT4, CDKN2A, and FYN as putative multi-cancer biomarkers. The combined framework of networks and spectral graph theory developed here can be applied to identify symmetrical patterns in other disease networks to predict proteins as potential disease biomarkers., Comment: 20 Pages, 5 Figures

Complex networks or graphs provide a powerful framework to understand importance of individuals and their interactions in real-world complex systems. Several graph theoretical measures have been introduced to access importance of the individual in systems represented by networks. Particularly, eigenvector centrality (EC) measure has been very popular due to its ability in measuring importance of the nodes based on not only number of interactions they acquire but also particular structural positions they have in the networks. Furthermore, the presence of certain structural features, such as the existence of high degree nodes in a network is recognized to induce localization transition of the principal eigenvector (PEV) of the network's adjacency matrix. Localization of PEV has been shown to cause difficulties in assigning centrality weights to the nodes based on the EC. We revisit PEV localization and its relation with failure of EC problem, and by using simple model networks demonstrate that in addition to the localization of the PEV, the delocalization of PEV may also create difficulties for using EC as a measure to rank the nodes. Our investigation while providing fundamental insight to the relation between PEV localization and centrality of nodes in networks, suggests that for the networks having delocalized PEVs, it is better to use degree centrality measure to rank the nodes., Comment: 13 pages, 7 figures

We present a technique to engineer solitary states by means of delayed links in a network of neural oscillators and in coupled chaotic maps. Solitary states are intriguing partial synchronization patterns, where a synchronized cluster coexists with solitary nodes displaced from this cluster and distributed randomly over the network. We induce solitary states in the originally synchronized network of identical nodes by introducing delays in the links for a certain number of selected network elements. It is shown that the extent of displacement and the position of solitary elements can be completely controlled by the choice (values) and positions (locations) of the incorporated delays, reshaping the delay engineered solitary states in the network., Comment: 15 pages, 10 figures, Accepted for publication in Chaos, Solitons & Fractals, Comments welcome

Chimera referring to a coexistence of coherent and incoherent states, is traditionally very difficult to control due to its peculiar nature. Here, we provide a recipe to construct chimera states in the multiplex networks with the aid of multiplexing-delays. The chimera state in multiplex networks is produced by introducing heterogeneous delays in a fraction of inter-layer links, referred as multiplexing-delay, in a sequence. Additionally, the emergence of the incoherence in the chimera state can be regulated by making appropriate choice of both inter- and intra-layer coupling strengths, whereas the extent and the position of the incoherence regime can be regulated by appropriate placing and {strength} of the multiplexing delays. The proposed technique to construct such {engineered} chimera equips us with multiplex network's structural parameters as tools in gaining both qualitative- and quantitative-control over the incoherent section of the chimera states and, in turn, the chimera. Our investigation can be of worth in controlling dynamics of multi-level delayed systems and attain desired chimeric patterns., Comment: 12 pages, 8 figures, Accepted for publication in Euro. Phys. Lett., Supplementary Material avaliable in https://figshare.com/s/6a931b8cfbd5c6f9ded6; Comments Welcome

To date, explosive synchronization (ES) is shown to be originated from either degree-frequency correlation or inertia of phase oscillators. Of late, it has been shown that ES can be induced in a network by adaptively controlled phase oscillators. Here we show that ES is a generic phenomenon and can occur in any network by appropriately multiplexing it with another layer. We devise an approach which leads to the occurrence of ES with hysteresis loop in a network upon its multiplexing with a negatively coupled (or inhibitory) layer. We discuss the impact of various structural properties of positively coupled (or excitatory) and inhibitory layer along with the strength of multiplexing in gaining control over the induced ES transition. This investigation is a step forward in highlighting the importance of multiplex framework not only in bringing novel phenomena which are not possible in an isolated network but also in providing more structural control over the induced phenomena., Comment: 10 pages, 13 figures

It is known that explosive synchronization (ES) in an isolated network of Kuramoto oscillators with inertia is significantly enhanced by the presence of time delay. Here we show that time delay in one layer of the multiplex network governs the transition to synchronization and ES in the other layers. We found that a single layer with time-delayed intra-layer coupling, depending on the values of time delay, experiences a different type of transition to synchronization, e.g., ES or continuous, and the same type of transition is incorporated simultaneously in other layer(s) as well. Hence, a suitable choice of time-delay in only one layer can lead to a desired (either ES or continuous) transition simultaneously in the delayed and other undelayed layers of multiplexed network. These results offer a platform for a better understanding of the dynamics of those complex systems which are represented by multilayered framework and contain time delays in the communication processes., Comment: 14 pages, 9 figures

This review presents an account of the major works done on spectra of adjacency matrices drawn on networks and the basic understanding attained so far. We have divided the review under three sections: (a) extremal eigenvalues, (b) bulk part of the spectrum and (c) degenerate eigenvalues, based on the intrinsic properties of eigenvalues and the phenomena they capture. We have reviewed the works done for spectra of various popular model networks, such as the Erd\H{o}s-R\'enyi random networks, scale-free networks, 1-d lattice, small-world networks, and various different real-world networks. Additionally, potential applications of spectral properties for natural processes have been reviewed., Comment: 29 pages, 18 figures

Understanding the localization properties of eigenvectors of complex networks is important to get insight into various structural and dynamical properties of the corresponding systems. Here, we analytically develop a scheme to construct a highly localized network for a given set of networks parameters that is the number of nodes and the number of interactions. We find that the localization behavior of the principal eigenvector (PEV) of such a network is sensitive against a single edge rewiring. We find evidences for eigenvalue crossing phenomena as a consequence of the single edge rewiring, in turn providing an origin to the sensitive behavior of the PEV localization. These insights were then used to analytically construct the highly localized network for a given set of networks parameters. The analysis provides fundamental insight into relationships between the structural and the spectral properties of networks for PEV localized networks. Further, we substantiate the existence of the eigenvalue crossing phenomenon by considering a linear-dynamical process, namely the ribonucleic acid (RNA) neutral network population dynamical model. The analysis presented here on model networks aids in understanding the steady-state behavior of a broad range of linear-dynamical processes, from epidemic spreading to biochemical dynamics associated with the adjacency matrices., Comment: 15 pages and 8 figures

We investigate spatio-temporal patterns occurring in a two-layer multiplex network of oscillatory FitzHugh-Nagumo neurons, where each layer is represented by a nonlocally coupled ring. We show that weak multiplexing, i.e., when the coupling between the layers is smaller than that within the layers, can have a significant impact on the dynamics of the neural network. We develop control strategies based on weak multiplexing and demonstrate how the desired state in one layer can be achieved without manipulating its parameters, but only by adjusting the other layer. We find that for coupling range mismatch weak multiplexing leads to the appearance of chimera states with different shapes of the mean velocity profile for parameter ranges where they do not exist in isolation. Moreover, we show that introducing a coupling strength mismatch between the layers can suppress chimera states with one incoherent domain (one-headed chimeras) and induce various other regimes such as in-phase synchronization or two-headed chimeras. Interestingly, small intra-layer coupling strength mismatch allows to achieve solitary states throughout the whole network.

The Integrated Information is a quantitative measure from information theory how tightly all parts of a system are interconnected in terms of information exchange. In this study we show that astrocyte, playing an important role in regulation of information transmission between neurons, may contribute to a generation of positive Integrated Information in neuronal ensembles. Analytically and numerically we show that the presence of astrocyte may be essential for this information attribute in neuro-astrocytic ensembles. Moreover, the proposed "spiking-bursting" mechanism of generating positive Integrated Information is shown to be generic and not limited to neuroglial networks, and is given a complete analytic description.

Symmetry breaking spatial patterns, referred to as chimera states, have recently been catapulted into the limelight due to their coexisting coherent and incoherent hybrid dynamics. Here, we present a method to engineer a chimera state by using an appropriate distribution of heterogeneous time delays on the edges of a network. The time delays in interactions, intrinsic to natural or artificial complex systems, are known to induce various modifications in spatiotemporal behaviors of the coupled dynamics on networks. Using a coupled chaotic map with the identical coupling environment, we demonstrate that control over the spatial location of the incoherent region of a chimera state in a network can be achieved by appropriately introducing time delays. This method allows for the engineering of tailor-made one cluster or multi-cluster chimera patterns. Furthermore, borrowing a measure of eigenvector localization from the spectral graph theory, we introduce a spatial inverse participation ratio, which provides a robust way for the identification of the chimera state. This report highlights the necessity to consider the heterogeneous time delays to develop applications for the chimera states in particular and understand coupled dynamical systems in general., Comment: 7 pages,6 figures. Supplemental Material available at https://figshare.com/s/1f836f334db39dc3cd58

Information of localization properties of eigenvectors of the complex network has applicability in many different areas which include networks centrality measures, spectral partitioning, development of approximation algorithms, and disease spreading phenomenon. For linear dynamical process localization of principal eigenvector (PEV) of adjacency matrices infers condensation of the information in the smaller section of the network. For a network, an eigenvector is said to be localized when most of its components are near to zero with few taking very high values. Here, we provide three different random-sampling-based algorithms which, by using the edge rewiring method, can evolve a random network having a delocalized PEV to a network having a highly localized PEV. In other words, we develop a learning framework to explore the localization of PEV through a random sampling-based optimization method. We discuss the drawbacks and advantages of these algorithms. Additionally, we show that the construction of such networks corresponding to the highly localized PEV is a non-convex optimization problem when the objective function is the inverse participation ratio. This framework is also relevant to construct a network structure for other lower-order eigenvectors., Comment: 15 pages, 7 figures

We study localization properties of principal eigenvector (PEV) of multilayer networks. Starting with a multilayer network corresponding to a delocalized PEV, we rewire the network edges using an optimization technique such that the PEV of the rewired multilayer network becomes more localized. The framework allows us to scrutinize structural and spectral properties of the networks at various localization points during the rewiring process. We show that rewiring only one-layer is enough to attain a multilayer network having a highly localized PEV. Our investigation reveals that a single edge rewiring of the optimized multilayer network can lead to the complete delocalization of a highly localized PEV. This sensitivity in the localization behavior of PEV is accompanied by a pair of almost degenerate eigenvalues. This observation opens an avenue to gain a deeper insight into the origin of PEV localization of networks. Furthermore, analysis of multilayer networks constructed using real-world social and biological data show that the localization properties of these real-world multilayer networks are in good agreement with the simulation results for the model multilayer network. The study is relevant to applications that require understanding propagation of perturbation in multilayer networks., Comment: 10 pages, 8 figures

Explosive synchronization (ES) of coupled oscillators on networks is shown to be originated from existence of correlation between natural frequencies of oscillators and degrees of corresponding nodes. Here, we demonstrate that ES is a generic feature of multiplex network of second-order Kuramoto oscillators and can exist in absence of a frequency-degree correlation. A monoplex network of second-order Kuramoto oscillators bearing homogeneous (heterogeneous) degree-distribution is known to display the first-order (second-order) transition to synchronization. We report that multiplexing of two such networks having homogeneous degree-distribution support the first-order transition in both the layers thereby facilitating ES. More interesting is the multiplexing of a layer bearing heterogeneous degree-distribution with another layer bearing homogeneous degree-distribution, which induces a first-order (ES) transition in the heterogeneous layer which was incapable of showing the same in the isolation. Further, we report that such induced ES transition in the heterogeneous layer of multiplex networks can be controlled by varying inter and intra-layer coupling strengths. Our findings emphasize on importance of multiplexing or impact of one layer on dynamical evolution of other layers of systems having inherent multiplex or multilevel architecture., Comment: 7 pages, 10 figures