Search

Many natural and human-made complex systems feature group interactions that adapt over time in response to their dynamic states. However, most of the existing adaptive network models fall short of capturing these group dynamics, as they focus solely on pairwise interactions. In this study, we employ adaptive higher-order networks to describe these systems by proposing a general framework incorporating both adaptivity and group interactions. We demonstrate that global synchronization can exist in those complex structures, and we provide the necessary conditions for the emergence of a stable synchronous state. Additionally, we analyzed some relevant settings, and we showed that the necessary condition is strongly related to the master stability equation, allowing to separate the dynamical and structural properties. We illustrate our theoretical findings through examples involving adaptive higher-order networks of coupled generalized Kuramoto oscillators with phase lag. We also show that the interplay of group interactions and adaptive connectivity results in the formation of stability regions that can induce transitions between synchronization and desynchronization

We report a new mechanism through which extreme events with a dragon king-like distribution emerge in a network of locally coupled Hindmarsh-Rose bursting neurons. We establish and substantiate the fact that depending on the choice of initial conditions, the neurons in the network are divided into clusters and whenever these clusters are phase synchronized intermittently, extreme events originate in the collective observable. This mechanism, which we name as intermittent cluster synchronization is proposed as the new precursor for the generation of extreme events in this system. These results are also true for electrical diffusive coupling. The distribution of the local maxima shows long tailed non-Gaussian while the interevent interval follows the Weibull distribution. The goodness of fit are corroborated using probability-probability plot and quantile-quantile plot. These extreme events become rarer and rarer with the increase in the number of different initial conditions., Comment: 11 pages, 12 figures

In this paper, we propose a network whose nodes are labeled by the composite numbers and two nodes are connected by an undirected link if they are relatively prime to each other. As the size of the network increases, the network will be connected whenever the largest possible node index $n\geq 49$. To investigate how the nodes are connected, we analytically describe that the link density saturates to $6/\pi^2$, whereas the average degree increases linearly with slope $6/\pi^2$ with the size of the network. To investigate how the neighbors of the nodes are connected to each other, we find the shortest path length will be at most 3 for $49\leq n\leq 288$ and it is at most 2 for $n\geq 289$. We also derive an analytic expression for the local clustering coefficients of the nodes, which quantifies how close the neighbors of a node to form a triangle. We also provide an expression for the number of $r$-length labeled cycles, which indicates the existence of a cycle of length at most $O(\log n)$. Finally, we show that this graph sequence is actually a sequence of weakly pseudo-random graphs. We numerically verify our observed analytical results. As a possible application, we have observed less synchronizability (the ratio of the largest and smallest positive eigenvalue of the Laplacian matrix is high) as compared to Erd\H{o}s-R\'{e}nyi random network and Barab\'{a}si-Albert network. This unusual observation is consistent with the prolonged transient behaviors of ecological and predator-prey networks which can easily avoid the global synchronization., Comment: 23 pages, 7 figures

Most complex systems are nonlinear, relying on emergent behavior from interacting subsystems, often characterized by oscillatory dynamics. Collective oscillatory behavior is essential for the proper functioning of many real world systems. Complex networks have proven efficient in elucidating the topological structures of both natural and artificial systems and describing diverse processes occurring within them. Recent advancements have significantly enhanced our understanding of emergent dynamics in complex networks. Among various processes, a substantial body of work explores the dynamical robustness of complex networks, their ability to withstand degradation in network constituents while maintaining collective oscillatory dynamics. Many physical and biological systems experience a decline in dynamic activities due to natural or environmental factors. The impact of such damages on network performance can be significant, and the system's robustness indicates its capability to maintain functionality despite dynamic changes, often termed aging. This review provides a comprehensive overview of notable research examining how networks sustain global oscillation despite increasing inactive dynamical units. We present contemporary research dedicated to the theoretical understanding and enhancement mechanisms of dynamical robustness in complex networks. Our focus includes various network structures and coupling functions, elucidating the persistence of networked systems. We cover system characteristics from heterogeneity in network connectivity to heterogeneity in dynamical units. Finally, we discuss challenges in this field and open areas for future studies., Comment: 52 pages, 33 figures

Networks incorporating higher-order interactions are increasingly recognized for their ability to introduce novel dynamics into various processes, including synchronization. Previous studies on synchronization within multilayer networks have often been limited to specific models, such as the Kuramoto model, or have focused solely on higher-order interactions within individual layers. Here, we present a comprehensive framework for investigating synchronization, particularly global synchronization, in multilayer networks with higher-order interactions. Our framework considers interactions beyond pairwise connections, both within and across layers. We demonstrate the existence of a stable global synchronous state, with a condition resembling the master stability function, contingent on the choice of coupling functions. Our theoretical findings are supported by simulations using Hindmarsh-Rose neuronal and R\"{o}ssler oscillators. These simulations illustrate how synchronization is facilitated by higher-order interactions, both within and across layers, highlighting the advantages over scenarios involving interactions within single layers.

Swarmalators are entities that swarm through space and sync in time and are potentially considered to replicate the complex dynamics of many real-world systems. So far, the internal dynamics of swarmalators have been taken as a phase oscillator inspired by the Kuramoto model. Here, for the first time, we examine the internal dynamics utilizing an amplitude oscillator capable of exhibiting periodic and chaotic behaviors. To incorporate the dual interplay between spatial and internal dynamics, we propose a general model that keeps the properties of swarmalators intact. This adaptation calls for a detailed study which we present in this paper. We establish our study with the Rossler oscillator by taking parameters from both the chaotic and periodic regions. While the periodic oscillator mimics most of the patterns in the previous phase oscillator model, the chaotic oscillator brings some new fascinating states., Comment: 12 pages, 9 figures; Accepted for publication in Physical Review E (2024)

