ANALYSIS OF A FUNCTIONAL RESPONSE WITH PREY-DENSITY DEPENDENT HANDLING TIME FROM AN EVOLUTIONARY PERSPECTIVE

The paper formulates and analyzes a Rosenzweig–MacArthur predator–prey model with prey-density–dependent handling time h(x)=1/(bx+c) and a handling-time–dependent conversion factor γ(h)=∫ρ(τ)e^{-τ/h}dτ (its equations (1)–(5)) and studies invasion and evolution via adaptive dynamics, with invasion fitness defined as the long-term average growth rate of a rare mutant in the resident environment (their eq. (31)). It presents numerical evidence that even weakly density-dependent handling (b>0) invades monomorphic and dimorphic fixed-handling residents and that increasing b eventually removes cycles via a (final) supercritical Hopf, eliminating coexistence (abstract and Sections 3–5) . The candidate solution correctly derives and exploits the invasion-fitness derivative at b_m=0 to show selection for b>0 when c_m is small, in line with the paper’s numerical findings. However, it then asserts a “type-I limit” as b→∞ with dy/dt=(εβx−δ)y based on ε=lim_{h→0}γ(h)/h; under the paper’s γ(h) this is false, because γ(h)→0 as h→0 so γ(h)f(x,h)→0, not εβx. Hence the Dulac-based “no cycles for large b” argument is invalid. The paper instead documents loss of cycles by bifurcation computations and phase portraits, including a subcritical and then supercritical Hopf in some regimes (Figures 5–6) . The paper’s support for its key claims is mostly numerical and scenario-based rather than general theorems; the model’s solution contains a flawed limiting argument and unproven global claims (e.g., uniform positivity of g'(x) and Hopf criticality). Therefore, both are partially correct but incomplete.