{"id":"26f73050-9990-4571-ada7-f4311fd58785","arxiv_id":"2607.20958","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Asymmetric initial connectivity in bipartite evolving networks enhances robustness, shifting the growth/extinction boundary and producing a re-entrant growing phase absent in unstructured systems.","lead":"This paper simulates evolving networks whose nodes split into two groups, with interactions only between groups, and shows that when new nodes in the two groups start with different numbers of links, the network can keep growing even though a standard calculation says it should collapse. This suggests a general robustness mechanism for bipartite systems such as plant-pollinator communities, microbial cross-feeding, and trade networks, and explains why asymmetry may be protect","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (24)-(27) mis-derive the shifted extinction cascades: the geometric-series denominator must change, so the claimed quantitative explanation of the re-entrant transition is unsupported.","rationale":"The reader's verdict was CONDITIONAL, and I agree with that overall assessment. However, the most load-bearing concern is not the empirical extinction law E(k) or the unproven link-balance condition (Eq. 19), but a concrete algebraic error in the cascade equations (24)-(27) that invalidates the claimed theoretical reproduction of the phase diagram. The central claim includes the statement that an extended mean-field analysis identifies the origin of the robustness enhancement and re-entrant transition. If Eqs. (24)-(27) are incorrect, then the explanatory mechanism—degree shift producing extinction suppression and low-degree extinction bias—is not quantitatively established. The simulation results may still be valid, but the paper's theoretical contribution is seriously weakened. The verdict remains CONDITIONAL because the empirical phase diagram could be correct, but the paper must correct the derivation and either provide corrected curves or openly mark those predictions as phenomenological. The concrete test is straightforward: re-derive the series and compare with Fig. 9. This is a higher-priority check than those proposed by the reader, though the reader's concerns about data availability and threshold choice remain relevant.","tokens_in":9528,"tokens_out":15017,"duration_ms":133943,"concrete_test":"Independently re-derive Ñ_αα, Ñ_αβ, Ñ_βα, Ñ_ββ by substituting Ĩα and Ĩβ into the geometric series (e.g., Ñ_αα = ĨαĨβ + Ĩα²Ĩβ² + ...). If the closed forms differ from Eqs. (24)-(27), recompute the dash-dotted curves in Fig. 9 using the corrected formulas. Then check whether the corrected n_αα still falls below the mean-field estimate; if the overestimation disappears, the separate E_i/E_e correction in §III.C.3 is unnecessary and the proposed 'strong reinforcement effect' is not supported by the theory.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theoretical claim—that the degree-shift mechanism quantitatively reproduces the phase boundary and re-entrant transition—rests on the cascade equations in §III.C.2. However, substituting the shifted impact rates Ĩα=(m_α/k_α)I_α and Ĩβ=(m_β/k_β)I_β (Eqs. 22-23) into the original geometric series does not yield Eqs. (24)-(27). For example, Ñ_αα should be ĨαĨβ/(1−ĨαĨβ) = [(m_α m_β/(k_α k_β)) I_α I_β] / [1 − (m_α m_β/(k_α k_β)) I_α I_β], not (m_α/k_α) N_αα(k_α,k_β). The denominator is modified by the same product factor, and the multiplicative factor is not a single m_χ/k_χ. The same error affects Eqs. (25)-(27). Consequently, the dash-dotted curves in Fig. 9 are not the correct solutions to the stated model; the reported overestimation of n_αα may be an artifact of this algebraic mistake, not a genuine physical effect requiring the separate E_i/E_e correction in §III.C.3. A secondary algebraic slip occurs in Eqs. (20)-(21): Eq. (18) implies k_β/k_α=(m_α/m_β)^{1/ν}, so the denominators should contain (m_α/m_β)^{1/ν}, not (m_α/m_β)^ν. This changes the predicted degrees quantitatively (e.g., m_α=45, m_β=16, ν=2/3 gives k_α≈10.7 rather than 20.4). Together, these errors undermine the paper's assertion that an extended mean-field analysis 'identifies the origin' of the enhanced robustness and re-entrant transition.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a directed bipartite extension of Shimada's evolving open-system model. New species of two types arrive with prescribed input degrees m_alpha and m_beta and random interaction weights; cascading extinctions are resolved before the next introduction. Simulations indicate that when m_alpha = m_beta the divergence-to-finite transition is unchanged from the unstructured case, whereas for asymmetric input degrees the system remains/grows in a much larger parameter region and exhibits a re-entrant transition back to the diverging phase (e.g., along m_beta = 16 with m_alpha near 40-45). The authors propose a mean-field explanation based on a shift of the emergent average degrees (k_alpha, k_beta) away from (m_alpha, m_beta), an empirical extinction law E(k) = c k^{-nu}, and modified cascade sums, and claim that this theory quantitatively reproduces the observed phase boundaries, including the re-entrant transition.","tokens_in":10020,"tokens_out":13861,"duration_ms":138305,"significance":"If correct, the finding that mere bipartite topology plus asymmetric input connectivity can enhance robustness and produce a re-entrant transition is a novel and potentially universal mechanism, relevant to plant-pollinator, microbial, trade, and other bipartite systems. The symmetric-case check against the known unstructured result is a useful internal validation, and the raw simulation evidence appears plausible. However, the theoretical part is not yet at the same standard: the central degree-shift argument rests on an unproven link-balance condition, an exponent error in the closed-form solution, and an empirically fitted extinction law. These issues must be fixed before the claim that the mean-field analysis 'identifies the origin' of the effect can be accepted.","major_comments":[{"comment":"Eq. (18) gives k_beta/k_alpha = (m_alpha/m_beta)^{1/nu}. Therefore k_alpha = (m_alpha+m_beta)/[1+(m_alpha/m_beta)^{1/nu}] and k_beta = (m_alpha+m_beta)/[1+(m_alpha/m_beta)^{-1/nu}], not the expressions with exponent nu in Eqs. (20)-(21). With nu=2/3 and (m_alpha,m_beta)=(45,16), the corrected formula gives k_alpha approximately 10.7 instead of 20.4. This materially changes the theoretical curves in Fig. 7 and the quantitative degree-shift prediction, so the error should be corrected and the figure re-examined.","section":"§III.C.1, Eqs. (18)-(21)"},{"comment":"The condition m_alpha + m_beta = k_alpha + k_beta is asserted without derivation. In a persistent state, each introduction adds m_chi links, while each extinction removes a number of links equal to the degree of the extinct node; the correct balance involves the average degree of extinct species, not the sum k_alpha+k_beta. In the diverging phase there is net link accumulation, so the equality is not generally expected. Since Eqs. (20)-(21) depend on Eq. (19), the predicted degree shift is not established. A proper link-budget derivation is needed.","section":"§III.C.1, Eq. (19)"},{"comment":"The inequality N_E(m_alpha,m_beta) < mu_I/(1-mu_I) may be correct, but the identification mu_I/(1-mu_I) = N_E(mu_m, mu_m) is not: mu_I is the arithmetic mean of I_alpha and I_beta evaluated at the actual m_alpha, m_beta, and this is generally not equal to mu_m E(mu_m)/2 for mu_m = (m_alpha+m_beta)/2. Therefore the claim that an asymmetric bipartite system is more robust than the symmetric system with the same average degree is not proven by this equation, even if it is supported by the simulations.","section":"§III.A, Eq. (14)"},{"comment":"The 'quantitative explanation' is largely an empirical consistency check. The extinction law E(k)=c k^{-nu} is fitted (Eq. 17), the boundary in Fig. 9 uses the observed average degrees k_chi and observed extinction out-degrees, and the separate functions E_i, E_e for introduction and deletion are introduced without specifying how they are obtained. Feeding simulation outputs back into the mean-field formulas can reproduce a phase boundary without independently explaining it. To support the explanatory claim, the authors should either derive E(k) or demonstrate its stability across parameters, and show that the corrected degree-shift formulas together with a parameter-free E predict the phase boundary without using the measured k_chi.","section":"§III.C.2–C.3, Eqs. (17), (29)-(33), Fig. 9"},{"comment":"The phase classification relies on the threshold lim_{t->inf} N_chi(t)/t >= 0.0006 (Fig. 4 caption), and the re-entrant transition is the headline claim. It should be checked that the detected growing/finite strips are robust to reasonable threshold choices (e.g., 0.0001 and 0.001) and to simulation duration. Without such a check, a slowly diverging phase could be misclassified as finite, affecting the reported re-entrant boundary.","section":"§III.B, Figs. 4–6"}],"minor_comments":[{"comment":"The notation here is ambiguous. If only the initial introduction event carries the m_chi/k_chi factor and all later cascade events use the resident degree k_chi, then the factored geometric-series forms are correct; the denominator need not change. The text should state this convention explicitly, because Eqs. (22)-(23) describe the direct impacts of 'an inclusion or extinction', which would suggest a different series.","section":"§III.C.2, Eqs. (24)-(27)"},{"comment":"The abstract states that the diverging phase persists even when both initial degrees individually exceed the critical point, but the illustrative sweep m_beta=16 has m_beta < m_c. Please point to a specific (m_alpha,m_beta) region in Fig. 6 where both m_alpha and m_beta exceed m_c and the system is in the diverging phase.","section":"Abstract and §III.B.2"},{"comment":"The text refers to dash-dotted lines in Fig. 9, but the caption describes solid lines for the improved estimate. Clarify which curves correspond to Eqs. (24)-(27) and which to Eqs. (29)-(32).","section":"Fig. 9 and surrounding text"},{"comment":"The quantities k_check_out_alpha and k_check_out_beta are used in Eq. (33), but §III.C.3 only discusses the extinction bias in group beta. Define both quantities and state how they are measured.","section":"Eq. (33) and Fig. 8"},{"comment":"No code or data availability statement is included; providing this would strengthen reproducibility. Minor wording: 'making it lying entirely outside' in the abstract should be revised.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The empirical phenomenon is interesting and likely publishable after substantial revision. The algebraic exponent error and the unproven link-balance condition are not cosmetic; they undercut the analytic explanation as written. The simulation evidence itself seems worth preserving, and the symmetric-case validation is a strength. I would not reject, but the theory sections need careful reworking before the claims of quantitative explanation can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core simulation result is genuinely new and I believe correct: for asymmetric bipartite topology the phase boundary shifts to higher connectivity, and there is a re-entrant transition back to the diverging phase. The symmetric case matches the known unstructured model, which is a good sanity check. The mechanism the authors propose—degree shift induced by the handshaking constraint and asymmetric extinction rates—is plausible and clearly explained.\n\nWhat the paper does well: the model is simple, the phase diagram is clean, and the authors do not oversell the mean-field theory. The observation that extinctions concentrate among newly introduced low-degree nodes is nice and supported by the degree distributions they show.\n\nThe soft spots are in the theory section, not in the simulation claim. E(k)=c k^{-\\nu} is fitted, not derived; Eq. (19) is assumed without a proper balance argument; and the final match in Fig. 9 feeds observed average degrees and extinction out-degrees back into the formulas, so it is partly a consistency check rather than an independent prediction. The divergence-speed threshold of 0.0006 is arbitrary, and there are no error bars or code/data. These are fixable.\n\nI checked the stress-test note on Eqs. (24)-(27). That specific criticism does not hold up: the series use \\tilde I only for the first event (the new node's m links) and I for subsequent extinctions of residents with average degree k. So the denominator should not be modified. But the secondary point is correct: Eq. (18) gives k_\\beta/k_\\alpha=(m_\\alpha/m_\\beta)^{1/\\nu}, so Eqs. (20)-(21) should have 1/\\nu in the exponent, not \\nu. For m_\\alpha=45, m_\\beta=16, this changes k_\\alpha from about 20 to about 11, so the degree-shift predictions in Fig. 7 are off until corrected.\n\nWho this is for: anyone interested in evolving networks, robustness in bipartite systems (ecology, microbiome, trade). The paper deserves a serious referee. I would send it to peer review with a request to tighten the theory and provide data/error bars, not desk reject.","headline":"Core simulation result is new and credible: asymmetric bipartite topology shifts the phase boundary and produces a re-entrant transition. The mean-field explanation is suggestive but relies on fits and observed simulation outputs, and Eqs. (20)-(21) contain a real exponent error.","tokens_in":10484,"tokens_out":9208,"would_cite":true,"duration_ms":81795,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37N25","05C82","92D40"],"pacs":["89.75.-k","87.23.-n"],"model":"deepseek-v4-flash","headline":"This paper claims that in evolving systems with a bipartite interaction structure, asymmetric initial connectivities between the two partitions markedly enhance robustness, shifting the collapse transition to higher connectivity and even pr","keywords":["bipartite networks","evolving open systems","robustness","phase transition","diverging phase","mean-field theory","degree shift","asymmetric connectivity"],"falsifier":"Measure the extinction probability E(k) directly within the evolving bipartite model at the emergent degrees k_alpha and k_beta predicted by Eqs. (20)-(21). If E(k) deviates measurably from c k^{-nu} with nu about 2/3, or if the link-balance condition does not hold in the growing phase, the degree-shift explanation collapses. Alternatively, run the same model with a different extinction law, such as an exponential in k, and check whether the phase-boundary shift and the re-entrant transition disappear; the claim predicts they would.","tokens_in":9392,"feed_emoji":"🕸️","tokens_out":5471,"duration_ms":48898,"temperature":0.7,"pith_summary":"The paper studies evolving open systems where new species enter and unfit ones go extinct, and asks whether the bipartite two-partition topology itself changes the known transition between an indefinitely growing phase and a bounded finite phase. It establishes that for symmetric initial degrees the transition is unchanged, but for asymmetric initial degrees the system becomes markedly more robust: the phase boundary moves to higher connectivity, and the diverging phase persists even when both initial degrees individually exceed the critical value of the unstructured system. The paper also reports a re-entrant transition, where increasing the asymmetry of one side while the other is fixed first drives the system to a finite phase and then back to the diverging phase. The origin is identified in a degree shift: the bipartite handshaking constraint combined with different extinction rates drives the emergent mean degrees well above the input connectivities, suppressing extinction probabilities while concentrating extinctions among newly introduced low-degree nodes. A sympathetic reader would care because this gives a simple and universal mechanism by which many real two-type systems, such as plant-pollinator or country-product networks, could be more robust than single-community models predict.","feed_headline":"Bipartite asymmetry shields evolving systems from collapse","feed_subtitle":"When the two partitions get unequal link counts, the growing phase persists even past the single-community tipping point.","key_machinery":"The argument is carried by three coupled ingredients: (i) the bipartite handshaking constraint N_alpha k_alpha = N_beta k_beta, which ties population ratios to degree ratios; (ii) the balance condition m_alpha + m_beta = k_alpha + k_beta, which closes the system and yields explicit formulas for the emergent degrees k_alpha and k_beta as functions of the input degrees; and (iii) the approximate extinction-probability law E(k) = c k^{-nu} with c about 1/2 and nu about 2/3, which converts higher degree into lower extinction probability. The corrected mean-field estimate further separates the extinction probability for species introduction from that for species deletion, capturing the concentrat","core_discovery":"The central discovery is that the robustness phase diagram of an evolving bipartite network is controlled not by the input degrees m_alpha and m_beta but by the emergent degrees k_alpha and k_beta, which can deviate strongly from the inputs when the inputs are asymmetric. Under asymmetry, the handshaking relation N_alpha k_alpha = N_beta k_beta, together with the link balance m_alpha + m_beta = k_alpha + k_beta and the empirical extinction law E(k) = c k^{-nu}, implies a self-consistent degree shift, given by Eqs. (20)-(21), in which the side with the smaller input degree gets a disproportionately large emergent degree. This degree elevation suppresses extinction probabilities across the com","pith_inferences":["Because the mechanism rests only on topology-induced degree elevation and extinction concentration, it should also apply when the two partitions have different intrinsic extinction functions, not just different insertion degrees: a testable extension the paper does not pursue.","In empirical bipartite communities with asymmetric interaction counts, one could look for signatures of the re-entrant transition: communities with strongly asymmetric interaction rates should show enhanced persistence compared with symmetric ones at the same mean interaction count.","The paper treats the extinction law as input; measuring E(k) directly inside the evolving bipartite system, rather than transferring it from the unstructured model, would either confirm or require modification of the predicted degree shift.","The re-entrant transition suggests that in engineered two-sided platforms, increasing activity asymmetry can act as a robustness lever, though the paper makes no engineering claim."],"forward_implications":["For any bipartite evolving system with asymmetric insertion rates, the collapse threshold is higher than for the corresponding unstructured system with the same average input degree; the two thresholds coincide only in the symmetric limit.","The diverging growing phase can persist even when both insertion degrees lie above the unstructured critical point, so robustness is not bounded by the single-community transition.","There exists a regime where increasing the input connectivity of one side drives the system from growing to finite and back to growing: a re-entrant transition absent in unstructured systems.","The mechanism is universal in that it depends only on the handshaking constraint, the balance of link gain and loss, and a monotonically decreasing extinction probability, not on the details of fitness or rewiring.","The extended mean-field framework with degree shift and separate introduction and deletion extinction probabilities can quantitatively predict the phase boundary, making it a usable tool for other bipartite assembly models."],"fun_headline_variants":["Asymmetric bipartite links shield evolving systems from collapse","Unequal connectivity delays phase transition in evolving networks","Bipartite asymmetry keeps systems growing past critical threshold","Why lopsided bipartite networks are more robust","Asymmetry fortifies evolving bipartite systems"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the extinction probability of a node obeys E(k) = c k^{-nu} with c about 1/2 and nu about 2/3, fitted on the same class of simulations and then assumed unchanged in the evolving bipartite system, and that the link-balance condition m_alpha + m_beta = k_alpha + k_beta closes the mean-field equations; if either fails, the predicted degree elevation and the robustness enhancement do not follow from the argument, even if the raw simulation phase d","fun_headline_variants_meta":{"raw":{"variants":["Asymmetric bipartite links shield evolving systems from collapse","Unequal connectivity delays phase transition in evolving networks","Bipartite asymmetry keeps systems growing past critical threshold","Why lopsided bipartite networks are more robust","Asymmetry fortifies evolving bipartite systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1287,"prompt_tokens":761,"completion_tokens":526,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":505,"completion_tokens_details":{"reasoning_tokens":452}},"tokens_in":505,"tokens_out":526,"duration_ms":6170,"temperature":1.0,"reasoning_tokens":452,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T08:52:15.381537+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the extinction probability E(k) directly within the evolving bipartite model at the emergent degrees k_alpha and k_beta predicted by Eqs. (20)-(21). If E(k) deviates measurably from c k^{-nu} with nu about 2/3, or if the link-balance condition does not hold in the growing phase, the degree-shift explanation collapses. Alternatively, run the same model with a different extinction law, such as an exponential in k, and check whether the phase-boundary shift and the re-entrant transition disappear; the claim predicts they would.","supporting_citations":[],"review_version":1}