# On an Intrinsic Stochastic Fitzhugh: Nagumo Model

^{*}

**Corresponding Author:**Elazab NS, Department of Mathematics, Majmaah University, Saudi Arabia, Tel: 966 16 404 4444, Email: [email protected]

*
Received Date: Sep 19, 2017 /
Accepted Date: Oct 13, 2017 /
Published Date: Oct 27, 2017 *

### Abstract

The Fitzhugh-Nagumo model for excitable systems with a high excitation parameter solves the question of selfoscillatory and self-adaptivity in these systems. This is not the case in systems with low excitation parameter. An intrinsic stochastic model that accounts for endogenous fluctuations is proposed. This model solves the question of self-oscillatory and self-adaptivity in systems with low excitation parameter.

**Keywords:**
Intrinsic stochasticity; Self adaptivity; Self-oscillatory; Fitzhugh-Nagumo model

#### Introduction

A first model that describes an excitable membrane was proposed by Hodgkin and Huxely (HH) [1]. This model solved the question of self-oscillatory in an excitable system that is oscillations between resting and ring membrane potentials, through external inputs from ion channels (or extrinsic noise).

Indeed the system undergoes intrinsic noises from randomness coherent with the processes of opening and closing the ion channels. This had been suggested in recent works [2,3]. A simple model that maintains the main aspects of the HH model equation had been proposed by Fitzhugh [4] and, Arimoto and Nagumo [5] (FHN). It reads

(1)

where u is the membrane potential and v is the recovery current [6-8]. In the equation (1),0<ε ≪ 1, a is the refractory parameter, 0<a<1, and b is the excitation parameter [9-13]. In case of a high excitation parameter b, *b*>(1−*a*)^{2}/4, the eqn. (1) shows an excitable system with a single equilibrium state which is a stable spiral [14-18]. In this case the phase portrait in the uv-plane shows spiraling trajectories. That is in an FHN system with one stable equilibrium state, the question of self-oscillatory was also solved as in the HH model. Consequently, high excitation is sufficient for self-oscillatory and self-adaptivity [3] in excitable systems.

Numerous studies of the effects of induced-noises in a stable FHN system, coherent to input resonances, had been carried out in the literature [6-18]. In these works induced-noises had been considered either in the activation potential or in the recovery current equations. The phase portrait for stochastic FHN systems shows an induced limit cycle solution. Further, intrinsic stochasticity had been introduced in FHN systems empirically apart from some works [18], where two mechanisms had been suggested. Also, everywhere in the literature it had been assumed that the stochastic noise is Gaussian. We think that, after a recent review in this area [13], the effect of intrinsic stochasticity on a bitable FHN system had not been carried out yet in the literature. This is the case that will be considered here. The mechanism suggested, accounts for endogenous fluctuations in both the activation potential and the recovery current in the absence of external resonances. Which is completely a new mechanism?

We shall present for an approach that an intrinsic stochasticity is induced by the fluctuations in the activation potential and in the recovery current due to the successive opening and closing of channels. Indeed these fluctuations enhance activation near the equilibrium states. This will be clarified later on in theorem 2.1.

In an excitable system with a low excitation parameter b, where *b*<(1−*a*)^{2}/4. The FHN equation describes a bistable medium where these two stable equilibrium states are and , . In this case the question of self-oscillatory is not evident by the simple equation (1). It needs further investigations different from those existing in the literature for a FHN system with one stable equilibrium state. The analysis of eqn. (1) shows that when the system starts from near the state of zero potential (resting state), then u evolutes towards the state (ring state) stimulated by the recovery current, with *v*a then the solution of (1) evolutes towards whatever the behavior of the recovery current. In this case, it was claimed in the literature that the system will return to the state through a long excursion [12]. We think that this do not hold due to the fact that; as the equilibrium states (0,0) and are hyperbolic then an FHN system attains these states asymptotically. On the other hand the numerical solution of the eqn. (1) by using Runge- Kuttamethod does not con rm this statement. In section 3, it will be shown that the solution of eqn. (1) evolutes towards and does not return to . Thus the FHN model with low excitation parameter is not self-adaptive or self-oscillatory.

We think that the system returns to if it is affected by a great stimulus that may arise from endogenous fluctuations (intrinsic noise). Indeed the duration of the potential components in different levels may depend on the strength of the stimulus for intensities near the threshold value. This is accompanied by a long duration of each level (or stage). The duration accounts for the latent period, ring, overshooting, depolarization and hyperpolarization periods. The successive repetition of this sequel may lead to fluctuations in the current. Alternatively, fluctuation in the potential may be argued to the random alteration of the nerve tissue from being a passive conductor to be an active one. Or, fluctuations may be argued to the low threshold of excitability of a nerve tissue. We may think that a model that describes the time evolution of an excitable medium may not be deterministic. Due to excitability, a FHN system may undergo fluctuations, so that we may write

(2)

where is an ensemble average over the space of all realizable fluctuations in FHN systems, namely

(3)

and *dm _{S}* is the measure endowed by this space. We mention that a similar analysis had been carried out for a discrete ensemble of FHN elements [13,14].

In eqn. (2) δu and δv are the fluctuation about the average with<*δu*>=<*δv*>=0. Hereafter, fluctuations are assumed to be smooth that is *δu*(*t*) and *δv*(*t*) are taken to be continuously differentiable functions.

#### The Model

By substituting from eqn. (2) into eqn. (1) and by conserving only terms quadratic in δu, we get [15]

(4)

(5)

By averaging both sides of (1) over the ensemble, we have

(6)

By using eqn. (2) into (3) and a direct calculation gives

(7)

(8)

In the eqns. (4) and (5) terms in (*δu*)^{2} and higher were neglected. By the same way the equations for u and v are given by

(9)

(10)

In the eqn. (6), we need to find<*u*^{2}>. To this end we con ne ourselves to the case when the fluctuations in the membrane potential and in the recovery current are decorrelated (decoupled)<*δvδu*>=0. By setting and by using the eqns. (4) and (5), we find closed form equations for , namely

(11)

and in the eqn. (8) initial conditions are taken σ_{i}(0)≪b, i=1,2, b is the activation parameter, practically *σ*_{i}(0) ≈ *b*/10. We mention that the case when <*δvδu*>6≠0 will be considered in section 4. The eqn. (8) integrates to

(12)

From the eqn. (9), we and that *σ*_{2}(*t*) → 0, when *t* → ∞. Also, at the equilibrium states *u*=0 and where is given in section 1, we and that *p*(0)=−*a*<0, and . Thus we have *σ*_{1}(*t*) → 0 when . Consequently the equilibrium states are unchanged due to fluctuations. Thus in the assumption made in the above, the FHN with intrinsic stochasticity is given by

(13),

where *p*(*u*) is defined in eqn. (6). We analyze the eqn. (7) and prove that, in 0 ≤ *t*<*T*, where , it describes a selfadaptive system. That is if the system is near the states then fluctuations temporate (increase or decrease) their instantaneous values so that the states are no longer stationary.

**Theorem 2.1**

The intrinsic FHN eqn. (10) describes a self-adaptive system in 0 ≤ *t*<*T*.

**Proof**

We mention that in the absence of the last term in the first eqn. (8), (namely when *σ*_{1}(0)=0) and if the system starts from the solution in eqn. 10 is respectively as they are the equilibrium states. Now in the presence of fluctuations, we assume the following:

(i) When in this case we find that , and for *t*> 0 increases so that. .

When then we find . When we find that so that for. Thus attains

(ii) When in this case we have , consequently for and decreases so that when then we find that and thus. When we find that and thus attains the state . This completes the proofs.

**Theorem 2.2**

The eqn. 10 describes a self-oscillatory system near the state , in 0<*t*<*T* where theorem 2.1, if the variance of the initial fluctuation satisfies .

**Proof**

By linearizing the eqn. (10) near , the last term in the first eqn. (10); *σ*_{1}(*t*) becomes , and the eqn. 10 becomes

(14)

or

(15)

According to when , respectively. We consider the in eqn. (11) where by using the Grown walls lemma it solves to

(16)

A similar result holds for in eqn. (12).

From the second eqn. (13) a periodic solution exists when

(17)

The above equation determines the initial variance of fluctuations in the membrane potential that induce an oscillatory behavior.

After this theorem, we find that an oscillatory solution holds for a sufficiently small initial value of the variance in fluctuations, namely *σ*_{1}(0).

In the next section we shall find numerical solutions of eqn. (10) and show that numerical results do con rm the above theorems.

#### Numerical Results

Our aim here is to solve the eqn. (10) for initial conditions and hereafter the bar on the variables will be omitted for simplicity. In the first eqn. (10) *v*(*t*) is replaced by the formal equation;

(18)

We will present for a method for finding approximate analytic solutions of in eqn. (10) [16]. A comparison between this method and some well-known ones is done in some cases. The reason for adopting this method is that it can be applied to find numerical solutions for equations with fluctuations in eqns. (18,19). It is based on using the following steps.

Inspecting the equilibrium points of equations. We have shown that in the case where the fluctuations in the membrane potential and the recovery current is decorrelated, the equilibrium points are not changed due to fluctuations. That is these equilibrium states are; and .

By dividing the first eqn. (10) by and then by integrating formally to get

(19)

where *f*(*u*,*v*), *q*(*u*) and *σ*_{1}(*t*) are given in eqn. (10).

In an analog to the discritization made for finding the fixed point numerically, the eqns. (1) and (2) are written in the form ( for *n*>1)

(20)

For *n*=0, . For more details [16].

Now we give some numerical solutions of eqn. (1) for initial conditions Numerical results for the membrane potential calculated by using Runge-Kutta method and by using the method presented in this section for the second approximation, namely *u*_{2}(*t*) ( when *σ*_{1}(0)=0) that is in the absence of fluctuations. The results are solid and dotted curves respectively. The specific values of the parameters are given in the legend. The two solutions show the same qualitative behavior for the potential. That is the potential when *t* → ∞ and *u*(*t*) does not return to the state *u*=0, which does not agree with that claimed [12] (namely the claim that *u*(*t*) reaches and returns to u=0 after a long excursion).

#### Fluctuations-Coupling Effects

Here, we consider the effects of coupling between the fluctuations in the membrane potential and the recovery current, namely when<*δvδu*>=*σ*_{12}(*t*)≠0. From the eqn. (4) and (5) the closed form equations for *σ*_{1}, *σ*_{2}, and *σ*_{12} are given by

(21)

where *p*(*x*)=−3*x*^{2}+2(1+*a*)*x*−*a*. It is worth noticing that, in this general case, the FHN intrinsic stochastic model is given the equations in eqn. (21) and equation

(22)

These five equations have to be solved with initial conditions namely for given *u*(0), *v*(0), *σ*_{1}(0), *σ*_{12}(0), and *σ*_{2}(0).

By iteration, the solution of eqn. (21) can be written as

(23)

Where . We notice that the matrices *H*(*t*_{1}) and *H*(*t*_{2})do not commute, that is the commutator

By introducing is the time ordering operator, namely

(24)

The eqn. (23) can be written in the form

(25)

Now as

(26)

where .

By using Zassenhaus formula [17] for non -commutative matrices

(27)

By considering the norm of the commutators where λ* _{ij}* are the eigenvalues of the matrix . When

*i*=2,3 we find that . To carry out numerical computations we use the eqn. (27) by neglecting C3 and higher limiting calculations for . Numerical results for the membrane potential, recovery current, mean square of the fluctuations in the potential and recovery current and the mean of the correlated fluctuations in both are displayed in the same initial conditions as respectively.

#### Conclusions

We have constructed an intrinsic stochastic FHN-model, for systems with low excitation parameter that accounts for endogenous fluctuations. A closed form for the set of equations for the ensemble averages of the membrane potential, recovery current and variances in their fluctuations had been given in eqns. (22) and (21). Theoretical proofs had shown that a system, which is described by this model, is self-adaptive and self-oscillatory. Numerical results had been carried out by including fluctuations effects and they confirmed the theoretical predictions. Consequently this model conserves the main features as in an excitable system with high excitation parameter. The model presented, accounts for fluctuations about the mean and it may be considered as a simple model for describing smooth-noisy systems.

#### References

- Hodgkin AL, Huxley AF (1952) A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology 117: 500-544.
- Deflice LJ (2000) Linear analysis in non-linear systems; signal to noise ratios in noisy neurons. Chaos Solitons & Fractals 11: 1823-1826.
- Vajtai R, Kiss LB, Gingle Z (2000) A self-adaptive stochastic resonator with logarithmic transfer. Chaos Solitons & Fractals 11: 1933-1935.
- Fitzhugh R (1961) Impulses and physiological states in theoretical models of nerve membrane. Biophysical Journal 1: 445-466.
- Arimoto S, Nagumo JS, Yoshzawa S (1962) An active pulse transmission line simulating nerve axon. Proceedings of the IEEE 50: 2061-2071.
- Lindner B, Schimansky-Geier L (2000) Physical Review E 61: 6103.
- Collins JJ, Chow CC, Imho TT (1995) A periodic stochastic resonance. Physical Review E52, R 3321-3324.
- Pei XA, Bachman K, Moss F (1995) The detection threshold, noise and stochastic resonance in the Fitzhugh-Nagumo neuron model. Physics Letters A 206: 61-65.
- Longtin A (1993) Stochastic resonance in neuron models. Journal Statistical Physics 70: 309-327.
- Chilov DR, Longtin A, Muller-Gerking J (1997) Stochastic resonances of neural ensembles. Physical Review E 55: 1798-1808.
- Longtin A, Chilov DR (1998) Stochastic and deterministic resonances in excitable systems. Physical Review Letter 82: 4012.
- Murray JD (1993) Mathematical Biology. (2nd edn.), Springer.
- Lindner B, Garcia-Ojalvo J, Neiman AB, Schimansky-Geier L (2004) Effects of noise in excitable systems. Physical Review 392: 321-424.
- Zaks MA, Neiman AB, Feisted S, Schimansky-Geier L (2003) Physical Review E 68: 196-206.
- Birmingham TJ, Bornatici M (1972) Resonant diffusion in the presence of strong plasma turbulence. Phys. Fluids. 4: 1778-1789.
- Abdel-Gawad HI (2008), On a finite memory transport in the Fisher equation. Applied Mathematical Model 32: 1883-1893.
- Brian CH (2003) Lie group Lie algebra and representation. An elementary introduction Springer.
- Tang J, Jia Y, Yi M, Ma J, Li J (2008) Multiplicative-noise-induced coherence resonance via two mechanisms in bistable neural models. Physical Review E 77: 601905.

Citation: Elazab NS (2017) On an Intrinsic Stochastic Fitzhugh: Nagumo Model. J Appl Computat Math 6: 376. Doi: 10.4172/2168-9679.1000376

Copyright: ©2017 Elazab NS. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.