Complex Dynamic Behaviour of Food Web Model with Generalized Fractional Operator

: We apply a new generalized Caputo operator to investigate the dynamical behaviour of the non-integer food web model (FWM). This dynamical model has three population species and is nonlinear. Three types of species are considered in this population: prey species, intermediate predators, and top predators, and the top predators are also divided into mature and immature predators. We calculated the uniqueness and existence of the solutions applying the ﬁxed-point hypothesis. Our study examines the possibility of obtaining new dynamical phase portraits with the new generalized Caputo operator and demonstrates the portraits for several values of fractional order. A generalized predictor–corrector (P-C) approach is utilized in numerically solving this food web model. In the case of the nonlinear equations system, the effectiveness of the used scheme is highly evident and easy to implement. In addition, stability analysis was conducted for this numerical scheme.


Introduction
Energy and materials follow one path between species in a food chain model, whereas food webs are more complex because they connect many food chains. Different trophic levels are found in a food web. There are various categories of organisms within the trophic levels, including producers, consumers, and decomposers. The structure of a food web is typically represented by a lattice arrangement. Using a system of differential equations, it is possible to design food chains and food webs. Food chains, in ecology, are a chain of organisms feeding on the organism next to them, while food webs are a collection of food chains joined together. It has been of interest to several researchers to analyse the dynamical behaviour of the food chain model and the web model [1][2][3][4]. A modular food web theory, which studies the structural and functional properties of low-species-based food webs, aims to determine how the structure and interactions mediate ecosystem stability [5,6]. Several species in nature have life cycles that are divided into at least two stages: mature and immature. These stages have different characteristics. Food web models (FWMs) depicting a stage structure have been extensively studied [7,8]. The impact of cannibalism on ecological systems has been studied extensively over the past few decades. Aquatic, as well as terrestrial food webs have cannibalistic populations. This subject has been addressed by several studies [9][10][11][12]. Stage-structured populations frequently engage in cannibalism, whether in the wild or in watery food webs. The cannibalism model was examined and investigated by Diekmann et al. [13]. A watery food chain in which a predator cannibalizes was studied by Bhattacharyya and Pal [14]. The dynamics of the system are therefore influenced by cannibalism in a very significant way. Fishes, birds, mammals, and others are among the animals that have cannibalistic natures.
There are over 300 years of development behind fractional calculus, and today, this is still an important concept of studying real-world problems [15][16][17][18][19][20]. The literature of fractional calculus has introduced a variety of fractional derivatives, including Caputo [21], Atangana-Baleanu [22], and Caputo-Fabrizio [23], which are the most widely used derivatives. Fractional differential equations can describe dynamic processes within biological and ecological systems with a higher degree of accuracy and reliability since most biological and ecological mathematical models have long-term memories. An understanding of fractional species systems can provide new possibilities for describing the dynamic behaviours of multi-species food web ecosystems, given the complexity and existence of nonlinear effects [24]. In addition, fractional-order forms have a number of advantages, such as a meticulous illustration and an accurate interpretation of operation rules. In order to further explore the dynamics of systems with competition, predation, and parasitism, classical integer differential equations of ecosystems are replaced with fractional differential equations. The literature contains a variety of nonlocal operators, which are used extensively in applied mathematics. An integral and fractional derivative introduced by Katugampola in 2014 generalizes both Riemann-Liouville and Hadamard integrals and derivatives [25,26]. The generalized Caputo operator was recently constructed by Odibat et al. [27]. In the literature, the generalized Caputo operator has been applied in various ways. A recent study by Rubayyi T. Alqahtani et al. [28] utilized a generalized Caputo operator to model bioethanol production. The new generalized Caputo operator was used to analyse the COVID-19 model [29]. In the paper [30], the author investigated irregular meshes with finite difference methods to determine the error estimates when the Caputo operator of the solution of the FDEs has a low smoothness. The paper [31] developed the asymptotic expansion formula for the trapezoidal approximation of the fractional integral, and the author applied the expansion formula to calculate approximations for fractional integrals of orders α, α + 1, α + 2, α + 3, and α + 4.
In this paper, we extend the classical integer-order food web model (FWM) to a noninteger food web model through a generalized Caputo operator. Moreover, we discuss a generalized predictor-corrector numerical solution that is a generalization of the P-C numerical scheme [32,33] to study the complexity of the food web model's behaviour, and we analyse the stability of this scheme. With the generalized Caputo operator, a nonuniform grid is used in the P-C scheme instead of the uniform grid in the Caputo operator. ϕ and ρ are the only parameters needed to generalize the Caputo integral operator, which provides a great deal of theoretical and numerical equipment for fractional mathematical modelling. There are numerous applications of this P-C technique in various fields of FDEs. In this study, we analysed the behaviour of the food web model on various different fractional orders and on another parameter of the derivative, which gives us different dynamic phase diagrams of this food web model. Due to Caputo derivatives describing better certain physical problems involving memory effects, we defined the generalized fractional derivatives in a Caputo form. This Caputo version of generalized fractional derivatives would prove useful for researchers interested in describing real-world phenomena using fractional operators. Finally, by noticing that lim ρ→0 x lead to Hadamard and Caputo-Hadamard results, we saw that the limiting case as ρ → 0 leads to those results. Furthermore, when ρ → 1, the fractional derivatives of Riemann-Liouville and of Caputo were obtained.

