Research article Special Issues

Innovative approach for developing solitary wave solutions for the fractional modified partial differential equations

  • Received: 20 August 2023 Revised: 10 September 2023 Accepted: 19 September 2023 Published: 09 October 2023
  • MSC : 33B15, 34A34, 35A20, 44A10, 35A22

  • The current work investigates solitary wave solutions for the fractional modified Degasperis-Procesi equation and the fractional gas dynamics equation with Caputo's derivative by using a modified extended direct algebraic method. This method transforms the targeted fractional partial differential equations (FPDEs) into more manageable nonlinear ordinary differential equations, which are then turned into systems of nonlinear algebraic equations with a series-based solution assumption. Using Maple 13, the solitary wave solutions are then obtained by solving the obtained systems. The method produces multiple innovative solitary wave solutions for both equations, which are graphically depicted as 3D and 2D graphs and provide important insights into their behaviors. These insights help us to comprehend wave behavior and the physical processes represented by these equations. Furthermore, the suggested technique exhibits dependability and efficacy in dealing with complicated FPDEs, which bodes well for future studies on the subject.

    Citation: Saima Noor, Azzh Saad Alshehry, Asfandyar Khan, Imran Khan. Innovative approach for developing solitary wave solutions for the fractional modified partial differential equations[J]. AIMS Mathematics, 2023, 8(11): 27775-27819. doi: 10.3934/math.20231422

    Related Papers:

    [1] Kimun Ryu, Wonlyul Ko . Stability and bifurcations in a delayed predator-prey system with prey-taxis and hunting cooperation functional response. AIMS Mathematics, 2025, 10(6): 12808-12840. doi: 10.3934/math.2025576
    [2] Teekam Singh, Ramu Dubey, Vishnu Narayan Mishra . Spatial dynamics of predator-prey system with hunting cooperation in predators and type I functional response. AIMS Mathematics, 2020, 5(1): 673-684. doi: 10.3934/math.2020045
    [3] Na Min, Hongyang Zhang, Xiaobin Gao, Pengyu Zeng . Impacts of hunting cooperation and prey harvesting in a Leslie-Gower prey-predator system with strong Allee effect. AIMS Mathematics, 2024, 9(12): 34618-34646. doi: 10.3934/math.20241649
    [4] Weili Kong, Yuanfu Shao . Bifurcations of a Leslie-Gower predator-prey model with fear, strong Allee effect and hunting cooperation. AIMS Mathematics, 2024, 9(11): 31607-31635. doi: 10.3934/math.20241520
    [5] Ali Yousef, Ashraf Adnan Thirthar, Abdesslem Larmani Alaoui, Prabir Panja, Thabet Abdeljawad . The hunting cooperation of a predator under two prey's competition and fear-effect in the prey-predator fractional-order model. AIMS Mathematics, 2022, 7(4): 5463-5479. doi: 10.3934/math.2022303
    [6] Vikas Kumar, Nitu Kumari . Controlling chaos in three species food chain model with fear effect. AIMS Mathematics, 2020, 5(2): 828-842. doi: 10.3934/math.2020056
    [7] Yalong Xue, Fengde Chen, Xiangdong Xie, Shengjiang Chen . An analysis of a predator-prey model in which fear reduces prey birth and death rates. AIMS Mathematics, 2024, 9(5): 12906-12927. doi: 10.3934/math.2024630
    [8] Xiaoming Su, Jiahui Wang, Adiya Bao . Stability analysis and chaos control in a discrete predator-prey system with Allee effect, fear effect, and refuge. AIMS Mathematics, 2024, 9(5): 13462-13491. doi: 10.3934/math.2024656
    [9] Wei Li, Qingkai Xu, Xingjian Wang, Chunrui Zhang . Dynamics analysis of spatiotemporal discrete predator-prey model based on coupled map lattices. AIMS Mathematics, 2025, 10(1): 1248-1299. doi: 10.3934/math.2025059
    [10] Weili Kong, Yuanfu Shao . The effects of fear and delay on a predator-prey model with Crowley-Martin functional response and stage structure for predator. AIMS Mathematics, 2023, 8(12): 29260-29289. doi: 10.3934/math.20231498
  • The current work investigates solitary wave solutions for the fractional modified Degasperis-Procesi equation and the fractional gas dynamics equation with Caputo's derivative by using a modified extended direct algebraic method. This method transforms the targeted fractional partial differential equations (FPDEs) into more manageable nonlinear ordinary differential equations, which are then turned into systems of nonlinear algebraic equations with a series-based solution assumption. Using Maple 13, the solitary wave solutions are then obtained by solving the obtained systems. The method produces multiple innovative solitary wave solutions for both equations, which are graphically depicted as 3D and 2D graphs and provide important insights into their behaviors. These insights help us to comprehend wave behavior and the physical processes represented by these equations. Furthermore, the suggested technique exhibits dependability and efficacy in dealing with complicated FPDEs, which bodes well for future studies on the subject.



    In recent ecological research, the intricate dynamics inherent in three-species food chain models have sparked substantial interest, prompting a deeper understanding of complex ecological relationships [2,8,9,12,23,25,33,43,44]. This fascination has driven investigations into omnivory, a phenomenon where species feed across multiple trophic levels, with particular attention directed towards intraguild predation (IGP) [19,48,50,51]. IGP, a specialized form of omnivory, involves intricate interactions between a resource (prey), an intermediate predator preying exclusively on the prey and a top predator (referred to as an intraguild predator or IG predator) that consumes both the prey and the intermediate predators [14,16]. Given its pronounced impact on species coexistence and the emergence of chaos, the dynamics of IGP have become a central topic in ecological studies. Panja [38] analyzed a fractional-order prey-predator model with intraguild predation and observed instability in the system for integer-order derivatives. In contrast, the system remained stable with fractional-order derivatives. In [20], the authors explored intraguild predation dynamics, focusing on the impact of interference and predation strength on the coexistence of intraguild predators and prey. Gómez-Hernández et al. [13] explored an intraguild predation model and found that interspecific competition consistently reduces infection rates in both species.

    Previous research has unveiled various mathematical models to capture the complex dynamics of IGP and its intricate interactions within predator-prey systems. In the realm of ecological modeling, diffusion stands as a fundamental factor closely tied to real environmental processes, influencing a myriad of dynamic behaviors. A wealth of research has delved into the dynamics of diffusive versions of ecological models, shedding light on their intricate mechanisms. Notably, Chen and Wu, in their work [4], illuminated the role of attractive predator-taxis in driving Turing bifurcation and the consequent emergence of spatiotemporal patterns within a predator-prey model characterized by nonlinear harvesting. Conversely, they found that the presence of repulsive predator-taxis or the complete absence thereof did not give rise to such phenomena. Moreover, this study provides valuable insights into how diffusion ratios and nonlinear harvesting practices impact the spatial distribution of species and the direction of Turing bifurcation.

    In [5], the authors undertook a study focusing on the stability of bifurcating solutions within a predator-prey diffusive model with prey-taxis, where they not only confirmed the occurrence of steady state bifurcations but also explored the influence of repulsive prey-taxis in promoting this phenomenon. This research further delved into the determination of bifurcating solution stability. These collective findings offer valuable insights into the realm of ecological dynamics and diffusion-based ecological modeling. Hastings and Powell [15] highlighted the potential for chaos within three-species food chain models. On the other hand, Holt and Polis [17] introduced a Lotka-Volterra model characterized by linear functional responses, revealing that increased intensity of IGP could induce destabilization within the system [26]. Research in the field of intraguild predation (IP) has led to diverse insights into predator-prey dynamics. Hsu et al. [18] examined IGP scenarios using Beddington-DeAngelis responses, studying coexistence and competitive exclusion. They highlighted the significance of the interference parameter in shaping stability or oscillation. Kang and Wedekin [21] explored type III functional responses for intermediate predators and IG predators, revealing distinctions between generalist and specialist IG predators, and identifying multiple attractors. Freeze et al. [10] introduced ratio-dependent functional responses, investigating persistence, permanence and stability through innovative methods. These studies have laid crucial groundwork for understanding predator-prey interactions and set the stage for our exploration of cooperative hunting effects in the context of IGP dynamics. Building on this foundation, Tanabe and Namba [47] explored the emergence of chaotic dynamics within IGP scenarios, omitting intraspecific competition among IG prey and predators. Recent advancements by Sen et al. [44] and Majumder et al. [25] extended these models by incorporating saturating functional responses and intraspecific competition among IG prey and predator, leading to enhanced insights into the rich dynamics such as chaos of predator-prey interactions. However, Majumder et al. [25] control chaotic behavior of the system through a memory effect.

    Furthermore, cooperative hunting behaviors have emerged as crucial factors shaping ecological stability. A range of adaptive characteristics have evolved among predators, including sharp teeth, claws, fangs and other tools, to facilitate efficient hunting. Predators have developed various strategies such as bait pursuit, collective hunting and other approaches. For instance, carnivores [24], lions [36], wolves [42] and African wild dogs [7] frequently collaborate to capture and subdue their prey. Similarly, ants [27] and spiders [49] exhibit cooperative behaviors in hunting, drawing considerable attention. The concept of hunting cooperation was initially introduced by Cosner et al. [6] in a predator-prey model. He proposed a functional response for predators that engage in spatially linear hunting formations and aggregate when encountering significant prey groups. Cooperative foraging has been explored in diverse predator-prey models, revealing scenarios of both destabilization and stabilization [1,3,6,11,31,32,34]. The incorporation of cooperative hunting behavior among predators in the prey-predator model reveals novel insights into the complex dynamics of ecological systems [45,46]. This growing understanding of cooperative behaviors' influence on ecological dynamics provides a strong motivation for our study.

    We undertake a comprehensive investigation into the interplay between hunting cooperation among IG predators and IG prey, aiming to elucidate its impact on the intricate dynamics of predator-prey interactions. Our objective is to decipher how this cooperative dimension influences the emergence of chaotic oscillations and its potential to shape the overall stability or instability of ecological communities. By incorporating this cooperative aspect into our model, we aspire to provide novel insights into the intricate interactions between predator and prey, specifically within the context of IGP and the cooperative behaviors that can arise in these ecological systems. The distinctive feature of our study, setting it apart from existing literature [25,41,44,47] addressing intraguild predation (IGP) scenarios, is the inclusion of cooperative hunting behavior within the model. Through this endeavor, we seek to contribute to the deeper comprehension of the complexities inherent in predator-prey interactions and the factors that govern their dynamics.

    We focus on a three-species Lotka-Volterra model that incorporates intraguild predation dynamics. The model involves three major species denoted as x, y and z, representing a basal prey, an intermediate consumer (IG prey) and an omnivorous top predator (IG predator), respectively. We adopt a linear functional response between these species, as discussed in [47]. Additionally, we account for intraspecific competition within the intraguild prey and predator populations, as elaborated in [44]. Given the limited availability of resources, competition becomes essential for their sustenance within the community. In our model, we introduce the concept of hunting cooperation, which confers benefits to both intraguild prey and predator. Moreover, the success in capturing a basal prey and intraguild prey increases with the density of intraguild prey and predator, respectively, as highlighted in [37]. The model can be expressed as follows:

    dxdt=r11xd11xα11x2(α12+c11y)xyα13xz=F1(x,y,z),dydt=e11(α12+c11y)xyα22y2(α23+c22z)yzd22y=F2(x,y,z),dzdt=e22(α23+c22z)yz+α31xzα33z2d33z=F3(x,y,z), (2.1)

    with non-negative initial conditions. The model parameters r11, d11, α11, α12, c11, α13, e11, α22, α23, c22, d22, e22, α31, α33 and d33 are all positive constants and their ecological descriptions are provided in Table 1. Our model is distinct from previous works like Kang and Wedekin [21] and Tanabe and Namba [47] due to the incorporation of intraspecific competitions and hunting cooperation behavior among intraguild prey and predator. Furthermore, compared to models studied by Sen et al. [44] and Majumder et al. [25], our model employs different functional response types between intraguild prey and predator, while also involving cooperative hunting.

    Table 1.  Details of parameters with their values and sources.
    Parameters Description Values Source
    r11 Birth rate of basal prey 5 [25]
    d11 Death rate of basal prey 0.1 [28]
    α11 Intraspecific competition within basal prey 0.38 [25]
    α12 Consumption rate of IG prey to basal prey 1 [25]
    c11 IG prey cooperation during hunting basal prey 0.15 [39]
    α13 Consumption rate of IG predator to basal prey 35 [25]
    e11 Conversion coefficient of basal prey to IG prey 0.9 [29]
    α22 Intraspecific competition within IG prey 0.2 [25]
    α23 Consumption rate of IG predator to IG prey 1 [25]
    c22 IG predator cooperation during hunting IG prey 0.9 [39]
    d22 Natural death rate of IG prey 1 [25]
    e22 Conversion coefficient of IG prey to IG predator 0.95 [30]
    α31 basal prey consumption into reproduction for IG predator 0.1 [25]
    α33 Intraspecific competition within IG predator 0.3 [25]
    d33 Natural death rate of IG predator 1.2 [25]

     | Show Table
    DownLoad: CSV

    The positivity of state variables is a fundamental concept in ecological modeling, indicating that the populations of the involved species remain non-negative within the system. Ensuring the positivity of solutions in the proposed model (2.1) is crucial for understanding the persistence of populations.

    To establish the positivity of the system (2.1), we integrate its equations and utilize the given initial conditions to obtain the following expressions:

    x(t)=x(0)exp(t0[r11d11α11x(s)(α12+c11y(s))y(s)α13z(s)]ds),y(t)=y(0)exp(t0[e11(α12+c11y(s))x(s)α22y(s)(α23+c22z(s))z(s)d22]ds),z(t)=z(0)exp(t0[e22(α23+c22z(s))y(s)+α31x(s)α33z(s)d33]ds).

    The positivity of these solutions follows from the fact that the right-hand side of each equation is non-negative given the non-negative initial conditions. Consequently, all solutions of the system (2.1) are guaranteed to be non-negative.

    While we have established the positivity of solutions, we omit the discussion of the boundedness of solutions, as the boundedness property can be similarly proven and is straightforward in this context.

    An equilibrium point in a dynamical system is a state or configuration where the system's variables remain constant over time. In other words, it's a point at which the system's behavior does not change, and the rates of change of the variables are all zero. Equilibrium points are also known as fixed points or steady states. Mathematically, for a dynamical system described by a set of equations, an equilibrium point is found by setting all the derivatives of the variables to zero. Equilibrium points play a crucial role in the analysis of dynamical systems. They serve as important reference states to understand the behavior of the system around those points.

    The proposed IGP system has the following feasible equilibria.

    (1) Population-free equilibrium E0(0,0,0) exists irrespective of any parametric restriction. From an ecological viewpoint, this equilibrium represents a theoretical scenario where all species involved in the model are entirely absent from the ecosystem. While such a scenario may not directly correspond to a real ecological state, it is a critical reference point.

    (2) Basal prey equilibrium E1(r11d11α11,0,0) exists under the parametric condition r11>d11. From an ecological viewpoint, this equilibrium reflects a scenario where the prey population is the only species present with intermediate and top predator both at zero. Moreover, the existence of E1 suggests that when the prey species experiences a birth rate greater than its death rate, it can maintain its population in the absence of predators.

    (3) Basal prey-IG prey equilibrium E2(x2,y2,0), where x2 and y2 are positive solution(s) of following isoclines:

    r11d11α11x(α12+c11y)y=0, (4.1)
    e11(α12+c11y)xα22yd22=0. (4.2)

    We can obtain x from first equation and then putting the value of x in the second equation we get a cubic equation for y. This is little bit difficult to find the root (i.e., the component of equilibrium point). Therefore, we numerically check the existence E2 by plotting the above isoclines in xy space as depicted in Figure 1. Figure exhibits unique E2 for the parameter values mentioned in the Table 1. The ecological significance of E2 lies in its representation of a complex ecological interaction where both basal prey and intraguild prey species coexist without the presence of intraguild predators. This equilibrium suggests the existence of resource partitioning mechanisms or the presence of ecological niches that allow both the basal prey and intraguild prey to thrive without direct predation pressure from higher trophic levels.

    Figure 1.  Isoclines (4.1) and (4.2) are plotted here. Red dot represents the intersection of these lines. Parameters are chosen from Table 1.

    (4) Basal prey-IG predator equilibrium E3(α13d33α33d11+α33r11α11α33+α13α31,0,α31r11α11d33α31d11α11α33+α13α31). E3 exists if the parametric conditions α13d33+α33r11>α33d11 and α31r11>α11d33+α31d11 hold. Ecologically, this equilibrium illustrates the mechanisms that allow basal prey and intraguild predator to coexist without one driving the intermediate predator to extinction. By understanding the conditions that lead to the coexistence of these two species, ecologists can make more informed decisions regarding the preservation of species and ecosystem health.

    (5) Coexistence equilibrium E4(x4,y4,z4), where the components x4, y4 and z4 can be obtained from following isoclines.

    r11d11α11x(α12+c11y)yα13z=0, (4.3)
    e11(α12+c11y)xα22y(α23+c22z)zd22=0, (4.4)
    e22(α23+c22z)y+α31xα33zd33=0. (4.5)

    Due to the difficulties of finding parametric values of the components of E4, we find it numerically by plotting above isoclines in xyz space as illustrated in Figure 2. It is apparent from the figure that for table parameter, E4 exists uniquely. The coexistence equilibrium E4 represents a critical ecological scenario within a model where all three interacting species coexist within an ecosystem. This equilibrium reflects the richness and stability of biodiversity, highlighting the ability of different species to thrive without the extinction of any species.

    Figure 2.  Equations (4.3)–(4.5) are plotted here. Red dot denotes intersecting point of the planes, which is the E4. Parameters values are same as in Table 1.

    Mathematically, local stability analysis involves linearizing the system's equations of motion around the equilibrium point. This is typically done by calculating the system's Jacobian matrix at the equilibrium, which describes how small changes in the variables affect the rates of change of the variables themselves. The eigenvalues of this Jacobian matrix are crucial indicators of the local stability properties of the system.

    Local stability analysis provides valuable insights into the qualitative behavior of a dynamical system around an equilibrium point. However, it is important to note that local stability analysis only considers behavior near the equilibrium and may not fully capture more complex behaviors, such as oscillations, bifurcations and global attractors, that can arise in nonlinear systems. Here, we check local stability of the proposed system around the system's equilibria. First, we evaluate the Jacobian matrix of the system (2.1) as illustrated below.

    J(x,y,z)=[A11c11xy(c11y+α12)xα13xe11(c11y+α12)yA22c22yz(c22z+α23)yα31ze22(c22z+α23)zA33],

    where

    A11=2α11xd11+r11(c11y+α12)yα13z,A22=e11c11xy+e11(c11y+α12)x2α22y(c22z+α23)zd22,A33=e22c22zy+e22(c22z+α23)y+α31x2α33zd33.

    (1) Eigenvalues of the Jacobian matrix J(E0) are d11+r11, d22 and d33. From this eigenvalue analysis, we can say that E0 is stable under the parametric condition d11>r11. This stability condition is opposite to the existence condition of E1. Thus, E0 and E1 are related via a transcritical bifurcation. Further, we verify the stability of E0 numerically as described in Figure 3. In the figure, we choose r11=0.05 and d11=0.1. It is apparent from this figure that E0 is stable as condition d11>r11 holds.

    (2) Eigenvalues of the Jacobian matrix J(E1) are d11r11, α12d11e11α12e11r11+α11d22α11 and α11d33+α31d11α31r11α11. Hence, E1 is locally asymptotically stable if the parametric conditions d11<r11, α12d11e11+α11d22>α12e11r11 and α11d33+α31d11>α31r11 hold. E1 loses its stability when either of the conditions fails. We verify its local stability numerically also as described in Figure 4. The parameters used in this figure satisfy stability conditions, as a result, it is stable as depicted in the figure.

    (3) Jacobian matrix evaluated at E2 is given by

    J(E2)=[2α11x2d11+r11(c11y2+α12)y2c11x2y2(c11y2+α12)x2α13x2e11(c11y2+α12)y2e11c11x2y2+e11(c11y2+α12)x22α22y2d22α23y200e22α23y2+α31x2d33].
    Figure 3.  Stability of E0. Here, the green dot denotes the stable equilibrium point. Parameters are taken from Table 1 except r11=0.05.
    Figure 4.  Stability of E1. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=0.5.

    One eigenvalue of J(E2) is e22α23y2+α31x2d33 and the other two eigenvalues are the roots of the characteristic polynomial of the matrix

    ˜J(E2)=[2α11x2d11+r11(c11y2+α12)y2c11x2y2(c11y2+α12)x2e11(c11y2+α12)y2e11c11x2y2+e11(c11y2+α12)x22α22y2d22].

    Now, eigenvalues of ˜J(E2) have negative real part if tr(˜J(E2))<0 and det(˜J(E2))>0. Thus, E2 is stable under the following conditions

    e22α23y2+α31x2<d33, tr(˜J(E2))<0  and  det(˜J(E2))>0.

    As it is difficult to follow these conditions for stability of E2, we numerically verify its stability, Figure 5. The parameters used for numerical verification satisfy above stability conditions of E2. As a result E2 is locally asymptotically stable as described in the figure. (5) The Jacobian matrix evaluated at E3 is given by

    J(E3)=[B11B12B13B21B22B23B31B32B33],
    Figure 5.  Stability of E2. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=0.75.

    where

    B11=2α11(α13d33α33d11+α33r11)α11a33+α13α31d11+r11+α13(α11d33+α31d11α31r11)(α11α33+α13α31),B12=α12(α13d33α33d11+α33r11)(α11α33+α13α31),  B13=α13(α13d33α33d11+α33r11)(α11α33+α13α31),  B21=0,B22=e11α12(α13d33α33d11+α33r11)(α11α33+α13α31)+(((d11+r11)c22+α13α23)α31+a11(α23α33c22d33))((d11r11)α31+d33α11)(α11α33+α13α31)2d22,B23=0,  B31=(α11d33+α31d11α31r11)α31(α11α33+α13α31),B32=e22(α11α23α33α11c22d33+α13α23α31α31c22d11+α31c22r11)(α11d33+α31d11α31r11)(α11α33+α13α31)2,B33=α31(α13d33α33d11+α33r11)(α11α33+α13α31)+2α33(α11d33+α31d11α31r11)(α11α33+α13α31)d33.

    The characteristic equation associated with the matrix J(E3) is given by

    λ3+Γ2λ2+Γ1λ+Γ0=0, (5.1)

    where

    Γ2=tr(J(E3)),  Γ1=(B11B22+B22B33+B33B11),  Γ0=det(J(E3)).

    By Routh-Hurwitz criteria, E3 is stable if Γ2>0, Γ0>0 and Γ2Γ1>Γ0. The expression of Γi, i=0,1,2, are large to say about the stability E3. So, we numerically investigate its stability as depicted in Figure 6. It is apparent from the figure that E3 is stable since the parameter used here holds its stability conditions.

    (6) The Jacobian matrix evaluated at a representative coexisting equilibrium point

    E4=(x4,y4,z4)
    Figure 6.  Stability of E3. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except e11=0.01.

    is denoted as

    J(E4)=[J4ij]3×3

    with J4ij representing the matrix components when J(x,y,z) is computed at

    E4=(x4,y4,z4).

    The associated characteristic equation for the matrix J(E4) is given by

    λ3+Δ2λ2+Δ1λ+Δ0=0, (5.2)

    where the coefficients are defined as follows:

    Δ2=tr(J(E4)),  Δ1=(J411J422+J422J433+J433J411),  Δ0=det(J(E4)).

    The stability of E4 can be determined using the Routh-Hurwitz criteria, which states that E4 is stable if Δ2>0, Δ0>0 and Δ2Δ1>Δ0. The explicit expressions for the components of E4 are not obtainable due to their complexity, leading to the inability to derive parametric restrictions explicitly for the stability of E4. To validate the stability of E4, a numerical example is employed. The chosen parameter values in Figure 7 adhere to the Routh-Hurwitz criteria for the stability of E4. Consequently, E4 is demonstrated to be stable, as depicted in Figure 7.

    Figure 7.  Stability of E4. Here, green and red dots denote the stable and unstable equilibrium points, respectively. Parameters are taken from Table 1 except r11=2.5.

    Bifurcations in a dynamical system are critical points where the qualitative behavior of the system changes as one or more parameters are varied. In other words, they mark a point at which the stability, periodicity, or overall behavior of the system undergoes a significant transformation. Bifurcations can lead to the emergence of new patterns, such as oscillations, chaos or the creation and destruction of equilibrium points. Bifurcations are fundamental to the study of nonlinear dynamical systems and play a pivotal role in understanding the complex behaviors that can arise in such systems. They are often associated with changes in the number or stability of equilibrium points or the appearance of periodic solutions. Bifurcation analysis involves studying how the system's behavior changes as the parameters are continuously varied. In our proposed IGP system, two types of bifurcations occur: Transcritical and Hopf, which are described in this section.

    Theorem 1. The system described by (2.1) undergoes a transcritical bifurcation at the parameter value

    r11=r[TB]11=d11

    concerning the equilibrium point E0(0,0,0).

    Proof. Let us consider the Jacobian matrix J(E0)|r[TB]11. One of its eigenvalues is zero, and the other two are negative. Since two eigenvalues are consistently negative and the remaining one is zero, concluding stability-instability is not feasible at this point.

    Let

    V=(1,0,0)TandW=(1,0,0)T

    be the eigenvectors of J(E0) and J(E0)T corresponding to the zero eigenvalue, respectively. Next, we define

    F=[F1, F2, F3]T.

    The conditions outlined by Sotomayor's theorem [40] for a transcritical bifurcation are as follows:

    (1) WTFr11(E0,r[TB]11)=0,

    (2) WT[DFr11(E0,r[TB]11)V]=0,

    (3) WT[D2F(E0,r[TB]11)(V,V)]=2α110.

    It's important to note that for a non-degenerate transcritical bifurcation, the second transversality condition should not be zero [40]. Hence, the system defined by (2.1) experiences a degenerate transcritical bifurcation [22] around E0 when

    r11=r[TB]11.

    Consequently, E0 undergoes a transcritical bifurcation as the parameter r11 crosses the critical value

    r11=r[TB]11=d11,

    and the conditions outlined in the theorem are satisfied. Next, we will establish the presence of a Hopf bifurcation concerning the model parameter r11. We have designated the critical value of r11 as r[HB]11.

    Theorem 2. The system defined by (2.1) exhibits a Hopf bifurcation if the condition

    (Δ13ω2)(dΔ0dr11ω2dΔ2dr11)+2Δ2ω2dΔ1dr110

    is met when the parameter r11 crosses the critical value r11=r[HB]11.

    Proof. Considering that the equilibrium point E4 is dependent on the parameter r11, we apply the Hopf bifurcation theorem [29]. The Hopf bifurcation occurs if the variation matrix J(E4) has one eigenvalue with a negative real part, where

    Re(dλdr11)r11=r[HB]110,

    and the other two eigenvalues are purely imaginary. For purely imaginary eigenvalues, the characteristic Eq (5.2) should satisfy Δ2Δ1Δ0=0. We assume that iω represents the purely imaginary eigenvalue corresponding to the parameter r[HB]11.

    By differentiating both sides of the characteristic Eq (5.2) with respect to r11, we obtain

    (3λ2+2Δ2λ+Δ1)dλdr11+λ2dΔ2dr11+λdΔ1dr11+dΔ0dr11=0,dλdr11=λ2dΔ2dr11+λdΔ1dr11+dΔ0dr113λ2+2Δ2λ+Δ1.

    Substituting λ=iω and considering the real part, we arrive at

    (Re(dλdr11))λ=iω=(Δ13ω2)(dΔ0dr11ω2dΔ2dr11)+2Δ2ω2dΔ1dr11(Δ13ω2)2+(2Δ2ω)2.

    Therefore, the system described by (2.1) experiences a Hopf bifurcation concerning the parameter r11 if the condition

    (Δ13ω2)(dΔ0dr11ω2dΔ2dr11)+2Δ2ω2dΔ1dr110

    is satisfied at r11=r[HB]11.

    In this section, we will delve into our bifurcation analysis results. We will start by examining the bifurcation structure concerning the parameter r11, as depicted in Figure 8.

    Figure 8.  Bifurcation structure of the proposed model. Parameters are adopted from Table 1.

    This figure uncovers dynamic behaviors as r11 undergoes variations. For lower values of r11, the system exclusively features the stable equilibrium E0. The presence of any other equilibria is absent within this parameter range. As r11 increases, a critical value is reached, leading to the emergence of the stable equilibrium E1. Simultaneously, the previously stable equilibrium E0 shifts from stability to instability, retaining this behavior for all subsequent higher values of r11. With further increments in r11, another transcritical point comes into view. At this juncture, the stable equilibrium E2 arises, causing the previously stable E1 to transition into an unstable state. This instability of E1 continues for higher values of r11. Finally, the system encounters an additional transcritical point, resulting in the emergence of the stable equilibrium E4, while the previously stable E2 transforms into an unstable state.

    Similarly, Figure 9 presents a transcritical bifurcation diagram involving the parameter e11. For lower values of e11, the system displays a stable equilibrium E2. Beyond a critical value of e11, the stable equilibrium E4 comes into play. Consequently, the previously stable E2 shifts into an unstable state. These analyses provide valuable insights into the system's behavior with varying parameters r11 and e11. They reveal critical points where equilibria undergo shifts in stability and new equilibria emerge, contributing to a deeper comprehension of predator-prey interactions within the context of changing parameter values. Transcritical bifurcation is ecologically significant as it marks a critical point in the dynamics of predator-prey relationships and can lead to shifts in population equilibrium, which, in turn, can impact the stability and biodiversity of ecosystems. Understanding these dynamics is crucial for effective ecological management and conservation efforts.

    Figure 9.  Transcritical bifurcation. Parameters are adopted from Table 1. Critical value of e11 is e[TC]11=0.083.

    Continuing our exploration, we now delve into the system's behavior surrounding the coexistence equilibrium E4. This investigation involves the examination of phase portraits and time series solutions for various values of r11, as showcased in Figure 10. As the parameter r11 undergoes increments, the system's dynamics around the coexistence equilibrium E4 progress through distinct phases, leading to diverse patterns of behavior. For lower values of r11, the system exhibits a stable behavior around the coexistence equilibrium, manifesting as confined trajectories in the phase portrait, Figure 10a. As r11 increases, the system's behavior shifts to one-period oscillations. Trajectories trace out regular, repetitive patterns in the phase portrait, Figure 10b. Further increments in r11 lead to the emergence of two-period oscillations, Figure 10c. Trajectories now complete two cycles before repeating themselves in the phase portrait. The system continues to evolve with higher values of r11, resulting in four-period oscillations, Figure 10d. Trajectories now complete four cycles within the same phase portrait. As r11 further increases, the system transitions into chaotic behavior, Figure 10e. Trajectories become highly sensitive to initial conditions, leading to unpredictable and complex patterns in the phase portrait. Moreover, we observed that the chaotic attractors size increases on increaing r11 (see Figure 11a). Finally, for larger values of r11 both the predator populations extinct making the axial equilibrium point E1 stable (see Figure 11b). In essence, Figure 10 vividly portrays the evolution of the system's dynamics around the coexistence equilibrium E4 in response to changing values of r11. This comprehensive visualization allows us to discern the progression from stable behavior to chaotic complexity and highlights the intricate nature of predator-prey interactions in the context of varying parameter settings. Chaotic dynamics in predator-prey systems can lead to irregular and unpredictable population cycles. This unpredictability can help maintain biodiversity by preventing any one species from dominating the ecosystem for extended periods, fostering coexistence and preventing extinctions.

    Figure 10.  Phase portraits and respective time series for the parameter r11=2.58 in (a); r11=2.85 in (b); r11=3.9 in (c); r11=4.35 in (d); r11=5 in (e). Rest are adopted from Table 1.
    Figure 11.  Time series for the parameter r11=8 in (a) and r11=15.11 in (b). Rest are adopted from Table 1.

    Moving forward, our analysis takes us to the exploration of species' maximum and minimum densities while systematically varying a key model parameter, as detailed in Figure 12. The insights gleaned from this visualization provide a comprehensive understanding of how the system's behavior evolves in response to specific parameter alterations. Figure 12a provides a revealing depiction of the system's response to changes in the parameter r11. Initially, the system sustains a stable coexistence state. This phase occurs after a transcritical bifurcation point where a switch occurs between equilibrium states. Then the system enters a phase of one-period oscillations, a behavior characterized by regular and predictable cycles. As r11 continues to increase, the system progresses into a phase of many-period oscillations. This pattern is marked by an increasing number of cycles within each oscillatory pattern, indicating a higher level of complexity in the population dynamics. Ultimately, the system experiences a transition into a chaotic state. In this phase, trajectories in the system become highly sensitive to initial conditions, leading to unpredictable and irregular behavior. This chaotic behavior underscores the intricate interactions that can arise within an ecosystem. Similar patterns of behavior are observed for the parameter α23 (Figure 12f), with the primary distinction being the absence of an initial transcritical bifurcation. Furthermore, we observe that parameter c22 (Figure 12b) influences the system's dynamics by inducing two cycles, leading to an increase in the number of cycles and finally chaos. The presence of chaos and transitions between stable states and periodic states can reflect the system's ability to adapt to changing environmental conditions. Ecosystems that exhibit such behavior may be more resilient, can support biodiversity, are capable of adjusting to disturbances and evolving in response to environmental shifts. On the other hand, parameters α11 (Figure 12c), α22 (Figure 12d), d11 (Figure 12e) and α12 (Figure 12g) exert control over the chaotic behavior by modulating the transition from many-period oscillations to one-period oscillations, bringing about a more ordered pattern and ultimately stable coexistence. The ability to control chaotic behavior through specific parameters is ecologically significant because it offers opportunities for managing, restoring and preserving ecosystems. Lastly, Figure 12h demonstrates that the parameter c11 is capable of driving the system towards more orderly behavior from an irregular behavior, as it produces a limit cycle through multiple cycles.

    Figure 12.  Bifurcation structure of the proposed system. Parameter values are adopted from Table 1 except r11=4.35 in {b, c, d, e, f, g, h}.

    Finally, we present the two-parameter bifurcation structure, varying the IG prey cooperation during hunting basal prey parameter (c11) and Consumption rate of IG predator to IG prey parameter (α23, as shown in Figure 13. This figure reveals the division of the entire parameter space into five subregions, each exhibiting distinct dynamics. The two-parameter bifurcation structure provide a comprehensive view of the model's behavior under varying two parameters, highlighting the occurrence of stable regions, limit cycles, two period, four period and chaotic dynamics. We note that for low values of consumption rate of IG predator, the cooperation of IG predator helps the system to keep in a stable state position. However, increasing in the consumption rate of IG predators push on the system into chaotic oscillations through periodic-doubling oscillations.

    Figure 13.  Two parameter bifurcation structure. Table 1 is used for parameter values.

    In conclusion, we delve into a three-species Lotka-Volterra model embodying intraguild predation dynamics. The model incorporates complex interactions between a basal prey, an intraguild prey and an omnivorous top predator. These interactions are defined by linear functional responses, alongside considerations for intraspecific competition and cooperative hunting, both of which play essential roles in shaping population dynamics within the ecological community. The proposed model diverges from prior works through its integration of intraspecific competition and cooperative hunting. We analyze equilibria and parametric conditions under which they are stable. Numerical stability of all steady states are justified. Our exploration encompasses a comprehensive bifurcation analysis that reveals intricate behaviors driven by parameter variations.

    The examination of bifurcation structures, particularly concerning parameters such as r11 and e11, showcases critical shifts in equilibria stability and the emergence of new equilibrium states. We uncover transitions from stable to unstable states and vice versa, offering insights into the interplay between predator-prey interactions and varying parameter values. Moreover, our investigation into the dynamics around the coexistence equilibrium E4 illuminates a diverse range of behaviors as parameter r11 evolves. These behaviors include stable coexistence, one-period oscillations, multi-period oscillations and ultimately chaotic dynamics. This detailed analysis provides a nuanced understanding of the system's response to changing parameter values.

    Furthering our comprehension, the exploration of species' maximum and minimum densities in response to parameter variations demonstrates distinct patterns. These patterns range from stable coexistence to chaotic behavior, each induced by different parameter influences. Parameters such as α12, c11 and d11 exert significant control over the system's chaotic dynamics, promoting stable or orderly behavior. Strong hunting cooperation behavior of intraguild predator produces chaos in the system whereas the same behavior of intraguild prey can control the chaotic nature of the system. In contrast to the findings of Sen et al. [44], our research demonstrates that a chaotic regime can emerge from a state of stable coexistence when considering the consumption rate of intraguild predators towards intraguild prey it may be due to the assumption of hunting cooperation. However, Vandermeer [50] observed chaotic oscillations in an intraguild predation (IGP) system in the absence of self-regulation for the same parameter. A recent study by Namba et al. [35] highlighted the stabilizing influence of strong self-regulation on a similar tri-trophic food web featuring type-I functional responses. Sen et al. [44] reported that robust self-regulation within both IG prey and predator populations can drive the dynamics of the system from a chaotic state to stable steady-state coexistence due to the saturated consumption of IG prey by IG predators. In our research, we extend this concept by revealing that strong self-regulation has the potential to guide the system's dynamics from a chaotic regime to stable coexistence. This stabilization is facilitated by the phenomenon of hunting cooperation. Also for low values of consumption rate of IG predator, the cooperation of IG predator helps the system to keep in a stable state position but increasing in the consumption rate of IG predators push on the system into chaotic oscillations through periodic-doubling oscillations.

    In summary, this study offers a comprehensive understanding of complex predator-prey interactions in the context of intraguild predation dynamics. The integrated effects of intraspecific competition and cooperative hunting, coupled with systematic parameter variations, provide valuable insights into the intricate behaviors exhibited by ecological systems. Our findings converge as well as diverge from the existing literature. By shedding light on the intricate interplay between species and parameters, this research contributes to the broader field of ecological modeling and enhances our understanding of the complex dynamics that characterize natural ecosystems.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This work has been supported by the Deanship of Scientific Research at King Khalid University through the Research Group program with Grant Number RGP2/202/44. The authors are very grateful to the editor and reviewers for their careful reading of the manuscript.

    The authors declare that they have no conflicts of interest.



    [1] H. Y. Jin, Z. A. Wang, L. Wu, Global dynamics of a three-species spatial food chain model, J. Differ. Equ., 333 (2022), 144–183. https://doi.org/10.1016/j.jde.2022.06.007 doi: 10.1016/j.jde.2022.06.007
    [2] P. Liu, J. Shi, Z. A. Wang, Pattern formation of the attraction-repulsion Keller-Segel system, Discrete Cont. Dyn. S., 18 (2013), 2597–2625. https://doi.org/10.3934/dcdsb.2013.18.2597 doi: 10.3934/dcdsb.2013.18.2597
    [3] H. Khan, R. Shah, P. Kumam, M. Arif, Analytical solutions of fractional-order heat and wave equations by the natural transform decomposition method, Entropy, 21 (2022), 597. https://doi.org/10.3390/e21060597 doi: 10.3390/e21060597
    [4] H. Y. Jin, Z. A. Wang, Global stabilization of the full attraction-repulsion Keller-Segel system, arXiv Preprint, 2019.
    [5] L. Liu, S. Zhang, L. Zhang, G. Pan, J. Yu, Multi-UUV maneuvering counter-game for dynamic target scenario based on fractional-order recurrent neural network, IEEE T. Cybernetics, 53 (2022), 4015–4028. https://doi.org/10.1109/TCYB.2022.3225106 doi: 10.1109/TCYB.2022.3225106
    [6] A. M. Zidan, A. Khan, M. K. Alaoui, W. Weera, Evaluation of time-fractional Fishers equations with the help of analytical methods, AIMS Math., 7 (2022), 18746–18766. https://doi.org/10.3934/math.20221031 doi: 10.3934/math.20221031
    [7] R. Luo, Z. Peng, J. Hu, On model identification based optimal control and it's applications to multi-agent learning and control, Mathematics, 11 (2023), 906. https://doi.org/10.3390/math11040906 doi: 10.3390/math11040906
    [8] X. Lyu, X. Wang, C. Qi, R. Sun, Characteristics of cavity dynamics, forces, and trajectories on vertical water entries with two spheres side-by-side, Phys. Fluids, 35 (2023). https://doi.org/10.1063/5.0166794 doi: 10.1063/5.0166794
    [9] E. M. Elsayed, R. Shah, K. Nonlaopon, The analysis of the fractional-order Navier-Stokes equations by a novel approach, J. Funct. Space., 2022 (2022), 1–18. https://doi.org/10.1155/2022/8979447 doi: 10.1155/2022/8979447
    [10] M. Ekici, A. Sonmezoglu, E. M. Zayed, A new fractional sub-equation method for solving the space-time fractional differential equations in mathematical physics, Comput. Methods Diffe., 2 (2014), 153–170. https://doi.org/10.9734/BJMCS/2013/2908 doi: 10.9734/BJMCS/2013/2908
    [11] H. Khan, R. Shah, J. F. G. Aguilar, D. Baleanu, P. Kumam, Travelling waves solution for fractional-order biological population model, Math. Model. Nat. Pheno., 16 (2021), 32. https://doi.org/10.1051/mmnp/2021016 doi: 10.1051/mmnp/2021016
    [12] A. Ara, N. A. Khan, O. A. Razzaq, T. Hameed, M. A. Z. Raja, Wavelets optimization method for evaluation of fractional partial differential equations: An application to financial modelling, Adv. Differ. Equ., 2018 (2018), 1–13. https://doi.org/10.1186/s13662-017-1461-2 doi: 10.1186/s13662-017-1461-2
    [13] H. Sun, Y. Zhang, D. Baleanu, W. Chen, Y. Chen, A new collection of real world applications of fractional calculus in science and engineering, Commun. Nonlinear Sci., 64 (2018), 213–231. https://doi.org/10.1016/j.cnsns.2018.04.019 doi: 10.1016/j.cnsns.2018.04.019
    [14] P. Sunthrayuth, N. H. Aljahdaly, A. Ali, I. Mahariq, A. M. Tchalla, ψ-Haar wavelet operational matrix method for fractional relaxation-oscillation equations containing ψ-Caputo fractional derivative, J. Funct. Space., 2021 (2021), 1–14. https://doi.org/10.1155/2021/7117064 doi: 10.1155/2021/7117064
    [15] M. K. Alaoui, K. Nonlaopon, A. M. Zidan, A. Khan, Analytical investigation of fractional-order cahn-hilliard and gardner equations using two novel techniques, Mathematics, 10 (2022), 1643. https://doi.org/10.3390/math10101643 doi: 10.3390/math10101643
    [16] N. J. Ford, J. Xiao, Y. Yan, A finite element method for time fractional partial differential equations, Fract. Calc. Appl. Anal., 14 (2011), 454–474. https://doi.org/10.2478/s13540-011-0028-2 doi: 10.2478/s13540-011-0028-2
    [17] D. Fulger, E. Scalas, G. Germano, Monte Carlo simulation of uncoupled continuous-time random walks yielding a stochastic solution of the space-time fractional diffusion equation, Phys. Rev. E, 77 (2008), 021122. https://doi.org/10.1103/PhysRevE.77.021122 doi: 10.1103/PhysRevE.77.021122
    [18] A. R. Alharbi, M. B. Almatrafi, Exact solitary wave and numerical solutions for geophysical KdV equation, J. King Saud Univ. Sci., 34 (2022), 102087. https://doi.org/10.1016/j.jksus.2022.102087 doi: 10.1016/j.jksus.2022.102087
    [19] M. B. Almatrafi, A. Alharbi, K. Lotfy, A. A. El-Bary, Exact and numerical solutions for the GBBM equation using an adaptive moving mesh method, Alex. Eng. J., 60 (2021), 4441–4450. https://doi.org/10.1016/j.aej.2021.03.023 doi: 10.1016/j.aej.2021.03.023
    [20] A. R. Alharbi, M. B. Almatrafi, M. A. Abdelrahman, Analytical and numerical investigation for Kadomtsev-Petviashvili equation arising in plasma physics, Phys. Scripta, 95 (2020), 045215. https://doi.org/10.1088/1402-4896/ab6ce4 doi: 10.1088/1402-4896/ab6ce4
    [21] M. A. Abdelrahman, M. B. Almatrafi, A. Alharbi, Fundamental solutions for the coupled KdV system and its stability, Symmetry, 12 (2020), 429. https://doi.org/10.3390/sym12030429 doi: 10.3390/sym12030429
    [22] A. R. Alharbi, M. B. Almatrafi, Analytical and numerical solutions for the variant Boussinseq equations, J. Taibah Univ. Sci., 14 (2020), 454–462. https://doi.org/10.1080/16583655.2020.1746575 doi: 10.1080/16583655.2020.1746575
    [23] A. M. Zidan, A. Khan, M. K. Alaoui, W. Weera, Evaluation of time-fractional Fishers equations with the help of analytical methods, AIMS Math., 7 (2022), 18746–18766. https://doi.org/10.3934/math.20221031 doi: 10.3934/math.20221031
    [24] K. Xu, Y. Guo, Y. Liu, X. Deng, Q. Chen, Z. Ma, 60-GHz compact dual-mode on-chip bandpass filter using GaAs technology, IEEE Electr. Device L., 42 (2021), 1120–1123. https://doi.org/10.1109/LED.2021.3091277 doi: 10.1109/LED.2021.3091277
    [25] Z. Li, K. Wang, W. Li, S. Yan, F. Chen, S. Peng, Analysis of surface pressure pulsation characteristics of centrifugal pump magnetic liquid sealing film, Front. Energy Res., 10 (2022), 937299. https://doi.org/10.3389/fenrg.2022.937299 doi: 10.3389/fenrg.2022.937299
    [26] H. Chen, W. Chen, X. Liu, X. Liu, Establishing the first hidden-charm pentaquark with strangeness, Eur. Phys. J. C, 81 (2021), 409. https://doi.org/10.1140/epjc/s10052-021-09196-4 doi: 10.1140/epjc/s10052-021-09196-4
    [27] M. Ayata, O. Ozkan, A new application of conformable Laplace decomposition method for fractional Newell-Whitehead-Segel equation, AIMS Math., 5 (2020), 7402–7412. https://doi.org/10.3934/math.2020474 doi: 10.3934/math.2020474
    [28] H. Khan, D. Baleanu, P. Kumam, J. F. Al-Zaidy, Families of travelling waves solutions for fractional-order extended shallow water wave equations, using an innovative analytical method, IEEE Access, 7 (2019), 107523–107532. https://doi.org/10.1109/ACCESS.2019.2933188 doi: 10.1109/ACCESS.2019.2933188
    [29] B. Q. Li, Y. L. Ma, Rich soliton structures for the Kraenkel-Manna-Merle (KMM) system in ferromagnetic materials, J. Supercond. Nov. Magn., 31 (2018), 1773–1778. https://doi.org/10.1007/s10948-017-4406-9 doi: 10.1007/s10948-017-4406-9
    [30] B. Q. Li, Y. L. Ma, The non-traveling wave solutions and novel fractal soliton for the (2+1)-dimensional Broer–Kaup equations with variable coefficients, Commun. Nonlinear Sci., 16 (2011), 144–149. https://doi.org/10.1016/j.cnsns.2010.02.011 doi: 10.1016/j.cnsns.2010.02.011
    [31] S. Mukhtar, S. Noor, The numerical investigation of a fractional-order multi-dimensional Model of Navier-Stokes equation via novel techniques, Symmetry, 14 (2022), 1102. https://doi.org/10.3390/sym14061102 doi: 10.3390/sym14061102
    [32] S. Duran, Exact solutions for time-fractional Ramani and Jimbo-Miwa equations by direct algebraic method, Adv. Sci. Eng. Medicine, 12 (2020), 982–988. https://doi.org/10.1166/asem.2020.2663 doi: 10.1166/asem.2020.2663
    [33] R. Y. Molliq, M. S. M. Noorani, I. Hashim, R. R. Ahmad, Approximate solutions of fractional Zakharov-Kuznetsov equations by VIM, J. Comput. Appl. Math., 233 (2009), 103–108. https://doi.org/10.1016/j.cam.2009.03.010 doi: 10.1016/j.cam.2009.03.010
    [34] B. Zheng, Exp-function method for solving fractional partial differential equations, The Scientific World J., 2013 (2013). https://doi.org/10.1186/1687-1847-2013-199 doi: 10.1186/1687-1847-2013-199
    [35] Y. L. Ma, B. Q. Li, Y. Y. Fu, A series of the solutions for the Heisenberg ferromagnetic spin chain equation, Math. Method. Appl. Sci., 41 (2018), 3316–3322. https://doi.org/10.1002/mma.4818 doi: 10.1002/mma.4818
    [36] Y. L. Ma, B. Q. Li, C. Wang, A series of abundant exact travelling wave solutions for a modified generalized Vakhnenko equation using auxiliary equation method, Appl. Math. Comput., 211 (2009), 102–107. https://doi.org/10.1016/j.amc.2009.01.036 doi: 10.1016/j.amc.2009.01.036
    [37] B. Q. Li, Y. L. Ma, Periodic solutions and solitons to two complex short pulse (CSP) equations in optical fiber, Optik, 144 (2017), 149–155. https://doi.org/10.1016/j.ijleo.2017.06.114 doi: 10.1016/j.ijleo.2017.06.114
    [38] M. Zhang, Y. L. Ma, B. Q. Li, Novel loop-like solitons for the generalized Vakhnenko equation, Chinese Phys. B, 22 (2013), 030511. https://doi.org/10.1088/1674-1056/22/3/030511 doi: 10.1088/1674-1056/22/3/030511
    [39] Y. L. Ma, B. Q. Li, Novel optical soliton structures for a defocusing Lakshmanan-Porsezian-Daniel optical system, Optik, 284 (2023), 170931. https://doi.org/10.1016/j.ijleo.2023.170931 doi: 10.1016/j.ijleo.2023.170931
    [40] B. Q. Li, Y. L. Ma, Hybrid soliton and breather waves, solution molecules and breather molecules of a (3+1)-dimensional Geng equation in shallow water waves, Phys. Lett. A, 463 (2023), 128672. https://doi.org/10.1016/j.physleta.2023.128672 doi: 10.1016/j.physleta.2023.128672
    [41] H. Yasmin, N. H. Aljahdaly, A. M. Saeed, R. Shah, Probing families of optical soliton solutions in fractional perturbed Radhakrishnan-Kundu-Lakshmanan model with improved versions of extended direct algebraic method, Fractal Fract., 7 (2023), 512. https://doi.org/10.3390/fractalfract7070512 doi: 10.3390/fractalfract7070512
    [42] Q. Al-Mdallal, K. A. Abro, I. Khan, Analytical solutions of fractional Walter's B fluid with applications, Complexity, 2018 (2018), 1–10. https://doi.org/10.1155/2018/8131329 doi: 10.1155/2018/8131329
    [43] H. Jiang, F. Liu, I. Turner, K. Burrage, Analytical solutions for the multi-term time-fractional diffusion-wave/diffusion equations in a finite domain, Comput. Math. Appl., 64 (2012), 3377–3388. https://doi.org/10.1016/j.camwa.2012.02.042 doi: 10.1016/j.camwa.2012.02.042
    [44] B. Zheng, Exp-function method for solving fractional partial differential equations, The Scientific World J., 2013 (2013). https://doi.org/10.1186/1687-1847-2013-199 doi: 10.1186/1687-1847-2013-199
    [45] H. Khan, S. Barak, P. Kumam, M. Arif, Analytical solutions of fractional Klein-Gordon and gas dynamics equations, via the (G'/G)-expansion method, Symmetry, 11 (2019), 566. https://doi.org/10.3390/sym11040566 doi: 10.3390/sym11040566
    [46] Y. Tian, J. Liu, Direct algebraic method for solving fractional Fokas equation, Thermal Sci., 25 (2021), 2235–2244. https://doi.org/10.2298/TSCI200306111T doi: 10.2298/TSCI200306111T
    [47] H. Gnerhan, F. S. Khodadad, H. Rezazadeh, M. M. Khater, Exact optical solutions of the (2+1) dimensions Kundu-Mukherjee-Naskar model via the new extended direct algebraic method, Mod. Phys. Lett. B, 34 (2020), 2050225. https://doi.org/10.1142/S0217984920502255 doi: 10.1142/S0217984920502255
    [48] M. Younis, M. Iftikhar, Computational examples of a class of fractional order nonlinear evolution equations using modified extended direct algebraic method, J. Comput. Methods Sci., 15 (2015), 359–365. https://doi.org/10.3233/JCM-150548 doi: 10.3233/JCM-150548
    [49] H. Yasmin, N. H. Aljahdaly, A. M. Saeed, R. Shah, Investigating symmetric soliton solutions for the fractional coupled Konno-Onno system using improved versions of a novel analytical technique, Mathematics, 11 (2023), 2686. https://doi.org/10.3390/math11122686 doi: 10.3390/math11122686
    [50] H. Yasmin, N. H. Aljahdaly, A. M. Saeed, R. Shah, Investigating families of soliton solutions for the complex structured coupled fractional Biswas-Arshed model in birefringent fibers using a novel analytical technique, Fractal Fract., 7 (2023), 491. https://doi.org/10.3390/fractalfract7070491 doi: 10.3390/fractalfract7070491
    [51] T. A. A. Ali, Z. Xiao, H. Jiang, B. Li, A class of digital integrators based on trigonometric quadrature rules, IEEE T. Ind. Electron., 2023. https://doi.org/10.1109/TIE.2023.3290247 doi: 10.1109/TIE.2023.3290247
    [52] Y. Fang, H. Min, X. Wu, W. Wang, X. Zhao, G. Mao, On-ramp merging strategies of connected and automated vehicles considering communication delay, IEEE Intel. Transp. Sy., 23 (2022), 15298–15312. https://doi.org/10.1109/TITS.2022.3140219 doi: 10.1109/TITS.2022.3140219
    [53] W. Kuang, H. Wang, X. Li, J. Zhang, Q. Zhou, Y. Zhao, Application of the thermodynamic extremal principle to diffusion-controlled phase transformations in Fe-C-X alloys: Modeling and applications, Acta Mater., 159, (2018), 16–30. https://doi.org/10.1016/j.actamat.2018.08.008 doi: 10.1016/j.actamat.2018.08.008
    [54] A. M. Wazwaz, New solitary wave solutions to the modified forms of Degasperis-Procesi and Camassa-Holm equations, Appl. Mathe. Comput., 186 (2007), 130–141. https://doi.org/10.1016/j.amc.2006.07.092 doi: 10.1016/j.amc.2006.07.092
    [55] D. Li, S. S. Ge, T. H. Lee, Fixed-Time-Synchronized consensus control of multiagent systems, IEEE T. Control Netw., 8 (2021), 89–98. https://doi.org/10.1109/TCNS.2020.3034523 doi: 10.1109/TCNS.2020.3034523
    [56] Q. Meng, Q. Ma, Y. Shi, Adaptive Fixed-time stabilization for a class of uncertain nonlinear systems, IEEE T. Automat. Contr., 2023 (2023). https://doi.org/10.1109/TAC.2023.3244151 doi: 10.1109/TAC.2023.3244151
    [57] H. Yuan, B. Yang, System dynamics approach for evaluating the interconnection performance of Cross-Border transport infrastructure, J. Manag. Eng., 38 (2022). https://doi.org/10.1061/(ASCE)ME.1943-5479.0001015 doi: 10.1061/(ASCE)ME.1943-5479.0001015
    [58] V. P. Dubey, R. Kumar, J. Singh, D. Kumar, An efficient computational technique for time-fractional modified Degasperis-Procesi equation arising in propagation of nonlinear dispersive waves, J. Ocean Eng. Sci., 6 (2021), 30–39. https://doi.org/10.1016/j.joes.2020.04.006 doi: 10.1016/j.joes.2020.04.006
    [59] K. Zhang, A. S. Alshehry, N. H. Aljahdaly, R. Shah, N. A. Shah, M. R. Ali, Efficient computational approaches for fractional-order Degasperis-Procesi and Camassa-Holm equations, Results Phys., 50 (2023), 106549. https://doi.org/10.1016/j.rinp.2023.106549 doi: 10.1016/j.rinp.2023.106549
    [60] S. Das, R. Kumar, Approximate analytical solutions of fractional gas dynamic equations, Appl. Math. Comput., 217 (2011), 9905–9915. https://doi.org/10.1016/j.amc.2011.03.144 doi: 10.1016/j.amc.2011.03.144
    [61] M. A. Dokuyucu, Caputo and Atangana-Baleanu-Caputo fractional derivative applied to garden equation, Turk. J. Sci., 5 (2020), 1–7.
    [62] S. M. M. Alizamini, H. Rezazadeh, K. Srinivasa, A. Bekir, New closed form solutions of the new coupled Konno-Oono equation using the new extended direct algebraic method, Pramana, 94 (2020), 1–12. https://doi.org/10.1007/s12043-020-1921-1 doi: 10.1007/s12043-020-1921-1
  • This article has been cited by:

    1. Parvaiz Ahmad Naik, Muhammad Amer, Rizwan Ahmed, Sania Qureshi, Zhengxin Huang, Stability and bifurcation analysis of a discrete predator-prey system of Ricker type with refuge effect, 2024, 21, 1551-0018, 4554, 10.3934/mbe.2024201
    2. Weirun Tao, Zhi-An Wang, Global well-posedness and Turing–Hopf bifurcation of prey-taxis systems with hunting cooperation, 2025, 0956-7925, 1, 10.1017/S0956792525000026
    3. Alyaa Hussain Naser, Dahlia Khaled Bahlool, Dynamics of an Intraguild Predation Food Web Cooperation Model Under the Influence of Fear and Hunting, 2025, 13, 2079-3197, 128, 10.3390/computation13060128
  • Reader Comments
  • © 2023 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(1318) PDF downloads(38) Cited by(1)

Figures and Tables

Figures(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog