Next Article in Journal
Novel Approaches for Differentiable Convex Functions via the Proportional Caputo-Hybrid Operators
Previous Article in Journal
An Averaging Principle for Stochastic Fractional Differential Equations Driven by fBm Involving Impulses
Previous Article in Special Issue
Global Exponential Stability of Fractional Order Complex-Valued Neural Networks with Leakage Delay and Mixed Time Varying Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model

1
Department of Mathematics and Computer Science, West University of Timişoara, 300223 Timişoara, Romania
2
Institute for Advanced Environmental Research, West University of Timişoara, 300223 Timişoara, Romania
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2022, 6(5), 257; https://doi.org/10.3390/fractalfract6050257
Submission received: 13 April 2022 / Revised: 4 May 2022 / Accepted: 5 May 2022 / Published: 7 May 2022
(This article belongs to the Special Issue Frontiers in Fractional-Order Neural Networks)

Abstract

:
The aim of this work is to describe the dynamics of a fractional-order coupled FitzHugh–Nagumo neuronal model. The equilibrium states are analyzed in terms of their stability properties, both dependently and independently of the fractional orders of the Caputo derivatives, based on recently established theoretical results. Numerical simulations are shown to clarify and exemplify the theoretical results.

1. Introduction

During the last few years, an increased number of scientific papers have examined the relevance of using fractional-order derivatives when modeling real-world phenomena. Thus, it has been suggested that fractional-order systems are efficiently providing more reasonable and realistic results in a substantial number of practical applications [1,2,3,4,5] in comparison with their associated integer-order counterpart. This is justified by the fact that fractional-order derivatives are provided with both memory and hereditary properties, as [6] suggests that a possible physical meaning for the order of a fractional derivative could be the index of memory.
It has been recently underlined [7,8] that fractional-order derivatives or differences could be involved in mathematically modeling neuronal dynamics. Due to their ability to introduce capacitive memory effects [9], fractional-order membrane potential dynamics have emphasized their advantage in reproducing the electrical activity of neurons observed from an experimental point of view. Very recently, ref. [10] proposed a novel mathematical model of neuronal electromechanics employing fractional-order derivatives of variable order to model multiple temporal scales, accounting for both local and nonlocal chemo-mechanical interactions observed experimentally [11]. Several types of fractional-order single-neuronal models have been investigated in recent years: leaky integrate-and-fire [12], Hindmarsh–Rose [13,14], Morris–Lecar [15,16,17], FitzHugh–Nagumo [18] and more general Hodgkin–Huxley models [9,19].
Coupled FitzHugh–Nagumo-type integer-order systems have recently been investigated. On a general note, a review on chimera states in neuronal networks has been presented in [20]. In particular, several new phenomena of spiral wave chimeras have been discovered in [21], where a two-dimensional nonlocally coupled FitzHugh–Nagumo system with an open boundary condition was considered. Chimera patterns in several types of two-dimensional networks of coupled neurons have been explored in [22], whereas a system of two identical FitzHugh–Nagumo units with a mutual linear coupling in the fast variables was investigated in [23]. Moreover, complex activities in neuronal systems can be better understood due to the results presented in [24], where a coupled network consisting of an arbitrary number of nonidentical FitzHugh–Nagumo neurons is investigated in terms of its stability properties and chaotic behavior. Coupled FitzHugh–Nagumo equations have also been considered in [25], whereas the firing activities of a fractional-order FitzHugh–Rinzel bursting neuron model and its coupled dynamics have been investigated in [26].
Regarding coupled FitzHugh–Nagumo neuronal models, ref. [27] presents novel patterns in fractional-in-space nonlinear coupled FitzHugh–Nagumo models with a Riesz fractional derivative; the chimera state in the network of fractional-order FitzHugh–Nagumo neurons is investigated in [28], whereas an analytical study of a fractional-order multiple chaotic FitzHugh–Nagumo neuron model using the multistep generalized differential transform method is presented in [29].
The main reference works on the topic of fractional calculus and the qualitative theory of fractional-order systems are [30,31,32,33]. A significant component in the qualitative theory of fractional-order systems is represented by stability analysis. Recently, main results were developed with respect to the stability properties of fractional-order systems [34,35]. Hence, ref. [36] presents a generalization of the well-known stability theorem of Matignon [37]. Furthermore, linearization theorems for fractional-order systems are presented [38].
The aim of this work is to investigate the stability of the equilibrium states of a fractional-order coupled FitzHugh–Nagumo neuronal model, by applying theoretical results recently obtained in [39,40]. Both fractional-order independent results and fractional-order dependent results are compactly enumerated, for their later applicability in the analysis of the neuronal model mentioned above. To the best of our knowledge, this is the first time that a complete theoretical stability analysis has been done for a coupled FitzHugh–Nagumo neuron model with different fractional orders.
The paper is organized as follows. Section 2 is designed for itemizing some preliminaries and the fractional-order dependent and independent results. The main results of the work are described in Section 3: the mathematical model is presented, the existence of equilibrium states is discussed, and the stability of the equilibrium states is later analyzed, both independent and dependent on the fractional orders of the considered neuronal model. The results are supported by illustrative numerical simulations, revealing rich spiking behavior. Several conclusions are drawn in Section 5.

2. Preliminaries

Let us consider the n-dimensional fractional-order system with Caputo derivatives [30,32,41]:
c D q x ( t ) = f ( t , x )
where q = ( q 1 , q 2 , . . . , q n ) ( 0 , 1 ) n and f : [ 0 , ) × R n R n is a continuous function on the whole domain of definition, Lipschitz-continuous with respect to the second variable, such that
f ( t , 0 ) = 0 for any t 0 .
The existence and uniqueness of solutions of the initial value problem associated with system (1) is guaranteed by the previously mentioned properties of the function f [31].
For the local stability analysis of a general nonlinear system (1), we rely on the linearization technique in a neighborhood of an equilibrium, which has been previously explored in [37]. It is important to emphasize that once the linearized system is established, stability and instability results may be obtained by analyzing the distribution of the roots of the corresponding characteristic equation [40].
For the stability analysis undertaken in the following sections, we recall some results recently obtained in [39,42], concerning the distribution of the roots of the following characteristic equation:
s q 1 + q 2 + β 1 s q 1 + β 2 s q 2 + γ = 0 ,
where s q 1 and s q 2 represent the principal values (first branches) of the corresponding complex power functions [43]. For completeness, we reformulate these previous results in the following propositions:
Proposition 1
(Fractional-order independent results).
  • All the roots of the characteristic Equation (2) are in the open left half-plane, regardless of the fractional orders q 1 and q 2 , if and only if the following inequalities are satisfied:
    γ > 0 , β 1 + β 2 > 0 and min { 1 , γ } + min { β 1 , β 2 } > 0 .
  • The characteristic Equation (2) has a root in the open right half-plane, regardless of the fractional orders q 1 and q 2 , if and only if either one of the following conditions is satisfied:
    i.
    γ < 0 ;
    ii.
    γ > 0 and β 1 + β 2 + γ + 1 0 ;
    iii.
    β 1 < 0 , β 2 < 0 and β 1 β 2 γ > 0 .
In what follows, we assume that the fractional orders q 1 and q 2 are arbitrarily fixed inside the domain
D = { ( q 1 , q 2 ) R 2 : 0 < q 1 < q 2 1 } .
Moreover, as γ < 0 implies that the system is unstable, for any choice of the fractional orders q 1 and q 2 , we will assume that γ > 0 .
As in [39,42], we define a family of smooth parametric curves in the ( β 1 , β 2 ) -plane by
Γ γ , q 1 , q 2 : β 1 = γ q 2 q 1 + q 2 h ω , q 1 , q 2 β 2 = γ q 1 q 1 + q 2 h ω , q 2 , q 1 , ω > 0 ,
where h : ( 0 , ) × D R is given by:
h ω , q 1 , q 2 = ω q 1 q 1 + q 2 ω ρ ( q 1 , q 2 q 1 ) ρ ( q 2 , q 2 q 1 )
with the function ρ defined as
ρ ( a , b ) = sin a π 2 sin b π 2 , a [ 0 , 1 ] , b [ 1 , 0 ) ( 0 , 1 ] .
It has been previously shown in [42] that the curve Γ ( γ , q 1 , q 2 ) is the graph of a smooth, decreasing, convex bijective function ϕ γ , q 1 , q 2 : R R in the ( β 1 , β 2 ) -plane. Moreover, the curve Γ ( γ , q 1 , q 2 ) lies outside the first quadrant of the ( β 1 , β 2 ) -plane.
With these notations, we have the following result [39]:
Proposition 2
(Fractional-order dependent results).
Let γ > 0 and 0 < q 1 < q 2 1 be arbitrarily fixed. Consider the curve Γ ( γ , q 1 , q 2 ) and the function ϕ γ , q 1 , q 2 : R R defined above.
i.
The characteristic Equation (2) has a pair of complex conjugated roots on the imaginary axis of the complex plane if and only if ( β 1 , β 2 ) Γ ( γ ; q 1 , q 2 ) .
ii.
All the roots of the characteristic Equation (2) are in the open left half-plane if and only if
β 2 > ϕ γ , q 1 , q 2 ( β 1 ) .
iii.
If β 2 < ϕ γ , q 1 , q 2 ( β 1 ) , the characteristic Equation (2) has at least one root in the open right half-plane.

3. Fractional-Order Coupled FithHugh–Nagumo-Type Neuronal Model

Consider the following fractional-order coupled FitzHugh–Nagumo neuronal model, which is a modified version of the model considered in [25], replacing the classical integer-order derivative with fractional-order Caputo derivatives:
c D q 1 v 1 ( t ) = v 1 ( v 1 a ) ( 1 v 1 ) w 1 + g ( v 1 v 2 ) c D q 2 w 1 ( t ) = ε ( v 1 β w 1 ) c D q 1 v 2 ( t ) = v 2 ( v 2 a ) ( 1 v 2 ) w 2 + g ( v 2 v 1 ) c D q 2 w 2 ( t ) = ε ( v 2 β w 2 )
where v 1 and v 2 represent the membrane potential of the two neurons, w 1 and w 2 are recovery variables, a , β , ε , g are positive constants, and q 1 and q 2 are the fractional orders of the Caputo derivatives, with 0 < q 1 q 2 1 .
In this paper, identical neurons are considered; therefore, the fractional orders of the Caputo derivatives in system (3) coincide for the membrane voltage equations and for the recovery variable equations, respectively. In many previous works concerning single-neuron models [9], q 2 = 1 has been chosen for the recovery variable, while q 1 ( 0 , 1 ] accounts for the incorporation of a non-ideal fractional-order capacitive element in the membrane potential equation [44]. However, it has been pointed out in [19] that an increase in the diversity of spike patterns and shapes is observed in a single-neuron model, when using a power-law behaving conductance, which emphasizes the importance of considering q 2 ( 0 , 1 ] .
We observe that
c D q 2 w 1 ( t ) = ε ( v 1 β w 1 ) = ε β 1 β v 1 w 1 = ϕ ( α v 1 w 1 ) c D q 2 w 2 ( t ) = ε ( v 2 β w 2 ) = ε β 1 β v 2 w 2 = ϕ ( α v 2 w 2 )
where ϕ : = ε β > 0 and α : = 1 β > 0 .
Considering the function I ( v i , w i ) = w i v i ( v i a ) ( 1 v i ) , with i { 1 , 2 } , system (3) is equivalently written as:
c D q 1 v 1 ( t ) = g ( v 1 v 2 ) I ( v 1 , w 1 ) c D q 1 v 2 ( t ) = g ( v 2 v 1 ) I ( v 2 , w 2 ) c D q 2 w 1 ( t ) = ϕ ( α v 1 w 1 ) c D q 2 w 2 ( t ) = ϕ ( α v 2 w 2 )

3.1. Existence of Equilibrium States

In this section, we investigate the existence of equilibrium points of the fractional-order neuronal model (4), which are the solutions ( v 1 * , v 2 * , w 1 * , w 2 * ) of the algebraic system:
g ( v 1 v 2 ) = I ( v 1 , w 1 ) g ( v 2 v 1 ) = I ( v 2 , w 2 ) α v 1 = w 1 α v 2 = w 2
which reduces to solving the following two-dimensional algebraic system:
g ( v 1 v 2 ) = I ( v 1 ) g ( v 2 v 1 ) = I ( v 2 )
where I ( v i ) = I ( v i , α v i ) = v i 3 ( 1 + a ) v i 2 + ( α + a ) v i , with i { 1 , 2 } .
Adding the two equations of system (5), we obtain:
I ( v 1 ) + I ( v 2 ) = 0 .
We distinguish the following two cases:
  • Case 1: v 1 * = v 2 *
In this case, we have I ( v 1 * ) = 0 , which is equivalent to
v 1 * ( ( v 1 * ) 2 ( 1 + a ) v 1 * + α + a ) = 0 .
It follows that v 1 * = 0 of ( v 1 * ) 2 ( a + 1 ) v 1 * + α + a = 0 .
The discriminant of the previous quadratic equations is Δ = ( 1 a ) 2 4 α . Therefore, we have the following situations:
  • If 4 α > ( 1 a ) 2 , then system (4) has a unique equilibrium point
    ( v 1 * , v 2 * , w 1 * , w 2 * ) = ( 0 , 0 , 0 , 0 ) .
  • If 4 α = ( 1 a ) 2 , then system (4) has two equilibrium states
    ( v 1 * , v 2 * , w 1 * , w 2 * ) ( 0 , 0 , 0 , 0 ) , 1 + a 2 , 1 + a 2 , α ( 1 + a ) 2 , α ( 1 + a ) 2 .
  • If 4 α < ( 1 a ) 2 , then system (4) has three equilibrium states
    ( v 1 * , v 2 * , w 1 * , w 2 * ) ( 0 , 0 , 0 , 0 ) , f ( a , Δ ) , f ( a , Δ ) , α f ( a , Δ ) , α f ( a , Δ ) ,
    where f ( a , Δ ) = 1 + a ± Δ 2 .
Case 2: v 1 * v 2 *
The equilibrium states ( v 1 * , v 2 * , w 1 * , w 2 * ) = ( v 1 * , v 2 * , α v 1 * , α v 2 * ) of system (4), with v 1 * v 2 * , are the solutions of the algebraic system
g ( v 1 v 2 ) = I ( v 1 ) g ( v 2 v 1 ) = I ( v 2 )
which is equivalent to the system
g ( v 1 v 2 ) = v 1 3 ( 1 + a ) v 1 2 + ( α + a ) v 1 g ( v 2 v 1 ) = v 2 3 ( 1 + a ) v 2 2 + ( α + a ) v 2 .
Subtracting the equations of the previous system, we obtain
2 g ( v 1 v 2 ) = v 1 3 v 2 3 ( 1 + a ) ( v 1 2 v 2 2 ) + ( α + a ) ( v 1 v 2 )
which, as we consider asymmetric equilibria, is equivalent to
2 g = v 1 2 + v 1 v 2 + v 2 2 ( 1 + a ) ( v 1 + v 2 ) + α + a .
Moreover, as I ( v 1 ) + I ( v 2 ) = 0 , it follows that
v 1 3 + v 2 3 ( 1 + a ) ( v 1 2 + v 2 2 ) + ( α + a ) ( v 1 + v 2 ) = 0 .
Denoting s = v 1 + v 2 and p = v 1 · v 2 , Equation (6) is equivalent to
p = s 2 ( 1 + a ) s + α + a 2 g
whereas Equation (7) becomes
s 3 3 s p ( 1 + a ) ( s 2 2 p ) + ( α + a ) s = 0 .
Combining the previous two relations, it results that
s 3 2 ( 1 + a ) s 2 + ( α + 3 a 3 g + a 2 + 1 ) s ( 1 + a ) ( α + a 2 g ) = 0 .
As we search for real distinct values of v 1 and v 2 , which are the roots of the quadratic equation v 2 s v + p = 0 , the inequality s 2 > 4 p must also be satisfied, or equivalently, from relation (8):
3 4 s 2 ( 1 + a ) s + α + a 2 g < 0 .
Therefore, determining the asymmetrical equilibrium states reduces to solving the cubic Equation (10), subject to the constraint (11). Hence, depending on the values of the parameters a, α and g, we may obtain zero to three pairs of asymmetrical equilibrium states of the form ( v 1 * , v 2 * , α v 1 * , α v 2 * ) and ( v 2 * , v 1 * , α v 2 * , α v 1 * ) .
Remark 1.
In Figure 1, one can deduce the number of equilibrium states of system (4) depending on the parameters of the system a, α and for fixed values of the parameter g. For values of the parameters a, α belonging to the green region, the system has three equilibrium states. If the parameters are situated in the lighter green region, the system has a total of five equilibrium states. In the orange region, the system has a total of seven equilibrium states, and in the red region, there is a total of nine equilibrium states of system (4).

3.2. Stability of Equilibrium States

In order to establish the stability of the equilibrium states previously determined, consider the Jacobian matrix associated with system (4) at an arbitrary equilibrium point ( v 1 * , v 2 * , w 1 * , w 2 * ) :
J = α + g I ( v 1 * ) g 1 0 g α + g I ( v 2 * ) 0 1 ϕ α 0 ϕ 0 0 ϕ α 0 ϕ
which can be written in the following block matrix form:
J = M I 2 ϕ α I 2 ϕ I 2 ,
where I 2 denotes a two-dimensional unit matrix and
M = M ( v 1 * , v 2 * ) = α + g I ( v 1 * ) g g α + g I ( v 2 * )
is a symmetric matrix with real eigenvalues denoted by:
μ ± = μ ± ( v 1 * , v 2 * ) = g + α I ( v 1 * ) + I ( v 2 * ) 2 ± g 2 + I ( v 1 * ) I ( v 2 * ) 2 2 .
The characteristic equation associated with the equilibrium ( v 1 * , v 2 * , w 1 * , w 2 * ) is
det ( J diag ( s q 1 , s q 1 , s q 2 , s q 2 ) ) = 0 ,
or equivalently:
det ( ϕ + s q 2 ) ( M s q 1 I 2 ) + ϕ α I 2 = 0
which can be expressed as:
det M s q 1 + ϕ α ϕ + s q 2 I 2 = 0 .
Hence, we obtain that s is a root of the characteristic Equation (13) if and only if s q 1 + ϕ α ϕ + s q 2 is an eigenvalue of the matrix M. It follows that the characteristic Equation (13) can be written as:
s q 1 + q 2 + ϕ s q 1 μ s q 2 + ϕ ( α μ ) · s q 1 + q 2 + ϕ s q 1 μ + s q 2 + ϕ ( α μ + ) = 0 .
Applying the theoretical results from Section 3, we obtain the following fractional-order-independent characterizations of the asymptotic stability and instability of the equilibrium states:
Proposition 3.
Denoting
μ s = min 1 , ϕ , ϕ α 1 + ϕ and μ u = min α , ϕ α 1 + ϕ + 1 ,
the inequality μ s < μ u holds. Moreover, the equilibrium state ( v 1 * , v 2 * , w 1 * , w 2 * ) of system (4) is:
  • asymptotically stable, regardless of the fractional orders q 1 and q 2 , if and only if μ + < μ s ;
  • unstable, regardless of the fractional orders q 1 and q 2 , if and only if μ + > μ u .
Proof. 
The proof of the inequality μ s < μ u is trivial.
On one hand, the equilibrium state ( v 1 * , v 2 * , w 1 * , w 2 * ) is asymptotically stable, regardless of the fractional order q 1 and q 2 , if and only if all the roots of the characteristic Equation (14) are in the open left half-plane, or equivalently, if and only if all the roots of both equations
s q 1 + q 2 + ϕ s q 1 μ s q 2 + ϕ ( α μ ) = 0 , with μ { μ , μ + }
have negative real parts. Considering β 1 = ϕ > 0 , β 2 = μ and γ = ϕ ( α μ ) , based on the first statement of Proposition 1, this holds if and only if
μ < min { α , ϕ } and min { 1 , ϕ ( α μ ) } + min { ϕ , μ } > 0 ,
which reduces to
μ < min 1 , ϕ , ϕ α 1 + ϕ = μ s .
Taking into account that μ < μ + , it is enough to require μ + < μ s .
On the other hand, the equilibrium state ( v 1 * , v 2 * , w 1 * , w 2 * ) is unstable, regardless of the fractional order q 1 and q 2 , if and only if the characteristic Equation (14) has at least one root in the open right half-plane, i.e., if either one of the equations from (15) has a root with a positive real part. From the second statement of Proposition 1, this holds if either γ < 0 , or γ > 0 and β 1 + β 2 + γ + 1 0 , which translates to
μ > min α , ϕ α 1 + ϕ + 1 = μ u .
As μ + > μ , it is enough to require μ + > μ u , and the desired conclusion is obtained. □
Remark 2.
Proposition 3 gives necessary and sufficient conditions for the fractional-order-independent stability and instability of the equilibrium states of system (4), formulated in terms of simple inequalities involving the eigenvalue μ + , which incorporates information about the equilibrium state, according to (15). It is important to remark that, based on the above result, when μ + ( μ s , μ u ) , the stability properties of the equilibrium state depend on the fractional orders q 1 and q 2 , as will be described in what follows.
Proposition 4.
Assume that the fractional orders ( q 1 , q 2 ) are arbitrarily fixed such that 0 < q 1 q 2 1 and μ + ( μ s , μ u ) . There exists a unique value μ ( q 1 , q 2 ) ( μ s , μ u ) such that the equilibrium state ( v 1 * , v 2 * , w 1 * , w 2 * ) of system (4) is asymptotically stable if μ + < μ * ( q 1 , q 2 ) and unstable if μ + > μ * ( q 1 , q 2 ) .
Proof. 
The critical value μ * ( q 1 , q 2 ) is the value for which Equation (15) has a pair of roots s = ± i ω on the imaginary axis. Hence, following Proposition 2, μ * ( q 1 , q 2 ) is given by the unique solution of the system (with the unknowns μ ( μ s , μ u ) and ω > 0 )
ϕ = [ ϕ ( α μ ) ] q 2 q 1 + q 2 h ω , q 1 , q 2 μ = [ ϕ ( α μ ) ] q 1 q 1 + q 2 h ω , q 2 , q 1
As the qualitative behavior in a neighborhood of the equilibrium changes only at the critical value μ * ( q 1 , q 2 ) , the rest of the proof is a consequence of Proposition 3. □
Remark 3.
Considering the trivial equilibrium ( 0 , 0 , 0 , 0 ) of system (4) (or system (3)), a simple computation shows that μ + = g + | g | a .
In the case of excitable coupling ( g < 0 ), it follows that μ + = a < μ s , and hence, the trivial equilibrium is always asymptotically stable, for any choice of the fractional orders, or the coupling coefficient g.
However, in the case of phase-repulsive coupling ( g > 0 ), we have μ + = 2 g a , and hence, based on Propositions 3 and 4, we have:
  • if g < μ s + a 2 , the trivial equilibrium state is asymptotically stable, regardless of the fractional orders q 1 and q 2 ;
  • if g > μ u + a 2 , the trivial equilibrium state is unstable, regardless of the fractional orders q 1 and q 2 ;
  • if g μ s + a 2 , μ u + a 2 , the stability of the trivial equilibrium depends on the fractional orders q 1 and q 2 , as seen from Proposition 4.

4. Numerical Simulations

For the numerical simulations presented in this paper, we have employed the Adams–Bashforth–Moulton method developed for fractional-order systems in [45]. For a review of existing numerical methods that can be successfully applied to fractional-order systems with multiple fractional orders, we refer to [46].

4.1. Case 1: A Unique Equilibrium State

For the first set of numerical simulations, we consider the parameter values used in [25]: a = 0.3 , ϵ = 0.01 and β = 0.1 . Hence, we have α = 10 and ϕ = 0.001 . It can be easily seen that for this combination of system parameters (also see Figure 1), the only equilibrium state of system (3) is the trivial one. Moreover, the parameters given by Proposition 3 are: μ s = 0.001 and μ u = 1.00999 . From Remark 3, we deduce that if g < 0.1505 , the trivial equilibrium is asymptotically stable, independently of q 1 and q 2 , while for g > 0.654995 , the trivial equilibrium is unstable, for any fractional orders q 1 and q 2 . However, when g ( 0.1505 , 0.654995 ) , the stability of the trivial equilibrium depends on the choice of fractional orders q 1 and q 2 .
Based on the above observations, we first consider g = 0.2 , and hence, μ + = 0.1 . Fixing the fractional order q 2 = 1 , we determine the critical value of the fractional order q 1 by numerically solving system (16) for q 1 and ω , considering μ = μ + in this system. We find q 1 * = 0.633408 , and consequently, the trivial equilibrium is asymptotically stable from q 1 < q 1 * and unstable for q 1 > q 1 * , which is in agreement with the numerical simulations presented in Figure 2. For q 1 > q 1 * , periodic oscillations are observed in system (3), and the frequency of oscillations increases for larger values of q 1 . It can be conjectured that a supercritical Hopf-type bifurcation takes place in a neighborhood of the trivial equilibrium when q 1 = q 1 * ; however, at this point, the bifurcation theory of fractional-order systems is not well studied.
In contrast with the classical integer-order case, it may be observed from Figure 2 that as the fractional order q 1 decreases from 1 to q 1 * , single high-amplitude spikes alternate with small-amplitude oscillations, which do not fully decay to a quiescent phase. This remarkable difference between the integer-order and fractional-order models provides insight into the complex dynamics exhibited by systems of coupled neurons, and justifies a thorough further study of the different oscillatory patterns induced by the fractional orders.

4.2. Case 2: Five Coexisting Equilibrium States

