{"id":"c47fc7a6-43a8-4c72-8391-3b9043591ae9","arxiv_id":"2501.02386","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"This paper analyzes stability and bifurcations of a discrete predator-prey map built with conformable fractional derivatives, but the central bifurcation theorem is unsupported and contains an impossible condition on a positive parameter.","lead":"The paper discretizes a Lotka-Volterra predator-prey model with a conformable fractional derivative and analyzes its stability and flip and Neimark-Sacker bifurcations. The main bifurcation theorem is not proven: it states a result for a negative bifurcation parameter even though the parameter is positive, and the Lyapunov coefficient is never computed.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 3.1's invariant-curve direction is unsupported and internally inconsistent: A_NS requires b>0, yet the theorem classifies negative b as attracting and positive b as repelling, and the first Lyapunov coefficient it quotes is never evaluated.","rationale":"Good-faith reading: the paper's central contribution is the claim that the discrete map (2.3) undergoes a Neimark-Sacker bifurcation with a specified stability type as b crosses the A_NS curve. For this to be true, Theorem 3.1 would need (i) eigenvalues crossing the unit circle transversally, (ii) non-resonance, and (iii) a sign-definite first Lyapunov coefficient. The manuscript checks a transversality condition d|\\lambda1|/d\\epsilon\\ne0 and some non-resonance exclusions, and it writes the standard formula for a. But it never evaluates a, then asserts that negative b yields an attracting curve and positive b a repelling curve. This is not a consensus disagreement; it is an internal gap and contradiction. A_NS explicitly restricts b>0, so the negative-b branch cannot occur. Even if the symbol were a typo for 'a', no computation of a appears anywhere, and no numerical example computes it. Thus the main claim is unsupported as written. The reader's weakest_assumption (no derivation of the map from the fractional system) is also serious but is partially separable: all theorems are about the map (2.3) itself. The Lyapunov-coefficient defect attacks the theorem even if one grants the discretization. Therefore I agree with the REJECT verdict and partially with the reader's diagnosis; no verdict change is needed.","tokens_in":12750,"tokens_out":5339,"duration_ms":53892,"concrete_test":"Take the Example 5.1 parameters, which lie on A_NS: \\alpha=0.7, h=0.3, d=2, r=1, b=3.62597. Substitute the coefficients of (3.2) into the displayed formula for the first Lyapunov exponent a at \\epsilon=0 and evaluate its sign, either symbolically with exact arithmetic or with high-precision numerics. If a<0, the invariant curve is attracting even though b>0, directly disproving Theorem 3.1's classification; if a=0, the nondegeneracy condition for the claimed Neimark-Sacker bifurcation is absent. This single computation settles whether the theorem's central stability claim holds at the paper's own illustrative parameter values.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central bifurcation claim fails at Theorem 3.1. The set A_NS is defined with b>0 and b=d+\\alpha/h^\\alpha; b is therefore strictly positive. Theorem 3.1 nonetheless concludes 'when b is negative, an attracting invariant curve emanates from E*' and 'when b is positive, a repelling invariant curve emanates'. The negative-b branch contradicts the theorem's own domain. Even reading 'b' as a typo for the first Lyapunov coefficient 'a', the manuscript never computes a; it only writes a formula in terms of \\xi20,\\xi11,\\xi02,\\xi21 and immediately asserts the conclusion. For a Neimark-Sacker bifurcation the direction of the emergent invariant curve is fixed by the sign of a, and existence of the curve requires a\\ne0; neither is established. The stated theorem therefore does not follow from the preceding derivation. The same unsupported pattern appears in the flip-bifurcation section, where \\beta1 and \\beta2 are declared nonzero without evaluation. But Theorem 3.1 is the clearest and most load-bearing failure because it is the paper's main advertised result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a discrete predator-prey map (2.3) as a conformable fractional-order discretization of the continuous Lotka-Volterra system (2.2). It computes the positive fixed point, gives a stability classification in Theorem 2.2, and then claims a Neimark-Sacker bifurcation in Theorem 3.1 and a period-doubling bifurcation in Section 3. A hybrid control strategy is added in Section 4, and numerical bifurcation diagrams and phase portraits are presented in Section 5. The central advertised result is Theorem 3.1, which asserts that the system undergoes a Neimark-Sacker bifurcation at the positive fixed point when the parameter b varies near the set A_NS, with attracting or repelling invariant curves depending on the sign of b.","tokens_in":12938,"tokens_out":6788,"duration_ms":65336,"significance":"If fully established, the paper would provide a useful discrete fractional-order bifurcation analysis of a predator-prey model. The manuscript does contain constructive components: an explicit fixed-point formula, the characteristic polynomial, a stability classification, a hybrid control law with Jury conditions, and several numerical illustrations. Those parts are standard and the numerics are suggestive. However, the main theoretical claims are not supported by the text: the first Lyapunov coefficient in Theorem 3.1 is never evaluated, the flip-bifurcation coefficients are asserted nonzero without computation, and the derivation of map (2.3) from the fractional system is omitted. As a result, the paper does not currently establish its advertised bifurcation results.","major_comments":[{"comment":"Theorem 3.1 is internally inconsistent. The set A_NS is defined by the conditions b > 0 and b = d + α/h^α, so any point in A_NS has b > 0. The theorem nevertheless concludes that 'when b is negative, an attracting invariant curve emanates from E*' and 'when b is positive, a repelling invariant curve emanates from E*'. Since b is strictly positive on A_NS, the negative-b branch is outside the hypotheses and the statement has no content on that branch. Moreover, the first Lyapunov coefficient a is only written as a generic formula in terms of ξ20, ξ11, ξ02, and ξ21; the manuscript never substitutes the computed C and D coefficients, never evaluates a, and never verifies a ≠ 0. The transversality condition d|λ|/dϵ ≠ 0 is shown, but existence of the invariant curve and the direction of bifurcation require the sign and non-vanishing of a. Therefore Theorem 3.1 does not follow from the preceding derivation.","section":"Section 3, Theorem 3.1 and the set A_NS"},{"comment":"The map (2.3) is presented as the conformable fractional discretization of the continuous fractional predator-prey system (2.2), but no derivation is given anywhere before the conclusion. The 'piecewise constant argument technique, conformal fractional case' is mentioned only in Section 6, and no theorem, lemma, or calculation connects (2.2) to (2.3). This is load-bearing because all stability and bifurcation statements in the paper concern the map (2.3); if that map is not actually obtained from the fractional system, the ecological conclusions of the paper have no basis. The authors need to state the discretization rule and show the intermediate steps, or clearly identify (2.3) as an independent discrete model.","section":"Section 2, Eq. (2.3)"},{"comment":"The period-doubling bifurcation is not proven. After the center-manifold reduction, the manuscript declares β1 = ∂²Φ/(∂e∂ϵ*) + (1/2)(∂Φ/∂ϵ*)(∂²Φ/∂e²) ≠ 0 and β2 = (1/6)(∂³Φ/∂e³) + (1/2)(∂²Φ/∂e²)² ≠ 0, but it never computes these expressions from the coefficients c_i, d_i, h_i, nor does it show that they are nonzero for any parameter range. The same section also has apparent algebra/notation slips in the center-manifold coefficient comparison (for example, the equation for the e_n ϵ* term mixes d4 and d5 and the h2 equation is not consistent with the preceding line). Since the sign and non-vanishing of β1 and β2 are precisely the nondegeneracy conditions for a flip bifurcation, the claim that the system 'observes a period-doubling bifurcation' within A_PD is unsupported.","section":"Section 3, flip bifurcation"}],"minor_comments":[{"comment":"The denominator in the second bullet contains an undefined symbol 'c' (the expression '4bα²/[c d(-bh^α + h^α d + 2α)]'), which makes the stated sink condition unreadable.","section":"Theorem 2.2, item 1"},{"comment":"The phrase 'The positive fixed points E∗ are complex' is unclear; the intended meaning is presumably that the eigenvalues are complex, not that the fixed point is complex.","section":"Theorem 2.2, item 5"},{"comment":"The set A_PD excludes h^α d r/(bα) ≠ 4, while Theorem 2.2 item 4 excludes both h^α d r/(bα) ≠ 2 and ≠ 4; the manuscript does not explain why the value 2 is no longer excluded in the flip-bifurcation set.","section":"Section 3, A_PD definition"},{"comment":"The parameter listing 'd = 2.5, d = 2.91667' uses d twice; the second value is evidently meant to be b, which makes the numerical flip example difficult to reproduce as written.","section":"Section 5, Example 5.2"},{"comment":"The caption of Figure 2 says 'Bifurcation Diagrams for the Controlled System (4.1)', but the text refers to 'figure (5)' for the controlled system; the figure numbering is inconsistent.","section":"Figure captions and text"}],"recommendation":"reject","confidential_remarks":"The paper's main advertised theorem is internally inconsistent because it draws conclusions for negative b on a set where b is strictly positive, and the crucial Lyapunov coefficient is never evaluated. The missing derivation of the discretization is a foundational gap, and the flip-bifurcation coefficients are asserted without proof. These are load-bearing errors that would require a substantial new computation and a rewrite, not a local revision, so I do not see a path to acceptance in the current scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this paper dresses up a known map in fractional clothing, and the central bifurcation theorem contains an internal contradiction. I would not send it to peer review.\n\nWhat's good: the map (2.3) is the standard exponential discrete Lotka-Volterra map with step h^alpha/alpha, so the stability conditions in Theorem 2.2 are classical and mostly correct. The numerical examples do reproduce the expected Neimark-Sacker and flip bifurcations, and the hybrid chaos control is a standard add-on. The paper is organized and follows a recognizable template; the equations are presented in full, not hidden.\n\nThe problems are decisive. There is no derivation of (2.3) from the conformable fractional system; the discretization is simply asserted. That is the paper's only claimed novelty, so it is load-bearing. Theorem 3.1 states that with b negative an attracting invariant curve appears, but the set A_NS is defined with b>0 and b=d+alpha/h^alpha; the negative branch is outside the theorem's own domain. Even if one reads 'b' as a typo for the Lyapunov coefficient a, a is never evaluated—only a formula in terms of xi20, xi11, xi02, xi21. So the direction of the invariant curve is unsupported. The flip bifurcation section has the same pattern: beta1 and beta2 are declared nonzero without computation, and the center manifold expansion contains undefined coefficients (c7 missing in the e equation, a 'b9' that should be d9). The conclusion's claim of 'more accurate than [23]' is never quantified or compared.\n\nI agree with the reader's verdict and the stress-test note. This is not a paper with a few typos; the advertised result is not proven, and the novelty reduces to a rescaling. I'd desk reject and suggest the authors verify their algebra and compare directly with [23] if they want to claim improved accuracy. Serious referee time would be wasted.","headline":"A known map in fractional disguise: the main theorem contradicts its own parameter domain and the Lyapunov computation that would fix it is never done; not worth referee time.","tokens_in":13524,"tokens_out":3266,"would_cite":false,"duration_ms":31306,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["39A28","37G15","34A08","92D25"],"pacs":[],"model":"deepseek-v4-flash","headline":"Using a conformable-fractional discretization, the paper locates the exact parameter curves where a predator-prey coexistence state undergoes Neimark-Sacker and flip bifurcations, and gives a control term that stabilizes it.","keywords":["conformable fractional derivative","Lotka-Volterra predator-prey model","Neimark-Sacker bifurcation","flip bifurcation","stability of fixed points","chaos control","discrete-time dynamical system","invariant curve"],"falsifier":"Take the parameter values of Example 5.1, continue the map (2.3) through $b=d+\\alpha/h^{\\alpha}\\approx 3.62597$, and compute the first Lyapunov exponent from the published $\\xi$-coefficients: Theorem 3.1 stands only if the exponent changes sign exactly at the crossing and a closed invariant curve is visible in the phase portraits at the stated values. Separately, deriving (2.3) from (2.2) by the piecewise-constant conformable argument would settle whether the map is the true discretization.","tokens_in":12516,"feed_emoji":"🦊","tokens_out":11753,"duration_ms":103099,"temperature":0.7,"pith_summary":"This paper sets out to show that a conformable-fractional discretization of the continuous Lotka-Volterra predator-prey system leads to a discrete map whose local dynamics can be classified completely. The positive coexistence equilibrium $E^*$ is a sink, saddle, source, or non-hyperbolic point according to sharp inequalities on the parameters, and crossing the set $b = d + \\alpha/h^{\\alpha}$ forces the equilibrium to undergo a Neimark-Sacker bifurcation, with an attracting invariant curve appearing on one side and a repelling one on the other. The paper also proves a flip-bifurcation surface and proposes a hybrid control term that leaves $E^*$ fixed while allowing the bifurcation threshold to be shifted. If these results are right, ecologically interesting oscillations and chaos in the fractional predator-prey model can be located, delayed, or suppressed by adjusting $b$ and the control parameter $\\beta$.","feed_headline":"Exact bifurcation curves found for fractional predator-prey map","feed_subtitle":"A conformable-fractional discretization gives provable stability and bifurcation thresholds for the coexistence state.","key_machinery":"The central object is the discretized map (2.3), obtained by replacing the fractional derivative in (2.2) with an exponential step whose fractional order enters through the factor $h^{\\alpha}/\\alpha$: $x_{n+1}=x_n e^{(r(1-x_n)-b y_n)h^{\\alpha}/\\alpha}$, $y_{n+1}=y_n e^{(b x_n-d)h^{\\alpha}/\\alpha}$. The argument then runs through the Jacobian at $E^*$, the characteristic polynomial (2.6), the transversality derivative of the eigenvalue modulus, and the normal-form coefficients $\\xi_{20},\\xi_{11},\\xi_{02},\\xi_{21}$ for the first Lyapunov exponent; for the flip case it uses a center-manifold computation with the coefficients $\\beta_1$ and $\\beta_2$. This machinery converts bifurcation detection into algebra on explicit parameter curves.","core_discovery":"The paper's central claim is that system (2.3) -- the map $x_{n+1}=x_n e^{(r(1-x_n)-b y_n)h^{\\alpha}/\\alpha}$, $y_{n+1}=y_n e^{(b x_n-d)h^{\\alpha}/\\alpha}$ -- is the conformable fractional discretization of (2.2), and that its coexistence fixed point is $E^*=(d/b,\\,(b-d)r/b^2)$ for $b>d$. Its local stability is governed by the characteristic polynomial $\\lambda^2+p\\lambda+q=0$ with $p=-2+drh^{\\alpha}/(b\\alpha)$ and $q=\\frac{-drh^{\\alpha}(dh^{\\alpha}+\\alpha)+b(h^{2\\alpha}dr+\\alpha^2)}{b\\alpha^2}$. Theorem 2.2 lists exactly when $E^*$ is a sink, saddle, source, or flip-unstable. Theorem 3.1 states that on $A_{NS}=\\{0<\\alpha\\le 1,\\, b=d+\\alpha/h^{\\alpha},\\, 0<drh^{\\alpha}/(b\\alpha)<4\\}$, the eigenvalue moduli cross the unit circle with nonzero speed $d|\\lambda|/d\\epsilon = h^{3\\alpha}dr/(2\\alpha^2(h^{\\alpha}d+\\alpha))\\ne 0$, so a Neimark-Sacker bifurcation occurs at $E^*$; when the perturbed parameter is negative, an attracting invariant curve emanates, and when positive, a repelling one. A separate set $A_{PD}$ gives a period-doubling bifurcation, and the hybrid control system (4.1) preserves $E^*$ while its Jury inequalities provide a stability window.","pith_inferences":["Editorial inference: since the paper declares $b>0$ throughout, the 'b negative' clause in Theorem 3.1 lies outside the stated parameter domain; reading it ecologically would flip the sign of the predation term, so the attracting-curve statement is better understood as a property of the mathematical continuation rather than a direct biological scenario.","Editorial inference: all thresholds depend on the product $h^{\\alpha}/\\alpha$, so changing the fractional order $\\alpha$ and the step size $h$ together moves the bifurcation curves; this gives a direct way to test the model by checking whether the onset of oscillations shifts as predicted when $\\alpha$ is varied.","Editorial inference: the same conformable discretization and center-manifold apparatus could be applied to fractional predator-prey models with Holling type II or III functional responses; this paper treats only the bilinear Lotka-Volterra interaction."],"forward_implications":["The coexistence equilibrium is stable only inside the parameter regions listed in Theorem 2.2; outside them, small perturbations produce saddle or source behavior, so the map can serve as a stability checklist for a fractional predator-prey interaction.","At $b = d + \\alpha/h^{\\alpha}$ with $0<drh^{\\alpha}/(b\\alpha)<4$, the system generically undergoes a Neimark-Sacker bifurcation: an invariant closed curve appears near the threshold, and the dynamics become quasi-periodic, with the sign of the perturbation deciding whether the curve attracts or repels.","On the flip surface $h^{\\alpha}>2\\alpha/r$ and $b = h^{\\alpha} d r (h^{\\alpha} d + 2\\alpha)/(d h^{2\\alpha} r + 4\\alpha^2)$, the fixed point period-doubles, starting a cascade that the numerical runs connect to chaotic behavior.","The hybrid controller (4.1) leaves the positive fixed point unchanged and, whenever the Jury conditions (4.3) hold, stabilizes it, so the control parameter $\\beta$ can postpone or suppress bifurcation-induced oscillations."],"supporting_citations":[{"why":"Provides the conformable fractional derivative definition on which the discretization of the continuous model is based.","marker":"[20]"},{"why":"Cited with [25] as the conformable fractional discretization technique used to turn system (2.2) into the discrete map (2.3).","marker":"[24]"},{"why":"Cited with [24] for the conformable fractional-order discretization procedure that motivates the map (2.3).","marker":"[25]"},{"why":"Supplies Lemma 1, the eigenvalue stability criterion for discrete maps, used to classify the fixed point E* in Theorem 2.2.","marker":"[26]"},{"why":"Cited together with [26] as the source of the eigenvalue stability Lemma 1 used in Theorem 2.2.","marker":"[13]"},{"why":"Provides the integer-order discrete predator-prey baseline that the paper says its fractional results improve on.","marker":"[23]"},{"why":"Supplies the bifurcation and normal-form theory behind the Neimark-Sacker and flip bifurcation derivations.","marker":"[22]"}],"fun_headline_variants":["Fractional discretization reveals exact bifurcation curves","Provable Neimark-Sacker and period-doubling in predator-prey","Conformable map pinpoints stability and bifurcation thresholds","Exact stability windows for fractional predator-prey system"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the map (2.3) really is the conformable-fractional discretization of the continuous fractional Lotka-Volterra system; the paper states this but supplies no derivation, and every stability and bifurcation result is a statement about that map.","fun_headline_variants_meta":{"raw":{"variants":["Fractional discretization reveals exact bifurcation curves","Provable Neimark-Sacker and period-doubling in predator-prey","Conformable map pinpoints stability and bifurcation thresholds","Exact stability windows for fractional predator-prey system"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000641,"raw_usage":{"total_tokens":2979,"prompt_tokens":1004,"completion_tokens":1975,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":1906}},"tokens_in":620,"tokens_out":1975,"duration_ms":14029,"temperature":1.0,"reasoning_tokens":1906,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:13:39.638048+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the parameter values of Example 5.1, continue the map (2.3) through $b=d+\\alpha/h^{\\alpha}\\approx 3.62597$, and compute the first Lyapunov exponent from the published $\\xi$-coefficients: Theorem 3.1 stands only if the exponent changes sign exactly at the crossing and a closed invariant curve is visible in the phase portraits at the stated values. Separately, deriving (2.3) from (2.2) by the piecewise-constant conformable argument would settle whether the map is the true discretization.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the conformable fractional derivative definition on which the discretization of the continuous model is based."},{"cited_title":"and Yang, XJ","cited_arxiv_id":null,"evidence_quote":"Cited with [25] as the conformable fractional discretization technique used to turn system (2.2) into the discrete map (2.3)."},{"cited_title":"and Yang, X.-J","cited_arxiv_id":null,"evidence_quote":"Cited with [24] for the conformable fractional-order discretization procedure that motivates the map (2.3)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Lemma 1, the eigenvalue stability criterion for discrete maps, used to classify the fixed point E* in Theorem 2.2."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Cited together with [26] as the source of the eigenvalue stability Lemma 1 used in Theorem 2.2."},{"cited_title":"Complex dynamic behaviors of a discrete time predator prey system","cited_arxiv_id":null,"evidence_quote":"Provides the integer-order discrete predator-prey baseline that the paper says its fractional results improve on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the bifurcation and normal-form theory behind the Neimark-Sacker and flip bifurcation derivations."}],"review_version":1}