Hopf bifurcations in dynamics of excitable systems

Monica De Angelis

DOI: https://doi.org/10.1007/s11587-022-00742-0

2023-06-13

Abstract:A general FitzHugh-Rinzel model, able to describe several neuronal phenomena, is considered. Linear stability and Hopf bifurcations are investigated by means of the spectral equation for the ternary autonomous dynamical system and the analysis is driven by both an admissible critical point and a parameter which characterizes the system.

Chaotic Dynamics,Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to analyze the linear stability and Hopf bifurcation phenomena in the dynamical system of the FitzHugh - Rinzel (FHR) model. Specifically, the paper focuses on:

1. **Linear stability analysis**: By studying the linear stability of the critical points of the FHR system, determine the stability of these critical points under different parameter conditions.

2. **Hopf bifurcation analysis**: Explore the conditions under which the system changes from a stable state to a periodic oscillation state under specific parameter changes, that is, the occurrence conditions of Hopf bifurcation.

3. **Influence of bifurcation parameters**: Analyze the influence of different bifurcation parameters (such as the critical point \(\bar{U}\) and the coefficient \(-\eta = -\varepsilon\beta\)) on the dynamical behavior of the system, especially how they lead to steady - state bifurcations and oscillatory bifurcations.

### Main research contents

- **Model introduction**: The paper first introduces the FHR model, which is an extended FitzHugh - Nagumo model used to describe the bursting oscillation phenomenon of neurons. The model equations are as follows:
\[
\begin{cases}
\frac{dU}{dt}=-aU + U^{2}(a + 1)-\frac{1}{k}U^{3}-W + Y+I\\
\frac{dW}{dt}=\varepsilon(-\beta W + c+U)\\
\frac{dY}{dt}=\delta(-U + h - dY)
\end{cases}
\]
where \(U\) represents the transmembrane potential, \(W\) is the recovery variable, and \(Y\) is the slow current in the dendrite.

- **Linear stability analysis**: By linearizing the system and performing stability analysis near the critical points, use the Lienard - Chipart criterion to judge the stability of the system.

- **Hopf bifurcation conditions**: By introducing the "instability coefficient power" (ICP) method, determine the conditions for the system to undergo Hopf bifurcation. Specifically, analyze the bifurcation situations driven by the critical point \(\bar{U}\) and the coefficient \(-\eta\).

### Key conclusions

- **Hopf bifurcation driven by the critical point \(\bar{U}\)**:
- When \(\bar{U}\) satisfies certain conditions, the system is linearly globally attractive and asymptotically exponentially stable.
- When \(\bar{U}\) does not satisfy these conditions, the system may undergo a simple oscillatory bifurcation or a steady - state plus oscillatory bifurcation.

- **Hopf bifurcation driven by the coefficient \(-\eta\)**:
- By choosing \(\eta=\varepsilon\beta\) as the bifurcation parameter, the paper proves that under certain conditions, the system can undergo Hopf bifurcation even when \(\Gamma\geq0\).

### Discussion

The research results of this paper are not only of great significance for understanding the bursting oscillation phenomenon of neurons, but may also be applied to other fields, such as superconducting Josephson junctions and viscoelastic problems. By analyzing the dynamical behavior of the FHR model, theoretical support and new insights can be provided for research in these fields.