{"id":"47510d4c-9478-4a8d-b083-b5821ffa061a","arxiv_id":"2504.18035","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"This paper presents a stability and bifurcation analysis plus time-optimal control for an additional-food predator-prey model, but the control formulation and two proofs contain load-bearing errors.","lead":"This paper analyzes a predator-prey model with Holling type III feeding, extra food for predators, and competition among predators, and claims to derive stability, bifurcations, and time-optimal pest control strategies. It aims to guide biological pest control, but several core proofs contain mathematical errors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The time-optimal control sections minimize the wrong objective: T in §9.1–9.2 equals ∫(1+αξ+x²)ds, so minimizing pseudo-time S does not minimize actual time T; the Pontryagin analysis and simulations solve a different problem.","rationale":"The reader's strongest claim and weakest assumption identify exactly this issue, and my reading confirms it. The time-scaling transformation is standard in mechanics, but it preserves a time-minimization objective only when the new cost includes the scaling factor dt/ds and one minimizes ∫(dt/ds)ds, not ∫ds. The paper minimizes S, i.e., ∫1 ds, so the state- and control-dependent weight 1+αξ+x² is absent from the objective. The numerical results in §9.3 report 'optimal times' of 2.1 and 5.05, but these are S-values; the paper never evaluates the actual elapsed time of the transformed trajectories. A direct comparison on the paper's own parameter set would settle the issue. I also noted a separate concern in the boundedness proof of Theorem 3.1, where the y-dependent terms are maximized without an explicit bound for the cross term, but the control-objective mismatch is the single most load-bearing issue because it invalidates the advertised optimal pest-management application. This is a correctness risk, not a novelty or stylistic issue; no amount of simulation of (13) can establish time-optimality for (12). I therefore agree with the reader's REJECT verdict and do not change it.","tokens_in":923,"tokens_out":1341,"duration_ms":77396,"concrete_test":"Use the paper's §9.3 parameters for the α-control problem: γ=7, ξ=0.1, δ=3, m=1, ε=0.3, from (x0,y0)=(5,2) to (x̄,ȳ)=(1,4). First solve (13) in CasADi to obtain the optimal S, control α*(s), and state path. Recover real time by integrating t(s)=∫_0^s (1+α*(σ)ξ+x(σ)²)dσ and record T_recovered. Then solve the original problem (12) by direct transcription with the same endpoints and control bounds, using the true time variable. If T_recovered exceeds the minimal T from (12) by any nontrivial amount, the transformation is not time-preserving and the claimed equivalence in §9.1 fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 9.1 (and identically in Section 9.2), the paper introduces s by dt = (1 + αξ + x²) ds and asserts that this converts the time-optimal problem (12) into (13), whose objective is to minimize S. This equivalence is false. The original cost is T = ∫_0^T dt, which in s-coordinates becomes T = ∫_0^S (1 + αξ(s) + x(s)²) ds. Since α is the control and x varies along the trajectory, this integrand is not constant. A control with smaller S can yield a larger T if it keeps x² + αξ large; the optimization in (13) entirely drops the factor 1 + αξ + x² from the cost. Consequently, the Hamiltonian, adjoint equations, bang-bang and singular switching conditions (14)–(17), and the CasADi simulations in §9.3 characterize optimal controls for the pseudo-time problem, not for the stated minimum-time problem. The paper's central pest-management claim—that solving (13)/(19) provides time-optimal strategies—is therefore unsupported. This is a load-bearing error in the control contribution.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes a two-species prey–predator model with Holling type-III functional response, additional food supplied to predators, and intra-specific competition among predators. After nondimensionalization, the authors prove positivity and attempt to prove boundedness and global existence of solutions, classify equilibria and their local stability, and derive transcritical, saddle-node, and Hopf bifurcation results. They also present numerical evidence of a hysteresis loop, study the global dynamics in the α–ξ parameter plane, and formulate time-optimal control problems in which the quality or the quantity of additional food is the control. The control problems are treated by introducing a new time variable via dt = (1 + αξ + x²) ds, which the authors claim converts each time-optimal problem into a linear problem for the pseudo-time S; bang-bang and singular control conditions are derived and simulated with CasADi for pest-management scenarios.","tokens_in":118,"tokens_out":9352,"duration_ms":159660,"significance":"The model and the intended application to biological pest control are topical, and the paper contains a fairly complete nullcline-based stability classification as well as an explicit numerical demonstration of hysteresis. The control section is the most distinctive contribution, and the authors are transparent that this is an initial attempt. However, the central claims are not reliable as written: the boundedness proof is invalid, the saddle-node bifurcation theorem violates the stated Sotomayor condition, and the time-rescaling in the control sections changes the objective so that the Pontryagin analysis and simulations solve a pseudo-time problem rather than the advertised minimum-time problem. Because the latter error is load-bearing for the main novelty, the results cannot be accepted in their present form.","major_comments":[{"comment":"The boundedness proof uses an invalid inequality. After the displayed computation, the authors bound the right-hand side by γ(1+K)²/4 + ξ/ε + (K−m)²/(4ε), but the y-dependent terms ξy/(1+x²+αξ) + ((K−m)/δ)y − (ε/δ)y² cannot be bounded in this way. Even using ξy/(1+x²+αξ) ≤ ξy, the maximum over y is δ(ξ + (K−m)/δ)²/(4ε), not ξ/ε + (K−m)²/(4ε); the factor δ is missing and the numerator of the quadratic term is not handled. A concrete counterexample to the displayed bound is obtained with α=0, m=0, K=0, x=0, ξ=10, ε=1, δ=1, and y=5, where the left-hand side is 25 while the proposed bound is γ/4+10, which is smaller for γ<60. Since this bound is the only argument for ultimate boundedness and hence for global existence of the solution, Theorem 3.1 is not proved as stated. The argument may be repairable, but the present proof is not valid.","section":"Section 3.2, Theorem 3.1"},{"comment":"The saddle-node bifurcation claim fails the authors' own Sotomayor condition. In the proof, they assert that WᵀH_ξ(E₂;ξ*) = δ(δξ−m(1+αξ))/(ε(1+αξ)³) ≠ 0. But at the stated bifurcation value ξ* = m/(δ−mα), the factor δξ−m(1+αξ) is exactly zero, so this quantity is zero, not nonzero. Thus the Sotomayor nondegeneracy condition required for a saddle-node bifurcation is not satisfied. Moreover, at this parameter value E₂ coincides with E₀=(0,0), and the exchange between these equilibria has the character of a transcritical bifurcation rather than a saddle-node bifurcation. The theorem, and the numerical discussion of Figure 3 that refers to it, therefore do not establish the claimed bifurcation.","section":"Section 6.2, Theorem 6.2"},{"comment":"The transformation dt = (1+αξ+x²)ds does not preserve the time-optimal objective. The original problem (12) minimizes T = ∫₀ᵀ dt. Under the change of variable, this becomes T = ∫₀^S (1+αξ(s)+x(s)²)ds, and since α is the control and x varies along the trajectory, the integrand is not constant. The linear problem (13) minimizes only S, dropping the factor 1+αξ+x² from the cost. Hence a control that is optimal for (13) need not be time-optimal for (12): a trajectory with smaller pseudo-time S can have a larger actual time T if it keeps 1+αξ+x² large. Consequently, the Hamiltonian, the adjoint equations, the switching conditions (14)–(17), and the singular-arc analysis characterize the pseudo-time problem, not the stated minimum-time problem. This is a load-bearing error for the paper's main pest-management application.","section":"Section 9.1, Eqs. (12)–(13)"},{"comment":"The same objective error affects the quantity-control problem: Eq. (19) minimizes S after the same transformation, but the actual time in problem (18) is still ∫₀^S (1+αξ+x²)ds. The numerical simulations in Section 9.3 solve (13) and (19) and then report the resulting 'optimal time' values (2.1 and 5.05 units), but those numbers are the pseudo-time S, not the physical time T. The control trajectories and pest-management conclusions in Figures 11 and 12 are therefore not supported as solutions of the time-optimal problems (12) and (18).","section":"Sections 9.2–9.3"}],"minor_comments":[{"comment":"There are several typographical errors, including 'Hysterisis' in the Section 6.4 heading and the abstract, 'guarentees' in Section 3.2, and 'atmost' in Section 4. The manuscript would benefit from a careful proofreading pass.","section":"Throughout"},{"comment":"The MSC codes 60H10, 60J65, and 60J70 are stochastic-analysis codes and appear irrelevant to this deterministic dynamical-systems paper; the authors should replace them with appropriate codes such as 34C23, 92D25, and 49J15.","section":"MSC codes"},{"comment":"In the displayed simplification of the determinant, the variables x and y are used without explicitly indicating that the expression is evaluated at (x*, y*); this makes the formula hard to follow.","section":"Section 5.1, Eq. (10)"},{"comment":"The control bounds α_min, α_max, ξ_min, and ξ_max are never specified, and the simulation section does not state the values of the bounds or the discretization parameters used in the CasADi multiple-shooting implementation, which limits reproducibility.","section":"Section 9.3"},{"comment":"The division of the α–ξ parameter space into regions is presented through representative phase portraits and verbal descriptions, but no theorem establishes that the listed curves exhaust the possible dynamics; this should be stated more cautiously.","section":"Section 7"}],"recommendation":"reject","confidential_remarks":"The time-optimal control part is the main claimed novelty, and the equivalence between the original and transformed problems is asserted without proof and is in fact false. This is not a local fix: the entire Pontryagin analysis and the numerical optimal-control results would need to be reformulated for the true minimum-time objective. The boundedness and saddle-node issues are also substantive. I therefore recommend rejection rather than major revision, although the dynamical-systems portions of the manuscript may contain salvageable material for a future submission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a mixed bag. The biological ODE work is solid and useful; the optimal-control part is not. The model—additional food with Holling type III and intraspecific competition—is a natural extension of the authors' earlier work, and they give a fairly complete qualitative study: up to five interior equilibria, stability conditions, nullcline configurations, a hysteresis loop from slow periodic forcing, and a lower bound on the pest population at equilibrium. That part reads competently and would be worth having in the literature, once the proof issues below are cleaned up.\n\nThe soft spots are real. First, the boundedness proof in Section 3.2 contains an invalid inequality: they bound ξ y/(1+x²+αξ) as if it were just ξ y, then combine it with the other y-terms in a way that doesn't produce the stated M. The result is probably true, but the proof as written is wrong. Second, the \"saddle-node\" bifurcation in Theorem 6.2 is misclassified. At ξ = m/(δ−mα), the equilibrium E2 collides with E0, and their own WᵀH_ξ evaluates to zero. That is a transcritical (or pitchfork-like) exchange, not a saddle-node per Sotomayor. The figure may be showing something else, but the analytic claim is wrong.\n\nThe load-bearing flaw is in Section 9. The transformation dt = (1+αξ+x²)ds is fine as a time-reparameterization, but then the paper minimizes S instead of T. The actual objective is T = ∫(1+αξ+x²)ds, which depends on the control α and the state x. Minimizing S does not minimize T, and the Hamiltonian, switching conditions, and CasADi simulations in §9.3 all characterize the pseudo-time problem, not the stated minimum-time pest-control problem. This is not a minor slip; it invalidates the main applied claim of the paper.\n\nI would not take the control results seriously until this is redone. A referee should push for either a correct treatment with the true time objective or an honest reframing as a problem in the transformed time variable. The rest of the paper could be salvaged with modest effort. It's a paper that deserves a careful referee, but as it stands the central control contribution should not be accepted.","headline":"A systematic ODE analysis of a new predator-prey model, undermined by a control section that solves the wrong objective.","tokens_in":17275,"tokens_out":4037,"would_cite":false,"duration_ms":38871,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37A50","60H10","60J65","60J70"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that adding food to predators in a Holling type-III prey-predator model with predator competition creates bistability, hysteresis, and bang-bang time-optimal control schedules.","keywords":["prey-predator system","Holling type-III functional response","additional food","intra-specific competition","bifurcation","hysteresis loop","time-optimal control","pest management"],"falsifier":"For the parameter values and endpoint states of the simulations in Section 9.3, take the reported optimal trajectory, compute the actual elapsed time $T=\\int_0^S (1+\\alpha\\xi+x^2(s))\\,ds$, and solve the same rendezvous by minimizing $T$ directly; if the reported trajectory's $T$ exceeds the directly minimized $T$, the time-optimality claim is refuted.","tokens_in":16404,"feed_emoji":"🐛","tokens_out":9014,"duration_ms":86945,"temperature":0.7,"pith_summary":"The paper builds a two-species model in which prey grow logistically, predators attack with a Holling type-III (sigmoidal) response, predators receive additional food whose quality and quantity are parameters, and predators compete among themselves. It tries to establish the complete qualitative picture of this system: global boundedness, the existence and stability of up to five interior equilibria, transcritical and saddle-node bifurcations, Hopf cycles, and an S-shaped equilibrium curve that produces a hysteresis loop when predator competition varies slowly. It then claims that the quality or quantity of additional food can serve as a time-optimal control, with the optimal strategy being bang-bang or singular, and applies this to pest management. A sympathetic reader would care because these predictions decide when supplemental food helps biocontrol versus when it creates bistable pest outbreaks.","feed_headline":"Extra predator food makes pest control bistable with hysteresis","feed_subtitle":"A predator-prey model shows when supplemental food holds pests at a low-density branch and when it can't.","key_machinery":"The argument runs on the model (3) with the Holling type-III denominators $1+x^2+\\alpha\\xi$, the quadratic predator crowding term $-\\epsilon y^2$, the nullcline equations, the quintic equilibrium equation determining $x^*$, and the time reparameterization $dt=(1+\\alpha\\xi+x^2)\\,ds$. The reparameterization is the load-bearing device for the control claims: it makes the control appear linearly in the Hamiltonian, turning the search for minimum time into a linear problem whose optimum is bang-bang or singular.","core_discovery":"The paper claims that system (3), which couples logistic prey growth to a Holling type-III predation term and to predators that receive additional food and compete among themselves, has a complete qualitative theory: solutions stay positive and bounded; up to five interior equilibria can occur; stability is governed by $x^{*2} < 1+\\alpha\\xi$; the system undergoes transcritical and saddle-node bifurcations as the quantity of additional food $\\xi$ varies, and Hopf bifurcation as predator crowding $\\epsilon$ varies; the interior-equilibrium curve is S-shaped, producing two saddle-node folds and a hysteresis loop under a slow periodic sweep of $\\epsilon$; and in the $(\\alpha,\\xi)$ plane the dynamics split into regions of no interior equilibrium, unique stable interior equilibrium, or bistability between the predator-free and interior equilibria. For pest management, the paper argues the pest-free equilibrium is always saddle when it exists, so eradication is unstable, and pest population can be held to a low-pest branch by slow releases; the minimum stable pest density is $\\frac{\\epsilon}{1+\\epsilon/\\gamma}$. For control, it formulates two minimum-time problems and claims the change of independent variable $dt=(1+\\alpha\\xi+x^2)\\,ds$ converts them into control-linear problems whose optimal policies are bang-bang or singular.","pith_inferences":["From the equations as written, problems (13) and (19) minimize $S$, not the actual elapsed time $T$, because the reparameterization factor $1+\\alpha\\xi+x^2$ is state-dependent; the reported optimal times of 2.1 and 5.05 units may be values of $S$, so the schedules may not be time-optimal for the original system.","A repair is available in principle: minimize $\\int_0^S (1+\\alpha\\xi+x^2)\\,ds$ over the same linearized dynamics, or derive the maximum-principle conditions in the original time variable; the switching curves in the paper would change under either correction.","The hysteresis structure suggests a testable management protocol: ramp natural-enemy additions slowly to hold the low-pest branch, and avoid pulse releases large enough to cross the second fold into a high-pest jump."],"forward_implications":["Pest eradication is not a stable outcome: whenever the predator-only equilibrium exists, it is a saddle, so the operational goal must be a low-pest interior equilibrium.","The S-shaped equilibrium curve plus two folds means slow changes in predator crowding produce hysteresis: the system stays on one branch past the fold and then jumps, so release history, not just release amount, decides which pest level is reached.","The minimum stable pest population is set by predator crowding as $\\frac{\\epsilon}{1+\\epsilon/\\gamma}$; below that floor the interior equilibrium cannot be pushed, so competition strength determines how low pest density can go.","In the $(\\alpha,\\xi)$ plane, additional food divides the dynamics into regions with no interior equilibrium, a unique stable interior equilibrium, and bistability between predator-free and interior states, so the effect of food supplementation depends on initial populations.","The optimal control problems with quality or quantity of additional food as control are reduced, according to the paper, to linear-in-control problems, so the optimal strategies are bang-bang with possible singular arcs."],"supporting_citations":[{"why":"Supplies the additional-food functional response and parameter definitions used to build model (2).","marker":"[3]"},{"why":"Introduces the intra-specific competition term in predator-prey models that the paper adds to the Holling type-III setting.","marker":"[6]"},{"why":"Provides prior bifurcation analysis of predator-prey models with predator competition that this paper extends.","marker":"[7]"},{"why":"Provides the comparison inequality used in the boundedness proof.","marker":"[8]"},{"why":"Gives the bifurcation criterion used to establish transcritical and saddle-node bifurcations.","marker":"[9]"},{"why":"Supplies the optimal-control theory used to characterize bang-bang and singular controls.","marker":"[10]"},{"why":"Supplies the numerical optimization framework used for the simulated time-optimal trajectories.","marker":"[11]"}],"fun_headline_variants":["Extra predator food turns pest control bistable with hysteresis","Hysteresis in predator-prey systems from extra food and competition","Supplemental food and predator crowding create pest-control hysteresis","Bistable pest control emerges when predators get extra food"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The shortest-time conclusion depends on the assumption that replacing the real clock $t$ by the new variable $s$ through $dt=(1+\\alpha\\xi+x^2)\\,ds$ leaves the meaning of \"minimum time\" intact; since $T=\\int_0^S (1+\\alpha\\xi+x^2)\\,ds$, minimizing $S$ is not generally the same as minimizing $T$.","fun_headline_variants_meta":{"raw":{"variants":["Extra predator food turns pest control bistable with hysteresis","Hysteresis in predator-prey systems from extra food and competition","Supplemental food and predator crowding create pest-control hysteresis","Bistable pest control emerges when predators get extra food"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000496,"raw_usage":{"total_tokens":2453,"prompt_tokens":984,"completion_tokens":1469,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":1402}},"tokens_in":600,"tokens_out":1469,"duration_ms":11242,"temperature":1.0,"reasoning_tokens":1402,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:26:30.075860+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the parameter values and endpoint states of the simulations in Section 9.3, take the reported optimal trajectory, compute the actual elapsed time $T=\\int_0^S (1+\\alpha\\xi+x^2(s))\\,ds$, and solve the same rendezvous by minimizing $T$ directly; if the reported trajectory's $T$ exceeds the directly minimized $T$, the time-optimality claim is refuted.","supporting_citations":[{"cited_title":"Structural and dynamic stability of model predator-prey systems,","cited_arxiv_id":null,"evidence_quote":"Introduces the intra-specific competition term in predator-prey models that the paper adds to the Holling type-III setting."},{"cited_title":"Bifurcations and hydra effects in bazykin’s predator–prey model,","cited_arxiv_id":null,"evidence_quote":"Provides prior bifurcation analysis of predator-prey models with predator competition that this paper extends."},{"cited_title":"The gronwall inequality,","cited_arxiv_id":null,"evidence_quote":"Provides the comparison inequality used in the boundedness proof."},{"cited_title":"Perko, Differential equations and dynamical systems","cited_arxiv_id":null,"evidence_quote":"Gives the bifurcation criterion used to establish transcritical and saddle-node bifurcations."},{"cited_title":"Cesari, Optimization—theory and applications: problems with ordinary differential equations","cited_arxiv_id":null,"evidence_quote":"Supplies the optimal-control theory used to characterize bang-bang and singular controls."}],"review_version":1}