Ensembles of coupled nonlinear oscillators are a popular paradigm and an ideal benchmark for analyzing complex collective behaviors. The onset of cluster synchronization is found to be at the core of various technological and biological processes. The current literature has investigated cluster synchronization by focusing mostly on the case of attractive coupling among the oscillators. However, the case of two coexisting competing interactions is of practical interest due to their relevance in diverse natural settings, including neuronal networks consisting of excitatory and inhibitory neurons, the coevolving social model with voters of opposite opinions, ecological plant communities with both facilitation and competition, to name a few. In the present article, we investigate the impact of repulsive spanning trees on cluster formation within a connected network of attractively coupled limit cycle oscillators. We successfully predict which nodes belong to each cluster and the emergent frustration of the connected networks independent of the particular local dynamics at the network nodes. We also determine local asymptotic stability of the cluster states using an approach based on the formulation of a master stability function. We additionally validate the emergence of solitary states and antisynchronization for some specific choices of spanning trees and networks., Comment: 19 pages, 18 figures, Accepted for publication in Physical Review E (2024)

Self-organized bistability (SOB) stands as a critical behavior for the systems delicately adjusting themselves to the brink of bistability, characterized by a first-order transition. Its essence lies in the inherent ability of the system to undergo enduring shifts between the coexisting states, achieved through the self-regulation of a controlling parameter. Recently, SOB has been established in a scale-free network as a recurrent transition to a short-living state of global synchronization. Here, we embark on a theoretical exploration that extends the boundaries of the SOB concept on a higher-order network (implicitly embedded microscopically within a simplicial complex) while considering the limitations imposed by coupling constraints. By applying Ott-Antonsen dimensionality reduction in the thermodynamic limit to the higher-order network, we derive SOB requirements under coupling limits that are in good agreement with numerical simulations on systems of finite size. We use continuous synchronization diagrams and statistical data from spontaneous synchronized events to demonstrate the crucial role SOB plays in initiating and terminating temporary synchronized events. We show that under weak coupling consumption, these spontaneous occurrences closely resemble the statistical traits of the epileptic brain functioning., Comment: 11 Pages, 8 Figures, Accepted for publication in Physical Review E (2024)

Swarmalators are oscillators that swarm through space as they synchronize in time. Introduced a few years ago to model many systems which mix synchrony with self-assembly, they remain poorly understood theoretically. Here we obtain the first analytic results on swarmalators moving in two-dimensional (2D) plane by enforcing periodic boundary conditions; this simpler topology allows expressions for order parameters, stabilities, and bifurcations to be derived exactly. We suggest some future directions for swarmalator research and point out some connections to the Kuramoto model and the Vicsek model from active matter; these are intended as a call-to-arms for the sync community and other researchers looking for new problems and puzzles to work on.

Synchronization is one of the emerging collective phenomena in interacting particle systems. Its ubiquitous presence in nature, science, and technology has fascinated the scientific community over the decades. Moreover, a great deal of research has been, and is still being, devoted to understand various physical aspects of the subject. In particular, the study of interacting \textit{active} particles has led to exotic phase transitions in such systems which have opened up a new research front-line. Motivated by this line of work, in this paper, we study the directional synchrony among self-propelled particles. These particles move inside a bounded region, and crucially their directions are also coupled with spatial degrees of freedom. We assume that the directional coupling between two particles is influenced by the relative spatial distance which changes over time. Furthermore, the nature of the influence is considered to be both short and long-ranged. We explore the phase transition scenario in both the cases and propose an approximation technique which enables us to analytically find the critical transition point. The results are further supported with numerical simulations. Our results have potential importance in the study of active systems like bird flocks, fish schools and swarming robots where spatial influence plays a pertinent role., Comment: Accepted for publication in Chaos (2023)

We present a case study of swarmalators (mobile oscillators) which move on a 1D ring and are subject to pinning. Previous work considered the special case where the pinning in space and the pinning in the phase dimension were correlated. Here we study the general case where the space and phase pinning are uncorrelated, both being chosen uniformly at random. This induces several new effects, such as pinned async, mixed states, and a first order phase transition. These phenomena may be found in real world swarmalators such as systems of vinegar eels, Janus matchsticks, electrorotated Quincke rollers or Japanese tree frogs., Comment: 9 pages, 6 figures

Over the past few decades, the research community has been interested in the study of multi-agent systems and their emerging collective dynamics. These systems are all around us in nature, like bacterial colonies, fish schools, bird flocks, as well as in technology, such as microswimmers and robotics, to name a few. Flocking and swarming are two key components of the collective behaviours of multi-agent systems. In flocking, the agents coordinate their direction of motion, but in swarming, they congregate in space to organise their spatial position. We investigate a minimal mathematical model of locally interacting multi-agent system where the agents simultaneously swarm in space and exhibit flocking behaviour. Various cluster structures are found, depending on the interaction range. When the coupling strength value exceeds a crucial threshold, flocking behaviour is observed. We do in-depth simulations and report the findings by changing the other parameters and with the incorporation of noise., Comment: Accepted for publication in Chaos (2023)

Extreme precipitation in the Arctic region plays a crucial role in global weather and climate patterns. Bear Island (Bj{\o}rn{\o}ya) is located in the Norwegian Svalbard archipelago, which is, therefore, selected for our study on extreme precipitation. The island occupies a unique geographic position at the intersection of the high and low Arctic, characterized by a flat and lake-filled northern region contrasting with mountainous terrain along its southern shores. Its maritime-polar climate is influenced by North Atlantic currents, resulting in relatively mild winter temperatures. An increase in precipitation level in Bear Island is a significant concern linked to climate change and has global implications. We have collected the amount of daily precipitation as well as daily maximum temperatures from the meteorological station of Bj{\o}rn{\o}ya located on the island, operated by the Norwegian Centre for Climate Services for a period spanning from January 1, 1960 to December 31, 2021. We observe that the trend of yearly mean precipitation during this period linearly increases. We analyze the recorded data to investigate the changing pattern of precipitation extremes over the climate scales. We employ the generalized extreme value distribution to model yearly and seasonal maxima of daily precipitation amount and determine the return levels and return period of precipitation extremes. We compare the variability of precipitation extremes between the two time periods: (i) 1960-1990 and (ii) 1991-2021. Our analysis reveals an increase in the frequency of precipitation extremes occurrences between 1991 and 2021. Our findings establish a better understanding of precipitation extremes in Bear Island from a statistical viewpoint, with an observation of seasonal and yearly variability, especially, during the period of the last 31 years.

