On solutions related to FitzHugh-Rinzel type model

TL;DR

This paper analyzes a FitzHugh-Rinzel type dynamical system, deriving explicit solutions for the initial value problem and incorporating control terms, with potential applications in neuroscience modeling.

Contribution

It provides explicit solutions for a nonlinear integro-differential equation derived from a FitzHugh-Rinzel model, including control modifications.

Findings

Explicit fundamental solution $H(x,t)$ determined.

Initial value problem solved in the whole space.

Explicit solutions obtained with added control term.

Abstract

A ternary autonomous dynamical system of FitzHugh-Rinzel type is analyzed. The system, at start, is reduced to a nonlinear integro differential equation. The fundamental solution $ H(x,t) $ is explicitly determined and the initial value problem is analyzed in the whole space. The solution is expressed by means of an integral equation involving $ H(x,t) $. Moreover, adding an extra control term, explicit solutions are achieved.

Equations

$$\left\{\begin{array}[]{lll}\displaystyle{\frac{\partial\,u}{\partial\,t}}=\,D\,\frac{\partial^{2}\,u}{\partial\,x^{2}}\,-\,w\,\,+f(u)\\ \\ \displaystyle{\frac{\partial\,w}{\partial\,t}}\,=\,\varepsilon(-\beta w+c+u),\end{array}\right.$$

$$f(u)=u(a-u)(u-1)(0<a<1).$$

$$\left\{\begin{array}[]{lll}\displaystyle{\frac{\partial\,u}{\partial\,t}}=\,u-u^{3}/3+I_{ext}\,-\,w\,\,+y\\ \\ \displaystyle{\frac{\partial\,w}{\partial\,t}}\,=\,\varepsilon(-\beta w+c+u)\\ \\ \displaystyle{\frac{\partial\,y}{\partial\,t}}\,=\,\delta(-u+h-dy)\end{array}\right.$$

$$\left\{\begin{array}[]{lll}\displaystyle{\frac{\partial\,u}{\partial\,t}}=\,D\,\frac{\partial^{2}\,u}{\partial\,x^{2}}\,-\,w\,\,+y\,\,+f(u)\\ \\ \displaystyle{\frac{\partial\,w}{\partial\,t}}\,=\,\varepsilon(-\beta w+c+u)\\ \\ \displaystyle{\frac{\partial\,y}{\partial\,t}}\,=\,\delta(-u+h-dy).\end{array}\right.$$

$$f(u)=-au+\phi (u)with\phi (u)=u^2(a+1-u)0<a<1$$

$$\left\{\begin{array}[]{lll}\displaystyle{\frac{\partial\,u}{\partial\,t}}=\,D\,\frac{\partial^{2}\,u}{\partial\,x^{2}}-au\,\,\,-\,w\,\,+y\,\,+\varphi(u)\\ \\ \displaystyle{\frac{\partial\,w}{\partial\,t}}\,=\,\varepsilon(-\beta w+c+u)\\ \\ \displaystyle{\frac{\partial\,y}{\partial\,t}}\,=\,\delta(-u+h-dy).\end{array}\right.$$

$$u(x,0)=u_0,w(x,0)=w_{0}y(x,0)=y_0,(x\in \Re )$$

$$\left\{\begin{array}[]{lll}\displaystyle w\,=\,w_{0}\,e^{\,-\,\varepsilon\beta\,t\,}\,+\,\frac{c}{\beta}\,(1-e^{-\,\varepsilon\,\beta\,t})\,+\varepsilon\int_{0}^{t}\,e^{\,-\,\varepsilon\,\beta\,(\,t-\tau\,)}\,u(x,\tau)\,d\tau\\ \\ \displaystyle y\,=\,y_{0}\,e^{\,-\,\delta\,d\,t\,}\,+\,\frac{h}{d}\,(1-e^{-\,\delta\,d\,t})\,-\delta\int_{0}^{t}\,e^{\,-\,\delta\,d\,(\,t-\tau\,)}\,u(x,\tau)\,d\tau.\end{array}\right.$$

$$F(x,t,u)=\phi (u)-w_0(x)e^{-\epsilon \beta t}+y_0(x)e^{-\delta dt}-\frac{c}{\beta }(1-e^{-\epsilon \beta t})+\frac{h}{d}(1-e^{-\delta dt}),$$

$$\left\{\begin{array}[]{lll}\displaystyle u_{t}-Du_{xx}+au+\int^{t}_{0}[\varepsilon e^{-\varepsilon\beta(t-\tau)}+\delta e^{-\delta d(t-\tau)}]u(x,\tau)d\tau=F(x,t,u)\\ \\ \displaystyle\,u(x,0)\,=u_{0}(x)\,\,\,\,\,\,x\,\in\,\Re.\end{array}\right.$$

$$Lu\equiv u_t-D{u}_{xx}+au+{\int }_0^t[\epsilon e^{-\epsilon \beta (t-\tau )}+\delta e^{-\delta d(t-\tau )}]u(x,\tau )d\tau ,$$

$$\hat{u}(x,s)={\int }_0^{\infty }{e}^{-st}u(x,t)dt,\hat{F}(x,s)={\int }_0^{\infty }{e}^{-st}F[x,t,u(x,t)]dt,$$

$$\hat{u}(x,s)={\int }_{\Re }\hat{H}(x-\xi ,s)[u_0(\xi )+\hat{F}(\xi ,s)]d\xi ,$$

$$u(x,t)={\int }_{\Re }H(x-\xi ,t)u_0(\xi )d\xi +$$

$$+{\int }_0^{t}d\tau {\int }_{\Re }H(x-\xi ,t-\tau )F[\xi ,\tau ,u(\xi ,\tau )]d\xi .$$

$$H_1(x,t)=\frac{e^{-\frac{x^2}{4Dt}}}{2\sqrt{\pi Dt}}{e}^{-at}+$$

$$-\frac{1}{2}{\int }_0^t\frac{e^{-\frac{x^2}{4Dy}-ay}}{\sqrt{t-y}}\frac{\sqrt{\epsilon }{e}^{-\beta \epsilon (t-y)}}{\sqrt{\pi D}}{J}_1(2\sqrt{\epsilon y(t-y)})\}dy,$$

$$H_2={\int }_0^t{H}_1(x,y)e^{-\delta d(t-y)}\sqrt{\frac{\delta y}{t-y}}{J}_1(2\sqrt{\delta y(t-y)}dy.$$

$${\sigma }^2=s+a+\frac{\delta }{s+\delta d}+\frac{\epsilon }{s+\beta \epsilon },$$

$$H=H_1-H_2,$$

$${\mathcal{L}}_{t}H\equiv \hat{H}={\int }_0^{\infty }{e}^{-st}H(x,t)dt=\frac{1}{\sqrt{D}}\frac{e^{-\frac{|x|}{\sqrt{D}}\sigma }}{2\sigma }.$$

$$\,\,L_{t}\bigg{(}\frac{1}{2}\,\,\,\,\,\int^{t}_{0}\frac{e^{-\frac{x^{2}}{4Dy}\,-a\,y}}{\sqrt{t-y}}\,\,\,\frac{\,\sqrt{\varepsilon}\,\,e^{-\beta\varepsilon\,(\,t\,-\,y\,)}}{\sqrt{\pi\,D\,}}J_{1}(\,2\,\sqrt{\,\varepsilon\,y\,(t-y)\,}\,\,\bigg{)}\,dy\,=$$

$$=\frac{\sqrt{\epsilon }}{2\sqrt{\pi D}}{\int }_0^{\infty }{e}^{-(s+a)y-\frac{x^2}{4Dy}}dy{\int }_0^{\infty }{e}^{-(s+\beta \epsilon )t}{J}_1(2\sqrt{\epsilon yt})\frac{dt}{\sqrt{t}}$$

$${\int }_0^{\infty }{e}^{-pt}\sqrt{\frac{c}{t}}{J}_1(2\sqrt{ct})dt=1-e^{-c/p}(\Re ep>0),$$

$${\int }_0^{\infty }\frac{e^{-x^2/4t}}{\sqrt{\pi t}}{e}^{-(s+\alpha )t}dt=\frac{e^{-x\sqrt{s+\alpha }}}{\sqrt{s+\alpha }}$$

$${\hat{H}}_1(x,s)=\frac{1}{2\sqrt{\pi D}}{\int }_0^{\infty }{e}^{-\frac{x^2}{4Dy}-(s+a+\frac{\epsilon }{s+\beta \epsilon })y}\frac{dy}{\sqrt{y}}=\frac{1}{2\sqrt{D}}\frac{e^{-\frac{|x|}{\sqrt{D}}r}}{r}$$

$$r^2=s+a+\frac{\epsilon }{s+\beta \epsilon }.$$

$${\hat{H}}_2={\hat{H}}_1-\frac{1}{2\sqrt{\pi D}}{\int }_0^{\infty }{e}^{-\frac{x^2}{4Dy}-(s+a+\frac{\epsilon }{s+\beta \epsilon }+\frac{\delta }{s+d\delta })y}\frac{dy}{\sqrt{y}}$$

$$\hat{H}(x,s)=\frac{1}{2\sqrt{D}}\frac{e^{-\frac{|x|}{\sqrt{D}}\sigma }}{\sigma }$$

$$|H_1|=\frac{e^{-\frac{{\eta }^2}{4Dt}}}{2\sqrt{\pi Dt}}{e}^{-at}+\sqrt{\frac{\epsilon t}{\pi D}};$$

$$|H_2|\le \sqrt{\frac{\delta t}{4\pi D}}+2\sqrt{\frac{\epsilon \delta }{\pi D}}t$$

Full text

1 Introduction

The FitzHugh-Rinzel (FHR) system [1, 2, 3, 4] is a three-dimensional model deriving from the FitzHugh-Nagumo (FHN) model [5, 6, 7, 8, 9, 10, 11, 12] developed to incorporate bursting phenomena of nerve cells. Indeed, a number of different cell types exhibit a behaviour characterized by brief bursts of oscillatory activity alternated by quiescient periods during which the membrane potential only changes slowly, and this behaviour is called bursting, see e.g. [13]. Accordingly, bursting oscillations are characterized by a variable of the system that changes periodically from an active phase of rapid spike oscillations to a silent phase. These phenomena are becoming increasingly important as they are being investigated in many scientific fields. Indeed, phenomena of bursting have been observed as electrical behaviours in many nerve and endocrine cells such as hippocampal and thalamic neurons, mammalian midbrain and pancreatic in $\beta-$ cells, see e.g. [1] and references therein. Also, in the cardiovascular system, bursting oscillations are generated by the electrical activity of cardiac cells that excite the heart membrane to produce the contraction of ventricles and auricles [14]. Furthermore, bursting oscillations represent a topic of potential interest in dynamics and bifurcation mechanisms of devices and structures and in the analysis of nonlinear problems in mechanics [15, 16, 17, 18, 19]. Recent studies proved that the development of this field helps in the studying of the restoration of synaptic connections. Indeed, it seems that nanoscale memristor devices have potential to reproduce the behaviour of a biological synapse [20, 21]. This would lead in the future, also in case of traumatic lesions, to the introduction of electronic synapses to connect neurons directly.

The paper is organized as follows. In section 1.1 the mathematical problem is defined and the state of the art with the aim of the paper are discussed. In section 2, the explicit expression of the fundamental solution $H(x,t)$ is achieved and some properties are proved. In section 3 the integral solution for the initial value problem is given. In section 4 the insertion of an extra term allows us to obtain explicit solutions for the model.

1.1 Mathematical considerations, state of the art and aim of the paper

Generally, denoting by $D,\varepsilon,\beta,c$ constant parameters, the (FHN) model is a p.d.e. system such that

[TABLE]

where function $\,f(u)\,$ depends on the reaction kinetics of the model. In the literature $f(u)$ can assume a piecewise linear form, see, e.g. [22] and reference therein, or $f(u)=u-u^{3}/3$ [12]. However, in general, one has [5, 13]:

[TABLE]

As for the FitzHugh-Rinzel model, most of the articles consider the following system characterized by three o.d.e.:

[TABLE]

where $I_{ext},\varepsilon,\beta,c,d,h,\beta,\delta$ indicate arbitrary constants.

In this paper, in order to evaluate the contribute of a diffusion term, the following FitzHugh-Rinzel type system is considered:

[TABLE]

Indeed, the second order term with $D>0$ represents just the diffusion contribution and it can be associated to the axial current in the axon. It derives from the Hodgkin-Huxley (HH) theory for nerve membranes where, if $d$ represents the axon diameter and $r_{i}$ is the resistivity, the spatial variation in the potential $V$ gives the term $(d/4r_{i})V_{xx}$ from which term $D\,u_{xx}$ derives [23].

Moreover it is also assumed $\,\beta,\,\,d,\,\varepsilon,\,\,\delta\,\,$ as positive constants that together with $c,\,\,h,\,$ characterize the model’s kinetic.

Model (1.4) can be considered as a two time-scale slow-fast system with two fast variables $(u,w)$ and one slow variable $(y)$. However, if for instance, $\varepsilon=\delta$ the system can be considered as a two time-scale with one fast variable $u$ and two slow variables $(w,y)$. Otherwise, if $\delta$ and $\varepsilon$ have significant difference, it can also be considered as a three-time-scale system with the fast variable $u,$ the intermediate variable and the slow variable [24].

As for function $f(u)$ one considers the non-linear form expressed in formula (1.2). As a consequence, it results

[TABLE]

Then, the system (1.4) becomes

[TABLE]

Indicating by means of

[TABLE]

the initial values, from $(\ref{15})_{2,3}$ it follows:

[TABLE]

Consequently, denoting the source term by

[TABLE]

problem (1.6)-(1.7) can be modified into the following initial value problem $\,{\cal P}$:

[TABLE]

As for the state of art, mathematical considerations allow to assert that the knowledge of the fundamental solution $H(x,t)$ related to the linear parabolic operator ${L}:$

[TABLE]

leads to determine the solution of $\,{\cal P}$. Indeed, if $\,F(\,x,\,t,\,u\,)\,$ verifies appropriate assumptions, through the fixed point theorem, solution can be expressed by means of an integral equation, see f.i. [25, 26].

Moreover, according to [26], when operator $L$ assumes a similar but simpler form, many properties and inequalities are achieved.

The aim of the paper is to explicitly determine the fundamental solution $H(x,t)$ which involves naturally the diffusion constant $D.$ Then, the initial value problem in all the space is analyzed and the solution is deduced by means of an integral equation. Moreover, using a method of travelling wave, solutions of a modified FitzHugh-Rinzel type system have been explicitly determined pointing out the influence of the diffusion parameter D.

2 Fundamental solution and its properties

Indicating by $\,T\,$ an arbitrary positive constant, let us consider the initial- value problem (1.10) defined in the whole space $\Omega_{T}:$

$\ \Omega_{T}=\{(x,t):x\in\Re,\ \ 0<t\leq T\},$

and let us denote by

[TABLE]

the Laplace transform with respect to $\,t.\,$ If $\,\hat{H}(x,s)\,$ expresses the ${\cal L}_{t}$ transform of the fundamental solution $\,H(x,t),\,$ from (1.10) it follows:

[TABLE]

and formally it follows that

[TABLE]

So that, denoting by $J_{1}(z)\,$ the Bessel function of first kind and order $\,1,\,$ let us consider the following functions:

[TABLE]

[TABLE]

Besides, by setting

[TABLE]

and by denoting

[TABLE]

the following theorem holds:

Theorem 1**.**

In the half-plane $\Re e\,s>\,max(\,-\,a,\,-\beta\varepsilon,-\delta d\,)\,$ the Laplace integral $\,{\cal L}_{t}\,H\,\,$ converges absolutely for all $\,x>0,\,$ and it results:

[TABLE]

Moreover, function H(x,t) satisfies some properties typical of the fundamental solution of heat equation, such as:

a)   $H(x,t)\,\,\in C^{\infty},\,\,\,\,$  $\,\,\,t>0,\,\,\,\,x\,\,\,\in\Re,$

b)   for fixed $\,t\,>\,0,\,\,\,H\,$ and its derivatives are vanishing esponentially fast as $\,|x|\,$ tends to infinity.

c)   In addition, it results $\displaystyle\lim_{t\,\rightarrow 0}\,\,H(x,t)\,=\,0,$ for any fixed $\,\eta\,>\,0,\,$ uniformly for all $\,|x|\,\geq\,\eta.\,$

Proof.

Since for all real $\,z\,$ one has $\,|J_{1}\,(z)|\,\leq\,1,$ the Fubini -Tonelli theorem assures that

[TABLE]

[TABLE]

and being

[TABLE]

[TABLE]

it results:

[TABLE]

where

[TABLE]

Besides, since Fubini -Tonelli theorem and (2.19) one has:

[TABLE]

from which, since (2.20),

[TABLE]

is deduced.

Besides, by means of property of convolution for which $f*g=g*f,$ since (2) and (2.15), property a) is evident. Moreover, properties b) and c) are proved following theorem 3.2.1 of [25]. In particular, as for property c), for $|x|\geq\eta$ and since $|J_{1}(z)|\leq 1,$ it results:

[TABLE]

[TABLE]

from which property follows.

  ∎

Now, introducing the following functions:

[TABLE]

it results:

[TABLE]

[TABLE]

Moreover, by denoting

[TABLE]

the convolution with respect to $t,$ for $\,t\,>\,0,\,$ as proved in [26] by means of formula (20),(21) and (24), it results:

[TABLE]

where $\,K_{\varepsilon}\,$ is given by

[TABLE]

Hence, the following theorem can be proved:

Theorem 2**.**

For $\,t\,>\,0,\,$ it results ${L}\,H=0,$ i.e.

[TABLE]

Proof.

Let us consider that:

[TABLE]

Accordingly, given relation (2.30), one has

[TABLE]

Besides, considering that:

[TABLE]

one has:

[TABLE]

where it results:

[TABLE]

So that, denoting by

[TABLE]

one has:

[TABLE]

Consequently, since for equation (2.30) one has $K_{\varepsilon}(x,t)=\,H_{1}(x,t)\ast e^{\,-\varepsilon\beta\,t},$ by means of Fubini -Tonelli theorem and (2.28), it is proved that:

[TABLE]

On the other hand, the convolution $e^{-\,\delta d\,t}*H(x,t)$ is given by

[TABLE]

with

[TABLE]

[TABLE]