For the second set of numerical simulations, we consider the following parameter values: a = 1.5 , ϵ = 0.032 , β = 2 and g = 0.8 . Hence, we have α = 0.5 and ϕ = 0.064 . For this combination of system parameters (also see Figure 1), five equilibrium states coexist in system (3), given in Table 1.
The parameters given by Proposition 3 are: μ s = 0.0300752 and μ u = 0.5 . Therefore, as μ + < μ s for the first pair of asymmetric equilibrium states, we deduce that they are asymptotically stable, regardless of the fractional orders chosen in system (3). On the other hand, as μ + > μ u for the second pair of asymmetric equilibrium states, it follows that they are unstable, regardless of the fractional orders chosen in system (3). However, as μ + ( μ s , μ u ) for the trivial equilibrium state, from Proposition 4, the stability of the trivial equilibrium depends on the choice of fractional orders.
For the trivial equilibrium state, fixing the fractional order q 2 = 1 , we determine the critical value of the fractional order q 1 by numerically solving system (16) for q 1 and ω , considering μ = μ + = 0.1 in this system. We find q 1 * = 0.911087 , and consequently, the trivial equilibrium is asymptotically stable from q 1 < q 1 * and unstable for q 1 > q 1 * . This is illustrated in Figure 3. Indeed, for q 1 = 0.8 , all trajectories initiated from a small neighborhood of the trivial equilibrium converge to origin. On the hand, for q 1 = 0.95 , the trajectories initiated from a small neighborhood of the trivial equilibrium converge to one of the asymptotically stable equilibria of the first asymmetric pair from Table 1.

5. Conclusions

Stability analysis of the equilibrium states of a fractional-order coupled FitzHugh–Nagumo neuronal model was explored. Firstly, some preliminary results were enumerated. Then, the mathematical model was described using fractional-order derivatives of Caputo type. For the mathematical model, equilibrium states were determined, which were later investigated in terms of their stability. Moreover, numerical methods were provided in order to illustrate the theoretical results.
Directions for future research are the generalization of the present investigation to a network on non-identical neurons, and to higher-dimensional networks of FitzHugh–Nagumo neurons.

Author Contributions

Conceptualization, O.B. and E.K.; methodology, O.B. and E.K.; software, O.B. and E.K.; validation, O.B. and E.K.; formal analysis, O.B. and E.K.; investigation, O.B. and E.K.; resources, O.B. and E.K.; data curation, O.B. and E.K.; writing—original draft preparation, O.B. and E.K.; writing—review and editing, O.B. and E.K.; visualization, O.B. and E.K.; supervision, O.B. and E.K.; project administration, O.B. and E.K.; funding acquisition, O.B. and E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CNCS-UEFISCDI, Project No. PN-III-P4-PCE2021-0204.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cottone, G.; Paola, M.D.; Santoro, R. A novel exact representation of stationary colored Gaussian processes (fractional differential approach). J. Phys. A Math. Theor. 2010, 43, 085002. [Google Scholar] [CrossRef]
  2. Engheia, N. On the role of fractional calculus in electromagnetic theory. IEEE Antennas Propag. Mag. 1997, 39, 35–46. [Google Scholar] [CrossRef]
  3. Henry, B.; Wearne, S. Existence of Turing instabilities in a two-species fractional reaction-diffusion system. SIAM J. Appl. Math. 2002, 62, 870–887. [Google Scholar] [CrossRef] [Green Version]
  4. Heymans, N.; Bauwens, J.C. Fractal rheological models and fractional differential equations for viscoelastic behavior. Rheologica Acta 1994, 33, 210–219. [Google Scholar] [CrossRef]
  5. Mainardi, F. Fractional Relaxation-Oscillation and Fractional Phenomena. Chaos Solitons Fractals 1996, 7, 1461–1477. [Google Scholar] [CrossRef]
  6. Du, M.; Wang, Z.; Hu, H. Measuring memory with the order of fractional derivative. Sci. Rep. 2013, 3, 3431. [Google Scholar] [CrossRef]
  7. Anastasio, T. The fractional-order dynamics of brainstem vestibulo-oculomotor neurons. Biol. Cybernet. 1994, 72, 69–79. [Google Scholar] [CrossRef]
  8. Lundstrom, B.; Higgs, M.; Spain, W.; Fairhall, A. Fractional differentiation by neocortical pyramidal neurons. Nat. Neurosci. 2008, 11, 1335–1342. [Google Scholar] [CrossRef]
  9. Weinberg, S.H. Membrane capacitive memory alters spiking in neurons described by the fractional-order Hodgkin-Huxley model. PLoS ONE 2015, 10, e0126629. [Google Scholar] [CrossRef]
  10. Drapaca, C. Fractional calculus in neuronal electromechanics. J. Mech. Mater. Struct. 2016, 12, 35–55. [Google Scholar] [CrossRef]
  11. Grevesse, T.; Dabiri, B.E.; Parker, K.K.; Gabriele, S. Opposite rheological properties of neuronal microcompartments predict axonal vulnerability in brain injury. Sci. Rep. 2015, 5, 9475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Teka, W.; Marinov, T.M.; Santamaria, F. Neuronal spike timing adaptation described with a fractional leaky integrate-and-fire model. PLoS Comput. Biol. 2014, 10, e1003526. [Google Scholar] [CrossRef] [PubMed]
  13. Jun, D.; Guang-Jun, Z.; Yong, X.; Hong, Y.; Jue, W. Dynamic behavior analysis of fractional-order Hindmarsh–Rose neuronal model. Cogn. Neurodyn. 2014, 8, 167–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kaslik, E. Analysis of two-and three-dimensional fractional-order Hindmarsh-Rose type neuronal models. Fract. Calc. Appl. Anal. 2017, 20, 623–645. [Google Scholar] [CrossRef] [Green Version]
  15. Brandibur, O.; Kaslik, E. Stability properties of a two-dimensional system involving one Caputo derivative and applications to the investigation of a fractional-order Morris-Lecar neuronal model. Nonlinear Dyn. 2017, 90, 2371–2386. [Google Scholar] [CrossRef] [Green Version]
  16. Shi, M.; Wang, Z. Abundant bursting patterns of a fractional-order Morris–Lecar neuron model. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 1956–1969. [Google Scholar] [CrossRef]
  17. Upadhyay, R.K.; Mondal, A.; Teka, W.W. Fractional-order excitable neural system with bidirectional coupling. Nonlinear Dyn. 2017, 87, 2219–2233. [Google Scholar] [CrossRef]
  18. Brandibur, O.; Kaslik, E. Stability of two-component incommensurate fractional-order systems and applications to the investigation of a FitzHugh-Nagumo neuronal model. Math. Methods Appl. Sci. 2018, 41, 7182–7194. [Google Scholar] [CrossRef] [Green Version]
  19. Teka, W.; Stockton, D.; Santamaria, F. Power-Law Dynamics of Membrane Conductances Increase Spiking Diversity in a Hodgkin-Huxley Model. PLoS Comput. Biol. 2016, 12, e1004776. [Google Scholar] [CrossRef]
  20. Majhi, S.; Bera, B.K.; Ghosh, D.; Perc, M. Chimera states in neuronal networks: A review. Phys. Life Rev. 2019, 28, 100–121. [Google Scholar] [CrossRef]
  21. Guo, S.; Dai, Q.; Cheng, H.; Li, H.; Xie, F.; Yang, J. Spiral wave chimera in two-dimensional nonlocally coupled Fitzhugh–Nagumo systems. Chaos Solitons Fractals 2018, 114, 394–399. [Google Scholar] [CrossRef]
  22. Schmidt, A.; Kasimatis, T.; Hizanidis, J.; Provata, A.; Hövel, P. Chimera patterns in two-dimensional networks of coupled neurons. Phys. Rev. E 2017, 95, 032224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Eydam, S.; Franović, I.; Wolfrum, M. Leap-frog patterns in systems of two coupled FitzHugh-Nagumo units. Phys. Rev. E 2019, 99, 042207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mao, X. Complicated dynamics of a ring of nonidentical FitzHugh–Nagumo neurons with delayed couplings. Nonlinear Dyn. 2017, 87, 2395–2406. [Google Scholar] [CrossRef]
  25. Lavrova, S.; Kudryashov, N.; Sinelshchikov, D. On some properties of the coupled Fitzhugh-Nagumo equations. J. Phys. Conf. Ser. 2019, 1205, 012035. [Google Scholar] [CrossRef] [Green Version]
  26. Mondal, A.; Sharma, S.K.; Upadhyay, R.K.; Mondal, A. Firing activities of a fractional-order FitzHugh-Rinzel bursting neuron model and its coupled dynamics. Sci. Rep. 2019, 9, 15721. [Google Scholar] [CrossRef]
  27. Li, X.; Han, C.; Wang, Y. Novel Patterns in Fractional-in-Space Nonlinear Coupled FitzHugh–Nagumo Models with Riesz Fractional Derivative. Fractal Fract. 2022, 6, 136. [Google Scholar] [CrossRef]
  28. Ramadoss, J.; Aghababaei, S.; Parastesh, F.; Rajagopal, K.; Jafari, S.; Hussain, I. Chimera state in the network of fractional-order fitzhugh–nagumo neurons. Complexity 2021, 2021, 2437737. [Google Scholar] [CrossRef]
  29. Momani, S.; Freihat, A.; Al-Smadi, M. Analytical study of fractional-order multiple chaotic FitzHugh-Nagumo neurons model using multistep generalized differential transform method. Abstract Appl. Anal. 2014, 2014, 276279. [Google Scholar] [CrossRef] [Green Version]
  30. Podlubny, I. Fractional Differential Equations; Academic Press: Cambridge, MA, USA, 1999. [Google Scholar]
  31. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  32. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  33. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and breach Science Publishers: Yverdon Yverdon-les-Bains, Switzerland, 1993; Volume 1. [Google Scholar]
  34. Li, C.; Zhang, F. A survey on the stability of fractional differential equations. Eur. Phys. J. Special Top. 2011, 193, 27–47. [Google Scholar] [CrossRef]
  35. Rivero, M.; Rogosin, S.V.; Tenreiro Machado, J.A.; Trujillo, J.J. Stability of fractional order systems. Math. Prob. Eng. 2013, 2013, 356215. [Google Scholar] [CrossRef]
  36. Sabatier, J.; Farges, C. On stability of commensurate fractional order systems. Int. J. Bifurc. Chaos 2012, 22, 1250084. [Google Scholar] [CrossRef]
  37. Matignon, D. Stability Results For Fractional Differential Equations With Applications To Control Processing. In Proceedings of the Computational Engineering in Systems Applications, Lille, France, 9–12 July 1996; pp. 963–968. [Google Scholar]
  38. Cong, N.; Tuan, H.; Trinh, H. On asymptotic properties of solutions to fractional differential equations. J. Math. Anal. Appl. 2020, 484, 123759. [Google Scholar] [CrossRef] [Green Version]
  39. Brandibur, O.; Kaslik, E. Exact stability and instability regions for two-dimensional linear autonomous systems of fractional-order differential equations. Fract. Calc. Appl. Anal. 2021, 24, 225–253. [Google Scholar] [CrossRef]
  40. Brandibur, O.; Garrappa, R.; Kaslik, E. Stability of Systems of Fractional-Order Differential Equations with Caputo Derivatives. Mathematics 2021, 9, 914. [Google Scholar] [CrossRef]
  41. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Scientific Publishers: Cambridge, MA, USA, 2009. [Google Scholar]
  42. Brandibur, O.; Kaslik, E. Stability analysis of multi-term fractional-differential equations with three fractional derivatives. J. Math. Anal. Appl. 2021, 495, 124751. [Google Scholar] [CrossRef]
  43. Doetsch, G. Introduction to the Theory and Application of the Laplace Transformation; Springer: Berlin/Heidelberg, Germany, 1974. [Google Scholar]
  44. Westerlund, S.; Ekstam, L. Capacitor theory. IEEE Trans. Dielectr. Electr. Insul. 1994, 1, 826–839. [Google Scholar] [CrossRef]
  45. Diethelm, K.; Ford, N.; Freed, A. A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dyn. 2002, 29, 3–22. [Google Scholar] [CrossRef]
  46. Garrappa, R. Numerical solution of fractional differential equations: A survey and a software tutorial. Mathematics 2018, 6, 16. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Number of equilibrium states in the parametric plane ( a , α ) for different values of the parameter g { ± 0.2 , ± 0.6 , ± 1 } .