In a predator-prey metapopulation, the two traits are adversely related: synchronization and persistence. A decrease in synchrony apparently leads to an increase in persistence and, therefore, necessitates the study of desynchrony in a metapopulation. In this article, we study predator-prey patches that communicate with one another while being interconnected through distinct dispersal structures in the layers of a three-layer multiplex network. We investigate the synchronization phenomenon among the patches of the outer layers by introducing higher-order interactions (specifically three-body interactions) in the middle layer. We observe a decrease in the synchronous behavior or, alternatively, an increase in desynchrony due to the inclusion of group interactions among the patches of the middle layer. The advancement of desynchrony becomes more prominent with increasing strength and numbers of three-way interactions in the middle layer. We analytically validated our numerical results by performing the stability analysis of the referred synchronous solution using the master stability function approach. Additionally, we verify our findings by taking into account two distinct predator-prey models and dispersal topologies, which ultimately assert that the findings are generalizable across various models and dispersal structures., Comment: 12 Pages, 10 figures. Accepted for publication in Physical Review E, 2023

Swarmalators are oscillatory systems endowed with a spatial component, whose spatial and phase dynamics affect each other. Such systems can demonstrate fascinating collective dynamics resembling many real-world processes. Through this work, we study a population of swarmalators where they are divided into different communities. The strengths of spatial attraction, repulsion as well as phase interaction differ from one group to another. Also, they vary from inter-community to intra-community. We encounter, as a result of variation in the phase coupling strength, different routes to achieve the static synchronization state by choosing several parameter combinations. We observe that when the inter-community phase coupling strength is sufficiently large, swarmalators settle in the static synchronization state. On the other hand, with a significant small phase coupling strength the state of anti-phase synchronization as well as chimera-like coexistence of sync and async are realized. Apart from rigorous numerical results, we have been successful to provide semi-analytical treatment for the existence and stability of global static sync and the anti-phase sync states., Comment: Accepted for publication in Physical Review E (2023)

Higher-order interactions shape collective dynamics, but how they affect transitions between different states in swarmalator systems is yet to be determined. To that effect, we here study an analytically tractable swarmalator model that incorporates both pairwise and higher-order interactions, resulting in four distinct collective states: async, phase wave, mixed, and sync states. We show that even a minute fraction of higher-order interactions induces abrupt transitions from the async state to the phase wave and the sync state. We also show that higher-order interactions facilitate an abrupt transition from the phase wave to the sync state by bypassing the intermediate mixed state. Moreover, elevated levels of higher-order interactions can sustain the presence of phase wave and sync state, even when pairwise interactions lean towards repulsion. The insights gained from these findings unveil self-organizing processes that hold the potential to explain sudden transitions between various collective states in numerous real-world systems., Comment: 14 pages, 8 figures

We present a formula for determining synchronizability in large, randomized and weighted simplicial complexes. This formula leverages eigenratios and costs to assess complete synchronizability under diverse network topologies and intensity distributions. We systematically vary coupling strengths (pairwise and three-body), degree and intensity distributions to identify the synchronizability of these simplicial complexes of the identical oscillators with natural coupling. We focus on randomized weighted connections with diffusive couplings and check synchronizability for different cases. For all these scenarios, eigenratios and costs reliably gauge synchronizability, eliminating the need for explicit connectivity matrices and eigenvalue calculations. This efficient approach offers a general formula for manipulating synchronizability in diffusively coupled identical systems with higher-order interactions simply by manipulating degrees, weights, and coupling strengths. We validate our findings with simplicial complexes of R\"ossler oscillators and confirm that the results are independent of the number of oscillators, connectivity components and distributions of degrees and intensities. Finally, we validate the theory by considering a real-world connection topology using chaotic R\"ossler oscillators., Comment: 11 pages, 8 figures, accepted for publication in Phys. Rev. E

Synchronization has received a lot of attention from the scientific community for systems evolving on static networks or higher-order structures, such as hypergraphs and simplicial complexes. In many relevant real world applications, the latter are not static but do evolve in time, in this paper we thus discuss the impact of the time-varying nature of high-order structures in the emergence of global synchronization. To achieve this goal we extend the master stability formalism to account, in a general way, for the additional contributions arising from the time evolution of the higher-order structure supporting the dynamical systems. The theory is successfully challenged against two illustrative examples, the Stuart-Landau nonlinear oscillator and the Lorenz chaotic oscillator.

We present a simple model of driven matter in a 1D medium with pinning impurities, applicable to magnetic domains walls, confined colloids, and other systems. We find rich dynamics, including hysteresis, reentrance, quasiperiodicity, and two distinct routes to chaos. In contrast to other minimal models of driven matter, the model is solvable: we derive the full phase diagram for small $N$, and for large $N$, derive expressions for order parameters and several bifurcation curves. The model is also realistic. Its collective states match those seen in the experiments of magnetic domain walls, and its force-velocity curve imitates those of superconductor vortices.

Since Galileo's time, the pendulum has evolved into one of the most exciting physical objects in mathematical modeling due to its vast range of applications for studying various oscillatory dynamics, including bifurcations and chaos, under various interests. This well-deserved focus aids in comprehending various oscillatory physical phenomena that can be reduced to the equations of the pendulum. The present article focuses on the rotational dynamics of the two-dimensional forced damped pendulum under the influence of the ac and dc torque. Interestingly, we are able to detect a range of the pendulum's length for which the angular velocity exhibits a few intermittent extreme rotational events that deviate significantly from a certain well-defined threshold. The statistics of the return intervals between these extreme rotational events are supported by our data to be spread exponentially. The numerical results show a sudden increase in the size of the chaotic attractor due to interior crisis which is the source of instability that is responsible for triggering large amplitude events in our system. We also notice the occurrence of phase slips with the appearance of extreme rotational events when phase difference between the instantaneous phase of the system and the externally applied ac torque is observed., Comment: 10 pages, 7 figures, Comments are welcome

