REVIEW 4 major objections 4 minor 31 references
Three-species predator-prey model with respect to Caputo and Caputo-Fabrizio fractional operators
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper claims that the type of fractional differential operator—Caputo versus Caputo-Fabrizio—and the value of the fractional order can flip which equilibria of a three-species predator-prey system are stable, and it supplies a…
desk verdict A well-intentioned but mathematically broken numerical scheme invalidates the central claims, though the stability comparison and examples are also shaky. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the Caputo-Fabrizio fractional derivative, $\mathrm{CF}D_0^\alpha f(t)=\frac{M(\alpha)}{1-\alpha}\int_0^t \exp(-\frac{\alpha}{1-\alpha}(t-\tau))f'(\tau)\,d\tau$, whose exponential kernel removes the singularity of the Caputo kernel $(t-\tau)^{-\alpha}$. Around this operator the paper builds two tools: the linear stability criterion for CF systems, which replaces the Caputo wedge with a bounded circular boundary in the complex plane, and a corrected predictor-corrector (Adams-Bashforth-type) scheme with weights $b_{i,k+1}$ and $d_{i,k+1}$ derived from the equivalent integral equation (14). The scheme is what allows CF trajectories to be simulated; the eigenvalue criterion decides which equilibrium those trajectories should approach.
What would settle it
Run the proposed scheme on a scalar linear CF equation with a known exact solution at decreasing step sizes $h$; if the numerical solution does not converge to the direct integral solution of the CF initial-value problem, the claimed stability gain is an artifact of equation (14).
Extended reading notes
Core claim
On the paper's own terms, the discovery is that the Caputo-Fabrizio derivative produces a different stability geometry for fractional Lotka-Volterra systems. For the linearized system $D^\alpha u = Au$, the Caputo stability condition is $|\arg(\mathrm{spec}(A))| > \alpha\pi/2$, an unbounded wedge in the complex plane; the Caputo-Fabrizio condition, taken from Theorem 2.2, is the disjunction $|\lambda| \geq 1/(1-\alpha)$ (with $\lambda\neq 1/(1-\alpha)$), $\operatorname{Re}\lambda > 1/(1-\alpha)$, $\operatorname{Re}\lambda < 0$, or $|\operatorname{Im}\lambda| > 1/(2(1-\alpha))$, and Fig. 1 draws the unstable region as a bounded circle centered at $(0, 1/(2(1-\alpha)))$ with radius $1/(2(1-\alpha))$. The paper classifies eigenvalues into four classes—stable under both operators, stable only under Caputo, unstable under both, and stable only under Caputo-Fabrizio—and constructs a corrected Adams-Bashforth predictor-corrector method for the CF operator. Simulating the three-species model with this scheme reproduces the theoretical stability classification at each equilibrium. The central message is that operator type and order jointly determine which equilibrium the system approaches.
Load-bearing premise
The whole numerical method rests on equation (14), the paper's integral-equation form of a Caputo-Fabrizio initial-value problem; for the case $0<\alpha<1$ used in the examples, that formula invokes an undefined factorial and does not match the known integral representation of the Caputo-Fabrizio derivative.
Editorial extensions
If this is right
- For the same three-species interaction parameters, switching from the Caputo to the Caputo-Fabrizio operator can turn a stable equilibrium into an unstable one, so the choice of derivative is a substantive modeling decision.
- The value of the fractional order $\alpha$ can determine whether a given equilibrium is asymptotically stable; in the reported examples, lowering $\alpha$ expands the stability region and can allow the system to settle at equilibria that are unstable at higher $\alpha$.
- A reliable numerical study of Caputo-Fabrizio Lotka-Volterra systems needs the corrected predictor-corrector scheme rather than the earlier CF scheme, which the paper identifies as flawed.
- The four-class eigenvalue classification supplies a practical rule before simulating: locate the Jacobian eigenvalues to see whether Caputo, CF, both, or neither will stabilize the system.
- Because local stability of a nonlinear system is determined by its Jacobian eigenvalues, the same stability conditions and the same corrected scheme apply to other polynomial-type CF systems, not only to the three-species Lotka-Volterra model.
Reading between the lines
- The bounded-circle geometry implies that eigenvalues with large imaginary parts can be stable under CF while unstable under Caputo, so the operator could be chosen according to the spectral character of the linearization.
- The examples' sensitivity to initial conditions suggests that domains of attraction, not only local stability, decide the fate of trajectories; mapping those domains for each equilibrium and each operator would be a natural next test.
- Because the CF operator is interpreted in the paper as modeling distributed lag rather than memory, the stability differences could be probed with real predator-prey time series: if lag-based models fit better, CF would be the appropriate operator.
- For incommensurate orders (different $\alpha$ for prey and predators), the stability boundaries should be intersections of the individual regions, producing new equilibrium classifications; the paper lists incommensurate orders as a future direction.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies a three-species Lotka-Volterra system with Caputo and Caputo-Fabrizio (CF) fractional derivatives. It proposes a predictor-corrector Adams-Bashforth scheme for the CF equation, derives local stability conditions for the equilibria, and reports numerical comparisons that are claimed to show that the type of differential operator and the value of the fractional order significantly influence stability and dynamics. The paper's central numerical claim rests on Eq. (14), which is presented as the equivalent integral form of the CF initial value problem; the stability analysis and the simulations in Section 5 are built on that representation.
Significance. The motivation is reasonable: a correct, stable numerical method for CF Lotka-Volterra systems would be a useful contribution, and the paper correctly identifies that earlier CF predictor-corrector schemes [30] are flawed. The paper also attempts a direct comparison of Caputo and CF stability regions, which is an interesting and potentially useful goal. However, the central derivation is invalidated by a fundamental error in Eq. (14), so the proposed numerical method does not solve the stated model and the numerical comparisons in Section 5 are not evidence about the model's behavior. The analytical stability section also lacks a valid justification for applying linear eigenvalue conditions to nonlinear CF systems. The contribution, as it stands, cannot support the paper's claims.
major comments (4)
- [Sec. 4, Eq. (14)] Eq. (14) is not the equivalent integral form of the CF initial value problem for 0<α<1. For n=1 the first term contains (n-2)! = (-1)!, and the kernel (x-t)^{-1} is singular. Directly inverting Eq. (2) gives f(t) = f(0) + ((1-α)/M(α)) g(t) + (α/M(α)) ∫_0^t g(s) ds, which includes a non-integral (1-α)/M(α) g(t) term that Eq. (14) omits. Since the predictor-corrector scheme in Eqs. (15)-(18) is derived from Eq. (14), it solves a different equation, and all numerical results in Section 5 are therefore not simulations of the stated CF Lotka-Volterra model.
- [Sec. 3.2, Theorem 3.2] Theorem 3.2 asserts local asymptotic stability from the eigenvalues of the Jacobian J(ε*) with the proof described as 'straightforward with Theorem 2.2 and [26]'. No linearization theorem for nonlinear systems with the CF operator is stated or proved, and eigenvalue conditions for linear fractional systems do not automatically transfer to nonlinear equilibria. Consequently, the stability classifications in Table 1 and the interpretations of Examples 1-3 are not supported by the arguments given.
- [Sec. 3.2.5 and Fig. 1] The stability region depicted in Fig. 1 is inconsistent with Theorem 2.2 of the paper. The figure labels the disk centered at (0, 1/(2(1-α))) with radius 1/(2(1-α)) as the CF unstable region, but points in the left half-plane inside this disk satisfy Re(λ)<0, which Theorem 2.2 condition 3 declares asymptotically stable. In Sec. 3.2.5 the discriminant cases are also misstated: the trigonometric formulas (9)-(11) apply to Δ<0, not Δ>0, and Eq. (7) ends with -q/3 rather than -a/3.
- [Table 1] Table 1 contains repeated and inconsistent conditions for ε1. The Caputo row and the CF row both list 'a1a2 < a2a3 - a2' twice, whereas the text in Sec. 3.2.2 uses a1a4 and a1a6; the CF row also mixes 'a1a2 < ...' with the conditions 'a1a4 - a2a3/a2 > α/(1-α)' and 'a1a6 - a2a5/a2 > α/(1-α)'. As a result, the table cannot be used to check the asserted stability regions, and some of its entries are internally inconsistent.
minor comments (4)
- [Sec. 2.1] The normalization function M(α) is only constrained by M(0)=M(1)=1 and is never specified for other α; since it appears in Eq. (14) and in the discretization, the numerical simulations are not reproducible without stating the choice made.
- [Sec. 3.2.2] The sentence 'Since a1 > 0, 1 - a3 < 0, 1 - a4 < 0' appears to contain a typo: the third condition should likely be 1 - a5 < 0, since the eigenvalue is 1 - a5.
- [Sec. 4] The notation 'Adam-Bashforth' in the keywords and text is a misspelling of 'Adams-Bashforth', and the paper contains several other typographical errors such as 'asymptoialy' and 'stablitity'.
- [Sec. 5] The axis labels and legends in Fig. 2 (right) and Fig. 5 (right) are difficult to read; the three-dimensional trajectories would benefit from clearer labeling of which curve corresponds to which species.
Circularity Check
No significant circularity found; the paper's central claims are not derived from their own outputs, and the numerical method is an independent construction, even if equation (14) may be mathematically questionable.
full rationale
The paper's central claim is a new predictor-corrector scheme for the Caputo-Fabrizio operator, built on the asserted equivalent integral equation (14), and a stability comparison between Caputo and Caputo-Fabrizio operators. The stability criteria are imported from external references: Theorem 2.2 is attributed to Ref. [27] (Li, Cheng, Li, and Zhong), not to the present authors, and the Caputo criteria come from Matignon [26] and Ref. [28]. The numerical scheme adapts the Adams-Bashforth-multon method from Ref. [29] and claims to correct Ref. [30]; neither is the authors' own prior work. The self-citations in the introduction (Refs. [3]-[6]) are background examples of fractional calculus applications and are not used to justify the numerical method or the stability conclusions. There is no fitting of parameters to a subset of data and then calling the result a prediction, and no definition is constructed so that the target result is true by definition. Equation (14) may indeed be invalid for 0<alpha<1, since it contains (n-2)! for n=1 and omits the non-integral (1-alpha)/M(alpha) term of the true CF inverse; however, an incorrect mathematical equivalence is a correctness flaw, not circularity. The paper does not reduce its claimed result to its own assumptions or to a self-citation chain. Accordingly, the circularity score is 0.
Assumptions & free parameters
free parameters (1)
- Normalization function M(α) in the CF derivative =
unspecified; typically set to 1
assumptions (4)
- domain assumption The Caputo-Fabrizio derivative definition (Eq. 2) is a valid fractional operator and the linear stability theorem from [27] is correctly stated as Theorem 2.2.
- domain assumption Local stability of the nonlinear system (4) is determined by the eigenvalues of its Jacobian using the linear stability conditions.
- domain assumption The function F in (6) is Lipschitz on L[0,t'] with a uniform constant L.
- ad hoc to paper The equivalent integral formulation (14) correctly represents the CF initial value problem for 0<α<1.
Cite this review
Pith. "Pith review of Three-species predator-prey model with respect to Caputo and Caputo-Fabrizio fractional operators." pith.science (2026). https://pith.science/paper/PHXLA5CK
@misc{pith2026190803685,
author = {Pith},
title = {Pith review of: Three-species predator-prey model with respect to Caputo and Caputo-Fabrizio fractional operators},
year = {2026},
howpublished = {\url{https://pith.science/paper/PHXLA5CK}},
note = {Machine review of arXiv:1908.03685}
}
read the original abstract
We study distributed lag effects in three-dimensional Lotka-Volterra systems by applying the concept of fractional calculus. We derive a new numerical method that provides enhanced stability for the Caputo-Fabrizio operator based on Adams-Bashforth method, considering non-singular kernel in the definition of Caputo-Fabrizio operator. We investigate the stability conditions of this system with comparisons to the Caputo fractional derivative. Numerical results show that the type of differential operators and the value of orders significantly influence the stability of the numerical solution, and dynamics of the Lotka-Volterra system.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[30]
Y. Teng Toh, C. Phang, and J. R. Loh, “New predictor-correc tor scheme for solving nonlinear differential equations with caputo-fabrizio ope rator,” Mathematical Methods in the Applied Sciences , vol. 42, 2018
work page 2018
-
[26]
Stability results for fractional differential equa tions with applications to control processing,
D. Matignon, “Stability results for fractional differential equa tions with applications to control processing,” pp. 963–968, 1996
work page 1996
-
[1]
A new collec tion of real world applications of fractional calculus in science and engineering,
H. Sun, Y. Zhang, D. Baleanu, W. Chen, and Y. Chen, “A new collec tion of real world applications of fractional calculus in science and engineering,” Communications in Nonlinear Science and Numerical Simulation , vol. 64, pp. 213–231, 2018
work page 2018
-
[2]
Modeling some re al phenomena by fractional differential equations,
R. Almeida, N. R. Bastos, and M. T. T. Monteiro, “Modeling some re al phenomena by fractional differential equations,” Mathematical Methods in the Applied Sciences , vol. 39, no. 16, pp. 4846–4855, 2016
work page 2016
-
[3]
Mem- ory effects on epidemic evolution: The susceptible-infected-recov ered epidemic model,
M. Saeedian, M. Khalighi, N. Azimi-Tafreshi, G. R. Jafari, and M. Au sloos, “Mem- ory effects on epidemic evolution: The susceptible-infected-recov ered epidemic model,” Physical Review E , vol. 95, p. 022409, 2017
work page 2017
-
[4]
Fractional dynamics of network growth constrained by aging node interaction s,
H. Safdari, M. Zare Kamali, A. Shirazi, M. Khalighi, G. Jafari, and M. Ausloos, “Fractional dynamics of network growth constrained by aging node interaction s,” PLOS ONE , vol. 11, no. 5, pp. 1–13, 2016
work page 2016
-
[5]
S. Hosseinpour and A. Nazemi, “Solving fractional optimal contr ol prob- lems with fixed or free final states by haar wavelet collocation metho d,” IMA Journal of Mathematical Control and Information , vol. 33, no. 2, 2015. 18
work page 2015
-
[6]
Nonlinear fractiona l optimal control problems with neural network and dynamic optimization schemes,
S. Ghasemi, A. Nazemi, and S. Hosseinpour, “Nonlinear fractiona l optimal control problems with neural network and dynamic optimization schemes,” Nonlinear Dynamics , vol. 89, no. 4, 2017
work page 2017
Show all 31 references
-
[7]
The concepts and applications of frac tional order differential calculus in modeling of viscoelastic systems: A primer,
M. A. Matlob and Y. Jamali, “The concepts and applications of frac tional order differential calculus in modeling of viscoelastic systems: A primer,” Critical Reviews ™ in Biomedical Engineering, vol. 47, no. 4, 2019
2019
-
[8]
A new definition of fractional derivat ive without singular ker- nel,
M. Caputo and M. Fabrizio, “A new definition of fractional derivat ive without singular ker- nel,” Progress in Fractional Differentiation and Applications , vol. 1, no. 2, pp. 1–13, 2015
2015
-
[9]
Properties of a new fractional derivat ive without singular kernel,
J. Losada and J. Nieto, “Properties of a new fractional derivat ive without singular kernel,” Progress in Fractional Differentiation and Applications , vol. 1, pp. 87–92, 04 2015
2015
-
[10]
Applications of new time and spatial f ractional derivatives with exponential kernels,
M. Caputo and M. Fabrizio, “Applications of new time and spatial f ractional derivatives with exponential kernels,” Progress in Fractional Differentiation and Applications , vol. 2, pp. 1–11, 01 2016
2016
-
[11]
Controlling the wave movement o n the sur- face of shallow water with the caputo–fabrizio derivative with fract ional order,
A. Atangana and B. Alkahtani, “Controlling the wave movement o n the sur- face of shallow water with the caputo–fabrizio derivative with fract ional order,” Chaos Solitons and Fractals , pp. 1–7, 04 2016
2016
-
[12]
Cance r treatment model with the caputo-fabrizio fractional derivative,
M. Ali Dokuyucu, E. Celik, H. Bulut, and H. M. Baskonus, “Cance r treatment model with the caputo-fabrizio fractional derivative,” The European Physical Journal Plus , vol. 133, 03 2018
2018
-
[13]
Mathematica l analysis of hiv/aids infec- tion model with caputo-fabrizio fractional derivative,
S. Bushnaq, S. A. Khan, K. Shah, and G. Zaman, “Mathematica l analysis of hiv/aids infec- tion model with caputo-fabrizio fractional derivative,” Cogent Mathematics and Statistics , vol. 5, no. 1, p. 1432521, 2018
2018
-
[14]
A numerical frame - work for the approximate solution of fractional tumor-obesity mo del,
S. Arshad, D. Baleanu, O. Defterli, and S. , “A numerical frame - work for the approximate solution of fractional tumor-obesity mo del,” International Journal of Modeling, Simulation, and Scient ific Computing , vol. 10, 11 2018
2018
-
[15]
Caputo–fabrizio operator in terms of integer derivatives: memory or dis- tributed lag?,
V. E. Tarasov, “Caputo–fabrizio operator in terms of integer derivatives: memory or dis- tributed lag?,” Computational and Applied Mathematics , vol. 38, p. 113, May 2019. 19
2019
-
[16]
A mathematical model on fractional lotk a–volterra equations,
S. Das and P. Gupta, “A mathematical model on fractional lotk a–volterra equations,” Journal of Theoretical Biology , vol. 277, no. 1, pp. 1 – 6, 2011
2011
-
[17]
Banking system: Three level lotka-volterra m odel,
C. Calin-Adrian, “Banking system: Three level lotka-volterra m odel,” Procedia Economics and Finance , vol. 3, p. 251–255, 12 2012
2012
-
[18]
Volterra, Variazioni e fluttuazioni del numero di individui in specie a nimali conviventi
V. Volterra, Variazioni e fluttuazioni del numero di individui in specie a nimali conviventi . Atti della R. Accademia Nazionale dei Lincei, C. Ferrari, 1927
1927
-
[19]
Equilibrium points, stab ility and numerical solutions of fractional-order predator-prey and r abies models,
E. Ahmed, A. El-Sayed, and H. El-Saka, “Equilibrium points, stab ility and numerical solutions of fractional-order predator-prey and r abies models,” Journal of Mathematical Analysis and Applications , vol. 325, pp. 542–553, 01 2007
2007
-
[20]
Explosive route to chaos thro ugh a fractal torus in a generalized lotka-volterra model,
N. Samardzija and L. D. Greller, “Explosive route to chaos thro ugh a fractal torus in a generalized lotka-volterra model,” Bulletin of Mathematical Biology , vol. 50, pp. 465–491, 1988
1988
-
[21]
Fractional-o rder model of two-prey one-predator system,
M. Elettreby, A. Abdullah Al-Raezah, and T. Nabil, “Fractional-o rder model of two-prey one-predator system,” Mathematical Problems in Engineering , vol. 2017, pp. 1–12, 2017
2017
-
[22]
Memor y and mutualism in species sustainability: A time-fractional lotka-volterra model with harves ting,
M. M. Amirian, I. Towers, Z. Jovanoski, and A. J. Irwin, “Memor y and mutualism in species sustainability: A time-fractional lotka-volterra model with harves ting,” Heliyon, vol. 6, no. 9, p. e04816, 2020
2020
-
[23]
Analytical results on the unsteady rotational flow of fractional-order non-newton ian fluids with shear stress on the boundary,
M. M. Ghalib, A. A. Zafar, Z. Hammouch, M. B. Riaz, and K. Shabb ir, “Analytical results on the unsteady rotational flow of fractional-order non-newton ian fluids with shear stress on the boundary,” Discrete & Continuous Dynamical Systems-S , vol. 13, no. 3, p. 683, 2020
2020
-
[24]
A new study on the mathemat- ical modelling of human liver with caputo–fabrizio fractional derivativ e,
D. Baleanu, A. Jajarmi, H. Mohammadi, and S. Rezapour, “A new study on the mathemat- ical modelling of human liver with caputo–fabrizio fractional derivativ e,” Chaos, Solitons & Fractals, vol. 134, p. 109705, 2020
2020
-
[25]
A. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and applications of fractional differential equations , vol. 204. Elsevier Science Limited , 2006
2006
-
[27]
Stability analysis of a f ractional-order linear system described by the caputo-fabrizio derivative,
H. Li, J. Cheng, H.-b. Li, and S.-M. Zhong, “Stability analysis of a f ractional-order linear system described by the caputo-fabrizio derivative,” Mathematics, vol. 7, p. 200, 2019. 20
2019
-
[28]
Analysis of a fractional order pre y-predator model (3-species),
G. Selvam A, D. . R, and A. .D, “Analysis of a fractional order pre y-predator model (3-species),” Global Journal of Computational Science and Mathematics , vol. Volume 5, pp. 95–102, 2015
2015
-
[29]
Control of a fractional-order econ omical system via slid- ing mode,
S. Dadras and H. Momeni, “Control of a fractional-order econ omical system via slid- ing mode,” Physica A: Statistical Mechanics and its Applications , vol. 389, pp. 2434–2442, 2010
2010
-
[31]
A fast sampling method f or estimating the domain of attraction,
E. Najafi, R. Babuska, and G. Lopes, “A fast sampling method f or estimating the domain of attraction,” Nonlinear Dynamics , vol. 86, 2016. 21
2016
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.