Figure 1. Number of equilibrium states in the parametric plane ( a , α ) for different values of the parameter g { ± 0.2 , ± 0.6 , ± 1 } .
Fractalfract 06 00257 g001
Figure 2. Evolution of the membrane potentials v 1 (blue) and v 2 (red) of the two neurons, for fixed parameter values a = 0.3 , ϵ = 0.01 and β = 0.1 , coupling coefficient g = 0.2 , fixed fractional order q 2 = 1 and q 1 { 0.6 , 0.7 , 0.8 , 0.9 , 1 } (top to bottom). Initial conditions for system (3) have been chosen in a neighborhood of the trivial equilibrium.
Figure 2. Evolution of the membrane potentials v 1 (blue) and v 2 (red) of the two neurons, for fixed parameter values a = 0.3 , ϵ = 0.01 and β = 0.1 , coupling coefficient g = 0.2 , fixed fractional order q 2 = 1 and q 1 { 0.6 , 0.7 , 0.8 , 0.9 , 1 } (top to bottom). Initial conditions for system (3) have been chosen in a neighborhood of the trivial equilibrium.
Fractalfract 06 00257 g002
Figure 3. Evolution of the membrane potential v 1 for fixed parameter values a = 1.5 , ϵ = 0.032 and β = 2 , coupling coefficient g = 0.8 , fixed fractional order q 2 = 1 and q 1 { 0.8 , 0.95 } , considering multiple initial conditions for system (3) in a neighborhood of the trivial equilibrium (corresponding solutions plotted with different colors).
Figure 3. Evolution of the membrane potential v 1 for fixed parameter values a = 1.5 , ϵ = 0.032 and β = 2 , coupling coefficient g = 0.8 , fixed fractional order q 2 = 1 and q 1 { 0.8 , 0.95 } , considering multiple initial conditions for system (3) in a neighborhood of the trivial equilibrium (corresponding solutions plotted with different colors).
Fractalfract 06 00257 g003
Table 1. Equilibrium states of (3) for the parameter values: a = 1.5 , ϵ = 0.032 and β = 2 .
Table 1. Equilibrium states of (3) for the parameter values: a = 1.5 , ϵ = 0.032 and β = 2 .
Equilibrium ( v 1 * , v 2 * , w 1 * , w 2 * ) μ + Stability
( 2.01375 , 0.555812 , 1.00687 , 0.277906 ) 2.46674 asympt. stable
( 0.555812 , 2.01375 , 0.277906 , 1.00687 ) 2.46674 asympt. stable
( 0.38957 , 0.183994 , 0.194785 , 0.0919969 ) 1.02553 unstable
( 0.183994 , 0.38957 , 0.0919969 , 0.194785 ) 1.02553 unstable
( 0 , 0 , 0 , 0 ) 0.1 depends on q 1 , q 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brandibur, O.; Kaslik, E. Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model. Fractal Fract. 2022, 6, 257. https://doi.org/10.3390/fractalfract6050257

AMA Style

Brandibur O, Kaslik E. Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model. Fractal and Fractional. 2022; 6(5):257. https://doi.org/10.3390/fractalfract6050257

Chicago/Turabian Style

Brandibur, Oana, and Eva Kaslik. 2022. "Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model" Fractal and Fractional 6, no. 5: 257. https://doi.org/10.3390/fractalfract6050257

APA Style

Brandibur, O., & Kaslik, E. (2022). Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model. Fractal and Fractional, 6(5), 257. https://doi.org/10.3390/fractalfract6050257

Article Metrics

Back to TopTop