{"id":"a6d64005-4fad-4027-a585-9109e52388cf","arxiv_id":"2607.28275","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"On heterogeneous graphs, feedback Ising networks oscillate above a Hopf line fixed by α=⟨k²⟩/⟨k⟩, and strong heterogeneity can stabilize a phase-separated fixed point that destroys the limit cycle.","lead":"Heterogeneous network structure controls when feedback-driven Ising neural networks start oscillating, via a closed-form Hopf threshold set by the degree moment ratio. The same heterogeneity creates degree-ordered activation waves and a low-temperature phase-separated state that can kill the synchronized cycle.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The reader correctly isolates the annealed HMF closure (Eqs. 5–10, App. A) as the weakest assumption and still assigns ACCEPT with medium correctness risk for quenched robustness and biology—precisely the right calibration. Re-reading the linear stability, discriminant simplification to (c+α)²−8c, PS third-order expansion, and bimodal Jacobian reduction turned up no algebraic or logical break in the load-bearing claim. MC phase extraction (App. B) is heuristic but directionally supportive; the paper itself notes dark modes and cavity/DMFT extensions. No verdict shift is warranted.","tokens_in":17839,"tokens_out":526,"duration_ms":52936,"concrete_test":"Integrate the deterministic HMF ODEs for the symmetric bimodal graph across the line φ(β)=β_Hopf(c) (Sec. II.C.6) and, at the same (β,c,k1,k2), run quenched MC; if a stable finite-amplitude LC coexists with the stable PS fixed points on the ODE side, or if quenched MC retains oscillations where HMF predicts PS capture, the abstract’s “dynamically destroys the synchronized limit cycle” claim needs narrowing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Within the paper’s stated Curie–Weiss/HMF scope, the central claims hold. The Hopf threshold (Eq. 19) follows cleanly from the (m,u,h) Jacobian and Routh–Hurwitz; the m–u decoupling, kinematic-wave phase ordering, and PS pitchfork at β_c2=⟨k⟩/Var(k) are internally consistent; the bimodal reduction φ=β sech²(β Δk u_s/2) giving exact PS stabilization when φ=β_Hopf(c) is a genuine closed-form result. Quenched MC (ER/SF/log-normal/bimodal) already tracks the annealed Hopf boundary and shows PS after quench, so the annealed rank-one closure—while formally the weakest structural premise and correctly flagged by the reader—is not, on present evidence, breaking the headline claims. Residual gaps (supercriticality Lyapunov coefficient not recomputed; global bifurcation that eliminates the LC when PS stabilizes not fully charted; dark-mode splitting on quenched graphs) are real but secondary and largely acknowledged as future work.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript studies a kinetic Ising model with homeostatic linear feedback from global magnetization to the external field on random graphs with heterogeneous degree distributions. Within a Curie–Weiss heterogeneous mean-field closure, the authors linearize the (m,u,h) dynamics about the paramagnetic fixed point, apply Routh–Hurwitz, and obtain a closed-form Andronov–Hopf threshold β_c(c) controlled by the moment ratio α=⟨k²⟩/⟨k⟩ (Eq. 19), with limits recovering the known ferromagnetic HMF critical point. Degree heterogeneity decouples the per-neuron rate m from the per-synapse rate u, producing (i) kinematic waves of degree-ordered activation, quantified via cross-correlations C_k(τ), and (ii) a low-temperature phase-separated fixed point (m=0, u≠0) born in a supercritical pitchfork at β_c2=⟨k⟩/Var(k). For bimodal graphs they reduce PS stability exactly to a rescaled paramagnetic condition φ=β sech²(·)=β_Hopf(c). Analytical phase boundaries are compared to Monte Carlo on ER, scale-free, log-normal, and bimodal graphs.","tokens_in":18166,"tokens_out":1407,"duration_ms":30401,"significance":"If the results hold, the paper supplies a rare analytically tractable out-of-equilibrium phase diagram for feedback-driven Ising dynamics on heterogeneous networks, with an explicit topology dependence through α and two phenomena (kinematic waves, spike-phase separation) absent on regular graphs. The closed-form Hopf line, the bimodal reduction φ→β_Hopf, and the van Kampen derivation of the HMF equations are concrete technical contributions that can be reused and falsified. Monte Carlo agreement on several degree ensembles strengthens the claim within the annealed scope. The work cleanly extends prior homogeneous feedback-Ising results and connects network heterogeneity to oscillatory onset and limit-cycle destruction in a way that is useful for statistical mechanics of neural and complex systems.","major_comments":[{"comment":"Abstract and §II.C.5–6 claim that when the phase-separated fixed point stabilizes it “dynamically destroys the synchronized limit cycle.” Stabilization of the PS point is controlled analytically (especially the exact bimodal reduction φ(β_stab)=β_Hopf(c)), but the global mechanism that eliminates the LC (e.g., collision, basin capture, or fold of cycles) is not charted. Fig. 4 shows a quench into PS, which is consistent with capture but does not establish how the LC ceases to exist as a stable object. A minimal addition—continuation of the LC in β for a bimodal or high-variance case, or a clear statement that destruction is inferred from MC basins rather than proven—would make the central claim load-bearing and precise.","section":"§II.C, Abstract"},{"comment":"The Hopf bifurcation is repeatedly called supercritical (Abstract, Introduction, §II.A), yet the first Lyapunov coefficient is not recomputed for the heterogeneous (m,u,h) Jacobian. Supercriticality is inherited from the homogeneous/lattice feedback-Ising literature. Because the phase diagram and the interpretation of a stable macroscopic limit cycle rest on this, either a short Lyapunov-coefficient calculation (or normal-form reduction) for Eqs. (11)–(13), or an explicit caveat that supercriticality is assumed by continuity with the homogeneous case and supported by MC, should be added.","section":"§II.A.2, Eqs. (11)–(19)"}],"minor_comments":[{"comment":"Appendix B phase-detection thresholds (peak position >0.05, p(0)<1) are ad hoc. State that they were fixed from visual inspection and, if possible, show that the extracted β_c(c) is stable under modest threshold changes, or mark the MC phase boundary as semi-quantitative.","section":"Appendix B"},{"comment":"Fig. 2 insets confirm oscillatory vs fixed-point regimes but do not overlay the analytical β_c(c) on a quantitative MC scan (e.g., oscillation amplitude or spectral peak vs β). A single panel with measured onset vs Eq. (19) would strengthen the validation claim.","section":"Figure 2"},{"comment":"Notation: α is introduced as ⟨k²⟩/⟨k⟩ in the text after Eq. (15), while z is used for ⟨k⟩; keep this consistent in the ER specialization (Eq. 22) and in the PS section where z reappears.","section":"§II.A.2–3"},{"comment":"In §II.B the kinematic-wave interpretation is clear, but the claim that the structure is “encoded in the Hopf eigenvector through q_u/q_h=β_H/(α+iω_0)” appears only in the Conclusion. Either move a one-line derivation into §II.B or soften the Conclusion wording.","section":"§II.B, §III"},{"comment":"Typos/style: “F eedback”, “W eiss”, “RESUL TS” (space in headings); “withe noise” in Appendix A; “as done done” in Appendix B; arXiv stamp date “30 Jul 2026” looks like a placeholder.","section":"Throughout / Appendix A–B"},{"comment":"Briefly clarify the feedback update convention Δt=1/N (after each spin flip) versus the continuous ˙h=−c m used in HMF, so that MC and theory timescales are unambiguously matched.","section":"§II, Eq. (4) and Eqs. (9)–(10)"}],"recommendation":"minor_revision","confidential_remarks":"The central derivations are sound within the stated HMF scope; the annealed rank-one closure is the main structural limitation but is acknowledged and already stress-tested by quenched MC. I would not block on quenched dark-mode splitting. Fit to cond-mat.stat-mech / interdisciplinary statistical physics is good. Prior self-citations supply the homogeneous mechanism but the heterogeneous Jacobian, α-dependence, and PS analysis are new. Minor revision is appropriate; I would accept after the two major points are clarified without requiring a full global-bifurcation atlas."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real payload is the closed-form Hopf threshold β_c(c) that depends on the network only through α=⟨k²⟩/⟨k⟩, plus the clean demonstration that m and u decouple and produce two effects you do not get on regular graphs: degree-ordered kinematic waves and a pitchfork phase-separated state (m=0, u≠0) that can kill the limit cycle when heterogeneity is large.\n\nThey do the linearization and Routh–Hurwitz properly. The discriminant collapses nicely, the c→0 and c→∞ limits recover the classical HMF ferromagnet, and the bimodal reduction (β → φ = β sech²) that pins PS stability exactly to the paramagnetic Hopf line is a genuine closed-form result. Monte Carlo on ER, scale-free, log-normal and bimodal graphs tracks the annealed boundary and shows the quench into PS; that is enough support for the claims inside the stated approximation. Prior self-citations supply the homogeneous feedback setup; the heterogeneous Jacobian, α-dependence, C_k(τ) waves and PS bifurcation are derived here.\n\nSoft spots are real but secondary and mostly flagged. Everything sits on the annealed Curie–Weiss closure (rank-one interaction, dark modes stuck at −1). Quenched graphs will split those modes and may reshape the wave structure; they say so. They do not recompute the Lyapunov coefficient for supercriticality, the global bifurcation that eliminates the LC when PS stabilizes is not fully charted, and the phase-detection KDE thresholds are a bit ad hoc. No code ships. None of that breaks the headline analytics.\n\nThis is for people who already care about mean-field spin models, network Ising, or nonequilibrium neural criticality. If that is your lane, read it; the formulas are usable. It deserves a serious referee. I would engage.","headline":"Clean closed-form Hopf line on heterogeneous graphs plus two real new phenomena (kinematic waves, m=0/u≠0 phase separation); annealed HMF is the limit, not a break.","tokens_in":18811,"tokens_out":490,"would_cite":true,"duration_ms":8205,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.70.Ln","64.60.aq","87.19.lj","05.45.-a"],"model":"grok-4.5","headline":"Degree heterogeneity sets when feedback Ising networks start oscillating, and extreme heterogeneity can kill the synchronized cycle by locking the system into a phase-separated state.","keywords":["feedback Ising model","heterogeneous mean-field","Andronov-Hopf bifurcation","kinematic waves","phase separation","degree heterogeneity","neural synchronization","random graphs"],"falsifier":"On a large quenched scale-free or bimodal graph, measure whether the analytic Hopf boundary still matches the onset of global oscillations and whether the phase-separated fixed point really captures and extinguishes the limit cycle once degree variance exceeds the predicted threshold.","tokens_in":18700,"feed_emoji":"🌊","tokens_out":1024,"duration_ms":21702,"temperature":0.7,"pith_summary":"This paper asks how uneven connectivity changes collective oscillations in a simple neural-network model: binary spins on a random graph whose excitability is homeostatically pushed against the population firing rate. Using a heterogeneous mean-field theory checked against Monte Carlo runs, the authors show that the silent state loses stability through a Hopf bifurcation whose threshold is fixed by the single number α = ⟨k²⟩/⟨k⟩. Because heterogeneity splits the ordinary magnetization m from a synapse-weighted magnetization u, two effects appear that never occur on regular lattices: degree-ordered kinematic waves that sweep from leaves to hubs, and a low-temperature phase-separated fixed point with m = 0 but u ≠ 0. For sufficiently broad degree distributions that fixed point stabilizes and swallows the limit cycle. The result supplies closed-form phase boundaries that let network topology, not just feedback gain, be read as a control knob for macroscopic neural rhythms.","feed_headline":"Heterogeneity can kill neural sync once it births a phase-separated state","feed_subtitle":"A single degree-moment ratio sets the Hopf line and decides when oscillations die on random graphs.","key_machinery":"The Curie–Weiss heterogeneous mean-field reduction that tracks two distinct order parameters—firing rate per neuron m and firing rate per synapse u—together with the closed-form Hopf line β_c(c) = (3α − c − √((c + α)² − 8c))/(2(α² − c(α − 1))) obtained from the Routh–Hurwitz conditions on their joint Jacobian.","core_discovery":"On heterogeneous random graphs the paramagnetic fixed point of the feedback Ising model loses stability via a supercritical Andronov–Hopf bifurcation at the explicit threshold β_c(c) determined solely by the moment ratio α = ⟨k²⟩/⟨k⟩ and the feedback strength c; when degree variance is large enough, a pitchfork-born phase-separated state (m = 0, u ≠ 0) later stabilizes and dynamically destroys the synchronized limit cycle.","pith_inferences":["If cortical degree distributions broaden or narrow with aging or disease, the same homeostatic feedback could flip a circuit from rhythmic to phase-separated silence without altering cellular excitability parameters.","The rank-one dark modes that sit at eigenvalue −1 under the annealed approximation should split on quenched graphs, potentially producing additional slow kinematic modes visible in degree-resolved spectra.","Extending the same m–u decoupling to non-reciprocal or multi-state spins would likely generate coexisting avalanche and oscillatory regimes whose boundaries still collapse onto α."],"forward_implications":["Scale-free networks with 2 < γ ≤ 3 sit permanently in the oscillatory phase at any nonzero temperature in the thermodynamic limit, independent of feedback strength.","Kinematic waves—low-degree nodes leading, hubs lagging—appear both on the limit cycle and as measurable peak delays in subcritical cross-correlations ⟨m_k(t+τ)h(t)⟩.","For bimodal graphs the phase-separated point stabilizes exactly when an effective inverse temperature ϕ = β sech²(β Δk u_s/2) drops below the paramagnetic Hopf line, giving a sharp topology-dependent death of synchrony.","Slow drift of α (for example through gradual loss of hubs) can push a network across the phase boundary without any change in the feedback gain c."],"fun_headline_variants":["Degree heterogeneity kills sync via phase-separated state","Moment ratio sets Hopf line then ends oscillations on random graphs","Phase separation stabilizes and destroys the neural limit cycle","Kinematic waves run periphery-to-hub before sync dies","High degree variance quenches synchronized oscillations"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The annealed mean-field closure that every node of the same degree shares one magnetization and feels only a single average neighbor field, turning the interaction into a rank-one object.","fun_headline_variants_meta":{"raw":{"variants":["Degree heterogeneity kills sync via phase-separated state","Moment ratio sets Hopf line then ends oscillations on random graphs","Phase separation stabilizes and destroys the neural limit cycle","Kinematic waves run periphery-to-hub before sync dies","High degree variance quenches synchronized oscillations"]},"model":"grok-4.5","effort":"low","cost_usd":0.004784,"raw_usage":{"total_tokens":1365,"prompt_tokens":797,"num_sources_used":0,"completion_tokens":59,"cost_in_usd_ticks":47844000,"prompt_tokens_details":{"text_tokens":797,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":509,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":797,"tokens_out":59,"duration_ms":12274,"temperature":1.0,"reasoning_tokens":509,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T12:30:27.934981+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"On a large quenched scale-free or bimodal graph, measure whether the analytic Hopf boundary still matches the onset of global oscillations and whether the phase-separated fixed point really captures and extinguishes the limit cycle once degree variance exceeds the predicted threshold.","supporting_citations":[],"review_version":1}