1
|
Dutta S, Mondal A, Kundu P, Khanra P, Pal P, Hens C. Impact of phase lag on synchronization in frustrated Kuramoto model with higher-order interactions. Phys Rev E 2023; 108:034208. [PMID: 37849147 DOI: 10.1103/physreve.108.034208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Accepted: 08/25/2023] [Indexed: 10/19/2023]
Abstract
The study of first order transition (explosive synchronization) in an ensemble (network) of coupled oscillators has been the topic of paramount interest among the researchers for more than one decade. Several frameworks have been proposed to induce explosive synchronization in a network and it has been reported that phase frustration in a network usually suppresses first order transition in the presence of pairwise interactions among the oscillators. However, on the contrary, by considering networks of phase frustrated coupled oscillators in the presence of higher-order interactions (up to 2-simplexes) we show here, under certain conditions, phase frustration can promote explosive synchronization in a network. A low-dimensional model of the network in the thermodynamic limit is derived using the Ott-Antonsen ansatz to explain this surprising result. Analytical treatment of the low-dimensional model, including bifurcation analysis, explains the apparent counter intuitive result quite clearly.
Collapse
Affiliation(s)
- Sangita Dutta
- Department of Mathematics, National Institute of Technology, Durgapur 713209, India
| | - Abhijit Mondal
- Department of Mathematics, National Institute of Technology, Durgapur 713209, India
| | - Prosenjit Kundu
- Dhirubhai Ambani Institute of Information and Communication Technology, Gandhinagar, Gujarat 382007, India
| | - Pitambar Khanra
- Department of Mathematics, State University of New York at Buffalo, Buffalo 14260, USA
| | - Pinaki Pal
- Department of Mathematics, National Institute of Technology, Durgapur 713209, India
| | - Chittaranjan Hens
- Center for Computational Natural Science and Bioinformatics, International Institute of Informational Technology, Gachibowli, Hyderabad 500032, India
| |
Collapse
|
2
|
Dutta S, Kundu P, Khanra P, Hens C, Pal P. Perfect synchronization in complex networks with higher-order interactions. Phys Rev E 2023; 108:024304. [PMID: 37723785 DOI: 10.1103/physreve.108.024304] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 07/11/2023] [Indexed: 09/20/2023]
Abstract
Achieving perfect synchronization in a complex network, specially in the presence of higher-order interactions (HOIs) at a targeted point in the parameter space, is an interesting, yet challenging task. Here we present a theoretical framework to achieve the same under the paradigm of the Sakaguchi-Kuramoto (SK) model. We analytically derive a frequency set to achieve perfect synchrony at some desired point in a complex network of SK oscillators with higher-order interactions. Considering the SK model with HOIs on top of the scale-free, random, and small world networks, we perform extensive numerical simulations to verify the proposed theory. Numerical simulations show that the analytically derived frequency set not only provides stable perfect synchronization in the network at a desired point but also proves to be very effective in achieving a high level of synchronization around it compared to the other choices of frequency sets. The stability and the robustness of the perfect synchronization state of the system are determined using the low-dimensional reduction of the network and by introducing a Gaussian noise around the derived frequency set, respectively.
Collapse
Affiliation(s)
- Sangita Dutta
- Department of Mathematics, National Institute of Technology, Durgapur 713209, India
| | - Prosenjit Kundu
- Dhirubhai Ambani Institute of Information and Communication Technology, Gandhinagar, Gujarat 382007, India
| | - Pitambar Khanra
- Department of Mathematics, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - Chittaranjan Hens
- Center for Computational Natural Science and Bioinformatics, International Institute of Informational Technology, Gachibowli, Hyderabad 500032, India
| | - Pinaki Pal
- Department of Mathematics, National Institute of Technology, Durgapur 713209, India
| |
Collapse
|
3
|
Ghosh S, Khanra P, Kundu P, Ji P, Ghosh D, Hens C. Dimension reduction in higher-order contagious phenomena. CHAOS (WOODBURY, N.Y.) 2023; 33:2893033. [PMID: 37229635 DOI: 10.1063/5.0152959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Accepted: 05/01/2023] [Indexed: 05/27/2023]
Abstract
We investigate epidemic spreading in a deterministic susceptible-infected-susceptible model on uncorrelated heterogeneous networks with higher-order interactions. We provide a recipe for the construction of one-dimensional reduced model (resilience function) of the N-dimensional susceptible-infected-susceptible dynamics in the presence of higher-order interactions. Utilizing this reduction process, we are able to capture the microscopic and macroscopic behavior of infectious networks. We find that the microscopic state of nodes (fraction of stable healthy individual of each node) inversely scales with their degree, and it becomes diminished due to the presence of higher-order interactions. In this case, we analytically obtain that the macroscopic state of the system (fraction of infectious or healthy population) undergoes abrupt transition. Additionally, we quantify the network's resilience, i.e., how the topological changes affect the stable infected population. Finally, we provide an alternative framework of dimension reduction based on the spectral analysis of the network, which can identify the critical onset of the disease in the presence or absence of higher-order interactions. Both reduction methods can be extended for a large class of dynamical models.
Collapse
Affiliation(s)
- Subrata Ghosh
- Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700108, India
| | - Pitambar Khanra
- Department of Mathematics, State University of New York at Buffalo, Buffalo, New York 14260, USA
| | - Prosenjit Kundu
- Dhirubhai Ambani Institute of Information and Communication Technology, Gandhinagar, Gujarat 382007, India
| | - Peng Ji
- Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai 200433, China
- Key Laboratory of Computational Neuroscience and Brain-Inspired Intelligence (Fudan University), Ministry of Education, Shanghai 200433, China
| | - Dibakar Ghosh
- Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700108, India
| | - Chittaranjan Hens
- Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700108, India
- International Institute of Information Technology, Hyderabad 500 032, India
| |
Collapse
|
4
|
de Aguiar MAM. Generalized frustration in the multidimensional Kuramoto model. Phys Rev E 2023; 107:044205. [PMID: 37198798 DOI: 10.1103/physreve.107.044205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Accepted: 03/24/2023] [Indexed: 05/19/2023]
Abstract
The Kuramoto model describes how coupled oscillators synchronize their phases as the intensity of the coupling increases past a threshold. The model was recently extended by reinterpreting the oscillators as particles moving on the surface of unit spheres in a D-dimensional space. Each particle is then represented by a D-dimensional unit vector; for D=2 the particles move on the unit circle and the vectors can be described by a single phase, recovering the original Kuramoto model. This multidimensional description can be further extended by promoting the coupling constant between the particles to a matrix K that acts on the unit vectors. As the coupling matrix changes the direction of the vectors, it can be interpreted as a generalized frustration that tends to hinder synchronization. In a recent paper we studied in detail the role of the coupling matrix for D=2. Here we extend this analysis to arbitrary dimensions. We show that, for identical particles, when the natural frequencies are set to zero, the system converges either to a stationary synchronized state, given by one of the real eigenvectors of K, or to an effective two-dimensional rotation, defined by one of the complex eigenvectors of K. The stability of these states depends on the set eigenvalues and eigenvectors of the coupling matrix, which controls the asymptotic behavior of the system, and therefore, can be used to manipulate these states. When the natural frequencies are not zero, synchronization depends on whether D is even or odd. In even dimensions the transition to synchronization is continuous and rotating states are replaced by active states, where the module of the order parameter oscillates while it rotates. If D is odd the phase transition is discontinuous and active states can be suppressed for some distributions of natural frequencies.
Collapse
Affiliation(s)
- Marcus A M de Aguiar
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, Unicamp 13083-970, Campinas, São Paolo, Brazil
| |
Collapse
|
5
|
Kundu P, MacLaren NG, Kori H, Masuda N. Mean-field theory for double-well systems on degree-heterogeneous networks. Proc Math Phys Eng Sci 2022. [DOI: 10.1098/rspa.2022.0350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Many complex dynamical systems in the real world, including ecological, climate, financial and power-grid systems, often show critical transitions, or tipping points, in which the system’s dynamics suddenly transit into a qualitatively different state. In mathematical models, tipping points happen as a control parameter gradually changes and crosses a certain threshold. Tipping elements in such systems may interact with each other as a network, and understanding the behaviour of interacting tipping elements is a challenge because of the high dimensionality originating from the network. Here, we develop a degree-based mean-field theory for a prototypical double-well system coupled on a network with the aim of understanding coupled tipping dynamics with a low-dimensional description. The method approximates both the onset of the tipping point and the position of equilibria with a reasonable accuracy. Based on the developed theory and numerical simulations, we also provide evidence for multistage tipping point transitions in networks of double-well systems.
Collapse
Affiliation(s)
- Prosenjit Kundu
- Department of Mathematics, State University of New York at Buffalo, Buffalo, NY 14260-2900, USA
| | - Neil G. MacLaren
- Department of Mathematics, State University of New York at Buffalo, Buffalo, NY 14260-2900, USA
| | - Hiroshi Kori
- Department of Complexity Science and Engineering, The University of Tokyo, Chiba 277-8561, Japan
| | - Naoki Masuda
- Department of Mathematics, State University of New York at Buffalo, Buffalo, NY 14260-2900, USA
- Computational and Data-Enabled Science and Engineering Program, State University of New York at Buffalo, Buffalo, NY 14260-5030, USA
- Faculty of Science and Engineering, Waseda University, Tokyo 169-8555, Japan
| |
Collapse
|
6
|
Roy M, Poria S, Hens C. Assortativity-induced explosive synchronization in a complex neuronal network. Phys Rev E 2021; 103:062307. [PMID: 34271687 DOI: 10.1103/physreve.103.062307] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 05/25/2021] [Indexed: 11/07/2022]
Abstract
In this study, we consider a scale-free network of nonidentical Chialvo neurons, coupled through electrical synapses. For sufficiently strong coupling, the system undergoes a transition from completely out of phase synchronized to phase synchronized state. The principal focus of this study is to investigate the effect of the degree of assortativity over the synchronization transition process. It is observed that, depending on assortativity, bistability between two asymptotically stable states allows one to develop a hysteresis loop which transforms the phase transition from second order to first order. An expansion in the area of hysteresis loop is noticeable with increasing degree-degree correlation in the network. Our study also reveals that effective frequencies of nodes simultaneously go through a continuous or sudden transition to the synchronized state with the corresponding phases. Further, we examine the robustness of the results under the effect of network size and average degree, as well as diverse frequency setup. Finally, we investigate the dynamical mechanism in the process of generating explosive synchronization. We observe a significant impact of lower degree nodes behind such phenomena: in a positive assortative network the low degree nodes delay the synchronization transition.
Collapse
Affiliation(s)
- Mousumi Roy
- Department of Applied Mathematics, University of Calcutta, 92, A.P.C. Road, Kolkata 700009, India
| | - Swarup Poria
- Department of Applied Mathematics, University of Calcutta, 92, A.P.C. Road, Kolkata 700009, India
| | - Chittaranjan Hens
- Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700108, India
| |
Collapse
|
7
|
Urdapilleta E. Transition to synchronization in heterogeneous inhibitory neural networks with structured synapses. CHAOS (WOODBURY, N.Y.) 2021; 31:033151. [PMID: 33810717 DOI: 10.1063/5.0038896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/27/2020] [Accepted: 03/05/2021] [Indexed: 06/12/2023]
Abstract
Inhibitory neurons form an extensive network involved in the development of different rhythms in the cerebral cortex. A transition from an incoherent state, where all inhibitory neurons fire unrelated to each other, to a synchronized or locked state, where all or most neurons define a tight firing pattern, is maybe the most salient process to analyze when considering neuronal rhythms. In this work, we analyzed whether different patterns of effective synaptic connectivity may support a first-order-like transition in this path to synchronization. Such an "explosive" phenomenon may be relevant in neural processes, as normal cognitive processing in different tasks and some neurological disorders manifest an increased power in many neuronal rhythms, supported by an extended concerted spiking activity and an abrupt change to this state. Furthermore, we built an adaptive mechanism that supports the generation of this kind of network, which rapidly creates the underlying structure based on the ongoing firing statistics.
Collapse
Affiliation(s)
- Eugenio Urdapilleta
- Centro Atómico Bariloche and Instituto Balseiro, Comisión Nacional de Energía Atómica (CNEA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Universidad Nacional de Cuyo, Av. E. Bustillo 9500, R8402AGP San Carlos de Bariloche, Río Negro, Argentina
| |
Collapse
|