Outline of the Paper
We divided this whole work into the following sections: Section 2 represents the description of the food web model. The preliminaries of fractional calculus (FC) are covered in Section 3. The existence of solutions is demonstrated in Section 4. Solution uniqueness is demonstrated in Section 5. Section 6 presents generalized predictor-corrector numerical algorithms for fractional-order food web models using the generalized Caputo operator. Section 7 contains simulations and discussions of the numerical results. A conclusion is given in Section 8.

Description of the Food Web Model
Our study proposes and analyses a three-species FWM that includes cannibalism and a stage structure within top predator species. Predators at the top are generally divided into immature and mature stages. Initial stage individuals are unable to hunt or reproduce, as they are dependent on their mature parents for survival. Additionally, we constructed an ecological model that includes stage cannibalism and structure in the top predators as part of a three-species food web model. A food web model can be constructed using the above considerations [34].
Prey density (lower level species) at time t is denoted by x(t); intermediate predator density (middle-level species) at time t is denoted by y(t); top predator density (mature and immature of higher-level species) at time t is denoted by z(t), u(t). With an intrinsic growth rate r and carrying capacity H, the prey grows logistically. Based on the Lotka-Volterra functional response, the intermediate predator consumes the prey at the lowest level, with an attack rate a 1 and a conversion rate e 1 . In the absence of their food source, it continues to decay exponentially as a result of natural mortality rate d 1 . There are two kinds of top predators: mature and immature. Immature populations are assumed to grow exponentially along with their parents denoted by the mature population with growth rate b, while a part grows up to become a mature population with growth rate c. Additionally, both the mature and immature populations face natural death with mortality rates of d 2 and d 3 , respectively. With maximum attack rate a 3 and conversion rate e 2 , the mature top predator attacks the intermediate predator using the Lotka-Volterra response functional. When the availability of their preferred food becomes rare, they cannibalise the immature top predator based on the Lotka-Volterra functional response with maximum attack rate a 3 and conversion rate e 3 .

Equilibrium Points of the Food Web Model
In the part of this section, we calculate the equilibrium points corresponding to the food web model. The steady-state conditions for the model are as follows: is the interior equilibrium point. The characteristic equation forẍ is as follows: where γ 1 = −a 2 1 a 2 a 3 e 2 1 e 2 r < 0.
Thus, a simple computation shows that this exists only if and only if the following is true: with one set condition (γ 2 < 0 and γ 4 < 0) or (γ 3 > 0 and γ 4 > 0).

Preliminaries
Definition 2. Consider f (t) ∈ H 1 to be a differentiable function in the interval (a, b), a < b, and ϕ ∈ [0, 1], then we define the Caputo non-classical operator as Gamma functions are represented by Γ(.). Here is the definition of the gamma function: Definition 3. Here, the order of derivative ϕ > 0 and ρ > 0, and the generalized non-classical Definition 4. Here, the order of derivative ϕ > 0 (m − 1 < ϕ < m) and ρ > 0. For a function f(t), the generalized Riemann-type non-classical derivative GRL D ϕ,ρ Definition 5. Here, the order of derivative ϕ > 0, ρ > 0 and m = ϕ . For a function f(t), the generalized Caputo-type non-classical derivative GC D ϕ,ρ a+ is defined as Definition 6. Here, the order of derivative ϕ > 0 (m − 1 < ϕ < m) and ρ > 0. For a function f(t), the generalized Caputo-type non-classical derivative GC D ϕ,ρ The relation between the Riemann-Liouville and the generalized non-classical integral from the substitution χ ρ → χ is as follows When the lower limit is zero a = 0, the relation is

Existence of Solutions
The fixed-point assumption is used to investigate the existence of a solution for the fractional food web mathematical model. Now, a non-integer food web mathematical model can be described as follows.
The initial conditions of a mathematical model of a food web are as follows: Using the generalized Caputo-type non-classical integral, we have In order to simplify, we determine

y(t) is the kernel and x(t)
and x 1 (t) are two function, we can find By , and u(t) ≤ q 4 are the bounded functions; furthermore, we have Therefore, the Lipschitz condition holds for 1 if the inequality 0 ≤ 1 < 1 is the contraction of 1 . As we apply the same procedure to kernels 2 , 3 , and 4 , the following results emerge: Kernels 1 , 2 , 3 , and 4 are determined by Equation (16). Afterwards, we determine the associated integrals: furthermore, we obtain and the initial condition is When we subtract consecutive terms, we obtain Take the following Equation (23) is found using the triangular and norm properties.
and under the Lipschitz condition, the Kernels will exhibit the following outcomes: Therefore, we obtain the following: We obtain the same results for ∆ n (t), ℵ n (t), and ‫ג‬ n (t) when we follow the same procedure: Following the above conclusion, we can prove the new theorem.
Theorem 2. The generalized Caputo-type non-classical-order food web mathematical model has a unique solution if t max fulfills the following criteria.
Proof of Theorem 2. By assuming that x(t), y(t), z(t), and u(t) are bounded functions and considering that we have already shown that the kernels possess the Lipschitz condition, the following relation is then given: Considering that all the above functions exist and are smooth, we prove that these functions are the solution to the food web mathematical model. Hence, we assume When n → 0 is taken as the limit in Equation (31), we obtain A recursive process leads to the following equation: then, for t max , we obtain: We can obtain A n (t) → 0 at n → ∞ by taking the limits of both sides of the above equation, and B n (t) → 0, C n (t) → 0 and D n (t) → 0 can also be obtained by taking the limits of both sides. Therefore, the proof is complete.

Find the Uniqueness of the Solution
In this segment of the food web mathematical model, unique solutions are presented. Consider x 1 (t), y 1 (t), z 1 (t), and u 1 (t) to be the other solutions of the proposed system, then we have When the norm is applied to each side of Equation (35), the following result is obtained: The Lipschitz condition applied to the kernel yields In addition, we obtain the following: According to the above, the first differential equation of the financial model has a unique solution. Similarly, we show that y(t), z(t), and u(t) have unique solutions.

Generalized Predictor-Corrector Technique
We converted the model into a fractional Volterra type in order to obtain numerical solutions. We propose a P-C scheme with a generalized Caputo operator to solve the food web system.
Consider the Volterra integral form of the first equation of the food web system: we can write the above equation as follows: In order to simplify, we write A 1 (ς, x(ς)) instead of A 1 (ς, x, y, z, u). The interval [0, T] is divided into N subintervals {[t r , t r+1 ], r = 0, 1, 2, . . . , N − 1} with the mesh points as follows: Here, h = T ρ N . The approximate solution x n+1 ≈ x(t n+1 ) of Equation (41) can be calculated as follows: Let k = ς ρ ; therefore, the above equation becomes: The integral can now be discretized as follows: Using the trapezoidal rule, the right-hand side of (45) is evaluated relative to the weight function (t ρ n+1 − k) ϕ−1 . We can replace A 1 (k 1 ρ , x(k 1 ρ )) with its piecewise linear interpolant by choosing nodes at t ρ r (r = 0, 1, 2, . . . , n + 1). Then, we have: The corrector expression for x(t n+1 ), n = 0, 1, 2, . . . , N − 1 is as follows if the above term is substituted into (45): where Using the Adams-Bashforth method, we determine the predictor value x p (t n+1 ) for integral (44). We replace A 1 (k 1 ρ , x(k 1 ρ )) with A 1 (t r , x(t r )) at each integral in Equation (45) to obtain the following: Thus, we can conclude that: We now approximate x(t n+1 ) ≈ x n+1 to develop the P-C algorithm by replacing x p (t n+1 ) with x(t n+1 ) in Equation (47), as follows: We developed a P-C scheme specified in (50) and (51). In this case, the P-C algorithm for the whole model can be written as follows: here, h = T ρ N and x p n+1 , y p n+1 , z p n+1 , and u p n+1 are defined as follows: where A 1 , A 2 and A 3 are defined as follows: A 3 (t, z) = a 2 e 2 y(t)z(t) + a 3 e 3 z(t)u(t) + cu(t) − d 2 z(t), Remark 1. The comparison of our adaptive P-C formula with that of [12], based on the product integration methods described in [30], shows that the error should behave in this way: where p = min{2, 1 + ϕ}.

