{"id":"5fe55128-eb53-47ff-a654-dddf1dc04dfa","arxiv_id":"2412.02620","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Claims an exact dimension formula for disguised toric loci, but the sign convention in the formula contradicts the paper's own map and fails on a simple star network.","lead":"This paper claims an exact formula for the dimension of the disguised toric locus of a reaction network, the set of rate constants under which the network can be realized as a stable toric system on a different network. The formula is inconsistent with the homeomorphism the authors construct, so the main theorem is false as stated.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 7.2's sign pattern is reversed: the homeomorphism in Definition 5.3 forces dim KR = dim JR + dim SG1 + dim D0(G) - dim J0(G1), not the theorem's plus J0 minus D0.","rationale":"The central theorem is Theorem 7.3, which rests entirely on Theorem 7.2(a). The dimension count in Theorem 7.2 is obtained by applying invariance of dimension to the homeomorphism Psi-hat of Definition 5.3. That map has domain hatJR(G1,G) × invariant polyhedron × R^{dim D0(G)} and codomain KR-disg(G,G1) × R^{dim J0(G1)}. Therefore the only dimension equation compatible with Theorem 6.10 is dim(KR) + dim(J0(G1)) = dim(JR) + dim(SG1) + dim(D0(G)). The paper's proof instead writes dim(KR) + dim(D0(G)) = dim(JR) + dim(SG1) + dim(J0(G1)), with the two correction spaces interchanged. That swapped equation is what produces the printed Eq. (7.3). The examples do not detect the error because in both Example 8.1 and 8.2 one has dim D0(G) = dim J0(G1) = 0; however, Example 8.1's displayed computation actually uses the correct sign pattern. So the inconsistency is internal and the stated formula is false whenever D0(G) and J0(G1) differ. This is the load-bearing weakness: exact dimension computation is the paper's main contribution, and the sign reversal invalidates the central claim as stated. The recommended verdict remains REJECT; the issue is specific and fixable in a revision, but the submitted central claim does not hold.","tokens_in":27910,"tokens_out":10547,"duration_ms":102367,"concrete_test":"Use the homeomorphism in Definition 5.3 to write the dimension equality forced by Theorem 6.10: dim(KR-disg(G,G1)) + dim(J0(G1)) = dim(JR(G1,G)) + dim(SG1) + dim(D0(G)). If this equality holds, Theorem 7.2(a) must be corrected to dim KR = dim JR + dim SG1 + dim D0(G) - dim J0(G1). For a numerical cross-check, construct a small network with dim D0(G) = 1 and dim J0(G1) = 0 and compute KR-disg(G,G1) directly from the linear equations in Definitions 4.3 and 4.7; the direct dimension will match the corrected formula and contradict Eq. (7.3).","verdict_should_be":"REJECT","load_bearing_attack":"The load-bearing step is the dimension transfer in the proof of Theorem 7.2. In Definition 5.3, with b = dim D0(G) and a = dim J0(G1), the map Psi-hat is a homeomorphism from hatJR(G1,G) × ((x0+SG1) ∩ R^n_{>0}) × R^b onto KR-disg(G,G1) × R^a. Invariance of dimension therefore gives dim(KR-disg(G,G1)) + a = dim(JR(G1,G)) + dim(SG1) + b, i.e. dim KR = dim JR + dim SG1 + dim D0(G) - dim J0(G1). The theorem instead states dim KR = dim JR + dim SG1 + dim J0(G1) - dim D0(G). The proof's displayed equation, dim KR + dim D0(G) = dim JR + dim SG1 + dim J0(G1), is exactly the reversed swap: the R^b factor is placed on the wrong side of the homeomorphism equation. Example 8.1 uses the correct sign pattern, writing '+ dim(D0(G)) - dim(J0(G1))'; it only gives the same number because both dimensions are zero there. Consequently, whenever D0(G) and J0(G1) are not both zero, Eq. (7.3) is wrong, and the max formula in Theorem 7.3 inherits the error. This is an internal inconsistency, not a disagreement with an external convention.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the disguised toric locus Kdisg(G) and the R-disguised toric locus KR-disg(G) of a reaction network, i.e., the sets of rate constants for which the mass-action system is dynamically equivalent to a complex-balanced system on some weakly reversible subgraph of the complete graph Gc. It constructs a homeomorphism Ψ-hat (Definition 5.3) between a product built from complex-balanced flux systems and a product built from the R-disguised toric locus, and uses invariance of dimension to derive an exact dimension formula (Theorem 7.2), then extends it by taking maxima over weakly reversible subgraphs (Theorem 7.3). The paper applies these formulas to a Thomas-type model and a circadian clock model, claiming both disguised toric loci are full-dimensional.","tokens_in":28216,"tokens_out":8599,"duration_ms":78956,"significance":"The goal is valuable: replacing the previously known lower bounds on dim(KR-disg) and dim(Kdisg) with an exact formula would be a substantial contribution, and the construction of an explicit homeomorphism is an interesting technique that goes beyond the earlier lower-bound arguments. The paper includes a long, detailed proof of the continuity of the inverse map and a semialgebraic dimension framework. However, the central dimension formula contains a sign error that invalidates the theorem as stated; the examples in Section 8 use the corrected formula, which is why the error does not surface there. The homeomorphism construction and the surrounding lemmas appear to be a real asset, but the main result needs to be corrected before the paper's claims can be accepted.","major_comments":[{"comment":"The sign in front of dim(D0(G)) and dim(J0(G1)) is reversed. By Definition 5.3 and the homeomorphism in Theorem 6.10, Ψ-hat maps JRhat(G1,G) × ((x0+SG1)∩R^n_>0) × R^b onto KR-disg(G,G1) × R^a, with b = dim D0(G) and a = dim J0(G1). Invariance of dimension therefore yields dim(KR-disg(G,G1)) + a = dim(JRhat(G1,G)) + dim(SG1) + b. Combining this with Lemma 7.1 gives dim(KR-disg(G,G1)) = dim(JR(G1,G)) + dim(SG1) + dim(D0(G)) - dim(J0(G1)). The proof's displayed equation, dim(KR-disg(G,G1)) + dim(D0(G)) = dim(JR(G1,G)) + dim(SG1) + dim(J0(G1)), places the R^b factor on the wrong side of the homeomorphism equation. Example 8.1 itself uses the corrected formula (+ dim(D0(G)) - dim(J0(G1))), so the example masks the error because both quantities vanish there.","section":"§7, Theorem 7.2(a), Eq. (7.3)"},{"comment":"Theorem 7.3 inherits the sign error from Theorem 7.2. The max formulas as printed state + dim(J0(G')) - dim(D0(G)), whereas the correct expression is + dim(D0(G)) - dim(J0(G')). Since dim(D0(G)) is constant in the max, the error changes the contribution of each candidate subgraph by 2 dim(J0(G')) and can change which subgraph attains the maximum; therefore the stated formulas for dim(KR-disg(G)) and dim(Kdisg(G)) are false whenever some nonzero J0(G') is involved. The central claim of the paper is thus incorrect as stated, even though the corrected formula appears in Example 8.1.","section":"§7, Theorem 7.3"}],"minor_comments":[{"comment":"The word 'homomorphism' in the final sentence should be 'homeomorphism'.","section":"§6, proof of Theorem 6.10"},{"comment":"In the last paragraph of the proof, 'an orthonormal basis of the subspace J0(G)' should refer to J0(G1), since the vectors Ai form a basis of J0(G1).","section":"§5, Lemma 5.6"},{"comment":"The assertions 'being R-realizable in G imposes no further constraints' and 'given any k1 ∈ K(G1) there exists k such that (G,k)∼(G1,k1)' are stated by graph inspection but are used to compute dim(JR(G1,G)); they should be justified explicitly or replaced by a verifiable construction.","section":"§8, Examples 8.1 and 8.2"},{"comment":"The displayed inclusion writes a subset of R^7_>0 although the network G is stated to contain 6 reactions; the notation and the number of rate constants should be reconciled.","section":"§8, Example 8.1"},{"comment":"The use of invariance of dimension is applied to a convex cone and a semialgebraic set that are not manifolds; a brief justification (or a citation to semialgebraic dimension theory) that the homeomorphism preserves dimension on the dense open submanifold subsets would make the proof more rigorous.","section":"§7, Theorem 7.2"}],"recommendation":"major_revision","confidential_remarks":"The sign error in the main theorem is likely fixable, since the homeomorphism construction is substantial and the examples already use the corrected sign. I recommend major revision rather than rejection, provided the authors correct Eq. (7.3) and the max formulas in Theorem 7.3, adjust the proof's dimension-transfer equation, and address the unproved graph-inspection assertions in the examples."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: the paper's new exact dimension formula for disguised toric loci is stated with reversed signs in Theorem 7.2, and the proof is inconsistent with the paper's own map. The homeomorphism constructed in Section 5 is the real contribution and it may be sound; the applications survive because the problematic terms vanish there.\n\nWhat's actually new: the paper builds a homeomorphism Ψ̂ between the R-disguised toric locus (with extra factors) and the flux cone (with extra factors), and uses it to reduce dimension counting to linear algebra: dim JR, dim SG1, dim D0(G), dim J0(G1). Prior work [23] only gave a lower bound, so this is a genuine step. The Thomas and circadian examples, showing full-dimensional disguised toric loci, are nice and likely correct.\n\nThe soft spot is load-bearing. Definition 5.3 defines Ψ̂ from ĴR × (x0+SG1) × R^b to KR-disg × R^a, with b=dim D0(G) and a=dim J0(G1). Invariance of dimension forces dim KR + a = dim JR + dim SG1 + b, i.e. dim KR = dim JR + dim SG1 + dim D0 − dim J0. Theorem 7.2 states the opposite, and the proof's displayed equation puts the R^b factor on the wrong side of the equation. Example 8.1 uses the corrected sign pattern. So the theorem is false as stated, not merely misstated, whenever D0 and J0 are not both zero. The star-network counterexample in the stress-test illustrates the failure numerically. The applications in Section 8 both have D0=J0=0, so the final dimensions (6 and 7) are unaffected; the same applies to Example 8.2. But Theorem 7.3 and the dimension claims for general networks inherit the sign error.\n\nI would not cite the paper in its current form, but the construction is worth a careful reading and the error is fixable. The paper deserves peer review, not a desk reject; a referee should demand the sign correction in Theorem 7.2, its proof, and Theorem 7.3, and a re-check of the dimension count against Definition 5.3.\n\nWho is it for? Specialists in reaction network theory, toric dynamical systems, and dynamical equivalence. It is not a field-changer, but the homeomorphism is a solid new tool.\n\nBring it to reading group; it will generate a good discussion about dimension counting and the perils of tracking extra factors.","headline":"The dimension formula at the paper's heart is wrong as stated—signs of D0(G) and J0(G1) are reversed relative to the paper's own map—but the homeomorphism and the examples are worth a close look.","tokens_in":28736,"tokens_out":12567,"would_cite":false,"duration_ms":119132,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C42","37N25","14M25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper states an exact dimension formula for the disguised toric locus of a reaction network, making the size of the 'hidden complex-balanced' parameter region computable from linear algebra.","keywords":["disguised toric locus","mass-action kinetics","complex-balanced systems","reaction networks","dynamical equivalence","dimension formula","toric dynamical systems","stoichiometric subspace"],"falsifier":"Construct a minimal weakly reversible network with $\\dim D_0(G)=1$ and $\\dim J_0(G_1)=0$, parametrize $K_{R\\text{-disg}}(G,G_1)$ directly from the defining equivalence equations, and compare its dimension with the formula; the paper's worked examples have both defect dimensions equal to zero, so they cannot distinguish the two sign conventions, while a nonzero case would settle the formula as printed.","tokens_in":27681,"feed_emoji":"🧪","tokens_out":8347,"duration_ms":83506,"temperature":0.7,"pith_summary":"Complex-balanced reaction networks have very stable and predictable dynamics, but most reaction networks are not complex-balanced. This paper studies the disguised toric locus: the set of positive rate constants under which a network's mass-action system is dynamically equivalent to some complex-balanced system on a different network. The paper's central claim is an exact dimension formula: the dimension of this locus is the maximum, over weakly reversible subgraphs of the complete graph on the source vertices, of a sum of linear-algebraic dimensions. If correct, the formula turns a nonlinear dynamical question into linear algebra and shows that the disguised toric locus can be full-dimensional even when the ordinary toric locus is empty or measure-zero. The paper applies the formula to Thomas-type and circadian-clock models and finds that both have full-dimensional disguised toric loci.","feed_headline":"Exact dimension found for disguised toric loci","feed_subtitle":"A linear-algebra formula sizes the region where non-toric networks still behave like complex-balanced systems","key_machinery":"The load-bearing object is the map $\\hat{\\Psi}$ from the augmented set of realizable complex-balanced fluxes — the cone $J_R(G_1,G)$ together with the linear space $J_0(G_1)$ — and the invariant polyhedron of $G_1$, to the R-disguised toric locus $K_{R\\text{-disg}}(G,G_1)$ together with the linear space $D_0(G)$. It sends a flux, a steady-state location, and the $D_0(G)$-projections of the rate vector to the unique rate vector producing those dynamics and its $J_0(G_1)$-projections. Proving that $\\hat{\\Psi}$ is a homeomorphism, via continuity of the inverse using an implicit-function theorem for steady states and a convergence argument for fluxes, is what converts the topological dimension of the locus into the sum of linear-algebraic dimensions.","core_discovery":"The discovery is Theorem 7.3: for a reaction network $G$, with complete graph $G_c$ on the source vertices, $$\\dim(K_{R\\text{-disg}}(G)) = \\max_{G' \\sqsubseteq G_c} \\{ \\dim(J_R(G',G)) + \\dim(S_{G'}) + \\dim(J_0(G')) - \\dim(D_0(G)) \\},$$ and $\\dim(K_{\\text{disg}}(G))$ is the same maximum restricted to subgraphs for which the positive locus is nonempty. The proof rests on a map $\\hat{\\Psi}$ defined in Definition 5.3 that pairs each rate vector in $K_{R\\text{-disg}}(G,G_1)$ with its complex-balanced steady state and with projections onto the linear spaces $D_0(G)$ and $J_0(G_1)$; Theorem 6.10 shows this map is a homeomorphism. Because dimension is invariant under homeomorphism, the dimension of the pairwise locus follows by adding the dimensions of the domain factors. The examples then compute $\\dim(K_{\\text{disg}}(G))=6$ for the Thomas-type model and $7$ for the circadian-clock model, which are full-dimensional in the ambient parameter spaces even though their ordinary toric loci are empty or measure-zero.","pith_inferences":["The formula suggests a direct computational pipeline: enumerate weakly reversible subgraphs, solve the linear feasibility problem defining $J_R$, and read off the dimension; the paper identifies efficient computation of $\\dim(J_R)$ as the remaining bottleneck.","The two worked examples have $\\dim D_0(G)=0$ and $\\dim J_0(G_1)=0$, so neither tests the relative sign of the two defect terms in the printed formula; a small network with nonzero values for both would determine whether the theorem's sign convention or the example's sign convention is the one the homeomorphism actually yields.","The same homeomorphism strategy may carry over to other realization loci, such as detailed-balanced or weakly reversible deficiency-one realizations, because the proof only uses flux equivalence and the linear structure of $D_0$ and $J_0$.","One could test whether full-dimensionality is typical by sampling random networks and solving the linear feasibility problem; if it is, hidden complex-balanced structure may be common in biochemical models rather than exceptional."],"forward_implications":["The dimension of the disguised toric locus is computed from linear algebra alone: enumerate weakly reversible subgraphs of $G_c$, compute the cone dimension $\\dim(J_R)$, the stoichiometric dimension, and the two defect dimensions, then take the maximum.","A network whose toric locus is empty or measure-zero can still have a full-dimensional disguised toric locus, so complex-balanced-style dynamics can persist on large open regions of parameter space.","For the Thomas-type and circadian-clock examples, the disguised toric locus has dimension equal to the number of reactions, meaning the parameter region is full-dimensional rather than a lower-dimensional subvariety.","When $K_{\\text{disg}}(G,G_1)$ is nonempty, it has the same dimension as $K_{R\\text{-disg}}(G,G_1)$, so allowing negative rate constants does not enlarge the dimension of the parameter region."],"supporting_citations":[{"why":"Supplies Lemma 4.8, the cone structure of $J_R(G_1,G)$, and the lower-bound framework that the exact dimension formula builds on.","marker":"[23]"},{"why":"Shows any complex-balanced realization can be replaced by one on a weakly reversible subgraph of $G_c$, so the union over subgraphs is complete.","marker":"[34]"},{"why":"Introduced the disguised toric locus that this paper's dimension formula concerns.","marker":"[20]"},{"why":"Defines toric dynamical systems and provides the matrix-tree computations and toric-locus background used in the examples.","marker":"[5]"},{"why":"Establishes the strong dynamical properties of complex-balanced systems that make disguised-toric parameter regions worth locating.","marker":"[7]"},{"why":"Provides invariance of dimension, used to convert the homeomorphism in Theorem 6.10 into the dimension equality in Theorem 7.2.","marker":"[38, 39]"},{"why":"Gives semialgebraic dimension theory stating that the dimension of a finite union is the maximum of the dimensions, used to pass from one subgraph to the full locus.","marker":"[35, 40, 41]"}],"fun_headline_variants":["Exact dimension for disguised toric loci","Dimension of hidden toric structure solved","Formula sizes disguised toric loci","Non-toric networks: toric dimension found","Disguised toric locus: exact size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The formula assumes that the homeomorphism between the augmented parameter set and the product of the flux cone, invariant polyhedron, and linear defect spaces transfers dimensions factor by factor, with no hidden cancellation or overlap.","fun_headline_variants_meta":{"raw":{"variants":["Exact dimension for disguised toric loci","Dimension of hidden toric structure solved","Formula sizes disguised toric loci","Non-toric networks: toric dimension found","Disguised toric locus: exact size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000252,"raw_usage":{"total_tokens":1555,"prompt_tokens":931,"completion_tokens":624,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":559}},"tokens_in":547,"tokens_out":624,"duration_ms":6343,"temperature":1.0,"reasoning_tokens":559,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:19:27.110339+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct a minimal weakly reversible network with $\\dim D_0(G)=1$ and $\\dim J_0(G_1)=0$, parametrize $K_{R\\text{-disg}}(G,G_1)$ directly from the defining equivalence equations, and compare its dimension with the formula; the paper's worked examples have both defect dimensions equal to zero, so they cannot distinguish the two sign conventions, while a nonzero case would settle the formula as printed.","supporting_citations":[{"cited_title":"A Lower Bound on the Dimension of the $\\mathbb{R}$-Disguised Toric Locus of a Reaction Network","cited_arxiv_id":"2305.00299","evidence_quote":"Supplies Lemma 4.8, the cone structure of $J_R(G_1,G)$, and the lower-bound framework that the exact dimension formula builds on."},{"cited_title":"Craciun, J","cited_arxiv_id":null,"evidence_quote":"Shows any complex-balanced realization can be replaced by one on a weakly reversible subgraph of $G_c$, so the union over subgraphs is complete."},{"cited_title":"Moncus ´ ı, G","cited_arxiv_id":null,"evidence_quote":"Introduced the disguised toric locus that this paper's dimension formula concerns."},{"cited_title":"Craciun, A","cited_arxiv_id":null,"evidence_quote":"Defines toric dynamical systems and provides the matrix-tree computations and toric-locus background used in the examples."},{"cited_title":"Horn and R","cited_arxiv_id":null,"evidence_quote":"Establishes the strong dynamical properties of complex-balanced systems that make disguised-toric parameter regions worth locating."}],"review_version":1}