Research article Special Issues

A two-state neuronal model with alternating exponential excitation

  • Received: 13 February 2019 Accepted: 04 April 2019 Published: 18 April 2019
  • We develop a stochastic neural model based on point excitatory inputs. The nerve cell depolarisation is determined by a two-state point process corresponding the two states of the cell. The model presumes state-dependent excitatory stimuli amplitudes and decay rates of membrane potential. The state switches at each stimulus time. We analyse the neural firing time distribution and the mean firing time. The limit of the firing time at a definitive scaling condition is also obtained. The results are based on an analysis of the first crossing time of the depolarisation process through the firing threshold. The Laplace transform technique is widely used.

    Citation: Nikita Ratanov. A two-state neuronal model with alternating exponential excitation[J]. Mathematical Biosciences and Engineering, 2019, 16(5): 3411-3434. doi: 10.3934/mbe.2019171

    Related Papers:

    [1] Shengyu Huang, Hengguo Yu, Chuanjun Dai, Zengling Ma, Qi Wang, Min Zhao . Dynamics of a harvested cyanobacteria-fish model with modified Holling type Ⅳ functional response. Mathematical Biosciences and Engineering, 2023, 20(7): 12599-12624. doi: 10.3934/mbe.2023561
    [2] Moitri Sen, Malay Banerjee, Yasuhiro Takeuchi . Influence of Allee effect in prey populations on the dynamics of two-prey-one-predator model. Mathematical Biosciences and Engineering, 2018, 15(4): 883-904. doi: 10.3934/mbe.2018040
    [3] Shengyu Huang, Hengguo Yu, Chuanjun Dai, Zengling Ma, Qi Wang, Min Zhao . Dynamic analysis of a modified algae and fish model with aggregation and Allee effect. Mathematical Biosciences and Engineering, 2022, 19(4): 3673-3700. doi: 10.3934/mbe.2022169
    [4] Md Nazmul Hassan, Kelsey Thompson, Gregory Mayer, Angela Peace . Effect of Excess Food Nutrient on Producer-Grazer Model under Stoichiometric and Toxicological Constraints. Mathematical Biosciences and Engineering, 2019, 16(1): 150-167. doi: 10.3934/mbe.2019008
    [5] Zhenliang Zhu, Yuming Chen, Zhong Li, Fengde Chen . Dynamic behaviors of a Leslie-Gower model with strong Allee effect and fear effect in prey. Mathematical Biosciences and Engineering, 2023, 20(6): 10977-10999. doi: 10.3934/mbe.2023486
    [6] Santanu Bhattacharya, Nandadulal Bairagi . Dynamic optimization of fishing tax and tourism fees for sustainable bioeconomic resource management. Mathematical Biosciences and Engineering, 2025, 22(7): 1751-1789. doi: 10.3934/mbe.2025064
    [7] Wenjie Yang, Qianqian Zheng, Jianwei Shen, Linan Guan . Bifurcation and pattern dynamics in the nutrient-plankton network. Mathematical Biosciences and Engineering, 2023, 20(12): 21337-21358. doi: 10.3934/mbe.2023944
    [8] Huanyi Liu, Hengguo Yu, Chuanjun Dai, Zengling Ma, Qi Wang, Min Zhao . Dynamical analysis of an aquatic amensalism model with non-selective harvesting and Allee effect. Mathematical Biosciences and Engineering, 2021, 18(6): 8857-8882. doi: 10.3934/mbe.2021437
    [9] Eric M. Takyi, Charles Ohanian, Margaret Cathcart, Nihal Kumar . Dynamical analysis of a predator-prey system with prey vigilance and hunting cooperation in predators. Mathematical Biosciences and Engineering, 2024, 21(2): 2768-2786. doi: 10.3934/mbe.2024123
    [10] Yong Yao . Dynamics of a delay turbidostat system with contois growth rate. Mathematical Biosciences and Engineering, 2019, 16(1): 56-77. doi: 10.3934/mbe.2019003
  • We develop a stochastic neural model based on point excitatory inputs. The nerve cell depolarisation is determined by a two-state point process corresponding the two states of the cell. The model presumes state-dependent excitatory stimuli amplitudes and decay rates of membrane potential. The state switches at each stimulus time. We analyse the neural firing time distribution and the mean firing time. The limit of the firing time at a definitive scaling condition is also obtained. The results are based on an analysis of the first crossing time of the depolarisation process through the firing threshold. The Laplace transform technique is widely used.


    In eutrophic lakes, cyanobacterial Microcystis colonies usually aggregate and float upwards to form nuisance mucilaginous Microcystis blooms, the development of mucilaginous cyanobacterial Microcystis blooms is a serious environmental and ecological problem, and information on the bloom-formation mechanism has been lacking until now [1]. Based on their long-term research work on cyanobacteria bloom in Taihu Lake, Qin et al. [2] believed that the outbreak of cyanobacteria bloom actually experienced four stages: 1) cyanobacteria cell proliferation stage; 2) cyanobacteria cell groups formation stage; 3) cyanobacteria cell groups floating stage; 4) cyanobacteria bloom outbreak stage. At the same time, they also deemed that in the process of cyanobacteria cell aggregation, the collision probability of cell groups could increase, and smaller cell masses were easier to bond, which could form larger groups and accelerate the emergence of cyanobacteria bloom. Furthermore, Chen et al. [3] pointed out that the aggregation and migration of Microcystis had many ecological advantages, which could promote the formation of Microcystis bloom. Kong et al. [4] proposed that the most direct reason for the formation of cyanobacteria bloom was the massive aggregation, growth, migration and floating to the water surface stage of Microcystis. Obviously, the formation of cyanobacteria cell aggregation is one of the keys to the outbreak of cyanobacteria bloom. Therefore, how to prevent and disperse cyanobacteria cell groups is one of the key control measures to prevent cyanobacteria bloom.

    Cyanobacterial blooms are commonly referred to as harmful algal blooms (HABs) due to their vast ecological impacts and ability to generate metabolites and taste, a variety of in-lake/reservoir control measures are implemented to reduce the abundance of nuisance cyanobacteria biomass[5]. Kibuye [5,6] gave a critical review on operation and performance of source water control strategies for cyanobacterial bloom: chemical control methods, mechanical/physical control methods and biological control methods. And they pointed out that each control method had site-specific strengths, limitations, and ecological impacts, which should hint that none of the reviewed control strategies can provide a comprehensive solution to cyanobacterial blooms. In view of the cyanobacteria blooms in Taihu Lake, the competent authorities mainly adopt the comprehensive method of combining mechanical control method and biological control method, this is because that Taihu Lake is a drinking water source. Although most mechanical and biological control strategies can offer long-term control, their application can be cost-prohibitive and treatment efficacy is influenced by source water geometry and continual nutrient inputs from external sources[6]. Therefore, how to coordinate mechanical control methods with biological control methods is very important, and it is also one of the problems worthy of our in-depth study.

    In order to curb the scale of cyanobacteria bloom in Taihu Lake, the management department of Taihu Lake should realize the transformation from focusing on source control and pollution interception to paying equal attention to source control, pollution interception and ecological regulation[7]. At the same time, there are many kinds of cyanobacteria in Taihu Lake, usually Microcystis[8], and allowing cyanobacteria cells to aggregate is one of the necessary conditions for cyanobacteria bloom [9]. Furthermore, once the cyanobacteria bloom breaks out, the regeneration from cell turnover and nutrient recycling by closely associated heterotrophic bacteria and microzooplankton grazers can help sustain bloom biomass, so the key biological factors to control cyanobacteria bloom must include grazing of zooplankton (possibly benthos and fish)[10]. However, some literatures have studied the dynamic relationship between cyanobacteria cell aggregation and zooplankton/fish grazing. Yang et al.[11] investigated the effect of grazing by different sorts of zooplankton on the induction of defensive morphology in the cyanobacterium Microcystis aeruginosa, and pointed out that cyanobacterium Microcystis aggregation could effectively block zooplankton predation and thus increase the survival of Microcystis aeruginosa. Burkert et al.[12] inquired into the interactions between the cyanobacteria and the potential grazer, pointed out that the predator came in direct contact with Microcystis, aggregates were formed. Chow-Fraser et al.[13] had evidence that grazers in oligotrophic lakes exert a greater impact on algae than those in eutrophic lakes. To summarise, cyanobacteria species can also benefit by their tendencies to congregate as large filamentous and colonial colonies, which reduces zooplankton predation.

    At present, the cyanobacteria bloom research team of Wenzhou University is trying to use a novel comprehensive method to deal with the problem of cyanobacteria bloom in Taihu Lake, which is a combination of mechanical control method and biological control method. Firstly, in the initial stage of cyanobacteria bloom, a part of cyanobacteria agglomerates were broken by physical jet treatment device (a cyanobacteria bloom treatment ship is developed by Wenzhou University), its purpose is to destroy the aggregation stage of cyanobacteria cells and improve the chance of individual cyanobacteria cells being grazed. Secondly, some streamline ecosystems and food web structures will be reconstructed to create unfavorable conditions for cyanobacteria, some filter-feeding organisms (such as zooplankton, mussels, bivalves and fish) can be added to feed on cyanobacteria [14,15,16]. Finally, with the help of the internal operation characteristics of aquatic ecosystem, we expect that aquatic ecosystem can quickly change from algae bloom ecosystem to algae-zooplankton-fish cycle ecosystem. However, whether or not cyanobacterial bloom can be successfully controlled, the novel comprehensive method can influence changes in competition, predation, and behavioral response to predation and competition in the restored habitat. Therefore, after the cyanobacteria bloom is treated by physical jet technology, the restoration and increase of zooplankton and filter-feeding fish how to affect the self purification ability and balance of aquatic ecosystem are worthy of our further study, whether algal population can coexist with zooplankton or filter-feeding fish in the range of relatively low algal population is also worth exploring.

    With the rapid development of numerical analysis and simulation technology, dynamic models have gradually become one of the powerful tools for studying natural science phenomena and problems, which can promote the relevant scientific problems to obtain some excellent research results, such as relevant literatures [17,18,19,20,21,22,23,24,25,26,27]. In references [17,18,19] applied reaction-diffusion equation to examine the calcium diffusion in the cells, investigated some nonlinear dynamics problems and got some interesting results. In reference [20] explored some dynamical behaviors of forced KdV equation to describe the free surface critical flow over a hole, and some meaningful results were given. In reference [21] inquired into chaotic dynamics of a fractional order HIV-1 model involving AIDS-related cancer cells, some meaningful results showed that order of the fractional derivative had a significant effect on the dynamics process. The paper [22] studied the dynamics of infectious diseases within a human host and in the population, and gave some meaningful theoretical and numerical simulation results. In reference [23] gave some good results to verify the effectiveness of linear feedback control for chaos. In reference [24] investigated a newly developed model for Hepatitis-B infection in sense of the Atangana-Baleanu Caputo fractional-order derivative, and obtained some good research conclusions. In references [25,26,27] looked into multiple bifurcations of discrete-time dynamic model, and some important theoretical and numerical results are obtained. In a word, the research on dynamic models has been developed rapidly, and a large number of excellent research results have emerged.

    Now, the disaster of aquatic ecosystem caused by cyanobacteria bloom is one of the hot issues in the world. After the cyanobacteria bloom is treated by chemical control methods, mechanical control methods and biological control methods, the essential characteristics of the internal self recovery and strengthening mechanism of aquatic ecosystem is one of the problems that need to be studied urgently. Furthermore, the algae bloom research team of Wenzhou University has built a physical spraying treatment ship, which mainly deals with the algae blooms in subtropical lakes and reservoirs, and reconstructs the biological food chain for ecological algae control after treatment. Thus, in order to understand the internal self recovery and strengthening mechanism of aquatic ecosystem, we must deeply explore the interactions between cyanobacteria (Microcystis aeruginosa) and potential grazer after cyanobacteria bloom treatment, and make clear the coexistence mode and change trend of cyanobacteria (Microcystis aeruginosa) and potential grazers based on water ecological model and dynamic simulation. Therefore, we believe that the originality of this paper is to illustrate the dynamic evolution relationship between cyanobacteria and potential grazers by means of the bifurcation dynamic evolution process of the model (2.1), and give the evolutionary characteristics of the coexistence mode of cyanobacteria and potential grazers.

    Bertalanffy [28] proposed a method to study the biological system using mathematical model in 1932, this approach is a combination of coordination, order and purpose, which can form three basic ideas of studying the biological system, namely: trophic-level analysis, system perspective and dynamic view. Since then, the research on ecological model and its related dynamics has developed unprecedentedly, and a large number of literature and works have appeared, such as fundamentals of ecological modelling [29], a nitrogen-phosphorus-blue-green algae model [30], some nitrogen and phosphorus dose-response models [31], a mathematical model including three types of zooplankton and two types of fish as well as two types of algae and nutrients [32], an aquatic bacteria-algae amensalism model [33], an ecological model considering two types of phytoplankton, three types of zooplankton, planktivorous fish, detritus and dissolved matters [34], an aquatic algae-fish model [35], one dimensional hydrodynamic-ecological model [36], a two-component model with nutrients and cyanobacteria [37], a four-component model including nutrient, unicellular algae, colonial algae and herbivorous zooplankton [38], an algal bloom mathematical model involving zooplankton-phytoplankton population [39]. Obviously, the ecological model has been developed rapidly in water eutrophication and algal bloom, and some meaningful theoretical results have been obtained. Although there is a certain understanding of the mechanism and control effect of cyanobacteria bloom in recent years, there are still many unsolved mysteries in the outbreak and control of cyanobacteria bloom.

    In order to make clear the coexistence mode and change trend of cyanobacteria (Microcystis aeruginosa) and potential grazers (mussels, bivalves, and silver carp) after cyanobacteria bloom is treated by comprehensive control strategy, which can help us explore the interactions between cyanobacteria and potential grazers, and explicit the effectiveness and feasibility of the comprehensive control strategy, we will construct a new aquatic ecological model based on the comprehensive control strategy and internal operation law of aquatic ecosystem. Therefore, some model assumptions and cyanobacteria bloom background will be given as follows:

    1) The growth environment of cyanobacteria and potential grazers is assumed to be a semi closed aquatic ecosystem, species biomass is always evenly distributed in space and changed instantaneously with time t. x(t) and y(t) represent respectively the biomass of total cyanobacteria (unicellular cyanobacteria and colonial cyanobacteria) and potential grazers. After cyanobacteria bloom is treated by mechanical control method (physical jet processing device), total cyanobacteria will be divided into colonial cyanobacteria α1x(t) and unicellular cyanobacteria α2x(t) with α1+α2=1. If α1=0 and α2=1, this shows that cyanobacteria population mainly exists in the form of unicellular cyanobacteria, and the cyanobacteria population has no aggregation behavior and can not form a group aggregation state. If α1=1 and α2=0, this indicates that cyanobacteria bloom has erupted in a large area, especially serious, and all cyanobacteria populations are clustered, which is an extreme state of existence and almost does not exist in reality.

    2) Because when unicellular cyanobacteria gather into groups, the light utilization rate of algae cells can be enhanced [40], the salt stress resistance and high light stress of algal cells can be improved [41,42], the ability to avoid the invasion of viruses, bacteria and algae phages can be strengthened [43], nutrient and inorganic carbon absorption capacity can be greatly promoted[44,45], the risk of predation can be greatly reduced [46]. In light of the above, cyanobacteria aggregation can find more suitable niche, obtain greater ecological advantages and enhance growth rate. Therefore, we assume that the growth kinetics function of cyanobacteria x(t) is f(α1α2)x(t)(1x(t)k) with intrinsic growth rate function f(α1α2) and maximum environmental capacity k,

    f(α1α2)={r1α1α2,0<α1<1,0<α2<1,rmin,α1=0,α2=1,rmax,α1=1,α2=0,

    where r1 is a median growth rate with α1=12 and α2=12, rmin is a growth rate without aggregation, and, rmax is a growth rate with overall aggregation state.

    3) During cyanobacteria blooms, cyanobacteria can induce morphological changes to colonial forms in order to avoid zooplankton grazing, this mechanism can be an adaptive reaction for survival that has developed through evolutionary processes[47,48]. For example, Microcystis exists as colonial forms, the relevant research results are in the papers[35,38,43]. However, unicellular cyanobacteria do not have this protective mechanism. Hence, we use Holling type II functional responses a2α2x(t)y(t)b+α2x(t) to describe the predation mechanism of potential grazers on unicellular cyanobacteria, this is because the predation process of potential grazers not only take time to capture, but also depend on the abundance of unicellular cyanobacteria, where a2 is a grazing coefficient and b is a half-saturation constant. At the same time, we use mathematical functional a1(α1x(t)g(α1α2))y(t) to describe the predation mechanism of potential grazers on colonial cyanobacteria, a1 is a grazing coefficient and predation function g(α1α2) avoided by colonial cyanobacteria, and

    g(α1α2)={mα1α2,0<α1<1,0<α2<1,0,α1=0,α2=1,mmax,α1=1,α2=0,

    where m is a median biomass of colonial cyanobacteria (which can never be hunted) with α1=12 and α2=12, 0 indicates that unicellular cyanobacteria do not have this ability, mmax is a maximum biomass of colonial cyanobacteria (which can never be hunted) with overall aggregation state.

    4) It is widely known that potential grazers can not eat only cyanobacteria, and must forage other species as alternate foods. Thus, the growth kinetics function of potential grazers is r2y(t)+e1a1(α1x(t)g(α1α2))y(t)+e2a2α2x(t)y(t)b+α2x(t) with alternate coefficient r2 and absorption conversion coefficients e1 and e2. Furthermore, cyanobacteria and potential grazers have their own other losses (such as natural death), we might as well set them to m1x(t) and m2y(t) with loss coefficient m1,m2.

    Based on the above assumptions and analysis, we can construct a new aquatic ecological control model and two special aquatic ecological models with non-negative initial conditions x(0)x00 and y(0)y00, which can be described as following:

    {dx(t)dt=r1α1α2x(1xk)a1(α1xα1α2m)ya2α2xyb+α2xm1x,dy(t)dt=r2y+e1a1(α1xα1α2m)y+e2a2α2xyb+α2xm2y, (2.1)
    {dx(t)dt=rminx(1xk)a2xyb+xm1x,dy(t)dt=r2y+e2a2xyb+xm2y, (2.2)
    {dx(t)dt=rmaxx(1xk)a1(xmmax)ym1x,dy(t)dt=r2y+e1a1(xmmax)ym2y. (2.3)

    The model (2.1) represents the mixed state of unicellular cyanobacteria and colonial cyanobacteria, the biomass of unicellular cyanobacteria and colonial cyanobacteria is controlled by critical parameter α1 or α2 because of α1+α2=1. If the value of α1 belongs to (0,0.2], the model (2.1) represents cyanobacteria cell proliferation stage. If the value of α1 belongs to (0.2,0.4], the model (2.1) represents cyanobacteria cell groups formation stage. If the value of α1 belongs to (0.4,0.6], the model (2.1) represents cyanobacteria cell groups floating stage. If the value of α1 belongs to (0.6,1), the model (2.1) represents cyanobacteria bloom outbreak stage. Overall, the model (2.1) can describe that the outbreak of cyanobacteria bloom in Taihu Lake actually experienced four stages [2]. The models (2.2) and (2.3) represent the stage of unicellular cyanobacteria and colonial cyanobacteria respectively, neither of these two states will exist in reality, therefore, we only study them as two limit states.

    With the help of bifurcation dynamic analysis and dynamic simulation of the model (2.1), the main purpose of this paper is to explore how the aquatic ecosystem can recover and strengthen itself after the comprehensive treatment of cyanobacteria bloom. To solve this problem, we firstly make a qualitative analysis of the model (2.1), deduce some critical threshold conditions to ensure the existence and stability of some equilibrium points, and give some critical conditions under which the model can induce Transcritical and Hopf bifurcation. Secondly, we will comprehensively investigate how comprehensive management measures affect the dynamic characteristics of cyanobacteria bloom and evaluate the effectiveness of comprehensive management measures in controlling cyanobacteria bloom. Finally, we carry out relevant dynamic experiments to verify the feasibility of the theoretical derivation results and identify the possible coexistence mode of cyanobacteria and potential grazers, then we can analyze interaction mechanism between cyanobacteria and potential grazers in the process of cyanobacteria bloom outbreak and control.

    In order to investigate the superiority of the model (2.1), the dynamic relationship between the biomass of cyanobacteria and potential grazers was qualitatively analyzed and numerically simulated with r1=0.56, r2=0.2, k=5, a1=0.6, a2=0.4, b=0.15, m1=0.18, m2=0.28, e1=0.6, e2=0.28, α1=0.35, α2=0.65, rmin=0.2, rmax=0.8, mmax=3, m=0.1.

    From the model (2.2), we can find that the relation expression of cyanobacteria x and potential grazers y is given by

    y=(rmin  (1xk)m1)(b+x)a2.

    It is easy to see from Figure 1(a) that the dynamic relationship of cyanobacteria x and potential grazers y is a concave function as the value of cyanobacteria x increases, and there is a unique maximum value, and the maximum value of cyanobacteria population is still low. Unfortunately, the value range of cyanobacteria x is relatively small, and the biomass of potential grazers y becomes negative with the increase of cyanobacteria x. Obviously, this dynamic result can not better reflect the running law of cyanobacteria bloom control test in the laboratory.

    Figure 1.  The dynamic evolution relationship between cyanobacteria and potential grazers.

    From the model (2.3), we can find that the relation expression of cyanobacteria x and potential grazers y is given by

    y=a1x+m1xrmax   (1xk)a1mmax.

    It can clearly find from Figure 1(b) that the dynamic relationship of cyanobacteria x and potential grazers y is a convex function as the value of cyanobacteria x increases, and there is a unique minimum value 0, which is not in line with the wild growth law of cyanobacteria population because cyanobacteria population will never die out. Furthermore, it is even more regrettable that the value of potential grazers y approaches positive infinity with the increase of cyanobacteria x value, which completely violates the dynamic relationship between cyanobacteria x and potential grazers y in the process of cyanobacteria bloom.

    From the model (2.1), we can find that the relation expression of cyanobacteria x and potential grazers y is given by

    y=r1α1xα2m1xa1(α1xα1mα2)+a2α2xb+α2x.

    It is worth pointing out that potential grazers y is monotonically increasing and has supremum r1a1α2m1a1α1 as the value of cyanobacteria x increases to positive infinity. Furthermore, the better result is that cyanobacteria x approaches a small fixed value when potential grazers y approaches positive infinity, this result is more in line with the natural law (see Figure 1(c)). Moreover, what is more gratifying is that there is a unique minimum value, through which we can regulate the recyclability of cyanobacteria bloom test, which is also our favorite routine manipulation.

    In a word, through the comparative analysis of numerical simulation results, it can be seen that the model (2.1) can more accurately describe the dynamic relationship between cyanobacteria and potential grazers under the mechanical and ecological integrated control, so we mainly study the relevant dynamic problems of the model (2.1).

    The existence and local stability of all possible equilibrium points will be investigated in detail. The equilibrium point is a kind of special solution of the model (2.1), which can determine some coexistence modes of cyanobacteria and potential grazers, and it is also the basis of subsequent bifurcation dynamics research. Thus, the existence and local stability of all possible equilibrium points will be futher investigated in this section.

    Now, we consider the algae isocline vertically and potential grazers isocline horizontally, respectively:

    { r1α1α2x(1xk)a1α1(xmα2)ya2α2xyb+α2xm1x=0, r2+e1a1α1(xmα2)+e2a2α2xb+α2xm2=0&y=0. (2.4)

    By solving the above equations with zero dash, there are two classes of equilibrium points: boundary equilibrium points E0 and E1, the coordinates of internal equilibria E1 and E2.

    It is obvious that the equilibrium points are the intersections of these nullclines. Thus, we easily see that the model (2.1) possesses two boundary equilibrium points E0(0,0) and E1(x1,0). For the possible positive equilibria, we only need consider the positive solutions of the following equations:

    { a1e1α1α2x2+(α2r2α2m2+a2e2α2+a1be1α1a1e1α1m)x+b(r2m2a1e1α1bmα2)=0, r1α1α2x(1xk)a1α1(xmα2)ya2α2xyb+α2xm1x=0. (2.5)

    For positive equilibrium points, the algae population x must satisfy mα2<x<k, Let Δ(x) denote the discriminant of the first equation of (2.5) and express Δ(x) in terms of, i.e,

    Δ(x)=B24AC,

    where, B=(r2m2+a2e2)α2+a1e1α1ba1e1α1m,

    A=a1e1α1α2, C=b(r2m2a1e1α1bmα2).

    Hence, we have

    Δ=((r2m2+a2e2)α2+a1e1α1ba1e1α2m)24a1e1α1α2b(r2m2+a1e1m).

    With respect to the conditions of the equilibrium points, we obtain the following theorem:

    Theorem 1. The equilibrium points of the model (2.1) are as follows: the model(2.1) always has a boundary equilibrium point E0(0,0). If α1α2>m1r1, then the model (2.1) has a boundary equilibrium point E1(x1,0) with x1=k(r1α1m1α2)r1α1.

    Theorem 2. For the positive equilibrium point, we have:

    1) When mI1, the model (2.1) has no positive equilibrium point E1(x1,y1),

    2) When r2m2+a2e2>0 and r2m2<0, if max{0,I1}<m<min{I2,I3,I4} or max{0,I1,I2}<m<min{I3,I4}, then the model (2, 1) has a unique positive equilibrium point E1(x1,y1),

    3) When r2m2+a2e2<0, if max{0,I1,I3}<m<min{I2,I4} or max{0,I1,I2,I3}<m<I4, then the model (2.1) has a unique positive equilibrium point E1(x1,y1),

    where

    x1=((r2m2+a2e2)α2+a1e1α1ba1e1α2m)+Δ2a1e1α1α2,y1=(r1α1α2m1)x1r1α1kα2x21a1(α1x1α1mα2)+a2α2x1b+α2x1,
    I1=(r2m2)α2a1e1α1,I2=α2(r2m2+a2e2)α1a1e1,I3=b(m2r2)r2m2+a2e2,I4=(kα22+bα2)(r2m2+a1e1α22ka1e1α1(α2k+b).

    Proof. Let the solution of the nullclines g(x)=Ax2+Bx+C=0, we know that the model (2.1) has the positive roots only if r2m2a1e1α1mα2<0, which means mI1. At the same time, we also show that if mI1, then Δ>0 and C<0, which means that g(x) has a unique positive root E1(x1,y1). Hence, if g(mα2)<0 and g(k)>0, there is a positive real root in mα2<x<k. Due to the fact that Δ>0, hence, if g(mα2)<0 and g(k)>0, we have the conditions with I3 and I4. If g(mα2)<0, it follows that r2m2+a2e2>0 for m<I3 or r2m2+a2e2<0 for m>I3. If m<I4, we have f(k)>0. Next, we just discuss the symbols of B.

    If m=I2, we have B=0, then g(x)=0 has a unique positive root x1=b(α2m2α2r2+a1e1α1m)a1e1α1α22. Similarly, if m<I2 or m>I2, then B>0 or B<0, we must know that g(x)=0 has a unique positive root x. From the above analysis, the model (2.1) has a unique positive equilibrium point if and only if one of the following situation:

    1) When m=I2, the model (2.1) has a unique equilibrium point E1(x1,y1).

    2) When r2m2+a2e2>0 and r2m2<0, if max{0,I1}<m<min{I2,I3,I4} or max{0,I1,I2}<m<min{I3,I4}, then the model (2.1) has a unique positive equilibrium point E1(x1,y1).

    3) When r2m2+a2e2<0, if max{0,I1,I3}<m<min{I2,I4} or max{0,I1,I2,I3}<m<I4, then the model (2.1) has a unique positive equilibrium point E1(x1,y1).

    Next, we will study the types and stability of the equilibrium points. Firstly, the Jacobian matrix of the model (2.1) is given as follows:

    J(x,y)=(r1α1α2(12xk)a1α1ya2α2by(b+α2x)2m1a1α1(xmα2)a2α2xb+α2xe1a1α1y+a2e2α2by(b+α2x)2r2m2+e1a1α1(xmα2)+e2a2α2xb+α2x).

    Theorem 3. With respect to the origin of this equilibrium point E0, we have:

    1) When α1α2<m1r1 and α1α2<m2r2e1a1m hold, the equilibrium point E0 is a stable node.

    2) When α1α2>m1r1 and α1α2<m2r2e1a1m hold, or α1α2<m1r1 and α1α2>m2r2e1a1m hold, the equilibrium point E0 is a saddle.

    Proof. The Jacobian matrix of the equilibrium point E0 is

    JE0=(r1α1α2m1α1a1mα20r2m2+α1e1a1mα2).

    Obviously, JE0 has two characteristic roots λ1=r1α1α2m1, λ2=r2m2+α1e1a1mα2. We applied the Routh-Hurwitz criterion to analyze the stability. The equilibrium point E0 is a stable node since λ1<0 and λ2<0, whiler1α1α2<m1 and r2+e1a1α1mα2<m2 hold. The equilibrium point E0 is a saddle since λ1<0 and λ2<0, while r1α1α2<m1 and r2+e1a1α1mα2<m2 hold. Similarly it can be proved that the equilibrium point E0 is a saddle if r1α1α2<m1 and r2+e1a1α1mα2>m2.

    Theorem 4. When α1α2>m1r1, then the equilibrium point E1(x1,0) is a stable node if r2m2α1e1a1mα2+e1η1+e2η2<0, and the equilibrium point E1(x1,0) is a saddle if r2m2α1e1a1mα2+e1η1+e2η2>0.

    Proof. The Jacobian matrix of E1 is

    JE1=(m1α1r1α2α1a1mα2η1η20r2m2α1e1a1mα2+e1η1+e2η2),

    where

    η1=a1k(α1r1α2m1)r1,η2=a2kα2(α1r1α2m1)α1r1bα2k(α1r1α2m1).

    Then the eigenvalues of JE1 are λ3=m1α1r1α2, λ4=r2m2α1e1a1mα2+e1η1+e2η2. It is obvious that λ3<0, hence the equilibrium point E1(x1,0) is a stable node if λ4<0, and the equilibrium point E1 is a saddle if λ4>0.

    Theorem 5. The conditions for the existence of the internal equilibrium point E1(x1,y1) are shown in Theorem 2. In this paper, it has been proved that there is and only one internal equilibrium point E1(x1,y1). If Tr(JE1(x1,y1))<0, then the internal equilibrium point E1(x1,y1) is locally asymptotically stable; otherwise, the internal equilibrium point E1(x1,y1) is unstable.

    Proof. The Jacobian matrix of the internal equilibrium point E1 is

    JE1=(r1α1α2(12x1k)a1α1y1a2α2by1(b+α2x1)2m1a1α1(x1mα2)a2α2x1b+α2x1e1a1α1y1+a2e2α2by1(b+α2x1)20).

    The determinant of the matrix JE1 is

    Det(JE1)=(e1a1α1y1+a2e2α2by1(b+α2x1)2)(a1α1(x1mα2)+a2α2x1b+α2x1).

    The trace of the matrix JE1 is

    Tr(JE1)=r1α1α2(12x1k)a1α1y1a2α2by1(b+α2x1)2m1.

    It is obvious that all terms of the product of determinants are positive, hence, we can judge that the determinants of the Jacobian matrix of E1 is always positive. Thus, when Tr(JE1)<0, the internal equilibrium point E1 has two negative eigenvalues, which means that the internal equilibrium point E1 is locally asymptotically stable. Otherwise, both characteristic roots are positive, and the internal equilibrium point E1 is unstable.

    In order to probe some possible coexistence evolution modes of cyanobacteria and potential grazers during the outbreak of cyanobacteria bloom, and clarify how the sheltering effect of colonial cyanobacteria affects the coexistence mode of cyanobacteria and potential grazers, we choose m as the control parameter to study the bifurcation dynamic behavior of the model (2.1), mainly including Transcritical bifurcation and Hopf bifurcation. If the model (2.1) can have Transcritical bifurcation, which implies that the cyanobacteria and potential grazers can coexist in a certain range at least, and there will be no extinction of at least one population with the advance of time. If the model (2.1) can have Hopf bifurcation, which shows that cyanobacteria and potential grazers can coexist in a periodic oscillation mode in a certain range, and the two populations can continue to survive with the advance of time. Therefore, bifurcation dynamics of the model (2.1) has important significance in population ecology and is worthy of our in-depth study.

    The existence of Transcritical bifurcation is studied at the boundary equilibrium point E1(x1,0), we should ensure that α1α2>m1r1 is satisfied. Thus, we state the following theorem.

    Theorem 6. The model (2.1) can undergo a Transcritical bifurcation at the boundary equilibrium poin E1(x1,0) when mmTC=α2α1e1a1(r2m2+e1η1+e2η2).

    Proof. Now, we use Sotomayor's Theorem to verify the transversality condition for the occurrence of Transcritical bifurcation at mmTC=α2α1e1a1(r2m2+e1η1+e2η2), it can be easily verified that Det(JE1)=0. Since λ3λ4=Det(JE1), which implies that JE1TC has a zero eigenvalue, that is to say λ4=0, thus the Jacobian matrix at E1 is

    JE1TC=(m1α1r1α2α1a1mα2η1η200).

    Let v and w be the two eigenvectors corresponding to the zero eigenvalue of the matrices JE1TC and JTE1TC respectively, and they are given by

    v=(v11),w=(01),

    where v1=α2(r2m2e1η2+e2η2)e1(α1r1α2m1).

    Furthermore, we can get

    Fm(E1,mTC)=(α1a1yα2α1e1a1yα2)(E1,mTC)=(00),
    DFm(E1,mTC)v=(0α1a1α20α1e1a1α2)(v11)(E1,mTC)=(α1a1α2α1e1a1α2),
    D2Fm(E1,mTC)(v,v)=(2F1x2v21+22F1xyv1v2+2F1y2v222F2x2v21+22F2xyv1v2+2F2y2v22)(E1,mTC)
    =(0α2(r2m2e1η2+e2η2)e1(α1(r1+m1)m1)(a1e1α1+a2e2bα2(b+α2(α1r1α2m1)r1α1)2)),

    clearly, they can satisfy the transversality conditions:

    wTFm(E1,mTC)=0,

    wT[DFm(E1,mTC)v]=α1e1a1α20,

    wT[D2Fm(E1,mTC)(v,v)]=2α2(r2m2e1η2+e2η2)e1(α1(r1+m1)m1)(a1e1α1+a2e2bα2(b+α2(α1r1α2m1)r1α1)2)0.

    Therefore, when mmTC=α2α1e1a1(r2m2+e1η1+e2η2), the model (2.1) can induce Transcritical bifurcation at the boundary equilibrium point E1.

    In the view of the proof of Theorem 5, we can know that Hopf bifurcation can occur at the positive equilibrium point E1. At the same time, the positive equilibrium point E1 can change its stability when the value of m passes through the critical magnitude. Thus, we summarize our findings in the following theorem.

    Theorem 7. Based on Theorem 5, the model (2.1) can induce a Hopf bifurcation when m=mHP.

    Proof. The characteristic equation of matrix JE1 is λ2Tr(JE1)λ+Det(JE1)=0, and three cross-sectional conditions of Hopf bifurcation occurs, i.e.,

    1) Tr(JE1)m=mHP=0,

    2) Det(JE1)m=mHP>0,

    3) dTr(JE1)dαm=mHP=0.

    Previously, we can show that Det(JE1)m=mHP>0. And when m=mHP, we have Tr(JE1)m=mHP=0. Thus we have dTr(JE1)dαm=mHP=0, where

    Tr(JE1)m=mHP=r1α1α2(12x1k)a1α1y1a2α2by1(b+α2x1)2m1=0.

    In conclusion, the model (2.1) can induce a Hopf bifurcation when m=mHP. Next, we will discuss the stability of the limit cycle by computing the first Lyapunov number. Translate the equilibrium point E1 to the origin by using the following transformation:

    x=xmx1,y=ymy1.

    Hence, we obtain

    {˙xm=a10xm+a01ym+a20x2m+a11xmym+a02y2m+a30x3m+a21x2mym+a12xmy2m+a03y3m+Q(xm,ym),˙ym=b10xm+b01ym+b20x2m+b11xmym+b02y2m+b30x3m+b21x2mym+b12xmy2m+b03y3m+P(xm,ym),

    where

        a10=a1α1y1+a2α2y1bx1α2+r1α1k+2r1α1x1α2k+a2α22x1y1(bx1α2)2m1,    a20=r1α1α2ka2α22y1(bx1α2)2a2α32x1y1(bx1α2)3,        a30=a2α32y1(bx1α2)3+a2α42x1y1(bx1α2)4,    a21=a2α22(bx1α2)2+a2α32x1(bx1α2)3,        a11=a1α1a2α2bx1α2a2α22x1(bx1α2)2,    a01=α1a1mα2+a1x1y1+a2α2x1bx1α2,        b10=ya1e1α1e2a2α2y1bx1α2e2a2α22x1y1(bx1α2)2,    b20=e2a2α22y1(bx1α2)2+e2a2α32x1y1(bx1α2)3,        b30=e2a2α32y1(bx1α2)3e2a2α42x1y1(bx1α2)4,    b21=e2a2α22(bx1α2)2e2a2α32x1(bx1α2)3,        b11=e2a2α2bx1α2+e2a2α22x1(bx1α2)2+a1e1α1,    b01=α1a1e1mα2m2+r2a1e1α1x1e2a2α2x1bx1α2,        a02=a12=a03=b02=b12=b03=0.

    And then Q(xm,ym), P(xm,ym) are power series in (xm,ym) with terms ximyjm, which can satisfy i+j4. Therefore, the first Lyapunov number to determine the stability of limit cycle is given by the formula

    l1=3π2a01Δ3/2{[a10b10(a211+a11b02+a02b11)+a10a01(b211+a20b11+a11b02)+b210(a11a02+2a02b02)2a10b10(b202a20a02)2a10a01(a220b20b02)a201(2a20b20+b11b20)+(a01b102a210)(b11b02a11a20)](a210+a01b10)[3(b10b03a01a30)+2a10(a21+b12)+(b10a12a01b21)]}. 

    It is well known that the limit cycle generated by the Hopf bifurcation is unstable if l1>0 and stable if l1<0. Nevertheless, the expression of the first Lyapunov number is too complex to judge the positive and negative. Therefore, we will give some numerical solutions in Section 6.

    Under the research framework of theoretical derivation, we will not only verify the feasibility and effectiveness of the theoretical results, but also explore how clustering affects the dynamic relationship between cyanobacteria and potential grazers, thus the parameter m and β=α1α2 are selected as control variables for numerical simulation with r1=0.56, r2=0.2, k=5, a1=0.6, a2=0.4, b=0.15, m1=0.18, m2=0.28, e1=0.6, e2=0.28, it should be noted that all parameter values are estimated according to biological significance and theoretical conditions. Furthermore, β=α1α2 represents the proportion of colonial cyanobacteria and unicellular cyanobacteria, which can also indicate the outbreak degree of cyanobacteria bloom. If β=α1α2>1 indicates that the manifestation of cyanobacteria is mainly colonial cyanobacteria, which is not conducive to the capture of potential grazers. If β=α1α2<1 indicates that the manifestation of cyanobacteria is mainly unicellular cyanobacteria, which is conducive to the capture of potential grazers. Moreover, it should be pointed out that the predation behavior of potential grazers on cyanobacteria is essentially the implementation of ecological control technology.

    When we take β=0.5385 with α1=0.35 and α2=0.65, which represents that the cyanobacteria in the whole ecosystem are in the bloom period, showing the phenomenon of aggregation. In other words, the current situation is not conducive to potential grazers to prey on cyanobacteria. Based on numerical simulation work, the bifurcation diagram of the model (2.1) with control parameter m has been shown in Figure 2(a). It is easy to find that the model (2.1) will experience a Hopf bifurcation and Transcritical bifurcation with the value of m changing within [0,1.8]. When the value of m gradually decreases from 1.8 to mTC=1.42282, the model (2.1) will occur a Transcritical bifurcation (Figure 3), the boundary equilibrium point E1 will change from a stable state to an unstable state, which suggests that potential grazers will not always be close to extinction.

    Figure 2.  The evolution process of coexistence mode of cyanobacteria and potential grazers. The red line represents the internal equilibrium point x1 changing by m, solid line represents the equilibrium point stable, dashed line represents the equilibrium point unstable and two vertical dashed lines represent the critical value of bifurcation. More detailed, Hp and TC are the critical values for the Hopf bifurcation, Transcritical bifurcation of the model, respectively. (a) mHP=0.2957018, mTC=1.42282; (b) mHP=0.1412414, mTC=1.74941.
    Figure 3.  Evolution process of steady-state coexistence mode with α1=0.35 and α2=0.65.

    And when the value of m is between (mHP,mTC), the model (2.1) has an asymptotically stable internal equilibrium point E1 (seeing Figure 3(b)), which indicates that cyanobacteria and potential grazers can coexist in a stable equilibrium state. If the value of m gradually decreases to mHP=0.2957018, the model (2.1) will occur Hopf bifurcation, the internal equilibrium point will change from a stable state to an unstable state, and a stable limit cycle will appear, the detailed numerical simulation results are shown in Figures 4 and 5. Furthermore, we can get that the Lyapunov exponent is 2.57010×107, which can further verify that the limit cycle is stable. Moreover, it should be explained that the dynamic characteristics of the model (2.1) have changed substantially, which means that the coexistence state of cyanobacteria and potential grazers changes from a stable equilibrium state to a stable periodic oscillation state.

    Figure 4.  Evolution process of periodic oscillation coexistence mode with α1=0.35 and α2=0.65. (a) We have stable periodic orbits bifurcate through Hopf bifurcation around E1(0.48682,0.16138) with m=mHP=0.2957018; (b) Local amplification of (a); (c) E1(0.50822,0.16493) is a locally asymptotically stable with m=0.315>mHP; (d) E1(0.46966,0.15840) is an unstable focus with m=0.28<mHP.
    Figure 5.  Dynamic evolution diagram of Hopf bifurcation with critical value m=mHP=0.2957018.

    When we take β=1.083 with α1=0.52 and α2=0.48, which represents that the cyanobacteria in the whole ecosystem are still in the bloom period, but the current situation is quite disadvantageous to potential grazers to prey on cyanobacteria because unicellular cyanobacteria account for 0.48. It is easy to discover from Figure 2(b) that the model (2.1) has almost the same dynamic evolution process. When the value of m gradually decreases from 2 to mTC=1.74941, a Transcritical bifurcation will occur (Figure 6). When the value of m continues to decrease to mHP=0.125, the instability of the equilibrium point can produce a stable limit cycle, the detailed dynamic results are shown in Figures 7 and 8.

    Figure 6.  Evolution process of steady-state coexistence mode with α1=0.52 and α2=0.48.
    Figure 7.  Evolution process of periodic oscillation coexistence mode with α1=0.35 and α2=0.65. (a) We have stable periodic orbits bifurcate through Hopf bifurcation around E1(0.38961,0.58723) with m=mHP=0.14124142; (b) Local amplification of (a); (c) E1(0.41055,0.60801) is locally asymptotically stable with m=0.155>mHP; (d) E1(0.36531,0.56250) is unstable focus with m=0.125<mHP.
    Figure 8.  Dynamic evolution diagram of Hopf bifurcation with critical value m=mHP=0.1412414.

    At the same time, the Lyapunov exponent is 5.15158×1014, which can show the stability of the limit cycle. From the essence of system dynamics evolution, whether the value of β is greater than 1 or less than 1, the model (2.1) has exactly the same dynamic evolution characteristics, that is, Transcritical bifurcation and Hopf bifurcation. That is to say, the ratio of colonial cyanobacteria to unicellular cyanobacteria does not affect the evolution characteristics of the dynamic state of the model (2.1). Moreover, it is also worth emphasizing that the value of parameter m seriously affects the dynamic relationship and coexistence mode between cyanobacteria and potential grazers.

    Based on the comparative analysis of numerical simulation results, a Transcritical bifurcation will occur at mTC=1.42282 when the value of β is less than 1, however a Transcritical bifurcation will occur at mTC=1.74941 when the value of β is greater than 1. This result means that during the outbreak of cyanobacteria bloom, the earlier the clustered cyanobacteria population is crushed by physical spray treatment technology, the more opportunities there are to change the dynamic relationship between cyanobacteria and potential grazers, and provide an entry point for the coexistence of cyanobacteria and potential grazers. However, the critical value of inducing Hopf bifurcation is just the opposite. When the value of β is less than 1, Hopf bifurcation can be induced at m=0.2957018, which suggests that potential grazers actively prey on unicellular cyanobacteria to obtain sufficient food, which can lead to a large increase in the biomass of potential grazers. This increased predation ability can lead to essential changes in the dynamic relationship between them, and then play a role in controlling the large-scale outbreak of cyanobacteria bloom. However, when the value of β is greater than 1, Hopf bifurcation can be induced at m=0.14124142, which means that the number of unicellular cyanobacteria can not basically maintain the growth needs of potential grazers, and a large number of colonial cyanobacteria need to be prey, which will cause colonial cyanobacteria to be gradually broken. This result also leads to a stable equilibrium between cyanobacteria and potential grazers in the range of m. Therefore, from the perspective of controlling cyanobacteria bloom by physical spray treatment technology, the colonial cyanobacteria turns into unicellular cyanobacteria after the cyanobacteria bloom is treated by physical spray, which is conducive to potential grazers to prey, and then lead to the coexistence mode of potential grazers and cyanobacteria in a stable state.

    Within the framework of physical and ecological integrated control and management of cyanobacteria bloom, according to the growth characteristics and comprehensive management characteristics of wild cyanobacteria population, we mainly focus on the impact of physical spraying treatment technology on the structural characteristics of cyanobacteria population and the predation effect of potential grazers, and give some ecological modeling assumptions. On this basis, a aquatic ecological dynamic new model is proposed to describe the dynamic relationship between cyanobacteria and potential grazers. At the same time, it is worth emphasizing from Figure 1 that the biggest advantage of this new model is that it integrates three kinds of dynamic action modes, of which the sub model (2.1) is the most distinctive and has not been studied by predecessors. Furthermore, the biggest advantage of the model (2.1) is that it can perfectly describe the dynamic relationship between cyanobacteria and potential grazers under the physical and ecological integrated control, most commendably, it can control cyanobacteria bloom from the perspective of minimum value, and the consumption of potential predators is limited.

    Underlying the analysis of mathematical theory, we first explore the existence and stability of the equilibrium point of the model (2.1), and give some critical conditions of some critical parameters. Then, we obtain some key threshold conditions to ensure that the model (2.1) has specific dynamic behaviors, such as Transcritical bifurcation and Hopf bifurcation. These theoretical results are the basis for our numerical simulation analysis, which provides theoretical support for in-depth exploration of the dynamic relationship between cyanobacteria and potential grazers.

    Based on the results of numerical analysis, the dynamic behavior type and dynamic evolution process of the model (2.1) are clarified, focusing on the Transcritical bifurcation and Hopf bifurcation. It is easy to find from Figure 2 that when cyanobacteria bloom break out and gather into clusters, we implement a physical and ecological comprehensive control strategy for cyanobacterial bloom, the dynamic relationship between cyanobacteria and potential grazers suddenly changes by Transcritical bifurcation, which can form coexistence mode in a stable equilibrium state, the detailed dynamic simulation results are shown in Figures 3 and 6. After then, the biomass of cyanobacteria and potential grazers gradually decrease with the reduction of agglomeration, finally, a stable periodic oscillation coexistence mode is formed to maintain a good cycle of ecological algae control strategy, these dynamic simulation results are shown in Figures 4, 5, 7 and 8. Furthermore, it is more meaningful to obtain from Figures 3 and 6 that Transcritical bifurcation represents evolution process of steady-state coexistence mode, and it is the most ecologically significant result form Figures 4 and 7 that Hopf bifurcation represents evolution process of periodic oscillation coexistence mode. In addition, it should be emphasized that the most important role of physical spraying treatment technology is to break up clumps of cyanobacteria in the process of controlling cyanobacteria bloom, which will not change the dynamic essential characteristics of cyanobacteria and potential grazers represented by the model (2.1).

    Based on theoretical and numerical research results, the main results of this paper are as follows: 1) The population density of cyanobacteria decreases gradually as the value of key parameter m decreases, that is to say, the implementation of physical spraying treatment technology can effectively control the population density of cyanobacteria. 2) Transcritical bifurcation can induce steate-state coexistence mode between cyanobacteria and potential grazers, and Hopf bifurcation can force cyanobacteria and potential grazers to form periodic oscillation coexistence mode. 3) The evolutionary characteristics of bifurcation dynamics suggest that physical and ecological integrated control technique can effectively inhibit the growth of cyanobacteria, and promote the formation of a stable coexistence mode between cyanobacteria and potential grazers within a controllable range.

    In general, some results have been obtained by proposing a new aquatic ecological model, however, the work of this paper still has some limitations, such as: 1) the characteristics of toxicity released by cyanobacteria population were not considered in the modeling, 2) cyanobacteria migration behavior is not mentioned, 3) modeling assumptions make that the application of the model (2.1) has certain limitations. Thus there is still much work to be further studied, for example, how does the release of cyanobacteria toxins affects the predatory ability of potential grazers during physical spray treatment, effect of cyanobacteria migration behavior on algae control effect of physical spraying technology, etc. In a word, it is hoped that the research results of this paper will provide some theoretical support for the application of physical and ecological comprehensive management of cyanobacteria bloom.

    This work was supported by the key International Cooperation Projects of China (Grant No: 2018YFE0103700), the National Natural Science Foundation of China (Grant No: 31570364, 61871293, 41876124, 61901303), by the Zhejiang Provincial Natural Science Foundation of China (Grant No: 2Q18C030002), by key research and development projects of Zhejiang Province (Grant No: 2021C03166), by the Science and Technology program of Cangnan, China (Grant No: 2018ZG29), by Science and Technology Major Program of Wenzhou, China (Grant No: 2018ZG002), by the Science and Technology Program of Wenzhou, China (Grant No.X20210097).

    The authors declare there is no conflict of interest.



    [1] A. L. Hodgkin and A. F. Huxley, A quantitative description of membrane current and its applica- tion to conduction and excitation in nerve, J. Physiol., 117 (1952), 500–544.
    [2] M. L. Lapicque, Recherches quantitatives sur l'excitation électrique des nerfs traitée comme une polarisation, (French) [Quantitative research on electrical excitation of nerves treated as a polarization], J. De Physiol. Pathol. Gen., 9 (1907), 620–635.
    [3] R. B. Stein, A theoretical analysis of neuronal variability, Biophys. J., 5 (1965), 173–194.
    [4] P. Lánsk´ y and V. Lánská, Diffusion approximation of the neuronal model with synaptic reversal potentials, Biol. Cybern., 56 (1987), 19–26.
    [5] V. Capasso and D. Bakstein, An Introduction to Continuous-Time Stochastic Processes. The- ory, Models, and Applications to Finance, Biology, and Medicine, Springer-Verlag, New York, Heidelberg, Dordrecht, London, 2012.
    [6] S. Olmi, D. Angulo-Garcia, A. Imparato, et al., Exact firing time statistics of neurons driven by discrete inhibitory noise, Sci. Rep., 7 (2017), 1577.
    [7] A. Di Crescenzo and B. Martinucci, Analysis of a stochastic neuronal model with excitatory inputs and state-dependent effects, Math. Biosci., 209 (2007), 547–563.
    [8] E. N. Brown, R. Barbieri, V. Ventura, et al., The time-rescaling theorem and its application to neural spike train data analysis, Neural Comput., 14 (2001), 325–346.
    [9] H. C. Tuckwel, Introduction to Theoretical Neurobiology, Volume 2, Nonlinear and Stochastic Theories, Cambridge University Press, 1988.
    [10] R. Rudnicki and M. Tyran-Kami´ nska, Piecewise Deterministic Processes in Biological Models, Spinger-Verlag, 2017.
    [11] H. U. Bauer and K. Pawelzik, Alternating oscillatory and stochastic dynamics in a model for a neuronal assembly, Phys. D Nonlin. Phenom., 69 (1993), 380–393.
    [12] M. Walczak and T. Błasiak, Midbrain dopaminergic neuron activity across alternating brain states of urethane anaesthetized rat. Eur. J. Neurosci., 45 (2017), 1068–1077.
    [13] A. D. Kolesnik and N. Ratanov, Telegraph Processes and Option Pricing, Springer-Verlag, Heidelberg-New York-Dordrecht-London, 2013.
    [14] S. Zacks, Sample Path Analysis and Distributions of Boundary Crossing Times, Lecture Notes in Mathematics, vol. 2203, Springer-Verlag, 2017.
    [15] N. Ratanov, First crossing times of telegraph processes with jumps, Methodol. Comput. Appl. Probab., (2019). https://doi.org/10.1007/s11009-019-09709-5
    [16] L. Beghin, L. Nieddu and E. Orsingher, Probabilistic analysis of the telegrapher's process with drift by means of relativistic transformations, J. Appl. Math. Stoch. Anal., 14 (2001), 11–25.
    [17] L. Bogachev and N. Ratanov, Occupation time distributions for the telegraph process, Stoch. Process. Appl., 121 (2011), 1816–1844.
    [18] O. López and N. Ratanov, On the asymmetric telegraph processes, J. Appl. Prob. 51 (2014), 569–589.
    [19] N. Ratanov, Self-exciting piecewise linear processes, ALEA Lat. Am. J. Probab. Math. Stat. 14 (2017), 445–471.
    [20] A. A. Pogorui, R. M. Rodrguez-Dagnino and T. Kolomiets, The first passage time and estimation of the number of level-crossings for a telegraph process, Ukr. Math. J. 67 (2015), 998–1007.
    [21] G. D'Onofrio, C. Macci and E. Pirozzi, Asymptotic results for first-passage times of some expo- nential processes, Methodol. Comput. Appl. Probab., 20 (2018), 1453–1476.
    [22] A. Di Crescenzo and A. Meoli, On a jump-telegraph process driven by an alternating fractional Poisson process, J. Appl. Probab., 55 (2018), 94–111.
    [23] M. Abundo, On the first hitting time of a one-dimensional diffusion and a compound Poisson process, Methodol. Comput. Appl. Probab., 12 (2010), 473–490.
    [24] N. Ratanov, Option pricing model based on a Markov-modulated diffusion with jumps, Braz. J. Probab. Stat. 24 (2010), 413–431.
    [25] L. Breuer, First passage times for Markov-additive processes with positive jumps of phase type, J. Appl. Prob. 45 (2008), 779–799.
    [26] V. Srivastava, S. F. Feng, J. D. Cohen, et al., A martingale analysis of first passage times of time-dependent Wiener diffusion models, J. Math. Psychol., 77 (2017), 94–110.
    [27] A. N. Shiryaev, On martingale methods in the boundary crossing problems for brownian motion, Sovrem. Probl. Mat., 8 (2007), 3–78.
    [28] M. T. Giraudo and L. Sacerdote, Jump-diffusion processes as models for neuronal activity, BioSys- tems, 40 (1997), 75–82.
    [29] N. L. Johnson, S. Kotz and N. Balakrishnan, Continuous Univariate Distributions, Vol. 1 (Wiley Series in Probability and Statistics) Wiley-Interscience, 1994.
    [30] K. S. Lomax, Business failures; another example of the analysis of failure data, J. Amer. Stat. Assoc., 49 (1954), 847–852.
    [31] H. C. Tuckwell, Introduction to theoretical neurobiology: Volume 1, Linear Cable Theory and Dendritic Structure, Cambridge University Press, 1988.
    [32] A. P. Prudnikov, Y. A. Brychkov and O. I. Marichev, Integrals and Series, Vol. 5. Inverse Laplace Transforms, Gordon and Breach Science Publ. 1992.
  • This article has been cited by:

    1. Mudasir Younis, Haroon Ahmad, Mahpeyker Ozturk, Deepak Singh, A novel approach to the convergence analysis of chaotic dynamics in fractional order Chua’s attractor model employing fixed points, 2025, 110, 11100168, 363, 10.1016/j.aej.2024.10.001
  • Reader Comments
  • © 2019 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(4660) PDF downloads(613) Cited by(5)

Figures and Tables

Figures(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog