Degree-Phase-Separation Monitor
Implementation & benchmark of arXiv:2607.28275 — Synchronization, Kinematic Waves and Spike-Phase-Separation in Feedback Ising Neural Networks on Heterogeneous Graphs
Source paper: Synchronization, Kinematic Waves and Spike-Phase-Separation in Feedback Ising Neural Networks on Heterogeneous Graphs arXiv:2607.28275 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Use the degree-resolved phase-separation mechanism as a diagnostic and regularizer for graph and recurrent networks. Penalize unintended divergence between peripheral-node and hub representations, or deliberately preserve bounded divergence when heterogeneous specialization is useful.
Formulas
Mathematical statement
Let \(m_k\) be the mean hidden activation of nodes with degree \(k\), let \(p_k\) be the degree distribution, and define edge-weighted activity \(u=\sum_kp_k k m_k/\langle k\rangle\). The degree-resolved mean-field dynamics are \(\dot m_k=-m_k+\tanh[\beta(Jku+h)]\). A symmetric fixed point can lose stability to a nonuniform degree mode, producing a pitchfork-like phase-separated state. Define the order parameter \(S=\sum_kp_k(m_k-m)^2\), where \(m=\sum_kp_km_k\). If \(F(\mathbf m)\) denotes the vector field, the separation transition is detected when the largest eigenvalue \(\lambda_{\mathrm{sep}}\) of its Jacobian restricted to zero-mean degree perturbations crosses zero. Penalizing \(S\) suppresses this instability.
Implementation notes
Add a degree-conditioned monitor to every message-passing or graph-recurrent layer. Partition nodes into logarithmic degree bins \(B_b\), compute the mean hidden vector \(m_b=|B_b|^{-1}\sum_{i\in B_b}x_i\), and use bin proportions \(p_b=|B_b|/N\). Compute \(S=\sum_bp_b\|m_b-m\|^2\), with \(m=\sum_bp_bm_b\), and add \(\lambda S\) to the task loss when synchronized representations are desired. Do not confuse this with ordinary global activation regularization: the penalty specifically detects degree-conditioned separation. Estimate the restricted instability by applying Jacobian-vector products to the recurrent layer and using 5–10 power iterations while projecting each vector to zero weighted mean across degree bins. This gives an empirical estimate of \(\lambda_{\mathrm{sep}}\) without constructing the full Jacobian. The degree moments and order parameter are computed exactly from each graph; the bifurcation eigenvalue is estimated numerically. First experiment: use a 3-layer GCN or GatedGCN for node classification on synthetic graphs with equal mean degree and increasing degree heterogeneity. Compare no penalty, global activation regularization, and the proposed degree-variance penalty. Sweep recurrent gain. The prediction is a sharp rise in \(S\) when \(\lambda_{\mathrm{sep}}\) crosses zero; the proposed penalty should shift or remove this transition while preserving global activation scale.
Verification
Stage 1 · Toy mechanism gate: Failed ✗
Stage 2 · Mechanism transferred to benchmark: Not tested
Stage 2 · Practical benchmark result: Not run
Methodology: Toy-system gate first; the benchmark stage runs only after a pass. How verification works
Stage 1 — Mechanism check agent confidence 9/10
Built the degree-resolved mean-field monitor, nonlinear simulator, Jacobian diagnostic, and regularizer comparison. The predicted collective thresholds g_c=<k>/<k^2> matched the nonlinear separation onset within the gain-grid resolution, and the finite-difference Jacobian check had error below 1e-11. However, the proposed zero-mean restricted eigenvalue remained approximately -1 rather than crossing zero, so the literal separation-eigenvalue claim was not validated; the experiment supports the phase-separation mechanism but exposes a mismatch between that diagnostic and this toy dynamics.
- Agent confidence
- 9/10
- Baseline
- At k=(1,5), gain=0.55: S=0.06985, mean activation=0.72145.
- Idea
- Degree penalty reg=0.2: S=0.03860, mean=0.70498; reg=0.5: S=0.02009, mean=0.69012. Global penalty reg=0.2 gave S=0.04773 but reduced mean to 0.44139.
- Thresholds
- For heterogeneity d=0,1,2,4, predicted gains were 0.3333, 0.3000, 0.2308, 0.1200; sampled nonlinear onsets were none, 0.32, 0.24, 0.12 respectively. At each predicted threshold, the largest full Jacobian eigenvalue was within 1.4e-16 of zero.
Limitations: Only a two-degree mean-field toy system was tested; no GCN, recurrent neural network, graph node-classification task, training/generalization evaluation, logarithmic degree bins, GPU experiment, or Jacobian-vector power iteration was benchmarked. The restricted-eigenvalue implementation and the paper's exact projection geometry may require refinement for a non-rank-one graph/recurrent layer.
How to run: python3 degree_phase_monitor.py && /home/maxwelhelp/main/bin/python3 verify_extra.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- degree_phase_monitor.py 3.8 KB View Raw
- report.md 1.9 KB View
- results.json 35.7 KB View Raw
- verify_extra.py 0.8 KB View Raw