Extreme events at the onset of epileptic-like chimeras in small-world networks of FitzHugh-Nagumo neurons

Javier Cubillos Cornejo,Miguel Escobar Mendoza,Ignacio Bordeu

2024-11-06

Abstract:In this work, we investigate the dynamics of complex networks of FitzHugh-Nagumo excitable oscillators, focusing on the impact of coupling strength, network size, and randomness on their collective dynamics. Considering Watts-Strogatz small-world network connectivities, the system exhibits three distinct dynamical phases: chaotic, intermittent, and synchronized, with the intermittent phase displaying transient, epileptic-like chimera states. We analyse the transition to synchronisation by means of the master stability function, and show that peaks in the proportion of extreme events of synchronisation, which correlate with the behaviour of the largest Lyapunov exponent of the system, precede the transitions between the distinct dynamical regimes and mark the onset of epileptic-like chimera states. Our findings contribute to a broader understanding of synchronisation in excitable systems real neural networks and offer insights into the conditions that may lead to pathological epileptic-like states. Furthermore, we discus the potential use of extreme events to study real neural data.

Adaptation and Self-Organizing Systems,Chaotic Dynamics,Biological Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the dynamic behavior of FitzHugh - Nagumo (FHN) oscillators in complex networks, especially the influence of coupling strength, network size and randomness on collective dynamics. Specifically, the researchers focus on three different dynamic phases exhibited by the system under Watts - Strogatz small - world network connectivity: the chaotic phase, the intermittent phase and the synchronous phase. Among them, the intermittent phase will show short - lived epileptic - like chimera states. The researchers analyzed the transition of synchronization through the master stability function (MSF) and showed that the peak proportion of extreme synchronization events is related to the behavior of the maximum Lyapunov index of the system. These peaks predict the transition between different dynamic phases and mark the emergence of epileptic - like chimera states.

### Main problem summary:
1. **Transition of dynamic phases**: The researchers explored how coupling strength, network size and connection randomness affect the dynamic phase transition of the FHN oscillator system.
2. **Extreme events and chimera states**: The researchers analyzed the frequency of extreme synchronization events and their relationship with the change of system dimension to predict the transition from chaos to intermittency and then to synchronization.
3. **Understanding of epileptic - like states**: The researchers paid special attention to the chimera states in the intermittent phase, which is similar to the instantaneous synchronization phenomenon of the brain during an epileptic seizure, aiming to provide a theoretical basis for understanding pathological epileptic - like states.

### Formulas and methods:
- **FitzHugh - Nagumo model**:
\[
\begin{pmatrix}
\dot{x}_i \\
\dot{y}_i
\end{pmatrix} =
\begin{pmatrix}
x_i-\frac{x_i^3}{3}-y_i \\
x_i + a
\end{pmatrix}+d\sum_{j = 1}^N A_{ij}R\begin{pmatrix}
x_j - x_i \\
y_j - y_i
\end{pmatrix}
\]
where \(\dot{x}_i\) and \(\dot{y}_i\) represent the time derivatives of the activation variable and the inhibition variable respectively, \(\epsilon= 0.05\), the parameter \(a\) controls the transition from the excitable state to the oscillatory state, \(d\) is the coupling strength, \(A_{ij}\) is the adjacency matrix, and \(R\) is the coupling matrix.

- **Kuramoto global order parameter**:
\[
R(t)=\left|\frac{1}{N}\sum_{j = 1}^N e^{i\theta_j(t)}\right|
\]
where \(\theta_j(t)=\arctan\left(\frac{y_j(t)}{x_j(t)}\right)\), which is used to measure the degree of synchronization of the system.

- **Master stability function**:
\[
\Lambda_{\text{max}}(\nu)
\]
where \(\nu = \gamma_2d\), \(\gamma_2\) is the second largest eigenvalue of the network Laplacian matrix, which is used to analyze the stability of the synchronous solution.

### Conclusion:
This study reveals the important influence of network randomness and coupling strength on the emergence of epileptic - like states, and proposes a method to predict dynamic phase transitions by using extreme synchronization events. These findings are helpful for a deeper understanding of the synchronization phenomenon in neural networks and its pathological mechanisms, and provide a new perspective for simulating and understanding pathological states such as epilepsy.