Abstract:Abstract Fractional-order dynamics of excitable systems can be physically described as a memory dependent phenomenon. It can produce diverse and fascinating oscillatory patterns for certain types of neuron models. To address these characteristics, we consider a nonlinear fast-slow FitzHugh-Rinzel (FH-R) model that exhibits elliptic bursting at a fixed set of parameters with a constant input current. The generalization of this classical order model provides a wide range of neuronal responses (regular spiking, fast-spiking, bursting, mixed-mode oscillations, etc.) in understanding the single neuron dynamics. So far, it is not completely understood to what extent the fractional-order dynamics may redesign the firing properties of excitable systems. We investigate how the classical order system changes its complex dynamics and how the bursting changes to different oscillations with stability and bifurcation analysis depending on the fractional exponent (0 < α ≤ 1). This occurs due to the memory trace of the fractional-order dynamics. The firing frequency of the fractional-order FH-R model is less than the classical order model, although the first spike latency exists there. Further, we investigate the responses of coupled FH-R neurons with small coupling strengths that synchronize at specific fractional-orders. The interesting dynamical characteristics suggest various neurocomputational features that can be induced in this fractional-order system which enriches the functional neuronal mechanisms.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to explore how fractional - order dynamics reshape the firing characteristics of neuronal systems, and analyze the complex dynamic changes of the fractional - order FitzHugh - Rinzel (FH - R) neuron model and its coupling dynamics. Specifically, the researchers focus on the following aspects:
1. **The influence of fractional - order dynamics on the dynamics of a single neuron**: By introducing fractional - order derivatives to describe the dynamic behavior of neurons, the researchers hope to more accurately simulate the memory - dependent phenomenon of neurons. They explored how fractional - order derivatives change the firing patterns of neurons, for example, from bursting to different oscillation patterns.
2. **Stability and bifurcation analysis**: Through stability analysis and bifurcation analysis, the researchers explored the dynamic changes of the system under different fractional - order exponents ($0 < \alpha \leq 1$). In particular, they studied how the fractional - order exponent affects the stability of the system and the location of bifurcation points.
3. **Synchronous activities of coupled neurons**: The researchers also explored how coupled FH - R neurons achieve synchronization under small coupling strengths, and analyzed the influence of different fractional - order exponents on synchronous activities.
4. **Biocomputational features of fractional - order dynamics**: The researchers pointed out that fractional - order dynamics can generate multiple neural computational features, enriching the understanding of functional neural mechanisms.
### Specific research contents
- **Definition of the fractional - order FH - R model**:
\[
\begin{cases}
\frac{d^\alpha v}{dt^\alpha}=v - \frac{v^3}{3}-w + I + y\\
\frac{d^\alpha w}{dt^\alpha}=a(v - bw - f(v, w, y))\\
\frac{d^\alpha y}{dt^\alpha}=c(y - df(v, w, y))
\end{cases}
\]
where $v$ represents the membrane voltage, $w$ represents the recovery variable, $y$ represents the slow modulation of current, $I$ represents the external stimulation current, $\alpha$ is the fractional - order exponent, and $a, b, c, d$ are system parameters.
- **Numerical solution method**:
The researchers used the fractional - order derivative defined by Caputo and adopted the L1 scheme for numerical solution. The specific formula is as follows:
\[
\frac{d^\alpha x(t)}{dt^\alpha}\approx\frac{1}{\Gamma(2 - \alpha)}\sum_{k = 0}^{N - 1}\left[\frac{x(t_{k + 1})-x(t_k)}{(t_{N - k - 1})^{1+\alpha}}-\frac{x(t_{k + 1})-x(t_k)}{(t_{N - k - 2})^{1+\alpha}}\right]
\]
- **Stability analysis**:
The researchers explored the stability conditions of the system by analyzing the Jacobian matrix of the system. For example, when the system is at the Hopf bifurcation point, the stability of the system can be judged by the following condition:
\[
\alpha < \frac{\pi}{2}\cdot\frac{\lambda_1+\lambda_2}{|\lambda_1\lambda_2|}
\]
- **Synchronous activities of coupled neurons**:
The researchers explored the synchronous behaviors of coupled FH - R neurons under different fractional - order exponents and found that fractional - order dynamics can significantly affect the nature of synchronous activities.
### Conclusion
Through the above research, the researchers revealed the important role of fractional - order dynamics in neuronal systems, especially its contributions in reshaping firing patterns, stability analysis, and synchronous activities. These findings not only enrich the understanding of neuronal dynamics but also provide a new perspective for future neural computing models.