{"id":"d009f80e-221b-45f1-aee7-1462b420dcc4","arxiv_id":"2505.00128","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Routing functions decompose the parameter space of polynomial ecological ODEs into connected regions of constant stable steady states, and applied to a coral-bacteria model they reveal bistability and no-stable-state regimes.","lead":"This paper applies a computational algebraic geometry method called routing functions to split a model's parameter space into regions where the same steady states are stable. It maps the stability landscape of a coral-bacteria symbiosis model and finds regions where no steady state is stable, which the authors interpret as limit cycles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The coral landscape depends on an unshown 31-factor boundary polynomial; without a certificate, all labels are conditional.","rationale":"The routing-function framework is a plausible application of the authors' earlier theory, and the Levins-Culver reproduction provides a useful sanity check. The central correctness risk is not the mathematical framework itself but the completeness and reproducibility of the symbolic boundary computation for the coral model. This is exactly the reader's weakest assumption, and it is load-bearing because every landscape label, including the white no-stable-state region, is a consequence of that boundary. The limit-cycle overclaim is real but secondary: it is an interpretation of the white region rather than a necessary part of the decomposition method. The paper's own conclusion already softens it, and the abstract should be corrected accordingly. Since the reader already marked the paper CONDITIONAL with moderate confidence, my concern does not move the verdict; it reinforces the need for the certificate and repository access before the coral claims can be independently checked.","tokens_in":14059,"tokens_out":6034,"duration_ms":73543,"concrete_test":"Release the Macaulay2 script with the exact generators of If and its five primary components. Independently recompute, for each component, the squarefree elimination ideal (If + <R>) ∩ C[b, βz, ˜b, ˜d] using a different term order or a different computer algebra system, and verify that the product of the 31 reported irreducible factors equals the squarefree generator up to a constant. Then, for each factor g_i, pick a generic point with g_i = 0 and positive coordinates and check that at least one Routh polynomial or stable-steady-state count changes sign across it; this certifies that no factor is spurious and that no real boundary factor is missing. If this certificate passes, the coral landscape labels stand; if it fails, the conclusions of Section 4 are unverified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step for the coral case study is the claim that the total boundary is exactly the zero set of a principal ideal <g> with 31 irreducible factors. This boundary is the input to Theorem 2.1: if g omits any hypersurface across which the number or type of stable steady states changes, the computed connected components can merge regions with different landscape labels, and labeling from one routing point per component will then be wrong. The paper does not display g, does not list the five primary-component elimination ideals, provides no certificate for the elimination or factorization, and cites a GitHub repository without a URL or commit. The reader therefore cannot distinguish a complete boundary computation from a partial one. A secondary overclaim is the inference from 'no stable steady state' to 'the model must have a limit cycle' in a three-dimensional system; without ruling out chaotic or other non-periodic attractors this does not follow, and the conclusion's softer 'limit cycles or more complex attractors' is the defensible version. The abstract's 'regions supporting limit cycles' is not established by the routing-function computation alone.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper applies the routing-function method of [7] to ecological ODE models, decomposing the positive parameter space minus a boundary hypersurface arrangement into connected components on which the number and type of stable steady states is constant. The boundaries considered are the singular, Routh-Hurwitz, and coordinate boundaries, each obtained via elimination ideals. The method is first demonstrated on the Levins-Culver competition-colonization model, reproducing the known coexistence-stability condition. It is then applied to a coral-bacteria model with parameters fixed to d = gamma_y = gamma_z = 1 and beta_y = 5; after computing a claimed total boundary with 31 irreducible factors, the authors slice at b = 2 and b = 0.5 and at eight values of beta_z, labeling each two-dimensional component by the stable steady states found at a routing point. The abstract claims that the method reveals regions supporting limit cycles, while Section 4.1.1 states that a white region with no stable steady states must have a limit cycle; the conclusion more cautiously says 'limit cycles or more complex attractors.'","tokens_in":14235,"tokens_out":7726,"duration_ms":80009,"significance":"The proposed pipeline is a worthwhile use of recent computational real algebraic geometry in mathematical biology: it gives a concrete, algorithmic way to obtain stability landscapes, and the Levins-Culver section is a clean sanity check in which the known condition beta_z > beta_y (beta_y + gamma_z - gamma_y)/gamma_y and beta_y > gamma_y is recovered without case-by-case analysis. The paper also provides a useful pedagogical summary of the three boundary types and of routing functions, and it names concrete software (HypersurfaceRegions.jl, HomotopyContinuation.jl, Macaulay2) that makes the computations reproducible in principle. However, the coral-landscape results are conditional on a boundary computation that is not verifiable from the manuscript, and one headline inference (no stable steady state implies a limit cycle) is not justified for a three-dimensional system. With added certificates and corrected language, the contribution would be useful to the q-bio audience.","major_comments":[{"comment":"The coral analysis rests entirely on the claim that eliminating x, y, z from the equilibrium ideal plus the Routh-Hurwitz conditions yields a principal ideal <g> with 31 irreducible factors, whose degrees are listed but whose factors are not shown. This is load-bearing for Figures 3 and 4 and for Theorem 2.1: if g omits any hypersurface across which the number or type of stable steady states changes, connected components of the complement can merge regions with different labels; a spurious factor would create artificial boundaries. The manuscript does not display g, does not list the five primary-component elimination ideals, and gives no certificate (e.g., a Macaulay2 script and transcript, or a repository URL with a commit hash). Please supply the elimination data in a supplement or make the repository identifiable, or explicitly mark the coral landscape results as conditional on this computation.","section":"Section 4, total boundary computation"},{"comment":"The white region in Figure 3(e) is described as a region where 'none of the steady states are stable and the model must have a limit cycle.' In a three-dimensional autonomous system, absence of stable steady states does not imply existence of a limit cycle; chaotic attractors, heteroclinic cycles, or unbounded dynamics are not excluded by the routing-function computation, which, even assuming the boundary is exact, only certifies that no stable steady state exists in that component and does not identify the attractor. The conclusion's phrase 'limit cycles or more complex attractors' is the defensible version. Please revise the abstract and Section 4.1.1 accordingly, or supply an additional argument (e.g., boundedness plus a Poincare-Bendixson-type analysis) that actually proves a periodic orbit.","section":"Section 4.1.1 and Abstract"},{"comment":"The manuscript does not state how the 31-factor boundary arrangement is restricted to each fixed (b, beta_z) slice before the two-dimensional routing function is constructed, and it does not list the active factors for the 16 slices. This information is needed to reproduce Figures 3 and 4 and to check that slice-dependent factors (for example, factors that become constant on a slice) are handled correctly. Please describe the restriction procedure and provide the restricted boundary data for each slice, either in the text or in a supplement.","section":"Section 4, slice decomposition"}],"minor_comments":[{"comment":"The symbol beta_y is overloaded: it denotes a Routh polynomial in the stability criterion and also a model parameter later in the paper. Please rename the Routh polynomials to avoid confusion.","section":"Section 2.1"},{"comment":"The primary decomposition notation 'If = T4i=1 J(i)f' is ambiguous; it should be an intersection, e.g., If = intersection_{i=1}^4 J_f^(i).","section":"Section 3"},{"comment":"There are several typos and infelicities: 'paramter landscape' in Section 2, 'for for' before 'the Levins-Culver model' in Section 4, and 'the large the number of parameters' in Section 4. These should be corrected.","section":"Various"},{"comment":"The statement that the computations 'can be found in a GitHub repository' should include a URL, a persistent identifier, and a commit hash or version tag so that the reported results are actually reproducible.","section":"Section 4"},{"comment":"The phrase 'as beta_z crosses four' and the description of an 'abrupt' landscape change are based on two sampled slices on either side of beta_z = 4; please phrase these observations as being observed at the sampled slices rather than as a proven characterization for all intermediate values.","section":"Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper's reliance on [7], co-authored by two of the present authors, is not a circularity problem in my view because [7] is a published algorithm with its own proofs and software. The main gate for the coral section is verification of the 31-factor boundary; if the authors cannot supply a certificate, the coral results should be reframed as exploratory. The limit-cycle language is an overclaim and can be fixed by rewording. I see no fundamental obstacle to publication after these revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Rough take: this is a faithful application of routing functions to ecology, and the coral-bacteria stability maps are new. The exposition of the three boundary types is clear, and the Levins-Culver reproduction is a good sanity check. But the abstract overclaims limit cycles, and the coral results rest on a boundary polynomial g that is never shown or certified, so the maps should be read as conditional on that computation being complete. What is actually new: the systematic package of singular, Routh-Hurwitz, and coordinate boundaries applied to a three-species ecological model, plus the first stability landscape decomposition for the tripartite symbiosis model. The 16 slices, including the bistability regions and the white no-stable-state region, are not in the cited literature. That is real output. Where it gets soft: Section 4.1.1 says a region with no stable steady states must have a limit cycle. In three dimensions that does not follow without ruling out chaotic or other non-periodic attractors. The conclusion uses the safer phrase \"limit cycles or more complex attractors,\" so the abstract's \"regions supporting limit cycles\" is not what the computation established. Fix the wording. The bigger issue is the unshown total boundary. The routing function graph is only as good as the hypersurface arrangement. If any of the 31 irreducible factors is missing or spurious, connected components merge or split and every region label can change. The paper gives degree counts but does not display g, list the elimination ideals for the five primary components, or provide a certificate. The GitHub repository is mentioned but no URL or commit hash is given. So a referee cannot verify the load-bearing step. This is fixable by making the code and data available, but until then the coral landscape is conditional, not final. The citation pattern is fine: [7] is the source of routing functions, published with independently available software, and the authors are explicit about the lineage. The self-citation is appropriate. No invented entities or fitted constants. Who this is for: mathematical biologists who want to map stability regions in polynomial or rational ODE models, and applied algebraic geometers looking for new application domains. It is not a theoretical advance beyond [7], but it is a well-motivated application with useful results. Recommendation: this deserves a serious referee. Conditional accept, with the expectation that the authors supply the missing boundary data and correct the limit-cycle language.","headline":"A genuinely useful application of routing functions to ecological stability landscapes, but the limit-cycle claim is too strong and the unshown 31-factor boundary polynomial makes the coral maps conditional.","tokens_in":706,"tokens_out":940,"would_cite":true,"duration_ms":31938,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92D40","37N25","14P10","34D20"],"pacs":[],"model":"deepseek-v4-flash","headline":"Routing functions decompose ecological parameter space into regions with a fixed set of stable steady states.","keywords":["stability landscape","routing function","real algebraic geometry","semi-algebraic sets","Routh-Hurwitz boundary","multistability","coral-bacteria symbiosis","limit cycles"],"falsifier":"Recompute the elimination ideal for the coral model with an independent symbolic computation and compare the 31 irreducible factors; then at representative parameter values inside each computed region, numerically integrate the ODEs and check that the predicted stable steady states are present and that in the white region no steady state is stable and orbits do not converge to equilibrium.","tokens_in":13789,"feed_emoji":"🗺️","tokens_out":9117,"duration_ms":84036,"temperature":0.7,"pith_summary":"The paper aims to turn the stability question for ecological ODE models into a symbolic geometry problem. It shows that outside a boundary built from singular, Routh-Hurwitz, and coordinate conditions, each connected region of parameter space has a constant number and type of stable steady states, and routing functions compute those regions. Applied to a coral-bacteria symbiosis model with several rates fixed, the method maps where coexistence, pathogen exclusion, and bistability are stable, and it identifies a region with no stable steady state, implying sustained oscillations. If the framework generalizes as claimed, ecologists could read collapse and recovery thresholds directly off such maps.","feed_headline":"Algebraic map exposes stable states and limit cycles of coral model","feed_subtitle":"Each region of parameter space holds the same stable states, so collapse and recovery thresholds can be read off a map.","key_machinery":"The key object is a routing function: for boundary polynomials $g_1,\\dots,g_m$, the rational map $r_c(a)=g_1(a)\\cdots g_m(a)/(1+\\sum_i(a_i-c_i)^2)^D$ with $D$ large enough, whose nondegenerate critical points (routing points) and gradient-ascent paths form a graph whose connected components match the connected components of $\\mathbb{R}^k_{>0}\\setminus B$. The boundary $B$ is assembled from singular, Routh-Hurwitz, and coordinate boundaries, each computed as an elimination ideal; this mechanism turns the stability question into a finite graph computation.","core_discovery":"The central claim is that routing functions—rational functions built from the product of the boundary polynomials—provide a computational route to the exact connected components of parameter space on which steady-state feasibility and stability are constant. The paper relies on a theorem stating that the graph of routing points and gradient connections is in bijection with the connected components of the complement of the total boundary, and it instantiates the method for the classical two-species competition-colonization model and for a three-species coral-bacteria model. In the coral model, after fixing death rates and one colonization rate, the total boundary is a principal ideal with 31 irreducible factors, and varying the remaining two parameters yields stability landscapes labeled by the stable equilibria, including bistable regions and a white region where no equilibrium is stable, so the dynamics there must be non-stationary.","pith_inferences":["I would expect the same routing-function decomposition to work when initial conditions are treated as parameters, which would turn questions of basin size and separatrix geometry into the same graph computation; the paper notes this possibility but does not develop it.","Because the no-steady-state region is identified by algebraic elimination rather than by simulation, the method offers a certificate-style indication of oscillations for parameter sets, something bifurcation continuation typically only suggests.","A natural stress test is to perturb the fixed rates ($d=\\gamma_y=\\gamma_z=1$, $\\beta_y=5$); the qualitative landscape, including the position of the limit-cycle window, may be sensitive to those choices, so the biological conclusions are conditional on the chosen slice."],"forward_implications":["The same boundary construction can be applied to any polynomial or rational ODE model, giving a complete stability-landscape map rather than a sampled approximation.","For the coral-bacteria model with the chosen fixed parameters, crossing any computed boundary changes the set of stable states, so the maps identify thresholds for pathogen exclusion versus coexistence.","The white region with no stable steady state implies sustained oscillations or more complex attractors, so the method can detect limit cycles without integrating trajectories.","Bistable regions mean that state perturbations, not just parameter changes, can switch the system between coral-only and coexistence, which is directly relevant to whether a coral population recovers after a shock."],"supporting_citations":[{"why":"Supplies the routing-function theorem and graph algorithm that justify the connected-components correspondence.","marker":"[7]"},{"why":"Introduces the coral-bacteria symbiosis model whose stability landscape is the paper's main case study.","marker":"[11]"},{"why":"Provides the Routh-Hurwitz criterion used to define the stability boundary.","marker":"[15]"},{"why":"Defines the classical competition-colonization model used as a validation example.","marker":"[18]"},{"why":"Origin of routing functions for connectivity queries on which the method builds.","marker":"[13]"},{"why":"Provides the computational routine for decomposing hypersurface arrangements used in the examples.","marker":"[3]"}],"fun_headline_variants":["Algebraic routing maps stability landscapes of ecological models","Coral symbiosis stability decoded via algebraic routing functions","Exact stability maps of ecosystems from polynomial algebra","Algebraic geometry reveals limit cycles in coral model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the computed 31-factor polynomial really is the complete total boundary for the coral model; if any boundary piece is missing or extraneous, every connected region and its stability label, including the claimed no-steady-state region, could change.","fun_headline_variants_meta":{"raw":{"variants":["Algebraic routing maps stability landscapes of ecological models","Coral symbiosis stability decoded via algebraic routing functions","Exact stability maps of ecosystems from polynomial algebra","Algebraic geometry reveals limit cycles in coral model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000545,"raw_usage":{"total_tokens":2589,"prompt_tokens":909,"completion_tokens":1680,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":1620}},"tokens_in":525,"tokens_out":1680,"duration_ms":13270,"temperature":1.0,"reasoning_tokens":1620,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:50:45.395593+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the elimination ideal for the coral model with an independent symbolic computation and compare the 31 irreducible factors; then at representative parameter values inside each computed region, numerically integrate the ODEs and check that the predicted stable steady states are present and that in the white region no steady state is stable and orbits do not converge to equilibrium.","supporting_citations":[{"cited_title":"Cummings, J","cited_arxiv_id":null,"evidence_quote":"Supplies the routing-function theorem and graph algorithm that justify the connected-components correspondence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the coral-bacteria symbiosis model whose stability landscape is the paper's main case study."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Routh-Hurwitz criterion used to define the stability boundary."},{"cited_title":"Levins and D","cited_arxiv_id":null,"evidence_quote":"Defines the classical competition-colonization model used as a validation example."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Origin of routing functions for connectivity queries on which the method builds."},{"cited_title":"Breiding, B","cited_arxiv_id":null,"evidence_quote":"Provides the computational routine for decomposing hypersurface arrangements used in the examples."}],"review_version":1}