Predator prey interactions are one of ecology's central research themes, but with many interdisciplinary implications across the social and natural sciences. Here we consider an often-overlooked species in these interactions, namely parasites. We first show that a simple predator prey parasite model, inspired by the classical Lotka Volterra equations, fails to produce a stable coexistence of all three species, thus failing to provide a biologically realistic outcome. To improve this, we introduce free space as a relevant eco-evolutionary component in a new mathematical model that uses a game-theoretical payoff matrix to describe a more realistic setup. We then show that the consideration of free space stabilizes the dynamics by means of cyclic dominance that emerges between the three species. We determine the parameter regions of coexistence as well as the types of bifurcations leading to it by means of analytical derivations as well as by means of numerical simulations. We conclude that the consideration of free space as a finite resource reveals the limits of biodiversity in predator prey parasite interactions, and it may also help us in the determination of factors that promote a healthy biota., Comment: 14 pages, 6 figures, Supplementary material related to this article can be found online at https://doi.org/10.1016/j.jtbi.2023.111446

With synchronization being one of nature's most ubiquitous collective behaviors, the field of network synchronization has experienced tremendous growth, leading to significant theoretical developments. However, most of these previous studies consider uniform connection weights and undirected networks with positive coupling. In the present article, we incorporate the asymmetry in a two-layer multiplex network by assigning the ratio of the adjacent nodes' degrees as the weights to the intralayer edges. Despite the presence of degree-biased weighting mechanism and attractive-repulsive coupling strengths, we are able to find the necessary conditions for intralayer synchronization and interlayer antisynchronization and test whether these two macroscopic states can withstand demultiplexing in a network. During the occurrence of these two states, we analytically calculate the oscillator's amplitude. In addition to deriving the local stability conditions for interlayer antisynchronization via the master stability function approach, we also construct a suitable Lyapunov function to determine a sufficient condition for global stability. We provide numerical evidence to show the necessity of negative interlayer coupling strength for the occurrence of antisynchronization, and such repulsive interlayer coupling coefficients can not destroy intralayer synchronization., Comment: 16 pages, 5 figures (Accepted for publication in the journal Physical Review E)

The stability analysis of synchronization in time-varying higher-order networked structures (simplicial complexes) is one of the challenging problem due to the presence of time-varying group interactions. In this context, most of the previous studies have been done either on temporal pairwise networks or on static simplicial complexes. Here, for the first time, we propose a general framework to study the synchronization phenomenon in temporal simplicial complexes. We show that the synchronous state exists as an invariant solution and obtain the necessary condition for it to be emerged as a stable state in fast switching regime. We prove that the time-averaged simplicial complex plays the role of synchronization indicator whenever the switching among simplicial topologies are adequately fast. We attempt to transform the stability problem into a master stability function form. Unfortunately, for the general circumstances, the dimension reduction of the master stability equation is cumbersome due to the presence of group interactions. However, we overcome this difficulty in two interesting situations based on either the functional forms of the coupling schemes or the connectivity structure of the simplicial complex, and demonstrate that the necessary condition mimics the form of a master stability function in these cases. We verify our analytical findings by applying them on synthetic and real-world networked systems. In addition, our results also reveal that with sufficient higher-order coupling and adequately fast rewiring, the temporal simplicial complex achieves synchrony even in a very low connectivity regime., Comment: 18 pages, 9 figures, any comments are welcome

A long-standing expectation is that two repulsively coupled oscillators tend to oscillate in opposite directions. It has been difficult to achieve complete synchrony in coupled identical oscillators with purely repulsive coupling. Here, we introduce a general coupling condition based on the linear matrix of dynamical systems for the emergence of the complete synchronization in pure repulsively coupled oscillators. The proposed coupling profiles (coupling matrices) define a bidirectional cross-coupling link that plays the role of indicator for the onset of complete synchrony between identical oscillators. We illustrate the proposed coupling scheme on several paradigmatic two-coupled chaotic oscillators and validate its effectiveness through the linear stability analysis of the synchronous solution based on the master stability function approach. We further demonstrate that the proposed general condition for the selection of coupling profiles to achieve synchronization even works perfectly for a large ensemble of oscillators., Comment: Accepted for publication in Physical Review E

The role of topological heterogeneity in the origin of extreme events in a network is investigated here. The dynamics of the oscillators associated with the nodes are assumed to be identical and influenced by mean-field repulsive interactions. An interplay of topological heterogeneity and the repulsive interaction between the dynamical units of the network triggers extreme events in the nodes when each node succumbs to such events for discretely different ranges of repulsive coupling. A high degree node is vulnerable to weaker repulsive interactions, while a low degree node is susceptible to stronger interactions. As a result, the formation of extreme events changes position with increasing strength of repulsive interaction from high to low degree nodes. Extreme events at any node are identified with the appearance of occasional large-amplitude events (amplitude of the temporal dynamics) that are larger than a threshold height and rare in occurrence, which we confirm by estimating the probability distribution of all events. Extreme events appear at any oscillator near the boundary of transition from rotation to libration at a critical value of the repulsive coupling strength. To explore the phenomenon, a paradigmatic second-order phase model is used to represent the dynamics of the oscillator associated with each node. We make an annealed network approximation to reduce our original model and thereby confirm the dual role of the repulsive interaction and the degree of a node in the origin of extreme events in any oscillator., Comment: 10 pages, 5 figures, Accepted for publication in Chaos: An Interdisciplinary Journal of Nonlinear Science

We study a population of swarmalators (swarming/mobile oscillators) which run on a ring and are subject to random pinning. The pinning represents the tendency of particles to stick to defects in the underlying medium which competes with the tendency to sync / swarm. The result is rich collective behavior. A highlight is low dimensional chaos which in systems of ordinary, Kuramoto-type oscillators is uncommon. Some of the states (the phase wave and split phase wave) resemble those seen in systems of Janus matchsticks or Japanese tree frogs. The others (such as the sync and unsteady states) may be observable in systems of vinegar eels, electrorotated Quincke rollers, or other swarmalators moving in disordered environments., Comment: Accepted for publication in Physical Review E (2023)

Understanding how the interplay between higher-order and multilayer structures of interconnections influences the synchronization behaviors of dynamical systems is a feasible problem of interest, with possible application in essential topics such as neuronal dynamics. Here, we provide a comprehensive approach for analyzing the stability of the complete synchronization state in simplicial complexes with numerous interaction layers. We show that the synchronization state exists as an invariant solution and derive the necessary condition for a stable synchronization state in presence of general coupling functions. It generalizes the well-known master stability function scheme to the higher-order structures with multiple interaction layers. We verify our theoretical results by employing them on networks of paradigmatic R\"{o}ssler oscillators and Sherman neuronal models, and demonstrate that the presence of group interactions considerably improves the synchronization phenomenon in the multilayer framework., Comment: 14 pages, 5 figures; Accepted for publication in Physical Review E (2022)

The persistence of biodiversity of species is a challenging proposition in ecological communities in the face of Darwinian selection. The present article investigates beyond the pairwise competitive interactions and provides a novel perspective for understanding the influence of higher-order interactions on the evolution of social phenotypes. Our simple model yields a prosperous outlook to demonstrate the impact of perturbations on intransitive competitive higher-order interactions. Using a mathematical technique, we show how alone the perturbed interaction network can quickly determine the coexistence equilibrium of competing species instead of solving a large system of ordinary differential equations. It is possible to split the system into multiple feasible cluster states depending on the number of perturbations. Our analysis also reveals the ratio between the unperturbed and perturbed species is inversely proportional to the amount of employed perturbation. Our results suggest that nonlinear dynamical systems and interaction topologies can be interplayed to comprehend species' coexistence under adverse conditions. Particularly our findings signify that less competition between two species increases their abundance and outperforms others., Comment: 17 pages, 10 figures

Designing navigation strategies for search time optimization remains of interest in various interdisciplinary branches in science. In here, we focus on microscopic self-propelled searchers namely active Brownian walkers in noisy and confined environment which are mediated by one such autonomous strategy namely resetting. As such, resetting stops the motion and compels the walkers to restart from the initial configuration intermittently according to an external timer that does not require control by the walkers. In particular, the resetting coordinates are either quenched (fixed) or annealed (fluctuating) over the entire topography. Although the strategy relies upon simple rules, it shows a significant ramification on the search time statistics in contrast to the original search. We show that the resetting driven protocols mitigate the performance of these active searchers based, robustly, on the inherent search time fluctuations. Notably, for the annealed condition, resetting is always found to expedite the search process. These features, as well as their applicability to more general optimization problems starting from queuing systems, computer science to living systems, make resetting based strategies universally promising.

Swarmalators have emerged as a new paradigm for dynamical collective behavior of multi-agent systems due to the interplay of synchronization and swarming that they inherently incorporate. Their dynamics have been explored with different coupling topologies, interaction functions, external forcing, noise, competitive interactions, and from other important viewpoints. Here we take a systematic approach and review the collective dynamics of swarmalators analytically and/or numerically. Long-term states of position aggregation and phase synchronization are revealed in this perspective with some future problems., Comment: 7 pages, 4 figures, comments welcome

Mounting evidence in recent years suggests that astrocytes, a sub-type of glial cells, not only serve metabolic and structural support for neurons and synapses but also play critical roles in regulation of proper functioning of the nervous system. In this work, we investigate the effect of astrocyte on the spontaneous firing activity of a neuron through a combined model which includes a neuron-astrocyte pair. First, we show that an astrocyte may provide a kind of multistability in neuron dynamics by inducing different firing modes such as random and bursty spiking. Then, we identify the underlying mechanism of this behavior and search for the astrocytic factors that may have regulatory roles in different firing regimes. More specifically, we explore how an astrocyte can participate in occurrence and control of spontaneous irregular spiking activity of a neuron in random spiking mode. Additionally, we systematically investigate the bursty firing regime dynamics of the neuron under the variation of biophysical facts related to the intracellular environment of the astrocyte. It is found that an astrocyte coupled to a neuron can provide a control mechanism for both spontaneous firing irregularity and burst firing statistics, i.e., burst regularity and size., Comment: 11 pages, 11 figures; accepted for publication in Chaos: An Interdisciplinary Journal of Nonlinear Science, 2022

Since the discovery of chimera states, the presence of a nonzero phase lag parameter turns out to be an essential attribute for the emergence of chimeras in a nonlocally coupled identical Kuramoto phase oscillators' network with pairwise interactions. In this letter, we report the emergence of chimeras without phase lag in nonlocally coupled identical Kuramoto network owing to the introduction of non-pairwise interactions. The influence of added nonlinearity in the coupled system dynamics in the form of simplicial complexes mitigates the requisite of a nonzero phase lag parameter for the emergence of chimera states. Chimera states stimulated by the reciprocity of the pairwise and non-pairwise interaction strengths and their multistable nature are characterized with appropriate measures and are demonstrated in the parameter spaces., Comment: 8 pages, 8 figures; accepted for publication in Physical Review E (Letters), 2022

Recent developments in complex systems have witnessed that many real-world scenarios, successfully represented as networks are not always restricted to binary interactions but often include higher-order interactions among the nodes. These beyond pairwise interactions are preferably modeled by hypergraphs, where hyperedges represent higher-order interactions between a set of nodes. In this work, we consider a multiplex network where the intralayer connections are represented by hypergraphs, called multiplex hypergraph. The hypergraph is constructed by mapping the maximal cliques of a scale-free network to hyperedges of suitable sizes. We investigate the intralayer and interlayer synchronization of such multiplex structures. Our study unveils that the intralayer synchronization appreciably enhances when the higher-order structure is taken into consideration inspite of only pairwise connections. We derive the necessary condition for stable synchronization states by the master stability function approach, which perfectly agrees with the numerical results. We also explore the robustness of interlayer synchronization and find that for the multiplex structures with many-body interaction, the interlayer synchronization is more persistent than multiplex networks with pairwise interaction., Comment: 10 pages, 7 figures; Accepted for publication in Chaos

Network science has evolved into an indispensable platform for studying complex systems. But recent research has identified limits of classical networks, where links connect pairs of nodes, to comprehensively describe group interactions. Higher-order networks, where a link can connect more than two nodes, have therefore emerged as a new frontier in network science. Since group interactions are common in social, biological, and technological systems, higher-order networks have recently led to important new discoveries across many fields of research. We here review these works, focusing in particular on the novel aspects of the dynamics that emerges on higher-order networks. We cover a variety of dynamical processes that have thus far been studied, including different synchronization phenomena, contagion processes, the evolution of cooperation, and consensus formation. We also outline open challenges and promising directions for future research., Comment: 15 pages, 6 figures; accepted for publication in Journal of the Royal Society Interface

Swarmalators are entities with the simultaneous presence of swarming and synchronization that reveal emergent collective behavior due to the fascinating bidirectional interplay between phase and spatial dynamics. Although different coupling topologies have already been considered, here we introduce time-varying competitive phase interaction among swarmalators where the underlying connectivity for attractive and repulsive coupling varies depending on the vision (sensing) radius. Apart from investigating some fundamental properties like conservation of center of position and collision avoidance, we also scrutinize the cases of extreme limits of vision radius. The concurrence of attractive-repulsive competitive phase coupling allows the exploration of diverse asymptotic states, like static $\pi$, and mixed phase wave states, and we explore the feasible routes of those states through a detailed numerical analysis. In sole presence of attractive local coupling, we reveal the occurrence of static cluster synchronization where the number of clusters depends crucially on the initial distribution of positions and phases of each swarmalator. In addition, we analytically calculate the sufficient condition for the emergence of the static synchronization state. We further report the appearance of the static ring phase wave state and evaluate its radius theoretically. Finally, we validate our findings using Stuart-Landau oscillators to describe the phase dynamics of swarmalators subject to attractive local coupling., Comment: 21 pages, 12 figures; accepted for publication in New Journal of Physics

The stability analysis of synchronization patterns on generalized network structures is of immense importance nowadays. In this article, we scrutinize the stability of intralayer synchronous state in temporal multilayer hypernetworks, where each dynamic units in a layer communicate with others through various independent time-varying connection mechanisms. Here, dynamical units within and between layers may be interconnected through arbitrary generic coupling functions. We show that intralayer synchronous state exists as an invariant solution. Using fast switching stability criteria, we derive the condition for stable coherent state in terms of associated time-averaged network structure, and in some instances we are able to separate the transverse subspace optimally. Using simultaneous block diagonalization of coupling matrices, we derive the synchronization stability condition without considering time-averaged network structure. Finally, we verify our analytically derived results through a series of numerical simulations on synthetic and real-world neuronal networked systems., Comment: 18 pages, 10 figures

Extreme events gain the attention of researchers due to their utmost importance in various contexts ranging from finance to climatology. This brings such recurrent events to the limelight of attention in interdisciplinary research. A comprehensive review of recent progress is provided to capture recent improvements in analyzing such very high-amplitude events from the point of view of dynamical systems and random walkers. We emphasize, in detail, the mechanisms responsible for the emergence of such events in complex systems. Several mechanisms that contribute to the occurrence of extreme events have been elaborated that investigate the sources of instabilities leading to them. In addition, we discuss the prediction of extreme events from two different contexts, using dynamical instabilities and data-based machine learning algorithms. Tracking of instabilities in the phase space is not always feasible and precise knowledge of the dynamics of extreme events does not necessarily help in forecasting extreme events. Moreover, in most studies on high-dimensional systems, only a few degrees of freedom participate in extreme events' formation. Thus, the notable inclusion of prediction through machine learning is of enormous significance, particularly for those cases where the governing equations of the model are explicitly unavailable. Besides, random walks on complex networks can represent several transport processes, and exceedances of the flux of walkers above a prescribed threshold may describe extreme events. We unveil the theoretical studies on random walkers with their enormous potential for applications in reducing extreme events. We cover the possible controlling strategies, which may be helpful to mitigate extreme events in physical situations like traffic jams, heavy load of web requests, competition for shared resources, floods in the network of rivers, and many more., Comment: 78 pages, 30 figures Accepted for publication in Physics Reports on 6th April 2022

Over the past two decades, complex network theory provided the ideal framework for investigating the intimate relationships between the topological properties characterizing the wiring of connections among a system's unitary components and its emergent synchronized functioning. An increased number of setups from the real world found therefore a representation in term of graphs, while more and more sophisticated methods were developed with the aim of furnishing a realistic description of the connectivity patterns under study. In particular, a significant number of systems in physics, biology and social science features a time-varying nature of the interactions among their units. We here give a comprehensive review of the major results obtained by contemporary studies on the emergence of synchronization in time-varying networks. In particular, two paradigmatic frameworks will be described in details. The first encompasses those systems where the time dependence of the nodes' connections is due to adaptation, external forces, or any other process affecting each of the links of the network. The second framework, instead, corresponds to the case in which the structural evolution of the graph is due to the movement of the nodes, or agents, in physical spaces and to the fact that interactions may be ruled by space-dependent laws in a way that connections are continuously switched on and off in the course of the time. Finally, our report ends with a short discussion on promising directions and open problems for future studies., Comment: 76 pages, 34 figures

