NANOSYSTEMS: PHYSICS, CHEMISTRY, MATHEMATICS, 2021, 12 (5), P. 553-562
Properties of an oriented ring of neurons with the FitzHugh-Nagumo model
E. G. Fedorov
ITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia [email protected]
PACS 87.19.La, 87.18.Sn
DOI 10.17586/2220-8054-2021-12-5-553-562
The transmission of an impulse through a neuron is provided by processes that occur at the nanoscale level. This paper will build a model for an oriented ring of connected neurons. To describe the process of impulse transmission through a neuron, the FitzHugh-Nagumo model is used, which allows one to set a higher abstraction level by simulating an impulse. In this case, when transmitting impulses between neurons, the delay is taken into account. For the constructed model, the dependence of the number of neurons on the dynamics of the network as a whole is studied, and local bifurcations are found. All results are verified numerically. It is shown that the period of self-oscillations of such a network depends on the number of neurons.
Keywords: Neural networks, Hopf bifurcations, FitzHugh-Nagumo system, communication delay, stability. Received: 15 July 2021 Revised: 19 September 2021
1. Introduction
The time delay can greatly affect the network of neurons and significantly expand the range of possible behaviors of the network. Because of this, models of networks are actively being investigated, where transmission between neurons occurs with a delay. There are works here on a network of connected neurons with an identical delay [1,2]. And works on different transmission lag times [3-6]. Most often, in such studies, a pair of connected neurons is considered.
The transmission of an impulse through a neuron is provided by processes that occur at the nanoscale level, for example, via the overflow of sodium, potassium, chlorine and calcium ions through ion channels. In this work, we study the effect of the delay in the transmission of an impulse between connected neurons on the dynamics of the system of neurons as a whole. The FitzHugh-Nagumo model will be used to describe one neuron. And the system as a whole will be an oriented ring of neurons connected one after another. In this case, any connection between neurons will be considered with a constant delay. The main goal is to determine the type of behavior depending on the delay time, as well as to test the hypothesis that there is no direct relationship between the number of neurons and the period of self-oscillation of the impulse in the system. For this, in the first part of the work, the analysis of equilibrium states and their stability is carried out. In the second part of the work, the results obtained are verified numerically.
2. Neuron network
We will consider a group of neurons, where each neuron can have unidirectional connections with other neurons. That is, the axon of a neuron can be called its output, and for another neuron it will be an input (through a bunch of synapses and dendrites). To describe one neuron, will use the FitzHugh-Nagumo model [7-9]:
where u(t) corresponds to the potential of the neuron at time t, and v(t) is the function of the internal state at time t, a, b, y are positive constant parameters of the neuron , and I is the external current to the neuron.
As an external current I for a neuron, we will consider the total impact of other neurons, and the output will be the value of the potential of a neuron with a delay t (delay for the transmission of an impulse through an axon). In this case, the effect between neurons will be sigmoidal (that is, it does not depend on the postsynaptic neuron). Thus, the model of the ¿-th neuron will be as follows:
—au + (a + 1)u2 — u3 — v + I, bu — yv,
(1)
—au + (a + 1)u2 — u3 — Vi + ctanh(uj),
bui — y Vi.
2.1. Ring of connected neurons
In this section, we will consider a specific network of neurons in the form of an oriented ring. That is, the input of the neuron i will be the output of the neuron i — 1 (i > 2), and the input of the neuron 1 will be the output of the neuron n. Moreover, all neurons are identical in their parameters. A schematic representation of the considered network in Fig. 1. We will write a model of the considered network of neurons:
/
S.
FIG. 1 . Schematic representation of an oriented ring of connected neurons
« i = — a«i + (a + l)«! — w3 — vi + c tanh(«n),
V i = 6«i — yvi,
« 2 = —a«2 + (a + 1)«2 — «2 — V2 + c tanh(w!),
v 2 = bu2 — Y-2, (3)
«n = — au„ + (a + 1)«^ — «n — Vn + c tanh(un_1), Vn = b«n — Y-n,
where «¿(t) corresponds to the potential of the ¿-th neuron at time t, and ¿¿(t) is the function of the internal state at time t, a, b, 7 are positive constant parameters of neurons, c positive constant, connection strength, «T (t) = «¿(t — t), i = 1, 2,.. .n.
2.1.1. Find equilibrium states. To search for equilibrium states of the system, it is necessary to solve:
0 = — a«i + (a + 1)«i — «3 — ¿1 + c tanh(«n), 0 = bwi — Y-i,
0 = —a«2 + (a + 1)«2 — «2 — ¿2 + c tanh(«i),
0 = b«2 — Y-2,
0 = -aun + (a + 1)u„ — un — + c tanh(un_i), 0 = 6u„ — Yv„,
(4)
b Ui
i € [1 . . .n],
c tanh(ui_i) = u3 — (a + 1)u
+ a + b
1 Y ,
i € [2 . . .n],
(5)
ctanh(un)
if — (a + 1)u1 + (a + b ) ui.
v
2
u
Some of the equations in the system can be represented as:
( f(ui-i) = g(ui), I f (u„) = g(ui),
(6)
where f (u) = ctanh(u), g(u) = u3 — (a + 1)u2 + (a +—) u. Let — > -(a2 — a +1). This assumption fits the
V 7/ Y 3
main parameters of the neurons under consideration [2,10], where a < 1, and — roughly the same 7. Throughout the rest of this section, we will assume everywhere that the assumption is fulfilled.
Proposition 1. In the described system, the equilibrium states of neurons are the same, that is u = uj, v = vj vi, j.
For the proof, we show that the function g(x) will be increasing. Derivative g(x): g'(u) = 3u2 —2(a+1)u+a+ —.
7
b
1
b
Then, if — > - (a2 — a +1), then (a + 1)2 — 3 a +— < 0, which coincides with the determinant g(u)/2. Hence,
y 3 \ yj
g(u) is an increasing function (the leading coefficient is positive). Now, for the system (6), the functions f (u) and g(u) are increasing, which means that the assumption that some two connected neurons have different values of the equilibrium state potential leads to contradiction. That is, ui = ujvi, j, hence vi = Vj vi, j. ■ Now, the system (5) can be reduced to an equivalent:
Vi Ui
b u
7
i G [1 . . .n], ui, i G [2 ... n],
ctanh(ui) = uf — (a + 1)u2 + ^a + Y^j u1.
(7)
Now the search for equilibrium states has been reduced to finding the roots of the equation c tanh(ui) = u3 — (a + 1)u2 + ^a + . Note that u3 =0 is always a root. The right side of the equation has only one inflection point and it is located on the right semiaxis. Now let's highlight a few cases:
1) a +— > c, then the considered equation has no negative solutions, but positive solutions can be 0, 1 or 2 (as c
Y
grows, the number of roots grows from 0 to 2, 1 root at the only critical value of c).
2) a +— < c, then the considered equation has only one negative and one positive solution.
Y
2.1.2. Define the stability of equilibrium states. To determine the stability of the equilibrium states of the system of neurons under consideration, we linearize the system in the general equilibrium state ui = u3, vi = v3, i G [2... n]. For this, we rewrite the system (3) as U(t) = AU(t) + BU(t — t), where U(t) = (ui(t), vi(t),..., u„(t), v„(t))T:
A =
—a —1 0 0 •• • 0 0 0 0 •• •0 0 c' 0
b —Y 0 0 •• • 0 0 0 0 •• •0 0 0 0
0 0 —a — 1 •• •0 0 c 0 •• •0 0 0 0
0 0 b —Y • • •0 0 , B = 0 0 •• •0 0 0 0
0 0 0 0 •• • —a —1 0 0 •• • c' 0 0 0
0 0 0 0 •• •b —Y 0 0 •• •0 0 0 0
(8)
where a = 3u2 — 2u1(a + 1) + a, C
cosh (ui)
Next, we find solutions to the characteristic equation det (A — A/ + Be Ar) =0 to determine stability:
= -A t
—a — A —1 b —7 — A 0
0 0
c'e-Ar
0 0
aA
0 0 —1 —7 — A
0 0
c'e-Ar
0
00 00 00
0 —a — A —1 0 b —y — A
(9)
0
To find the determinant, let's write it down on the top line. The first two elements will reduce the problem to a similar one with a lower dimension and this part of the determinant is equal to ((—a — A)(—7 — A) + b)n. The third element in the top line is unique and sets a set of elements in the matrix and the part corresponding to it is equal to
— (c'e-Ar (—7 — A))n. As a result:
((—a — A)(—7 — A) + b)n — (c'e-Ar (—7 — A))n = 0, (10)
(a + A)(y + A) + b = c' e-Ar + 2nk ¿¿(7 + A), k g [0 ...n — 1] . (11)
To determine stability in the system in the absence of delay, take t = 0. Then:
A2 + A(a + 7 — c'e 2nk ¿) + a7 + b — cV 2nk 4 = 0, k g [0 ...n — 1]. (12)
For even odd n, these equations can be split into pairs (k; — k) and there remains one equation for k = 0. For
n
even n these equations can be split into pairs (k; — k) and there will remain two equations for k = 0 and k = —. Next, we define the stability of the selected groups. For k = 0:
A2 + A(a + 7 — c') + a7 + b — c'7 = 0. (13)
This pair of roots will be stable if and only if:
(14)
I a + 7 — c' > 0, \ a + b — c' > 0.
That is a + min ^7, > c'. And, for the zero state of equilibrium:
( b \
a + min 7, — > c. (15)
7
n
For k = — (n:2) : 2 v 7
A2 + A(a + 7 + c') + a7 + b + c'7 = 0. (16)
This pair of roots will be stable if and only if:
i a + 7 + c' > 0,
| a + 7 + c' > <">
That is a + min ^7, > —c'. For a zero equilibrium state, this condition is automatically satisfied, since c > 0.
n
For couple (k; —k) (k = 0, k = ^) :
^A2 + A(a + 7 — c'e 2nki) + «7 + b — c'7e 2nk ^A2 + A(a + 7 — c'e + «7 + b — c'7e j
(18)
(A2 + A(a + 7) + «7 + b)2 — 2 (A2 + A(a + 7) + «7 + b) c'(7 + A) cos ^ + c'2(7 + A)2 = 0, (19)
A4 + 2 (a + 7 — c' cos A3 + ((a + 7)2 + 2(a7 + b) — 2(a + 27)c' cos ^ + c'2) A2 +
2 ((a + 7)(a7 + b) — (a + 7)c'7 cos ^ — (a7 + b)c' cos ^ + c'2^) A+ (20)
(a7 + b)2 — 2(a7 + b)c'7cos ^ + c'272 = 0. An equation of the 4-th degree with real coefficients is obtained, its stability can be checked by the Hurwitz or Mikhailov criterion.
Now, we will find candidates for bifurcation from stable to unstable equilibrium or vice versa (when there is a purely imaginary eigenvalue). To do this, we take A = iw:
(a + iw)(7 + iw) + b = c'e_iwT + 2nr ^7 + iw), (21)
{a7 — w2 + b + iw(a + 7) = c'(7 cos t ' + w sin t ' + i(w cos t ' — 7 sin t ')),
1 2nfr (22)
t ' = wt — ^,
n '
0
aY — w2 + b = c'(y cos t ' + w sin t '),
w(a + y) = c'(w cos t' — y sin t'), (23)
t' = wt — ^,
n '
c' (w2 + y2) cos t' = a (w2 + y2) + bY,
c' (w2 + y2) sint' = —w (w2 + y2) + bw, (24)
t ' = wt — ^,
n
n
let t' = - + nk, k g z, then:
tan (wt_) = ~(((2+72)+b(
V n / a(i2+72 )+6y
sign (cos (wt — ^)) = sign (a (w2 + y2) + bY), (25)
c'2 (w2 + y2)2 = (a (w2 + y2) + bY)2 + (—w (w2 + y2) + bw)2 ,
tan (wT - ) = -(((2+72)+b(
V n / a(w2 + 72) + &7 '
sign (cos (wt — 2t) ) = sign (a (w2 + y2) + bY), (26)
c'2 (w2 + y2) = a2 (w2 + n2) + 2baY + b2 + w2 (w2 + y2) — 2bw2, tan (wT - ) = -"("2+y2)+6"
V n / a(w2+72) + &7 1
sign (cos (wt — ) = sign (a (w2 + y2) + bY), (27)
w4 + (a2 + y2 — c'2 — 2b) w2 + a2Y2 + 2baY + b2 — c'2y2 = 0.
This means that the system can have 0, 1, or 2 series of decisions. Under the condition a (w2 + y2) + bY > 0 (which is always true for the zero equilibrium state), we obtain (if the condition is not true, then the series for t' must be offset by n):
w = ^2 (a2 + Y2 — c'2 — 2b) ± ^(a2 — y2 — c'2)2 — 4b ((a + y)2 — c'2) t = 1 arctanf ^fW)^ + ^ + , k, m g z,
w \ a((2+72)+07 I ( ni ' ' '
w = ^2 (a2 + Y2 — c'2 — 2b) ± ^(a2 — y2 — c'2)2 — 4b ((a + y)2 — c'2) t = 1 arctanf -((((22++722)++6^ + , k g z.
( I a((2+72)+07 I ni '
(28)
(29)
Next, we define the conditions for the presence of candidates for a local bifurcation. That is, the conditions under which (27) has a real or positive root for w2. If:
1) c' > |a +— |, in for w = 0 (27) is negative, therefore, with respect to w2 there is exactly one positive solution.
7
2) c' < |a +— |, (27) will be positive at 0, then it is necessary and sufficient that the roots are (the discriminant is
7
positive) and the coefficient at w2 was negative. That is:
I (a2 + y2 — c'2 — 2b)2 — 4 (a2y2 + 2baY + b2 — c'V) > 0, \ a2 + y2 — c'2 — 2b < 0,
j c'4 — 2 (a2 — y2 — 2b) c'2 + a4 — 2a2Y2 + Y4 — 4a2b — 8abY — 4bY2 > 0, \ c'2 > a2 + y2 — 2b.
For 2abY + 2bY2 + b2 < 0, the system is equivalent to c'2 > a2 + y2 — 2b, otherwise:
(30)
(31)
c'2 < a2 — y2 — 2b — ^2abY + 2bY2 + b2,
c'2 > a2 — y2 — 2b + 2^2abY + 2bY2 + b2, (32)
c'2 > a2 + y2 — 2b,
c'2 > a2 — y2 — 2b + 2v/2abY + 2bY2 + b2. (33)
2
2
Next, it is necessary to determine for each candidate for bifurcation its effect on stability (increases or decreases the number of eigenvalues with a positive real part). Let t0 be a candidate for bifurcation, wo, Ao = ¿wo correspond
to it. Define the sign of the expression equation (10):
<9k(A)
dr
<9k(A) :
dr T=To
dD / dD dW 9Ä
, where SR(x) is the real part of the number x, D - characteristic
n (c'e-A°T0 (y + Aq))" Ao
where
(34)
di = n((a + Aq)(y + Ao) + b)n-1(Y + a + 2Aq) - nc'"e-"A°T0 (y + Aq)"-1(-tq(y + Aq) + 1). (35)
Since t0 and A0 satisfy the characteristic equation, then
(C'e-A0T0(—y - Aq))" = ((-a - Aq)(-y - Aq) + b)"
then
<9k(A)
dr where
((a + Aq )(y + Aq)+ 6)"Ac d2
((a + Aq)(y + Aq) + 6)Aq(y + Aq)
d3
(36)
(37)
d2 = ((a + Ao)(y + Aq) + b)"-1(Y + a + 2Aq) - ((a + Aq)(y + Aq) + b)n -tq +
Y + Aq /
(38)
d3 = toAO + (aTo + 27T0 + 1)A0 + (2a7To + 7 2to + bTo + 27 )Ao + + a7 2to + b7To + 72 — b. (39)
Multiply in (37) the numerator and denominator by the conjugate of d3 and multiply the whole expression by the square of the modulus d3:
f dift(A)
sign
\ dT
j = sign ((w2 + Y2)2 - (b2 + 2abY + 2bY2)) .
(40)
Now, we can see that the condition for increasing/decreasing the number of roots with a positive real part does not depend on a specific representative of the t series, but is a property of the entire series or depends only on w. w0 is the root of (27):
(41)
in another view:
+ (a2 + y2 - c'2 - 2b) + a2Y2 + 2baY + b2 - c'V = 0, (w2 + y2)2 + (a2 - c'2 - y2 - 2b) (w2 + y2) + b2 + 2abY + 2bY2 = 0.
(42)
Lemma 2. Let x2 + ax + b = 0 have two real roots xi, x2 (xi < x2) and b > 0, then xi — b and x2 — b will have different signs. And if a > 0, then x2 — b < 0, and if a < 0, then xi — b < 0.
Let x2 + ax + b = 0 have two real roots xi, x2 (xi < x2) and b < 0, then xi — b and x2 — b will always be positive.
Substitute the values of the roots in x2 - b:
(-a ± Va2 - 4b) Y ; a2 - 4b T aV«2 - 4b
2
b
2
(43)
then sign(x2 — b) = sig^ Va2 — 4b ^ ^ . This implies the assertions of the lemma. ■
As previously received, if |a + b| < c', that is, one solution (27) for w (up to sign). |a + b | < c', so, a < c',
77
and therefore a2 — c'2 — 72 — 2b < 0. That is, among the two roots (27) for w2 ours is larger in absolute value. If
|a + -1 > c', then a2 + 72 — c'2 — 2b < 0, hence a2 — 72 — c'2 — 2b < 0.
Y
/d^(A) 1
Applying lemma (2), the expression for sign —-— I , the inequalities obtained above for wo and the
\ dT lr=T0
equation (let there be at least one solution for wo) obtain that sign
id k(A) 1 V dT '
-1, if and only if there are two
roots (up to a sign) for w and b2 + 2ab7 + 2b72 > 0 for all t, corresponding to a smaller w (larger period of the series).
d
1
T = T0
T=T0
T=T0
1
T = T0
T =T0
Then, we can conclude that:
1) If c > a +—, that is, there are three equilibrium states (« < 0,« = 0,« > 0), while the zero equilibrium state
7
will be unstable for t = 0, and will remain unstable with increasing t (there is one series of values that changes the number of roots with a positive real part, but it only increases this number).
b
2) If c < a +—, that is, either one equilibrium state (« = 0), or three (« = 0,« > 0,« > 0).
7
For equilibrium states, the following results were obtained:
bb
1) If c' > a + 7 or c' > a +—, then the equilibrium state is unstable at t = 0. If c' < a + 7 and c' < a +—,
77
then it is asymptotically stable.
2) If c' > |a +— |, then (the equilibrium state is unstable at t = 0) the equilibrium state will be unstable for any
7
t (there is one series of values, which changes the number of roots with a positive real part, but it only increases this number).
3) Otherwise, if 2a7 + 272 + b < 0 and c'2 > a2 + 72 — 2b, that is, two series of changing the number of roots with positive real part, both of which increase this number. If 2a7 + 272 + b > 0 and c'2 > a2 — 72 — 2b + 2v/2ab7 + 2b72 + b2, that is, two series of changing the number of roots with a positive real part, one of which decreases the number of roots, the other increases (the one with a shorter period). Otherwise, there are no succession series, the asymptotic stability of the equilibrium state does not depend on t.
As a special case, for the zero equilibrium state:
1) If c > a +—, then the equilibrium state will be unstable for any t and there are two other equilibrium states.
7
2) If a +— > c > a + 7, then the equilibrium state is unstable at t = 0 and there are two series of changing the
7
number of roots with a positive real part, one of which decreases the number roots, the other increases.
3) If c < a + 7, then the equilibrium state is asymptotically stable for t = 0, for c2 < a2 — 72 — 2b +
2v/2ab7 + 2b72 + b2 it will remain so for any t, otherwise there are two series of changing the number of roots with a positive real part, one of which decreases the number of roots, the other increases.
In the results obtained, it can be noted that the number and the states of equilibrium themselves do not depend on the number of neurons in the ring. Also, the main criteria for determining the stability of a series of candidates
2n
do not depend on the number of neurons. Only the period in this series depends: the period is equal to —. Thus,
nw
the direct relationship between the number of neurons in the ring and the pulse propagation period was not confirmed analytically.
2.1.3. Numerical verification of conclusion. Let us introduce several definitions: a sequence of the first kind is a sequence {Tk} such that for the characteristic equation (10) passing through the value of this sequence leads to an increase by two the number of roots with a positive real part. Similarly, introduce a sequence of the second kind - one that reduces by two the number of roots with a positive real part.
To obtain numerical solutions of the system, the Wolfram Mathematica package and the NDSolve method were used. The initial values were zero functions and a jump at the time t = 0 for the first neuron by 0.01. That is, wi(t|i < 0) = vi(t|t < 0) = 0, «4(0) = 0.01, Wi(t|t < 0) = Vj(t|t < 0) = 0, i g [2, ... n]. Consider the following values of neuron parameters: a = 0.15, b = 7 = 0.02, c = 0.18. The value c = 0.18 lies between the value c = 0.890512, which separates one state of equilibrium from three, and between c = a + 7 = 0.17. That is, from the results obtained earlier, we can conclude that, regardless of the number of neurons, there is one equilibrium state (« = 0), which is unstable at t = 0. For this state of equilibrium, there is a series of both the first and second kind.
Next, we consider the different number of neurons in the ring:
1) A ring of two neurons. Then two series of candidates are Ti k « 1.70691 + 25.715009k of the second kind, and T2,k « 14.431569 + 16.895513k of the first kind. Then, for this unique equilibrium state, we obtain four Hopf bifurcations: at t « 1.70691 - a transition from an unstable equilibrium to an asymptotically stable one (Fig. 2, 3), then for t « 14.431569, t « 27.42192 and t « 31.327082 (Fig. 4, 5).
2) A ring of three neurons. Then two series of candidates are Tijfc « 1.70691 + 17.143339k of the second kind, and T2,k « 8.799731 + 11.263675k of the first kind. Then, for this unique equilibrium state, we obtain four Hopf bifurcations: at t « 1.70691 - a transition from an unstable equilibrium to an asymptotically stable one (Fig. 6, 7), then for t « 8.799731, t « 18.850249 and t « 20.063406 (Fig. 8, 9).
2 neurons in a ring
FIG. 2. т = О
FIG. 3. т = 1О
liiiii
1000 2000 3000
Uff*
FIG. 4. т = 2О
Fig. 5. т = 29
3 neurons in a ring
FIG. 6. т = О
FIG. 7. т = З
"J—1 500 '—'—1—1 10'00 1
0.010
0.008
0.006
0.004
0.002
500
1000
1500
2000
0.002
0.010
0.8
0.008
0.6
0.006
0.4
0.004
0.2
0.002
1000
2000
3000
4000
5000
4000
5000
0.2
0.002
0.4
0.010
0.008
0.006
0.004
0.002
500
1000
1500
2000
0.002
-0.4
FIG. 8. т = 12
FIG. 9. т = 19.З
3) A ring of four neurons. Then two series of candidates are T1jfc « 1.70691 + 12.857505k of the second kind, and « 5.983812 + 8.447756k of the first kind. Then, for this unique equilibrium state, obtain two Hopf bifurcations: at t « 1.70691 - a transition from an unstable equilibrium to an asymptotically stable one (Fig. 10, 11), then at
t « 5.983812 (Fig. 12).
4 neurons in a ring
FIG. 10. t = 0
FIG. 11. t = 2
FIG. 12. t = 7
3. Conclusion
In this work, a system in the form of an oriented ring of neurons was considered. Connection between the two neurons was delayed. To describe the behavior of one neuron, the FitzHugh-Nagumo model was used, and on its basis, the final model for the system in the form of a ring of neurons was built. For this system, the equilibrium states were determined, criteria for determining the stability of equilibrium states and, as a consequence, criteria for the presence of periodic solutions in the system were obtained. It has been shown that the period of self-oscillations in the ring of neurons does not have a linear dependence on the number of neurons in the ring. In this case, for a fixed number of neurons, the dynamics will be affected only by the total delay in the transmission of an impulse along the ring, and not by individual pairs of neurons. All results are then verified numerically.
Acknowledgements
This work was supported by the Russian Foundation for Basic Research (RFBR) project number 20-31-90036.
References
[1] Schöll E., Hiller G., Hövel P., Dahlem M. Time-delayed feedback in neurosystems. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2009. 1891 (367). 1079-1096.
[2] BuriC N., TodoroviC D. Dynamics of FitzHugh-Nagumo excitable systems with delayed coupling. Physical Review E. 2003. 6 (67). 066222 p.
[3] Song Z., Xu J., Zhen B. Multitype activity coexistence in an inertial two-neuron system with multiple delays. International Journal of Bifurcation and Chaos. 2015. 13 (25). 1530040 p.
[4] Xu C., Zhang Q., Wu Y. Bifurcation analysis in a three-neuron artificial neural network model with distributed delays. Neural Processing Letters. 2016. 2 (44). 343-373.
[5] Karaoglu E., Yilmaz E., Merdan H. Stability and bifurcation analysis of two-neuron network with discrete and distributed delays. Neurocomputing. 2016.182. 102-110.
[6] Xu C. Local and global Hopf bifurcation analysis on simplified bidirectional associative memory neural networks with multiple delays. Mathematics and Computers in Simulation. 2018. 149. 69-90.
[7] FitzHugh R. Impulses and physiological states in theoretical models of nerve membrane. Biophysical journal. 1961. 6 (1). 445-466.
[8] Nagumo J., Arimoto S., Yoshizawa S. An active pulse transmission line simulating nerve axon. Proceedings of the IRE. 1962. 10 (50). 20612070.
[9] Fedorov E.G., Popov A.I., Popov I.Y. Metric graph version of the FitzHugh-Nagumo model. Nanosystems: Physics, Chemistry, Mathematics. 2019. 10 (6). 623-626.
[10] Buric N., Grozdanovic I., Vasovic N. Type I vs. type II excitable systems with delayed coupling. Chaos, Solitons and Fractals. 2005. 5 (23). 1221-1233.