As a consequence, it results:

[TABLE]

Therefore, given relations (2.30), (2.39), (2.41), theorem holds.

∎

3 Solution related to the (FHR) problem

To provide the solution by means of the integral expression (2), some convolutions need to be determined.

In order to evaluate $\int_{0}^{t}\,d\tau\,\int_{\Re}\,H(\xi,\,\tau)\,\,d\xi,$ let us start to observe that

[TABLE]

with

[TABLE]

Consequently, for (2.38), one has:

[TABLE]

So that, according to (2.17), it results:

[TABLE]

Now, let us evaluate

[TABLE]

Considering $(\ref{A22})_{3}$, after an integration by parts, one obtains:

[TABLE]

and, for (2.38), it results:

[TABLE]

Moreover, since (2.17) and (3.46), one deduces:

[TABLE]

Now, let us denote by

[TABLE]

the convolution with respect to the space $x,$ and let

[TABLE]

Considering (2.41) and (3.47) one has:

[TABLE]

Moreover, it results:

[TABLE]

where

[TABLE]

Consequently, given (1.9) and (2), for (3.44),(3.50),(3.51) and (3.52), it results:

[TABLE]

and this formula, together with relations (1.8), allow us to determine also $\,v(x,t)\,$ and $\,y(x,t)\,$ in terms of the data.