Understanding and predicting uncertain things are the central themes of scientific evolution. Human beings revolve around these fears of uncertainties concerning various aspects like a global pandemic, health, finances, to name but a few. Dealing with this unavoidable part of life is far tougher due to the chaotic nature of these unpredictable activities. In the present article, we consider a global network of identical chaotic maps, which splits into two different clusters, despite the interaction between all nodes are uniform. The stability analysis of the spatially homogeneous chaotic solutions provides a critical coupling strength, before which we anticipate such partial synchronization. The distance between these two chaotic synchronized populations often deviates more than eight times of standard deviation from its long-term average. The probability density function of these highly deviated values fits well with the Generalized Extreme Value distribution. Meanwhile, the distribution of recurrence time intervals between extreme events resembles the Weibull distribution. The existing literature helps us to characterize such events as extreme events using the significant height. These extremely high fluctuations are less frequent in terms of their occurrence. We determine numerically a range of coupling strength for these extremely large but recurrent events. On-off intermittency is the responsible mechanism underlying the formation of such extreme events. Besides understanding the generation of such extreme events and their statistical signature, we furnish forecasting these events using the powerful deep learning algorithms of an artificial recurrent neural network. This Long Short-Term Memory (LSTM) can offer handy one-step forecasting of these chaotic intermittent bursts. We also ensure the robustness of this forecasting model with two hundred hidden cells in each LSTM layer., Comment: Accepted for publication in Journal of Physics: Complexity (22 pages, 6 figures)

In the present article, we demonstrate the emergence and existence of the spiral wave chimera-like transient pattern in coupled ecological systems, composed of prey-predator patches, where the patches are connected in a three-dimensional medium through local diffusion. We explore the transition scenarios among the several collective dynamical behaviors together with transient spiral wave chimera-like states and investigate the long time behavior of these states. The transition from the transient spiral chimera-like pattern to the long time synchronized or desynchronized pattern appears through the deformation of the incoherent region of the spiral core. We discuss the transient dynamics under the influence of the species diffusion at different time instants. By calculating the instantaneous strength of incoherence of the populations, we estimate the duration of the transient dynamics characterized by the persistence of the chimera-like spatial coexistence of coherent and incoherent patterns over the spatial domain. We generalize our observations on the transient dynamics in three-dimensional grid of diffusive ecological systems by considering two different prey-predator systems., Comment: Accepted for publication in Chaos: An Interdisciplinary Journal of Nonlinear Science

Emerging collective behavior in complex dynamical networks depends on both coupling function and underlying coupling topology. Through this perspective, we provide a brief yet profound excerpt of recent research efforts that explore how the synergy of attractive and repulsive interactions influence the destiny of ensembles of interacting dynamical systems. We review the incarnation of collective states ranging from chimera or solitary states to extreme events and oscillation quenching arising as a result of different network arrangements. Though the existing literature demonstrates that many of the crucial developments have been made, nonetheless, we come up with significant routes of further research in this field of study.

Ecology and evolution are inherently linked, and studying a mathematical model that considers both holds promise of insightful discoveries related to the dynamics of cooperation. In the present article, we use the prisoner's dilemma (PD) game as a basis for long-term apprehension of the essential social dilemma related to cooperation among unrelated individuals. We upgrade the contemporary PD game with an inclusion of evolution-induced act of punishment as a third competing strategy in addition to the traditional cooperators and defectors. In a population structure, the abundance of ecologically-viable free space often regulates the reproductive opportunities of the constituents. Hence, additionally, we consider the availability of free space as an ecological footprint, thus arriving at a simple eco-evolutionary model, which displays fascinating complex dynamics. As possible outcomes, we report the individual dominance of cooperators and defectors as well as a plethora of mixed states, where different strategies coexist followed by maintaining the diversity in a socio-ecological framework. These states can either be steady or oscillating, whereby oscillations are sustained by cyclic dominance among different combinations of cooperators, defectors, and punishers. We also observe a novel route to cyclic dominance where cooperators, punishers, and defectors enter a coexistence via an inverse Hopf bifurcation that is followed by an inverse period doubling route.

The remarkable flexibility and adaptability of both deep learning models and ensemble methods have led to the proliferation for their application in understanding many physical phenomena. Traditionally, these two techniques have largely been treated as independent methodologies in practical applications. This study develops an optimized ensemble deep learning (OEDL) framework wherein these two machine learning techniques are jointly used to achieve synergistic improvements in model accuracy, stability, scalability, and reproducibility prompting a new wave of applications in the forecasting of dynamics. Unpredictability is considered as one of the key features of chaotic dynamics, so forecasting such dynamics of nonlinear systems is a relevant issue in the scientific community. It becomes more challenging when the prediction of extreme events is the focus issue for us. In this circumstance, the proposed OEDL model based on a best convex combination of feed-forward neural networks, reservoir computing, and long short-term memory can play a key role in advancing predictions of dynamics consisting of extreme events. The combined framework can generate the best out-of-sample performance than the individual deep learners and standard ensemble framework for both numerically simulated and real world data sets. We exhibit the outstanding performance of the OEDL framework for forecasting extreme events generated from Lienard-type system, prediction of COVID-19 cases in Brazil, dengue cases in San Juan, and sea surface temperature in Nino 3.4 region., Comment: 14 pages, 8 figures, any comments are welcome

We study the emergence of the traveling chimera state in a two-dimensional network of Hindmarsh-Rose burst neurons with the mutual presence of local and non-local couplings. We show that in the unique presence of the non-local chemical coupling modeled by a nonlinear function, the traveling chimera phenomenon occurs with a displacement in both directions of the plane of the grid. The introduction of local electrical coupling shows that the mutual influence of the two types of coupling can, for certain values, generate traveling chimera, imperfect-traveling, traveling multi-clusters, and alternating traveling chimera, ie the presence in the network under study, of patterns of coherent elements interspersed by other incoherent elements in movement and alternately changing their position over time. The confirmation of the states of coherence is done by introducing the parameter of instantaneous local order parameter in two dimensions. We extend our analysis through mathematical tools such as the Hamilton energy function to determine the direction of propagation of patterns in two dimensions., Comment: 8 pages, 8 figures, Accepted for publication in Physics Letters A (2021)