Numerical Results and Discussion
For the numerical simulation of the food web mathematical model, we propose a predictor-corrector (P-C) algorithm involving a generalized Caputo operator. Our numerical solution for the food web generalized Caputo derivative model illustrates the applicability and efficiency of the proposed algorithm. MATLAB was used to perform the simulations. The proposed algorithms should be beneficial for the simulation of non-integer models. The dynamical behaviours of the food web model were examined in our analysis. We considered the following parameter values and initial values in Table 1. Generalized Caputo-type fractional derivatives also possess the same properties as Caputo-type derivatives. In order to solve the fractional IVP efficiently, consistently, and accurately, the predictor-corrector (PC) scheme is one of the best available. We solved the projected model using the modified PC scheme in the current study. According to the generalized Caputo algorithm, adaptive PC schemes use a nonuniform grid, which differs from the derivative Caputo PC algorithm. In fractional calculus applications, the generalized fractional integral operator is a valuable tool for controlling and building mathematical models due to the effect of its parameters ϕ and ρ. This new generalized Caputo fractional derivative has extra features over the other fractional derivatives such as Caputo, Caputo-Fabrizio, and Atangana-Baleanu. There is another parameter ρ that is very helpful to graphical simulations when it comes to true data, in addition to the fractional-order parameter ϕ. Changing the parameter value ρ allows us to see more kinds of graphs. Figures 1 and 2 illustrate the three-dimensional and two-dimensional dynamic phase portrait of the fractional food web system with the generalized Caputo derivative, respectively, when ϕ = 1 and ρ = 1.1. Figure 3 exhibits the state variables x(t), y(t) z(t), and u(t) of the proposed model when ϕ = 1 and ρ = 1.1. It can be seen that the value of ρ strongly influences the characteristics of the fractional derivative, and this provides a different way of approaching control applications. Figures 4 and 5 illustrate the three-dimensional and two-dimensional dynamic phase portrait of the fractional food web system with the generalized Caputo derivative, respectively, when ϕ = 1 and ρ = 1.1. Figure 6 exhibits the state variables x(t), y(t) z(t), and u(t) of the proposed model when ϕ = 1 and ρ = 1.2. At fixed ρ, depending on the fractional-order value, our fractional food web system displays the complexity of the chaotic phase portrait. Hence, Figures 7-9 are the graphical illustrations of the proposed system at different fractional-order ϕ = 1, 0.95, 0.90 and fixed ρ. Further, we took the different values of ρ and ϕ to be fixed, then the fractional food web system exhibits different phase portraits, which are shown in Figures 10-12. During the simulation of models with two fractional parameters, we observed chaos, and we noticed that the dynamics became more complex. Intermediate predator y(t) 4 Prey population x(t) Immature predator u(t) 1.5

Conclusions
Here, we examined a three-species food web model. According to this model, top predators are stage-structured, with a mature predator having a cannibalism trait. In the absence of the predator, the prey grows logistically at the first level. Food consumption at different levels of the food web is described by Lotka-Volterra functional responses. The proposed mathematical model of the food web was examined using the generalized Caputo fractional derivative. Using a fixed-point hypothesis, this study presented an investigation of the existence and uniqueness of the fractional food web system. The algorithm described in this study is based on a numerical technique called 'predictorcorrector', which allows the approximate solution of the fractional food web model to be found. We demonstrated the stability of this numerical method. The fractional food web model was geometrically presented under the generalized Caputo operator for different choices of ϕ and ρ. The new dynamical behaviour and phase portrait were demonstrated for various fractional orders (ϕ) and the value of ρ. This graphical illustration showed how the order of derivatives and the system parameters greatly affect the system. The new generalized fractional derivative will be used in future efforts to model other biological systems with memory or with hereditary properties, as well as to identify other important properties of this new generalized derivative.

Conflicts of Interest:
The authors declare no conflict of interest.