4 Explicit solutions

Several methods have been developed to find exact solutions related to partial differential equations [27, 28, 29, 30, 31]. In this case, by referring to [7], an extra term is added in order to achieve some solutions. Accordingly, let us consider

[TABLE]

where $k\neq 0$ and let us assume $\varepsilon\beta=\delta\,d$ and $f(u)=2\,u\,(a-u)\,(u-1).$

Under these conditions, problem (1.10) turns into:

[TABLE]

where, by denoting $\displaystyle\varphi_{1}=2\,u^{2}\,(\,a+1\,-u\,)+k\int_{0}^{t}\,\,e^{\,-\,\varepsilon\,\beta\,(\,t-\tau\,)}\,\,u^{2}(x,\tau)\,\,d\tau,\,$ it results:

[TABLE]

and

[TABLE]

In order to find explicit solutions, let us introduce

[TABLE]

obtaining, from system (4.54), the following equation:

[TABLE]

Now, let us consider

[TABLE]

that is a solution of Riccati type equation:

[TABLE]

and let us assume

[TABLE]

Since

[TABLE]

in order to satisfy equation (4), constant $b$ needs to assume the following expression:

[TABLE]

and, moreover, it has to be

[TABLE]

and

[TABLE]

So, if for instance, it is assumed $\varepsilon\beta=1,$ $z_{0}=0$, $\varepsilon+\delta=0.04,k=0.01$ and $C=1$ with $D>0.019,$ for $A=\sqrt{D},$ by introducing

[TABLE]

equation (4.60) gives

[TABLE]

In Fig.1, solutions u(z) expressed by means of formula (4.64) are illustrated for different values of $D,$ by showing that the amplitude increases as $0<D<1$ increases.

Remark When fast variable $u$ simulates the membrane potential of a nerve cell, while slow variable $w$ and superslow variable $y$ determine the corresponding ion number densities, model (4.54) with its solutions can be of interest in applications to understand how impulses are propagated from one neuron to another. Moreover, as similarly underlined in [4], the knowledge of exact solutions can help in testing different applications of models in neuroscience.

Acknowledgements

The present work has been developed with the economic support of MIUR (Italian Ministry of University and Research) performing the activities of the project ARS$01_{-}00861$ “Integrated collaborative systems for smart factory - ICOSAF”.

The paper has been performed under the auspices of G.N.F.M. of INdAM.

The authors are grateful to anonymous referees for their helpful comments.