Formation of diverse patterns in spatially extended reaction-diffusion systems is an important aspect of study which is pertinent to many chemical and biological processes. Of special interest is the peculiar phenomenon of chimera state having spatial coexistence of coherent and incoherent dynamics in a system of identically interacting individuals. In the present article, we report the emergence of various collective dynamical patterns while considering a system of prey-predator dynamics in presence of a two-dimensional diffusive environment. Particularly, we explore the observance of four distinct categories of spatial arrangements among the species, namely spiral wave, spiral chimera, completely synchronized oscillations, and oscillation death states in a broad region of the diffusion-driven parameter space. Emergence of amplitude mediated spiral chimera states displaying drifted amplitudes and phases in the incoherent subpopulation is detected for parameter values beyond both Turing and Hopf bifurcations. Transition scenarios among all these distinguishable patterns are numerically demonstrated for a wide range of the diffusion coefficients which reveal that the chimera states arise during the transition from oscillatory to steady state dynamics. Furthermore, we characterize the occurrence of each of the recognizable patterns by estimating the strength of incoherent subpopulations in the two-dimensional space., Comment: 15 pages, 13 figures, Accepted for publication in Phys. Rev. E (2021)

Most previous studies on coupled dynamical systems assume that all interactions between oscillators take place uniformly in time, but in reality, this does not necessarily reflect the usual scenario. The heterogeneity in the timings of such interactions strongly influences the dynamical processes. Here, we introduce a time-evolving state-space dependent coupling among an ensemble of identical coupled oscillators, where individual units are interacting only when the mean state of the system lies within a certain proximity of the phase space. They interact globally with mean-field diffusive coupling in a certain vicinity and behave like uncoupled oscillators with self-feedback in the remaining complementary subspace. Interestingly due to this occasional interaction, we find that the system shows an abrupt explosive transition from oscillatory to death state. Further, in the explosive death transitions, the oscillatory state and the death state coexist over a range of coupling strengths near the transition point. We explore our claim using Van der pol, FitzHughNagumo and Lorenz oscillators with dynamic mean field interaction. The dynamic interaction mechanism can explain sudden suppression of oscillations and concurrence of oscillatory and steady state in biological as well as technical systems.

The role of a new form of dynamic interaction is explored in a network of generic identical oscillators. The proposed design of dynamic coupling facilitates the onset of a plethora of asymptotic states including synchronous states, amplitude death states, oscillation death states, a mixed state (complete synchronized cluster and small amplitude unsynchronized domain), and bistable states (coexistence of two attractors). The dynamical transitions from the oscillatory to death state are characterized using an average temporal interaction approximation, which agrees with the numerical results in temporal interaction. A first-order phase transition behavior may change into a second-order transition in spatial dynamic interaction solely depending on the choice of initial conditions in the bistable regime. However, this possible abrupt first-order like transition is completely non-existent in the case of temporal dynamic interaction. Besides the study on periodic Stuart-Landau systems, we present results for paradigmatic chaotic model of R\"ossler oscillators and Mac-arthur ecological model., Comment: Accepted for publication in Chaos: An Interdisciplinary Journal of Nonlinear Science, 2021

We propose a deterministic compartmental model of infectious disease which considers the test-kits as an important ingredient for the suppression and mitigation of epidemics. A rigorous simulation (with analytical argument) is provided to reveal the effective reduction of final outbreak size and peak of infection as a function of basic reproduction number in a single patch. Further, to study the impact of long and short-distance human migration among the patches, we have considered heterogeneous networks where the linear diffusive connectivity is determined by the network link structure. We numerically confirm that implementation of test-kits in the fraction of nodes (patches) having larger degrees or betweenness centralities can reduce the peak of infection (as well as final outbreak size) significantly. A next-generation matrix based analytical treatment is provided to find out the critical transmission probability in the entire network for the onset of epidemics. Finally, the optimal intervention strategy is validated in two real networks: global airport networks and transportation networks of Kolkata, India., Comment: 22 pages, 9 Figures

Interaction within an ensemble of coupled nonlinear oscillators induces a variety of collective behaviors. One of the most fascinating is a chimera state which manifests the coexistence of spatially distinct populations of coherent and incoherent elements. Understanding of the emergent chimera behavior in controlled experiments or real systems requires a focus on the consideration of heterogeneous network models. In this study, we explore the transitions in a heterogeneous Kuramoto model under the monotonical increase of the coupling strength and specifically find that this system exhibits a frequency-modulated chimera-like pattern during the explosive transition to synchronization. We demonstrate that this specific dynamical regime originates from the interplay between (the evolved) attractively and repulsively coupled subpopulations. We also show that the above mentioned chimera-like state is induced under weakly non-local, small-world and sparse scale-free coupling and suppressed in globally coupled, strongly rewired and dense scale-free networks due to the emergence of the large-scale connections.

Diverse collective dynamics emerge in dynamical systems interacting on top of complex network architectures. Along this line of research, temporal network has come out to be one of the most promising network platforms to investigate. Especially, such network with spatially moving agents has been established to be capable of modelling a number of practical instances. In this paper, we examine the dynamical outcomes of moving agents interacting based upon their physical proximity. For this, we particularly emphasize on the impact of competing interactions among the agents depending on their physical distance. We specifically assume attractive coupling between agents which are staying apart from each other, whereas we adopt repulsive interaction for agents that are sufficiently close in space. With this set-up, we consider two types of coupling configurations, symmetry-breaking and symmetry-preserving couplings. We encounter variants of collective dynamics ranging from synchronization, inhomogeneous small oscillation to cluster state and extreme events while changing the attractive and repulsive coupling strengths. We have been able to map all these dynamical behaviors in the coupling parameter space. Complete synchronization being the most desired state in absence of repulsive coupling, we present an analytical study for this scenario that agrees well with the numerical results., Comment: Accepted for publication in IEEE Transactions on Network Science and Engineering, 2020