Loading [MathJax]/jax/output/SVG/jax.js
Research article

Some results on transcendental entire solutions to certain nonlinear differential-difference equations

  • In this paper, we study the transcendental entire solutions for the nonlinear differential-difference equations of the forms:

    f2(z)+˜ωf(z)f(z)+q(z)eQ(z)f(z+c)=u(z)ev(z),

    and

    fn(z)+ωfn1(z)f(z)+q(z)eQ(z)f(z+c)=p1eλ1z+p2eλ2z,n3,

    where ω is a constant, ˜ω,c,λ1,λ2,p1,p2 are non-zero constants, q,Q,u,v are polynomials such that Q,v are not constants and q,u0. Our results are improvements and complements of some previous results.

    Citation: Nan Li, Jiachuan Geng, Lianzhong Yang. Some results on transcendental entire solutions to certain nonlinear differential-difference equations[J]. AIMS Mathematics, 2021, 6(8): 8107-8126. doi: 10.3934/math.2021470

    Related Papers:

    [1] Müzeyyen Sangurlu Sezen . Interpolative best proximity point results via γ-contraction with applications. AIMS Mathematics, 2025, 10(1): 1350-1366. doi: 10.3934/math.2025062
    [2] Naeem Saleem, Hüseyin Işık, Sana Khaleeq, Choonkil Park . Interpolative Ćirić-Reich-Rus-type best proximity point results with applications. AIMS Mathematics, 2022, 7(6): 9731-9747. doi: 10.3934/math.2022542
    [3] Khalil Javed, Muhammad Nazam, Fahad Jahangeer, Muhammad Arshad, Manuel De La Sen . A new approach to generalized interpolative proximal contractions in non archimedean fuzzy metric spaces. AIMS Mathematics, 2023, 8(2): 2891-2909. doi: 10.3934/math.2023151
    [4] Raju Gopi, Veerasamy Pragadeeswarar, Choonkil Park, Dong Yun Shin . Coupled common best proximity point theorems for nonlinear contractions in partially ordered metric spaces. AIMS Mathematics, 2020, 5(6): 6913-6928. doi: 10.3934/math.2020443
    [5] Umar Ishtiaq, Fahad Jahangeer, Doha A. Kattan, Manuel De la Sen . Generalized common best proximity point results in fuzzy multiplicative metric spaces. AIMS Mathematics, 2023, 8(11): 25454-25476. doi: 10.3934/math.20231299
    [6] Yao Yu, Chaobo Li, Dong Ji . Best proximity point of α-η-type generalized F-proximal contractions in modular metric spaces. AIMS Mathematics, 2024, 9(4): 8940-8960. doi: 10.3934/math.2024436
    [7] A. Sreelakshmi Unni, V. Pragadeeswarar, Manuel De la Sen . Common best proximity point theorems for proximally weak reciprocal continuous mappings. AIMS Mathematics, 2023, 8(12): 28176-28187. doi: 10.3934/math.20231442
    [8] Basit Ali, Muzammil Ali, Azhar Hussain, Reny George, Talat Nazir . Best proximity points in non-Archimedean fuzzy metric spaces with application to domain of words. AIMS Mathematics, 2022, 7(9): 16590-16611. doi: 10.3934/math.2022909
    [9] Arshad Ali Khan, Basit Ali, Reny George . On semi best proximity points for multivalued mappings in quasi metric spaces. AIMS Mathematics, 2023, 8(10): 23835-23849. doi: 10.3934/math.20231215
    [10] Mustafa Aslantas, Hakan Sahin, Raghad Jabbar Sabir Al-Okbi . Some best proximity point results on best orbitally complete quasi metric spaces. AIMS Mathematics, 2023, 8(4): 7967-7980. doi: 10.3934/math.2023401
  • In this paper, we study the transcendental entire solutions for the nonlinear differential-difference equations of the forms:

    f2(z)+˜ωf(z)f(z)+q(z)eQ(z)f(z+c)=u(z)ev(z),

    and

    fn(z)+ωfn1(z)f(z)+q(z)eQ(z)f(z+c)=p1eλ1z+p2eλ2z,n3,

    where ω is a constant, ˜ω,c,λ1,λ2,p1,p2 are non-zero constants, q,Q,u,v are polynomials such that Q,v are not constants and q,u0. Our results are improvements and complements of some previous results.



    Individual neurons are the basic units of the biological nervous system, they process information by working together in neuronal circuits with specific synaptic connectivity. In response to an applied current, the signal can be propagated from one neuron to the next neuron by generating electrical impulses or spikes[1]. Hodgkin and Huxley (HH) firstly proposed the mathematical description of the neural model based on a large number of biological experiments[2], which built a new field of neuronal electrophysiology and laid the foundation for theoretical neuroscience. However, the HH model is a complex nonlinear system, and is difficult to be analyzed theoretically. Subsequently, researchers have successively proposed many simplified neuron models such as the Fitzhugh-Nagumo (FN) model [3], Morris-Lecar (ML) model[4], Hindmarsh-Rose (HR) model[5], and the modified forms of them[6,7,8], to explore the physiological mechanisms of neurons.

    Previous studies have shown that many neurological dynamic diseases such as Epilepsy, Parkinson's disease, Alzheimer's disease and Schizophrenia may be induced by bifurcations caused by the changes of the regulatory parameters in the neural system, which leads to the disorders in the firing patterns of nervous systems[9,10,11,12]. In addition, the changes of firing patterns in neuron systems may cause neurons to produce a kind of plasticity, which finally affects learning, memory and other neural activities that rely on neuron-plasticity[13]. Therefore, a method is expected to be proposed to control the bifurcation mechanism of the neuron model, which can change the firing patterns from abnormal to normal under certain conditions.

    Bifurcation control has attracted increasing attention due to its potential applications in many fields, such as engineering, biomedicine, and meteorology systems etc[14]. The main idea of bifurcation control is to change the inherent bifurcation characteristics of the original system by designing the controller without changing the parameters, thereby obtaining some desirable bifurcation behaviors. This includes delaying or advancing the onset of the bifurcation points, modifying the shape or type of bifurcations, and so on. Many bifurcation control methods have been provided in previous studies[15,16,17,18,19]. Xie et al. [10] studied the dynamic behaviors of Hopf bifurcation under the control of the washout filter in HH model. The results show that the linear term of the controller can control the position of the Hopf bifurcation point, and the cubic term can change the criticality of the Hopf bifurcation point. The washout filter dynamic feedback controller is introduced into the two-dimensional HR model to change the firing pattern from type I to type II by delaying the location of its inherent Hopf bifurcation[12]. Huang et al. [19] studied the Hopf bifurcation control of the ML model with type I under the control of washout filter-aided dynamic feedback controller. However, the firing pattern cannot be changed from type I to type II by using the washout filter-aided dynamic feedback controller in ML model. Therefore, a new controller should be introduced to address this issue. In addition, the analytical expression for determining the stability index of Hopf bifurcation criticality has not been obtained. Moreover, none of these studies involved the bifurcation control of the ML model with Hc bifurcation type.

    Motivated by above analysis, in the present paper, we introduce a new controller proposed by Nguyen et al.[20] into the ML model to study the influence of the controller on the ML model with Hc bifurcation type. In addition, we give the analytical expression of stability index of the system to determine the criticality of Hopf bifurcation.

    The structure of the paper is organized as follows. In Sect. 2, we reviewed the bifurcation dynamics and the conditions for the emergence of Hopf bifurcations of the ML model. In Sect. 3, the dynamic state-feedback control law was applied to the ML neuron model with Hc bifurcation type. Also, we derived the analytical expression of the bifurcation stability index for the controlled ML model with Hc bifurcation type. Conclusions are given in Sect. 4.

    The ML neuron model is a simplified form of HH model and its mathematical description can be expressed by the following nonlinear differential equations:

    dVdt=1C{IgL(VVL)gCaMss(VVCa)gKN(VVK)},dNdt=ϕ(NssN)τN. (1)

    with

    Mss=12(1+tanh(VV1V2)),Nss=12(1+tanh(VV3V4)),τN=1/cosh(VV32V4). (2)

    where V is the membrane potential, and N is the recovery variable. I is the current stimulation and is treated as the bifurcation parameter. C is the membrane capacitance. Vi (i=K,Ca,L) represent the reversal potential, and gi (i=K,Ca,L) represent the maximum conductance. Vi (i=1,2,3,4) are constant potentials, and ϕ represents the temperature factor. Here, we named system (1) as the original system.

    There are three types of bifurcations in ML model by choosing different parameter sets, including Hopf bifurcation, saddle node bifurcation on invariant cycle (SNIC) and homoclinic (Hc) bifurcation (see Table 1). In this paper, we mainly focus on the Hopf bifurcation control of the ML model with Hc bifurcation type. The bifurcation diagram of the ML model with Hc bifurcation type is shown in Figure 1. From Figure 1 we can see that as the bifurcation parameter I increases, the ML model with Hc bifurcation type is divided into six types of dynamic regions. In region I (when I<9.949), the system has only one stable equilibrium point in the lower branch; In region II (when 9.949<I<34.94), the system contains three equilibrium points, the stable equilibrium point in the upper branch, the unstable equilibrium point in the middle branch and the stable equilibrium point in the lower branch; In region III (when 34.94<I<36.32), the system contains a stable limit cycle and three equilibrium points; In region IV (when 36.32<I<39.96), the system contains a stable limit cycle, an unstable limit cycle and three equilibrium points; In region V(when 39.96<I<40.56), the system contains a stable limit cycle, an unstable limit cycle and a stable equilibrium point; In region VI(when I>40.56), the system contains a stable equilibrium point. The corresponding phase diagrams of six typical dynamical behaviors are shown in Figure 2. Compared it with [19], the ML model with Hc bifurcation type has more dynamical behaviors than that with SNIC bifurcation type. In particular, Hc bifurcation occurs in the system when I=34.94, and the corresponding phase diagram is shown in Figure 3.

    Table 1.  Parameter setting in the simulation and experiment.
    Parameters Hopf type SNIC type Hc type
    ϕ 0.04 0.067 0.23
    gCa 4.4 4 4
    V3 2 12 12
    V4 30 17.4 17.4
    VCa 120 120 120
    VK -84 -84 -84
    VL -60 -60 -60
    gK 8 8 8
    gL 2 2 2
    V1 -1.2 -1.2 -1.2
    V2 18 18 18
    C 20 20 20

     | Show Table
    DownLoad: CSV
    Figure 1.  (a) Bifurcation diagram of the ML model with Hc bifurcation type. (b) The enlarged representation of (a). Here, the red and black lines represent stable equilibrium point and unstable equilibrium point, respectively. The green and blue circles represent stable limit cycle and unstable limit cycle, respectively. SNi (i=1,2) represent saddle node bifurcation points, H represents Hopf bifurcation point, and Hc represents homoclinic bifurcation point. IVI represent six types of dynamical behavior regions, respectively. SN1: I=9.949μA/cm2, Hc: I=34.94μA/cm2, H: I=36.32μA/cm2, SN2: I=39.96μA/cm2.
    Figure 2.  Six types of phase diagrams of the ML model with Hc bifurcation type. Here, the blue and red curves represent N-nullcline and V-nullcline, respectively. The black solid line and the black dashed line represent stable limit cycle (SLC) and unstable limit cycle (ULC), respectively. The red circle and the black circle represent the stable equilibrium point (SEP) and the unstable equilibrium point (UEP), respectively. (a) I=30μA/cm2, (b) I=30μA/cm2, (c) I=35.5μA/cm2, (d) I=37.5μA/cm2, (e) I=40μA/cm2, (f) I=45μA/cm2.
    Figure 3.  Homoclinic orbit of the ML model with Hc bifurcation type. I=34.94μA/cm2.

    Some researchers have proposed that the firing patterns can be determined by the types of bifurcation[21,22]. Specifically, type I firing pattern occurs when the system undergoes SNIC bifurcation or Hc bifurcation, while type II firing pattern occurs when the system undergoes Hopf bifurcation. For type I, the system can fire at any low frequency; for type II, the system fires from a positive minimum frequency. Figure 4 shows the firing patterns of the ML model under the three types of bifurcations. From Figure 4 we can see that the ML model with Hopf bifurcation type shows the type II firing pattern (see Figure 4b), and the corresponding time evolutions are shown in Figure 4a; the ML model with SNIC bifurcation or Hc bifurcation type both show the type I firing pattern (see Figure 4d and 4f), and the corresponding time evolutions are shown in Figure 4c and 4e, respectively.

    Figure 4.  (a) (c) (e) Time evolution diagrams of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively. (b) (d) (f) Frequency-I curves of the ML model with Hopf bifurcation, SNIC bifurcation and Hc bifurcation, respectively.

    Consider the following n-dimensional nonlinear system:

    dxdt=f(x,μ), (3)

    where xRn represents the state variables, μR represents the bifurcation parameter. Assume that the system has a fixed point x0 when the bifurcation parameter μ=μ0 and the corresponding eigenvalue λμ0=αμ0+iβμ0. The Hopf bifurcation occurs at μ=μ0 when the system satisfies the following conditions: βμ00, αμ0=0, and ddμ(αμ)μ=μ0>0. However, the eigenvalues of the system cannot be obtained easily when the system is a high-dimensional system. Therefore, the following Routh-Hurwitz stability criterion can be used to determine whether the Hopf bifurcation occurs[23]. Assume that the corresponding characteristic polynomial of system (3) is as follows:

    p(λ,μ)=p0(μ)λn+p1(μ)λn1++pn(μ). (4)

    Then, we obtain

    Hn(μ)=(p1(μ)p0(μ)00p3(μ)p2(μ)p1(μ)0p5(μ)p4(μ)p3(μ)0p2n1(μ)p2n2(μ)p2n3(μ)pn(μ)). (5)

    where pi(μ)=0, if i<0 or i>n. Then, we define

    D1(μ)=Det(H1(μ))=p1(μ),D2(μ)=Det(H2(μ))=Det(p1(μ)p0(μ)p3(μ)p2(μ)),Dn(μ)=Det(Hn(μ)). (6)

    Based on the Routh-Hurwitz stability criterion, the Hopf bifurcation occurs at μ=μ0 when the system satisfies the following conditions:

    (I) Eigenvalue crossing condition:

    pn(μ0)>0,Di(μ0)>0,         i=1,,n2,Dn1(μ0)=0. (7)

    (II) Transversality condition:

    dDn1(μ)dμμ=μ00. (8)

    In this section, we provide a dynamic state-feedback controller to control the Hopf bifurcation of the original system. Consider the following nonlinear system:

    u=u(x,y),˙y=g(x,y). (9)

    where yRm(1mn) represents the state variables of the controller, u(x,y) represents the feedback control and g(x,y) is a smooth function. In this paper, we utilize the following feedback control law proposed by Nguyen and Hong[20].

    ui(xi,yi)=k1ixi+k3i(xix0i)3liyi,˙yi=ui(xi,yi). (10)

    where x0i(i=1,2,...,m) are the equilibrium points of the system at the Hopf bifurcation point, k1i and k3i represent the linear term control gain and the cubic term control gain, respectively. li are constant parameters. Then we obtain the following controlled system:

    ˙x=f(x,μ)+u(x,y),˙y=g(x,y). (11)

    where

    u(x,y)=[u1(x1,y1),,um(xm,ym),0,,0]T,g(x,y)=[u1(x1,y1),,um(xm,ym)]T. (12)

    According to Eq (10), we have ui=0 when ˙y=ui=0. So, if x0 is the equilibrium point of the original system, (x0,y0) is the equilibrium point of the controlled system, where y0=(y01,y02,,y0m) and then y0i=(k1ix0i+k3i(x0ix01i)3)/li(i=1,2,,m). That is, the equilibrium points of the original system are the equilibrium points of the controlled model. The reason for using state feedback control to conduct the Hopf bifurcation control is that the equilibrium structure of the original system and the controlled system remains unchanged during the control process. In the present paper, we add the controller to the membrane potential V, and then we obtain the following controlled system:

    dVdt=1C{IgL(VVL)gCaMss(VVCa)gKN(VVK)}+k1V+k3(VV0)3ly,dNdt=ϕ(NssN)τN,dydt=k1V+k3(VV0)3ly. (13)

    where V0 is the equilibrium membrane potential of the original system at the Hopf bifurcation point.

    The Hopf bifurcation control of neuron models has been widely studied, including HH model, HR model, ML model with Hopf bifurcation type or SNIC bifurcation type[10,11,12,19,20]. However, the Hopf bifurcation control of the ML model with Hc bifurcation type has not been studied. In addition, the controllers used in previous researches only studied the system with single Hopf bifurcation point. Therefore, in the present paper, we utilize the dynamic state-feedback control proposed by Nguyen et al. to make the Hopf bifurcation point change from one to two. The parameter values of the model are ϕ=0.23,gCa=4,V3=12,V4=17.4,VCa=120,VK=84,VL=60,gK=8,gL=2,V1=1.2,V2=18,C=20,l=0.1, and I is treated as the bifurcation parameter.

    The Hopf bifurcation point is located at I=36.32 in the original ML model with Hc bifurcation type (Figure 1). To change the model from the Hc bifurcation type to the SNIC bifurcation type, the Hopf bifurcation point needs to be moved forward. Here, we change the Hopf bifurcation point from I=36.32 to I0=70. The equilibrium point of the original system is (V0=6.7697,N0=0.35407). Therefore, the equilibrium point of the controlled system is (V0=6.7697,N0=0.35407,y0=67.697k1), and then we obtain the Jacobian matrix of the controlled system at the equilibrium point:

    J(I0)=(0.1370+k136.30790.10.00610.23260k100.1). (14)

    Then, we obtain the characteristic polynomial of the Jacobian matrix J(I0):

    p(λ,I0)=p0λ3+p1λ2+p2λ+p3=0. (15)

    where p0=1,p1=0.1956k1,p2=0.19920.2326k1,p3=0.01896. Substituting pi (i=0,1,,m) into the crossing condition (I), we have

    p3=0.01896>0,D1(I0)=p1=0.1956k1>0,D2(I0)=p1p2p0p3=0.2326k210.2447k1+0.02=0, (16)

    Solving Eq (16), we have k1=0.0893. Then, we have:

    dD2(I)dII=I0=0.2032×1040. (17)

    which satisfy the transversal condition (II).

    The bifurcation diagram of the controlled system when k1=0.0893,k3=0 is shown in Figure 5a. As expected, the Hopf bifurcation point of the ML model with Hc bifurcation type moves forward from I=36.32 to I=70. Figure 5b shows the two-parameter bifurcation diagram of bifurcation parameter I with respect to k1, which further verifies that Hopf bifurcation occurs at I=70 when k1=0.0893. Now, the model shows a change from the Hc bifurcation type to the SNIC bifurcation type (see Figure 5a), but does not change the firing pattern of the system (see Figure 5c and 5d).

    Figure 5.  (a) Bifurcation diagram of the controlled system when k1=0.0893,k3=0,I0=70. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled model with different bifurcation parameter I when k1=0.0893,k3=0. (d) Frequency-I curves of the controlled system when k1=0.0893,k3=0.

    The results of this section show that the ML model can be transformed from the Hc bifurcation type to the SNIC bifurcation type by introducing a dynamic state-feedback control. In addition, there is only one Hopf bifurcation point in the system when k10,k3=0 (Figure 5a), that is, k1 can change the location of the Hopf bifurcation point of the ML neuron model with Hc bifurcation type, resulting in the change of bifurcation type from Hc type to SNIC type.

    Previous study has investigated the transition of ML model from SNIC bifurcation type to Hopf bifurcation type[18]. Here, we will study the transition of the ML model from Hc bifurcation type to Hopf bifurcation. First, the inherent Hopf bifurcation point in the Hc bifurcation type needs to be moved forward, and then a new Hopf bifurcation point is constructed in front of the moved inherent Hopf bifurcation point. Here, we expect to change the inherent Hopf bifurcation point from I=36.32 to I1=200, and then a new Hopf bifurcation point is constructed at I2=60.

    When I1=200, the equilibrium point of the controlled system is (V0=12.94,N0=0.52697,y0=129.4k1), and then we obtain the Jacobian matrix of the controlled system at the equilibrium point:

    J(I1)=(0.1375+k138.77600.10.00660.23010k100.1). (18)

    Then, we obtain the characteristic polynomial of the Jacobian matrix J(I1):

    p(λ,I1)=p0λ3+p1λ2+p2λ+p3=0. (19)

    where p0=1,p1=0.4676k1,p2=0.32430.2301k1,p3=0.028756. Substituting pi (i=0,1,,m) into the crossing condition (I), we obtain:

    p3=0.028756>0,D1(I0)=p1=0.4676k1>0,D2(I0)=p1p2p0p3=0.2301k210.4319k1+0.12288668=0. (20)

    Solving Eq (20), we have k1=0.34966383. Then, substituting k1=0.34966383 into the controlled system when I2=60, we obtain the equilibrium point of the controlled system (V0=6.1368,N0=0.33762,y0=21.46+3148.761k3). The Jacobian matrix of the controlled system at the equilibrium point is as follows:

    J(I2)=(0.5141+138.8506k336.05470.10.00600.233300.3497+138.8506k300.1). (21)

    Then, we obtain the characteristic polynomial of the Jacobian matrix J(I2):

    p(λ,I2)=p0λ3+p1λ2+p2λ+p3=0. (22)

    where p0=1,p1=0.1808138.8506k3,p2=0.103332.3938k3,p3=0.0178. Substituting pi (i=0,1,,m) into the crossing condition (I), we obtain:

    p3=0.0178>0,D1(I0)=p1=0.1808138.8506k3>0,D2(I0)=p1p2p0p3=4497.89857k238.4865k30.0365=0. (23)

    Solving Eq (23), we have k3=0.00205656. Figure 6 shows the dynamic behaviors of the controlled system when k1=0.34966383,k3=0.00205656. From Figure 6a we can see that the Hopf bifurcation point of the ML neuron model with Hc bifurcation type is moved forward from I=36.32 to I1=200, and a new Hopf bifurcation point is constructed at I2=60. Figure 6b shows the two-parameter bifurcation diagram of bifurcation parameter I with respect to k1, which further verifies that the model undergoes the Hopf bifurcation at I1=200,I2=60 when k1=0.0893. Figure 6c shows the time evolutions of the controlled system with different bifurcation parameter I. From Figure 6c we can see that the system fires from a positive minimum frequency, which is further verified in Figure 6d. Now, the model shows a change from Hc bifurcation type to Hopf bifurcation type (see Figure 6a), and the firing pattern changes from type I to type II (see Figure 6d).

    Figure 6.  (a) Bifurcation diagram of the controlled system with k1=0.34966383 and k3=0.00205656. (b) Two-parameter diagram of I with respect to k1. (c) Time evolution diagrams of the controlled system with different bifurcation parameter I when k1=0.34966383 and k3=0.00205656. (d) Frequency-I curves of the controlled system when k1=0.34966383 and k3=0.00205656.

    The results in this section show that the ML model can be transformed from the Hc bifurcation type to the Hopf bifurcation type by introducing dynamic state-feedback control. In addition, the model performs type II firing pattern when k10,k30, that is, the combination of k1 and k3 can change the firing patterns of the ML neuron model with Hc bifurcation type.

    As mentioned above, the criticality of the Hopf bifurcation can be changed by a dynamic state-feedback controller. In this section, we derive the stability coefficient β2 of the Hopf bifurcation by utilizing the Poincare Birkhoff (PB) normal form method. The Hopf bifurcation is subcritical when β2>0, resulting in an unstable limit cycle; the Hopf bifurcation is supercritical when β2<0, resulting in a stable limit cycle.

    Consider the following controlled system:

    ˙x=f(x,μ), (24)

    where x=(V,N,y)T, f=(f1,f2,f3)T. The equilibrium point of the system at μ0: x0(μ0)=(V0(μ0),N0(μ0),y0(μ0))T. Rewrite the controlled system as: ˙x=f(x,μ)=Jx+g(x), where J is the Jocabian matrix of the system at the equilibrium point, g(x) is the nonlinear term:

    J=fxx=x0=[f1Vf1Nf1yf2Vf2Nf3yf3Vf3Nf3y]V=V0,N=N0,y=y0g(x)=f(x,μ)Jx. (25)

    Then, we conduct the linear transformation of the controlled system: x=Py+x0(μ0), where P=[Re(v1),Im(v1),v3], Re(v1) and Im(v1) are the real and imaginary parts of the eigenvector v1 corresponding to eigenvalues λ1, and v3 is the eigenvector corresponding to the eigenvalue λ3. Then, we have:

    ˙x=P˙y=J(Py+x0(μ0))+g(Py+x0(μ0)), (26)

    Substituting Eq (25) into Eq (26), we have:

    P˙y=J(Py+x0(μ0))+f(Py+x0(μ0))J(Py+x0(μ0)), (27)

    where

    ˙y=P1J(Py+x0(μ0))+P1[f(Py+x0(μ0))J(Py+x0(μ0))]=P1JPy+P1f(Py+x0(μ0))P1JPy=Ay+P1f(Py+x0(μ0))Ay=Ay+F. (28)

    A=P1JP is the Jocabian matrix of the system, F is the nonlinear term:

    F=P1f(Py+x0(μ0))Ay. (29)

    That is,

    [F1F2F3]=P1[f1(Py+x0(μ0))f2(Py+x0(μ0))f3(Py+x0(μ0))]A[y1y2y3]. (30)

    Based on the PB normal form theory, we obtain the expression of the bifurcation stability coefficient β2:

    β2=2Re{g20g112|g11|213|g02|22ω0i+g212}. (31)

    where

    g20=14[2F1y212F1y21+22F2y1y2+i(2F2y212F2y2122F1y1y2)], (32)
    g11=14[2F1y21+2F1y22+i(2F2y21+2F2y22)], (33)
    g02=14[2F1y212F1y2122F2y1y2+i(2F2y212F2y21+22F1y1y2)], (34)
    g21=G21+2G1110w111+G1101w120, (35)
    G21=18[3F1y31+3F1y1y22+3F2y21y2+3F2y32+i(3F2y31+3F2y1y223F1y21y23F1y32)], (36)
    G1110=12[2F1y1y3+2F2y2y3+i(2F2y1y32F1y2y3)], (37)
    G1101=12[2F1y1y32F2y2y3+i(2F2y1y3+2F1y2y3)]. (38)

    Here, the vectors w111 and w120 satisfy the following linear equations:

    λ3w111=h11,   (λ32iω0)w120=h20. (39)
    h11=14(2F3y21+2F3y22),h20=14(2F3y212F3y222i2F3y1y2). (40)

    where λ3 is the real characteristic root of the Jacobian matrix, ω0 is the imaginary part of the pure virtual characteristic root of the Jacobian matrix, and ω0>0. Solving Eqs (39) and (40) we obtain:

    w111=14λ3(2F3y21+2F3y22), (41)
    w120=14(λ23+4ω20)[λ3(2F3y212F3y22)+4ω02F3y1y2]i12(λ23+4ω20)[ω0(2F3y212F3y22)λ32F3y1y2]. (42)

    Then, we obtain the nonlinear terms of the controlled system:

    f1=1c[Igl(Vl+V0+P11y1+P12y2+P13y3)12gca(1+tanh(V0+P11y1+P12y2+P13y3V1V2))(Vca+V0+P11y1+P12y2+P13y3)gk(N0+P21y1+P22y2+P23y3)(Vk+V0+P11y1+P12y2+P13y3)]+k1(V0+P11y1+P12y2+P13y3)+k3(V0+V0+P11y1+P12y2+P13y3)3l(y0+P31y1+P32y2+P33y3), (43)
    f2=ϕ[12(1+tanh(V0+P11y1+P12y2+P13y3V1V4))(N0+P21y1+P22y2+P23y3)]cosh(V0+P11y1+P12y2+P13y3V32V4), (44)
    f3=k1(V0+P11y1+P12y2+P13y3)+k3(V0+V0+P11y1+P12y2+P13y3)3l(y0+P31y1+P32y2+P33y3). (45)

    Because Ay is the first-order function with respect to y and the second or more derivative of y is required in the following calculations, the derivative of y is omitted. We only calculate the following equation: F=P1f(Py+x0(μ0)).

    Next, we calculate the stability coefficient β2 of the controlled system when k1=0.3497, k3=0.0021, μ0=I1=200. The equilibrium point of the controlled system is (V0=12.94,N0=0.52697,y0=45.2465), and the corresponding Jacobian matrix is:

    J=(0.212238.77600.10.00660.230100.349700.1). (46)

    Then, we obtain the eigenvalues and eigenvectors: λ1=0.4938i, λ2=0.4938i, λ3=0.1179.

    v1=(0.82150.00420.0090i0.11310.5588i),v2=(0.82150.0042+0.0090i0.1131+0.5588i),v3=(0.05120.00300.9987). (47)

    So, the matrix P=[Re(v1),Im(v1),v3] can be rewrite as:

    P=(0.821500.05120.00420.00900.00300.11310.55880.9987). (48)

    Then, we obtain:

    P1=(1.22683.29120.05300.521695.04460.25880.152952.80730.8505). (49)

    Substituting Eqs (48) and (49) into F=P1f(Py+x0(μ0)), we have:

    [F1F2F3]=[1.2268f13.2912f2+0.0530f30.5216f195.0446f20.2588f30.1529f152.8073f2+0.8505f3]. (50)

    Finally, substituting the parameters into the stability coefficient formula we have β2=9.9169×104<0, that is, the controlled system undergoes the supercritical Hopf bifurcation at (V0,N0)=(12.94,0.52697) when I1=200. The membrane potential curve and phase trajectory of the controlled system at the bifurcation point are shown in Figure 7a and 7b. From Figure 7b we can see that the phase trajectory of the controlled system converges to a stable limit cycle when I1=200. However, the phase trajectory of the original system converges to a stable equilibrium point when I1=200 (see Figure 7c and 7d). The similar results can be obtained when I2=60, so we omit them here.

    Figure 7.  Time evolutions and phase diagrams of the controlled system (a-b) and original system (c-d) when I1 = 200μA/cm2.

    In this paper, we investigate the Hopf bifurcation control of the ML model with Hc bifurcation type by utilizing the method of dynamic state-feedback control. The results show that (1) the location of the Hopf bifurcation point can be changed by introducing the state-feedback control without changing the model parameters; (2) the control gain k1 transforms the ML model from the Hc bifurcation type to the SNIC bifurcation type, but does not change the firing pattern of the model; (3) the cooperation of control gain k1 and k3 transforms the ML model from the Hc bifurcation type to the Hopf bifurcation type, and the firing pattern shows a change from Type I to Type II. In addition, we find that the criticality of the Hopf bifurcation of the system can be changed by introducing the control gain k3, that is, the criticality of the Hopf bifurcation is changed from the subcritical Hopf bifurcation to the supercritical Hopf bifurcation. Many neurodegenerative diseases are caused by the change of firing patterns. Therefore, the results of this paper may have potential applications for the treatment of these diseases.

    Compared to previous studies, the conclusions of this paper further enrich the research results. Some scholars have studied the Hopf bifurcation control based on the ML model with SNIC bifurcation type by introducing a washout filter-aided dynamic feedback controller[19]. In their research, they only studied how the location and criticality of the Hopf bifurcation changed under different control conditions. However, they cannot perform the change of the Hopf bifurcation point from one to two. Therefore, it is impossible to study how to control the firing pattern of the system changes from type I to type II. Here, we utilize a new dynamic state feedback controller to control the Hopf bifurcation of the ML model with Hc bifurcation type. This controller not only reproduces the previous results, but also realizes the control of two Hopf bifurcation points. We further study the effects of the controller on firing patterns, providing potential application value for the treatment of some dynamic diseases. In addition, the stability index to determine the criticality of the Hopf bifurcation is given in our paper. The mutual verification of theoretical analysis and numerical simulation improves the integrity of this research.

    This work was supported by the national Natural Science Foundation of China under Grant 11972115 and the Fundamental Research Funds for the Central Universities. The authors would like to thank the reviewers for their comments.

    The authors declare there is no conflicts of interest.



    [1] Z. X. Chen, Complex Differences and Difference Equations, Science Press, 2014.
    [2] M. F. Chen, Z. S. Gao, J. L.Zhang, Entire solutions of certain type of non-linear difference equations, Comput. Methods Funct. Theory, 19 (2019), 17–36.
    [3] Y. M. Chiang, S. J. Feng, On the growth of logarithmic differences, difference quotients and logarithmic derivatives of meromorphic functions, Trans. Amer. Math. Soc., 361 (2009), 3767–3791.
    [4] W. Chen, P. C. Hu, Q. Wang, Entire Solutions of Two Certain Types of Non-linear Differential-Difference Equations, Comput. Methods Funct. Theory., 21 (2021), 199–218.
    [5] J. B. Conway, Functions of One Complex Variable, I, 2Eds., springer.com: Springer, 1995.
    [6] G. G. Gundersen, Questions on meromorphic functions and complex differential equations, Comput. Methods Funct. Theory, 17 (2017), 195–209.
    [7] R. G. Halburd, R. J. Korhonen, Difference analogue of the lemma on the logarithmic derivative with applications to difference equations, J. Math. Anal. Appl., 314 (2006), 477–487.
    [8] R. G. Halburd, R. J. Korhonen, Nevanlinna theory for the difference operator, Ann. Acad. Sci. Fenn. Math., 31 (2006), 463–478.
    [9] R. G. Halburd, R. J. Korhonen, K. Tohge, Holomorphic curves with shift-invariant hyperplane preimages, Trans. Amer.Math. Soc., 366 (2014), 4267–4298.
    [10] W. K. Hayman, Meromorphic Functions, Oxford Mathematical Monographs, Clarendon Press, Oxford, 1964.
    [11] R. Korhonen, An extension of Picards theorem for meromorphic functions of small hyper-order, J. Math. Anal. Appl., 357 (2009), 244–253.
    [12] I. Laine, Nevanlinna theory and complex differential equations, W. de Gruyter, Berlin, 1993.
    [13] P. Li, C. C. Yang, On the nonexistence of entire solutions of certain type of nonlinear differential equations, J. Math. Anal. Appl., 320 (2006), 827–835.
    [14] L. W. Liao, Z. Ye, On solutions to nonhomogeneous algebraic differential equations and their application, J. Aust. Math. Soc., 97 (2014), 391–403.
    [15] J. R. Long, J. Heittokangas, Z. Ye, On the relationship between the lower order of coefficients and the growth of solutions of differential equations, J. Math. Anal. Appl., 444 (2016), 153–166.
    [16] E. Mues, N. Steinmetz, The theorem of Tumura-Clunie for meromorphic functions, J. Lond. Math.Soc., 23 (1981), 113–122.
    [17] N. Steinmetz, Zur Wertverteilung von Exponentialpolynomen, Manuscr. Math., 26 (1978/79), 155–167.
    [18] Z. T. Wen, J. Heittokangas, I. Laine, Exponential polynomials as solutions of certain nonlinear difference equations, Acta Math. Sci. Ser. B., 28 (2012), 1295–1306.
    [19] C. C. Yang, Applications of the Tumura–Clunie theorem, Trans. Am. Math. Soc., 151 (1970), 659–662.
    [20] H. X. Yi, On a theorem of Tumura–Clunie for a differential polynomial, Bull. Lond. Math. Soc., 20 (1988), 593–596.
    [21] C. C. Yang, H. X. Yi, Uniqueness theory of meromorphic functions, Mathematics and its Applications, 557. Kluwer Academic Publishers Group, Dordrecht, 2003.
  • This article has been cited by:

    1. Mi Zhou, Naeem Saleem, Basit Ali, Misha Mohsin, Antonio Francisco Roldán López de Hierro, Common Best Proximity Points and Completeness of ℱ−Metric Spaces, 2023, 11, 2227-7390, 281, 10.3390/math11020281
    2. Basit Ali, Arshad Ali Khan, Manuel De la Sen, Optimum Solutions of Systems of Differential Equations via Best Proximity Points in b-Metric Spaces, 2023, 11, 2227-7390, 574, 10.3390/math11030574
    3. Amer Hassan Albargi, Jamshaid Ahmad, Integral Equations: New Solutions via Generalized Best Proximity Methods, 2024, 13, 2075-1680, 467, 10.3390/axioms13070467
    4. Muzammil Ali, Basit Ali, Best Proximity Point Results for Fuzzy Proximal Quasi Contractions with Applications, 2024, 12, 2227-7390, 2295, 10.3390/math12142295
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2787) PDF downloads(109) Cited by(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog