Abstract
Simplicial complexes are mathematical constructions that describe higher-order interactions within the interconnecting elements of a network. Such higher-order interactions become increasingly significant in neuronal networks since biological backgrounds and previous outcomes back them. In light of this, the current research explores a higher-order network of the memristive Rulkov model. To that end, the master stability functions are used to evaluate the synchronization of a network with pure pairwise hybrid (electrical and chemical) synapses alongside a network with two-node electrical and multi-node chemical connections. The findings provide good insight into the impact of incorporating higher-order interaction in a network. Compared to two-node chemical synapses, higher-order interactions adjust the synchronization patterns to lower multi-node chemical coupling parameter values. Furthermore, the effect of altering higher-order coupling parameter value on the dynamics of neurons in the synchronization state is researched. It is also shown how increasing network size can enhance synchronization by lowering the value of coupling parameters whereby synchronization occurs. Except for complete synchronization, cluster synchronization is detected for higher electrical coupling strength values wherein the neurons are out of the completed synchronization state.
1. Introduction
Understanding the complicated functions of the brain of an individual has been fascinating and challenging for scientists and academics (Changeux and Dehaene, 1989). One important component of this pursuit is the development of neural models capable of capturing the activity and functionality of individual neurons as well as the networks they compose (Aihara et al., 1990; Boccaletti et al., 2006; Majhi et al., 2019). Neuronal models provide a framework for investigating the brain's computational capacities, information processing, and cognitive function emergence (Ibarz et al., 2011). Over the years, several neuronal models have been proposed to describe and simulate the behavior of neurons, such as Hodgkin and Huxley (1952), FitzHugh (1961), Morris and Lecar (1981), and Hindmarsh and Rose (1984) models defined by differential equations and Chialvo (1995), Rulkov (2002), Izhikevich and Hoppensteadt (2004), and Zandi-Mehran et al. (2020) described by difference equations. The memristor has recently developed as a novel electrical component with enormous potential for neuronal modeling. In recent years, the memristor has emerged as a revolutionary electronic component with great potential for neuronal modeling (Lin et al., 2020; Ding et al., 2023; Xu et al., 2023). The potential of memristors to imitate crucial synaptic actions makes them significant for neuronal modeling. Synapses or connections amongst neurons, are critical in the brain's information analysis and learning process. They exhibit plasticity, indicating that their strength and efficiency can alter depending on the activity and timing of neuronal inputs. Neuronal models based on memristors, including memristive Hodgkin-Huxley (Hu and Liu, 2019), memristive FitzHugh-Nagumo (Chen et al., 2020), memristive Morris-Lecar (Bao et al., 2022), memristive Hindmarsh-Rose (Bao et al., 2018, 2021), memristive Rulkov (mRulkov) (Li et al., 2021), and memristive Chialvo (Vivekanandhan et al., 2023) neuron models, have enormous potential for various applications. They can, for instance, be used to study neurological disorders and replicate brain activity.
The study of complex systems has yielded essential insights into the structures and dynamics of numerous interconnected systems, ranging from social networks (Shahal et al., 2020) and biological systems (Ma and Tang, 2017) to network marketing and trading networks (Kim et al., 2006; Feng et al., 2014; Cho et al., 2023), in the discipline of network science. Networks have traditionally been depicted as graphs, with nodes and edges capturing entities and their pairwise connections (Burgio et al., 2020; Ghosh et al., 2022). Real-world systems, on the other hand, frequently demonstrate interactions and linkages that extend beyond simple pairwise connections (Ince et al., 2009; Alvarez-Rodriguez et al., 2021; Battiston et al., 2021). To capture the rich connectivity patterns in complex systems, researchers have turned to higher-order network representations (Majhi et al., 2022). Higher-order networks consider interactions among groups of nodes rather than just pairwise connections, enabling a more comprehensive understanding of complex systems (Carletti et al., 2020; Lotito et al., 2022). One powerful mathematical framework for modeling higher-order networks is simplicial complexes (Skardal and Arenas, 2020; Gambuzza et al., 2021). A simplicial complex is a mathematical structure that reflects the interactions between nodes at various granularities (Gambuzza et al., 2021). It expands the graph concept by integrating higher-order interactions such as triangles, tetrahedra, and higher-dimensional simplices (Ghorbanchian et al., 2021). Each simplex in a simplicial complex represents a group of nodes connected in a particular way. A triangle, for example, depicts a group of three nodes, each connected to the other two. By leveraging simplicial complexes, researchers can analyze and characterize complex systems more nuancedly. Higher-order networks derived from simplicial complexes allow for exploring intricate relationships that might not be evident when considering only pairwise connections (Battiston et al., 2021). These higher-order structures provide a more detailed description of the system's organization and dynamics, leading to a deeper understanding of complex systems.
Synchronization is a phenomenon where elements in a system coordinate their behavior to achieve coherence (Boccaletti et al., 2002). This fundamental phenomenon can be witnessed in a variety of complex systems, including biological (Li et al., 2022), social (Sorrentino and Ott, 2007), and technological (Sivrikaya and Yener, 2004) networks. While much study has concentrated on the synchronization in conventional pairwise interactions (Lin et al., 2021; Fan et al., 2022, 2023), there is rising interest in understanding synchronization in higher-order networks that capture interactions beyond pairwise links (Battiston et al., 2020; Bick et al., 2023; Boccaletti et al., 2023). The goal of studying synchronization in higher-order networks is to decipher the intricate dynamics, emergent behaviors, and collective phenomena resulting from higher-order interactions. This emerging subject studies synchronization's emergence, development, and evolution in complex network structures such as simplicial complexes and hypergraphs. Researchers have been investigating synchronization in higher-order networks to unravel the fundamental principles driving collective dynamics and information processing in various fields spanning from neuroscience to social networks and beyond. For instance, Skardal et al. (2021) addressed how higher-order interactions effectively help to achieve optimized synchronization. Employing a random network of 500 phase oscillators optimized in order to get strongly synchronized, they found that strengthening the higher-order interactions resulted in improving the optimal synchronization. On the other hand, according to Gallo et al. (2022), directed higher-order interactions were reported as an impediment to attaining synchrony; otherwise, such interactions may stabilize unstable synchronized states. More clearly, studying a higher-order network of 8 Rössler oscillators structures in a random graph, the authors noticed that applying directed higher-order interactions could ruin the synchronization state. In contrast, as the asymmetry varied, synchronization was achievable. Multiplex higher-order network was the subject that Anwar and Ghosh (2022) pursued. They looked for synchronization in a two-layer network within each layer, the 500 Rössler oscillators interacting via diffusive pairwise and non-pairwise connections and structured in a scale-free configuration. Their results portrayed that higher-order interactions enhance intra-layer synchronization, and inter-layer synchronization becomes more robust compared to the pure pairwise case. Focusing on higher-order neuronal networks, Parastesh et al. (2022) analyzed the synchronization of a fully connected network of 20 Hindmarsh-Rose neurons interacting through two- and three-node interactions. Assuming different non-pairwise (electrical and chemical) interactions in combination with diffusive (electrical) pairwise connections, they showed that even weak second-order interactions could facilitate synchronization by reducing the pairwise strengths needed for achieving synchrony. In this regard, Mirzaei et al. (2022) pictured that scaling the synchronization patterns to lower coupling strength values was feasible via higher-dimensional interactions. To demonstrate the impact of adding three-, four-, and five-node chemical interactions stepwise on the synchronization state, they took into account a fully connected network of five conventional Rulkov maps with two-node diffusive inner linking and chemical connections. Mehrabbeik et al. (2023) recently investigated the synchronization of different higher-order networks made up of 10 Hindmarsh-Rose maps. The authors intended to figure out the impact of higher-order synaptic functions applied to two-node and three-node communication on network global synchronization. As a result, it was discovered that chemical synapses greatly improved synchronization compared to diffusively defined synapses, which include electrical and inner-linking functions.
Through the use of electrical and chemical synapses on two- and three-node interactions, the present study addresses synchronization in a higher-order network of mRulkov maps. The acquired results are compared to the case where electrical and chemical synapses are devoted simultaneously to two-node connections called hybrid synapses. The consequence of increasing the network size on synchronization state and synchronization manifold dynamics is additionally considered. Furthermore, the largest network (50 nodes) is sought for other synchronization patterns. The following is how the paper is set up: The mRulkov model and its dynamics are fully clarified in Section 2. In Sections 3 and 4, the subjected higher-order network model's linear stability analysis is covered in detail. Results are presented in Section 5, and Section 6 recaps the paper by emphasizing the findings.
2. The mRulkov neuron map
The mRulkov map is a mathematical neuron model that combines the dynamics of the 2D Rulkov map with the memristive nonlinearity proposed by Li et al. (2021). The original Rulkov map is a discrete-time dynamical system that exhibits complicated behavior, such as chaos, and is frequently used to resemble biological systems. The mRulkov model is a piecewise nonlinear model described by three dependent nonlinear difference equations as follows:
where x(n) and y(n) are the membrane potential and recovery state variables at time step n, and α, and β are parameters that determine the system's behavior. Also, the parameter ρ involves the effect of external factors on the model.
The addition of memristive nonlinearity to the Rulkov map results in a system that exhibits even more complex behavior, including hyperchaotic behavior and the emergence of extreme multistability (Li et al., 2021). The mathematical expression of the mRulkov map proposed in Li et al. (2021), to which the flux-controlled memristor is applied, is given as follows:
where φ is the flux variable, and μ is the magnetic induction strength induced by the membrane potential. For simplicity, ρ = 0 is considered in System (2). Considering α = 5, β = ε = 0.05, and μ = 0.55, the Rulkov model exhibits a chaotic dynamic as represented in Figure 1.
Figure 1
3. Discrete network model
Simplicial complexes, geometric objects that generalize triangles and tetrahedra to higher dimensions, are a powerful tool for representing higher-order interactions between groups of elements in a network. In a simplicial complex, each element is represented by a vertex, and higher-order interactions are represented by higher-dimensional simplices, such as edges (1-simplex), triangles (2-simplex), tetrahedra (3-simplex), pentahedra (4-simplex), and other polyhedra. As indicated by Figure 2, simplicial complexes allow for the representation of more complex relationships between elements than is possible with traditional network representations. Only links color-coded in gray have the key role in defining the structure of a network, as seen in Figure 2 (left panel). On the other hand, as shown in the right panel, a d-dimensional simplicial complex is formed when a set of nodes join together to form a morphological object with d faces. These linked neurons can also establish a novel network topology in which they interact in novel ways. For instance, a 2-simplex structure is formed when three nodes are joined together to create a triangle, and interactions involving three nodes must be taken into account.
Figure 2
The mathematical definition of a d-dimensional simplicial complex in the discrete domain is expressed as:
Here, , where X∈R3, contains the system's state variables belonging to the i-the system, F[Xi(n)] =[f(Xi(n)), g(Xi(n)), h(Xi(n))], where F:R3→R3, determines the nodes' dynamics, N is the number of nodes, and σ1, ..., σd are the coupling strength of interactions defined by 1 to d dimensional simplices. The network structure is defined by adjacency tensors A(1), ..., A(d), where (Wei and Ding, 2016; Gambuzza et al., 2021). Adjacency tensors are a generalization of adjacency matrices to higher-order networks. In an adjacency matrix, each element represents the presence or absence of a connection between a pair of nodes in a network, while in an adjacency tensor, each element represents the presence or absence of a connection between a set of nodes of arbitrary size (Lucas et al., 2020). For example, indicates that the i-th and j1-th nodes are connected through a link while shows that nodes with indices i, j1, ..., jd combine to form a polyhedron (d-simplex; Mirzaei et al., 2022; Parastesh et al., 2022; Mehrabbeik et al., 2023). Furthermore, G(1), ..., G(d) are the coupling functions describing the type of interactions among 2, ..., d+1 units building simplicial complexes.
A traditional network of N mRulkov maps coupled via hybrid (electrical and chemical) pairwise interactions and in a global configuration is mathematically depicted as follows:
where v is the reversal potential, is the chemical synaptic function that describes the kinetics of the neurotransmitter release and modulates the synaptic transmission using a sigmoidal function with a slope of r and a firing threshold of θ. The expression [v−xi(n)]Γ[xj(n)] represents the postsynaptic response at a synapse, where xi(n) and xj(n) are, respectively, the postsynaptic and presynaptic membrane potentials. The reversal potential v represents the neurotransmitter's equilibrium potential at the synapse. There is no net flow of ions when v = xi(n), and the synapse is said to be at its reversal potential. Moreover, Γ[xj(n)] denotes synaptic conductance and is influenced by several variables, including presynaptic activity. It influences how strongly neurotransmitters bind to receptors, affecting synaptic transmission and the postsynaptic membrane potential (Yamakou et al., 2020). If v−xi(n)>0, it results in depolarization (excitatory response), while if v−xi(n) < 0, it leads to hyperpolarization (inhibitory response; Shafiei et al., 2020). This response determines whether an action potential is generated a nd influences information transmission and processing within the neuronal network. Nonetheless, based on the general definition in Equation (3), considering a higher-order network of N globally coupled mRulkov neuron maps with first-order electrical connections and second-order chemical interactions (), Network (4) changes into the following equations:
In Networks (4) and (5), σ1 and σ2 denote the strength of the electrical and chemical synapses, respectively. Here, xj(n) and xk(n) are both considered as the postsynaptic neurons responsible for the release of the same neurotransmitters into the synaptic cleft of the postsynaptic neuron xi(n). For the subsequent investigations, the chemical synaptic parameters are set at v = θ = −1.4 and r = 50. The present study focuses on the collective dynamics of Network (5) to find the impact of higher-order interactions.
4. Stability analysis
The master stability function (MSF) is considered to assess the neurons' stability in the synchronization state (Pecora and Carroll, 1998). The MSF, which is influenced by the topology of the network, the degree of coupling between nodes, and the dynamics of each node, characterizes the linear stability of a synchronized state in a network. This function provides the necessary criteria for network synchronization by linearizing the dynamics of the network around the synchronized state and analyzing the resulting eigenvalues of the linearized system (Pecora and Carroll, 1998). To find the linearized system, a negligible perturbation δX(n) = [δx(n), δy(n), δφ(n)]T is locally added to the neurons in their synchronous state Xs(n) = [xs(n), ys(n), φs(n)]T. Thus, . Assuming V as a matrix with columns that are the eigenvectors of the Laplacian matrix generated from the adjacency matrix, the transformation leads to the desired linearized system.
In the first case, wherein no non-pairwise interactions are involved, and neurons are interacting through pairwise hybrid synapses (Network 4), the linearized system becomes:
where JF[Xs(n)] is the 3 × 3 Jacobian matrix of F[X(n)] evaluated at Xs(n) and is the derivative of Γ(x) in the synchronization state xs(n). When all neurons evolve synchronously, , and thus, and . Also, for a global network of N nodes, we have . Therefore, the dynamics of the synchronous neurons with hybrid synapses obey the following system:
which differs from the dynamics of an individual mRulkov neuron.
In the second case, wherein the non-pairwise chemical interactions, as well as the pairwise electrical connections, are involved (Network 5), the linearized system change into:
In the synchronization state, since , G(1)[Xs(n), Xs(n)]≡0 and G(2)[Xs(n), Xs(n), Xs(n)]≡[(v−xs(n))Γ2(xs(n)), 0, 0]. Moreover, due to the global coupling scheme, we have . Thus, the neurons' dynamics in the synchronization state can be described as follows:
which is different from System (7) or even an isolated neuron behavior.
The negative maximum Lyapunov exponent (Λ) of the linearized systems described in Systems (6) and (8) shows the stability of the synchronization state since the locally injected perturbations achieve zero, which means that the neurons remain in their synchronous state. The Appendix provides more information on how to conduct MSF analysis for both traditional and higher-order networks.
5. Results
Using the MSF analysis detailed in Section 4, here, the effect of applying higher-order interaction is investigated by examining Networks (4) and (5). The averaged synchronization error is also numerically determined in addition to the MSF technique, which analytically provides the essential synchronization criteria for the assessed networks to confirm the outcomes of the MSF method. The following definition applies to the averaged synchronization error:
where 〈...〉 and ||...|| are the functions calculating the averaged value over discrete time steps n and the Euclidean norm, respectively.
Figure 2 depicts the regions in the parameter plane σ1-σ2 (0 ≤ σ1 ≤ 0.24 and 0 ≤ σ2 ≤ 0.03) wherein N = 5 mRulkov neurons interacting through the hybrid synapses [Network (4)] are completely synchronous. Such region can be detected for Λ < 0 (Figure 3A) and E = 0 (Figure 3B). Similarly, Figure 4 provides the results of Lyapunov analysis of System (8) and synchronization error calculation of Network (5) in which the chemical synapse are considered as the second-order interactions alongside the first-order electrical connections for 0 ≤ σ1 ≤ 0.24 and 0 ≤ σ2 ≤ 0.01. Figure 4 reveals that when second-order chemical interactions are involved with the pairwise electrical connections, weaker strength of chemical couplings is needed to have the same synchronization patterns as shown in Figure 3. More clearly, higher-order interactions enhance synchronization by scaling the patterns to lower values of σ2.
Figure 3
Figure 4
The dynamics of synchronous neurons do not follow the equation for a single isolated neuron, as demonstrated by System (9); nonetheless, they are dependent on the network settings as well as the model parameters, specifically the number of neurons N and the chemical (second-order) coupling parameter σ2. Figure 5 illustrates how the dynamics of neurons in the synchronization state, or in other words, the dynamics of System (9), vary according to the coupling parameter σ2 by performing a simple dynamical analysis using the bifurcation and the Lyapunov exponents (LEs) diagrams. It can be seen that if first- and second-order coupling parameter values are selected in the synchronous regions demonstrated in Figures 4A, C, the mRulkov neurons are able to exhibit chaotic and periodic behaviors according to the dynamical analysis performed in Figure 5. More precisely, it can be recognized that the neurons synchronize with the dynamics shown in Figure 5 if the first-order coupling strength σ1 is chosen in the synchronous zone indicated in Figure 4. For instance, Figures 6A, B shows that N = 5 mRulkov neurons in a higher-order network defined in Network (5) achieve synchrony with periodic dynamics (LE1 = −0.2472, LE2 = −0.0656, and LE3 = 0) for σ1 = 0.1 and σ2 = 0.002; however, if σ1 = 0.1 and σ2 = 0.01 are selected, as indicated by Figures 6C, D, they behave chaotically (LE1 = −0.2065, LE2 = 0, and LE3 = 0.0499) in the synchronization state.
Figure 5
Figure 6
To further examine the impact of network size on the synchronization state of the higher-order network composed of mRulkov models, larger networks with more interacting neurons are taken into consideration. Figures 7A, B demonstrates the synchronous and asynchronous regions of Network (5) for 0 ≤ σ1 ≤ 0.06 and 0 ≤ σ2 ≤ 0.00035 when N = 20 neurons are involved. Similarly, the stability regions of Network (5) with N = 50 neurons are 0 ≤ σ1 ≤ 0.024 and 0 ≤ σ2 ≤ 0.0000505 shown in Figures 7C, D. More precisely, the results show that as the network size increases, the synchronization patterns are scaled to the lower values of both first-order and second-order coupling parameters. However, the amount of this decrease is not the same for σ1 and σ2. Figure 8 illustrates how the number of participating neurons affects the network's synchronization according to the variation of σ1 while σ2 = 0.0001 (Figure 8A), and σ2 while σ1 = 0.001 (Figure 8B).
Figure 7
Figure 8
Focusing on the asynchronous regions in Figures 7C, D, wherein N = 50 mRulkov neurons are configured in a higher-order network described by Network (5), cluster synchronization patterns can be detected. More precisely, it is found that the neurons evolve asynchronously in lower values of σ1 (in asynchronous regions) while in higher values, they tend to participate in forming synchronous clusters. For instance, as shown in Figures 9A, B, two-cluster synchronization is found for σ1 = 0.024 and σ2 = 0.00004, in which the synchronous neurons behave chaotically. Based on Figures 9C, D, the same two-cluster synchronization pattern is also identified for σ1 = 0.023 and σ2 = 0.00004; however, the neurons evolve periodically synchronously in each cluster.
Figure 9
6. Conclusions
The consequences of applying higher-order interactions on the synchronization of the mRulkov network were well-explained in this research. In order to achieve this goal, a network of globally connected mRulkov neurons was considered, in which electrical and chemical synapses were respectively applied to the two- and three-node interactions. Thereafter, the regions of the coupling parameter space wherein the neurons were completely synchronous were detected through the MSF approach, which was then verified by calculating the network's synchronization error. The same study was carried out on a network with pure pairwise hybrid interactions, wherein electrical and chemical pathways were considered active simultaneously. Our findings suggest that not neglecting the non-pairwise or multi-node interactions can improve global synchronization by scaling synchronization patterns to lower chemical coupling parameter values while leaving the electrical coupling strength the same. Furthermore, the higher-order coupling strength and network size were demonstrated to alter the behavior of neurons in the synchronization state. Therefore, through the bifurcation analysis, the behaviors of the synchronous neurons, regardless of the stability of the synchronization manifold, were investigated concerning the variation of the second-order chemical coupling parameter. Additionally, it was demonstrated that the neurons could exhibit periodic or chaotic behaviors in the synchronous zone of the parameter space. Through the Lyapunov analysis of the linearized system developed with the MSF formalism, the influence of the engaged neurons' number on the network synchronization was also investigated. The results showed that the synchronous patterns scale to smaller values of the coupling parameters as the network size grows. Furthermore, looking more closely at the asynchronous regions, cluster synchronization patterns were detected. It was shown that the synchronous neurons in the cluster have periodic or chaotic dynamics.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
MM, SJ, and MP designed and performed the research as well as wrote the paper. All authors contributed to the article and approved the submitted version.
Funding
MP was supported by the Slovenian Research and Innovation Agency (Javna agencija za znanstvenoraziskovalno in inovacijsko dejavnost Republike Slovenije) (Grant Nos. P1-0403, J1-2457, and N1-0232).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2023.1248976/full#supplementary-material
References
1
AiharaK.TakabeT.ToyodaM. (1990). Chaotic neural networks. Phys. Lett. A144, 333–340. 10.1016/0375-9601(90)90136-C
2
Alvarez-RodriguezU.BattistonF.de ArrudaG. F.MorenoY.PercM.LatoraV. (2021). Evolutionary dynamics of higher-order interactions in social networks. Nat. Hum. Behav. 5, 586–595. 10.1038/s41562-020-01024-1
3
AnwarM. S.GhoshD. (2022). Intralayer and interlayer synchronization in multiplex network with higher-order interactions. Chaos32:033125. 10.1063/5.0074641
4
BaoB.HuA.BaoH.XuQ.ChenM.WuH. (2018). Three-dimensional memristive Hindmarsh-Rose neuron model with hidden coexisting asymmetric behaviors. Complexity2018:3872573. 10.1155/2018/3872573
5
BaoH.HuaZ.LiuW.BaoB. (2021). Discrete memristive neuron model and its interspike interval-encoded application in image encryption. Sci. China Technol. Sci. 64, 2281–2291. 10.1007/s11431-021-1845-x
6
BaoH.YuX.XuQ.WuH.BaoB. (2022). Three-dimensional memristive Morris-Lecar model with magnetic induction effects and its FPGA implementation. Cogn. Neurodyn. 17, 1079–1092. 10.1007/s11571-022-09871-6
7
BattistonF.AmicoE.BarratA.BianconiG.Ferraz de ArrudaG.FranceschielloB.et al. (2021). The physics of higher-order interactions in complex systems. Nat. Phys. 17, 1093–1098. 10.1038/s41567-021-01371-4
8
BattistonF.CencettiG.IacopiniI.LatoraV.LucasM.PataniaA.et al. (2020). Networks beyond pairwise interactions: structure and dynamics. Phys. Rep. 874, 1–92. 10.1016/j.physrep.2020.05.004
9
BickC.GrossE.HarringtonH. A.SchaubM. T. (2023). What are higher-order networks?SIAM Rev. 65, 686–731. 10.1137/21M1414024
10
BoccalettiS.De LellisP.del GenioC. I.Alfaro-BittnerK.CriadoR.JalanS.et al. (2023). The structure and dynamics of networks with higher order interactions. Phys. Rep. 1018, 1–64. 10.1016/j.physrep.2023.04.002
11
BoccalettiS.KurthsJ.OsipovG.ValladaresD. L.ZhouC. S. (2002). The synchronization of chaotic systems. Phys. Rep. 366, 1–101. 10.1016/S0370-1573(02)00137-0
12
BoccalettiS.LatoraV.MorenoY.ChavezM.HwangD. U. (2006). Complex networks: structure and dynamics. Phys. Rep. 424, 175–308. 10.1016/j.physrep.2005.10.009
13
BurgioG.MatamalasJ. T.GómezS.ArenasA. (2020). Evolution of cooperation in the presence of higher-order interactions: from networks to hypergraphs. Entropy22:744. 10.3390/e22070744
14
CarlettiT.FanelliD.NicolettiS. (2020). Dynamical systems on hypergraphs. J. Phys. 1:035006. 10.1088/2632-072X/aba8e1
15
ChangeuxJ.-P.DehaeneS. (1989). Neuronal models of cognitive functions. Cognition33, 63–109. 10.1016/0010-0277(89)90006-1
16
ChenM.QiJ.WuH.XuQ.BaoB. (2020). Bifurcation analyses and hardware experiments for bursting dynamics in non-autonomous memristive FitzHugh-Nagumo circuit. Sci. China Technol. Sci. 63, 1035–1044. 10.1007/s11431-019-1458-5
17
ChialvoD. R. (1995). Generic excitable dynamics on a two-dimensional map. Chaos Solitons Fractals5, 461–479. 10.1016/0960-0779(93)E0056-H
18
ChoW.LeeD.KimB. J. (2023). A multiresolution framework for the analysis of community structure in international trade networks. Sci. Rep. 13:5721. 10.1038/s41598-023-32686-2
19
DingS.WangN.BaoH.ChenB.WuH.XuQ. (2023). Memristor synapse-coupled piecewise-linear simplified Hopfield neural network: dynamics analysis and circuit implementation. Chaos Solitons Fractals166:112899. 10.1016/j.chaos.2022.112899
20
FanW.ChenX.WuH.LiZ.XuQ. (2023). Firing patterns and synchronization of morris-lecar neuron model with memristive autapse. Int. J. Electron. Commun. 158:154454. 10.1016/j.aeue.2022.154454
21
FanW.WuH.LiZ.XuQ. (2022). Synchronization and chimera in a multiplex network of Hindmarsh-Rose neuron map with flux-controlled memristor. Eur. Phys. J. Spec. Top. 231, 4131–4141. 10.1140/epjs/s11734-022-00720-5
22
FengX.JoW. S.KimB. J. (2014). International transmission of shocks and fragility of a bank network. Phys. A403, 120–129. 10.1016/j.physa.2014.02.030
23
FitzHughR. (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophys. J. 1, 445–466. 10.1016/S0006-3495(61)86902-6
24
GalloL.MuoloR.GambuzzaL. V.LatoraV.FrascaM.CarlettiT. (2022). Synchronization induced by directed higher-order interactions. Commun. Phys. 5:263. 10.1038/s42005-022-01040-9
25
GambuzzaL. V.Di PattiF.GalloL.LepriS.RomanceM.CriadoR.et al. (2021). Stability of synchronization in simplicial complexes. Nat. Commun. 12:1255. 10.1038/s41467-021-21486-9
26
GhorbanchianR.RestrepoJ. G.TorresJ. J.BianconiG. (2021). Higher-order simplicial synchronization of coupled topological signals. Commun. Phys. 4:120. 10.1038/s42005-021-00605-4
27
GhoshD.FrascaM.RizzoA.MajhiS.RakshitS.Alfaro-BittnerK.et al. (2022). The synchronized dynamics of time-varying networks. Phys. Rep. 949, 1–63. 10.1016/j.physrep.2021.10.006
28
HindmarshJ.RoseR. (1984). A model of neuronal bursting using three coupled first order differential equations. Proc. R. Soc. B: Biol. Sci. 221, 87–102. 10.1098/rspb.1984.0024
29
HodgkinA. L.HuxleyA. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117:500. 10.1113/jphysiol.1952.sp004764
30
HuX.LiuC. (2019). Dynamic property analysis and circuit implementation of simplified memristive Hodgkin-Huxley neuron model. Nonlin. Dyn. 97, 1721–1733. 10.1007/s11071-019-05100-8
31
IbarzB.CasadoJ. M.SanjuánM. A. (2011). Map-based models in neuronal dynamics. Phys. Rep. 501, 1–74. 10.1016/j.physrep.2010.12.003
32
InceR. A. A.MontaniF.ArabzadehE.DiamondM. E.PanzeriS. (2009). On the presence of high-order interactions among somatosensory neurons and their effect on information transmission. J. Phys. Conf. Ser. 197:012013. 10.1088/1742-6596/197/1/012013
33
IzhikevichE. M.HoppensteadtF. (2004). Classification of bursting mappings. Int. J. Bifurcat. Chaos14, 3847–3854. 10.1142/S0218127404011739
34
KimB. J.JunT.KimJ.-Y.ChoiM. Y. (2006). Network marketing on a small-world network. Phys. A360, 493–504. 10.1016/j.physa.2005.06.059
35
LiK.BaoB.MaJ.ChenM.BaoH. (2022). Synchronization transitions in a discrete memristor-coupled bi-neuron model. Chaos Solitons Fractals165:112861. 10.1016/j.chaos.2022.112861
36
LiK.BaoH.LiH.MaJ.HuaZ.BaoB. (2021). Memristive Rulkov neuron model with magnetic induction effects. IEEE Trans. Ind. Inf. 18, 1726–1736. 10.1109/TII.2021.3086819
37
LinH.WangC.ChenC.SunY.ZhouC.XuC.HongQ. (2021). Neural bursting and synchronization emulated by neural networks and circuits. IEEE Trans. Circuits Syst. I Regul. Pap. 68, 3397–3410. 10.1109/TCSI.2021.3081150
38
LinH.WangC.HongQ.SunY. (2020). A multi-stable memristor and its application in a neural network. IEEE Trans. Circuits Syst. II Express Briefs67, 3472–3476. 10.1109/TCSII.2020.3000492
39
LotitoQ. F.MusciottoF.MontresorA.BattistonF. (2022). Higher-order motif analysis in hypergraphs. Commun. Phys. 5:79. 10.1038/s42005-022-00858-7
40
LucasM.CencettiG.BattistonF. (2020). Multiorder Laplacian for synchronization in higher-order networks. Phys. Rev. Res. 2:033410. 10.1103/PhysRevResearch.2.033410
41
MaJ.TangJ. (2017). A review for dynamics in neuron and neuronal network. Nonlinear Dyn. 89, 1569–1578. 10.1007/s11071-017-3565-3
42
MajhiS.BeraB. K.GhoshD.PercM. (2019). Chimera states in neuronal networks: a review. Phys. Life Rev. 28, 100–121. 10.1016/j.plrev.2018.09.003
43
MajhiS.PercM.GhoshD. (2022). Dynamics on higher-order networks: a review. J. R. Soc. Interface19:20220043. 10.1098/rsif.2022.0043
44
MehrabbeikM.AhmadiA.BakouieF.JafariA. H.JafariS.GhoshD. (2023). The impact of higher-order interactions on the synchronization of Hindmarsh-Rose neuron maps under different coupling functions. Mathematics11:2811. 10.3390/math11132811
45
MirzaeiS.MehrabbeikM.RajagopalK.JafariS.ChenG. (2022). Synchronization of a higher-order network of Rulkov maps. Chaos32:123133. 10.1063/5.0117473
46
MorrisC.LecarH. (1981). Voltage oscillations in the barnacle giant muscle fiber. Biophys. J. 35, 193–213. 10.1016/S0006-3495(81)84782-0
47
ParasteshF.MehrabbeikM.RajagopalK.JafariS.PercM. (2022). Synchronization in Hindmarsh-Rose neurons subject to higher-order interactions. Chaos32:013125. 10.1063/5.0079834
48
PecoraL. M.CarrollT. L. (1998). Master stability functions for synchronized coupled systems. Phys. Rev. Lett. 80, 2109–2112. 10.1103/PhysRevLett.80.2109
49
RulkovN. F. (2002). Modeling of spiking-bursting neural behavior using two-dimensional map. Phys. Rev. E65:041922. 10.1103/PhysRevE.65.041922
50
ShafieiM.JafariS.ParasteshF.OzerM.KapitaniakT.PercM. (2020). Time delayed chemical synapses and synchronization in multilayer neuronal networks with ephaptic inter-layer coupling. Commun. Nonlinear Sci. Numer. Simul. 84:105175. 10.1016/j.cnsns.2020.105175
51
ShahalS.WurzbergA.SibonyI.DuadiH.ShnidermanE.WeymouthD.et al. (2020). Synchronization of complex human networks. Nat. Commun. 11:3854. 10.1038/s41467-020-17540-7
52
SivrikayaF.YenerB. (2004). Time synchronization in sensor networks: a survey. IEEE Netw. 18, 45–50. 10.1109/MNET.2004.1316761
53
SkardalP. S.ArenasA. (2020). Higher order interactions in complex networks of phase oscillators promote abrupt synchronization switching. Commun. Phys. 3, 1–6. 10.1038/s42005-020-00485-0
54
SkardalP. S.Arola-FernándezL.TaylorD.ArenasA. (2021). Higher-order interactions can better optimize network synchronization. Phys. Rev. Res. 3:043193. 10.1103/PhysRevResearch.3.043193
55
SorrentinoF.OttE. (2007). Network synchronization of groups. Phys. Rev. E76:056114. 10.1103/PhysRevE.76.056114
56
VivekanandhanG.NatiqH.MerrikhiY.RajagopalK.JafariS. (2023). Dynamical analysis and synchronization of a new memristive Chialvo neuron model. Electronics12:545. 10.3390/electronics12030545
57
WeiY.DingW. (2016). Theory and Computation of Tensors: Multi-Dimensional Arrays. London: Academic Press.
58
XuQ.LiuT.DingS.BaoH.LiZ.ChenB. (2023). Extreme multistability and phase synchronization in a heterogeneous bi-neuron Rulkov network with memristive electromagnetic induction. Cogn. Neurodyn. 17, 755–766. 10.1007/s11571-022-09866-3
59
YamakouM. E.HjorthP. G.MartensE. A. (2020). Optimal self-induced stochastic resonance in multiplex neural networks: Electrical vs. chemical synapses. Front. Comput. Neurosci. 14:62. 10.3389/fncom.2020.00062
60
Zandi-MehranN.PanahiS.HosseiniZ.Hashemi GolpayeganiS. M. R.JafariS. (2020). One dimensional map-based neuron model: a phase space interpretation. Chaos Solitons Fractals132:109558. 10.1016/j.chaos.2019.109558
Summary
Keywords
simplicial complex, higher-order network, memristive Rulkov, synchronization, cluster synchronization
Citation
Mehrabbeik M, Jafari S and Perc M (2023) Synchronization in simplicial complexes of memristive Rulkov neurons. Front. Comput. Neurosci. 17:1248976. doi: 10.3389/fncom.2023.1248976
Received
27 June 2023
Accepted
11 August 2023
Published
31 August 2023
Volume
17 - 2023
Edited by
Thanos Manos, Université de Cergy-Pontoise, France
Reviewed by
Sorinel A. Oprisan, College of Charleston, United States; Kun Li, Hebei University of Technology, China
Updates
Copyright
© 2023 Mehrabbeik, Jafari and Perc.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Matjaž Perc matjaz.perc@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.