{"id":"5abf5f84-3fbb-4da5-a882-9b52b2780f08","arxiv_id":"2607.26213","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The conditional energy landscape of phase-based Ising machines has an astroid spinodal boundary separating monostable and bistable regimes, with barrier scaling ΔE ∝ μ^{3/2} at generic points and ΔE ∝ μ^2 at the longitudinal cusp.","lead":"A star-shaped boundary (an astroid) determines when an oscillator in a phase-based Ising machine can get trapped in a wrong spin state. The paper maps this boundary and the energy barrier growth, connecting it to classic Stoner–Wohlfarth switching.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The astroid boundary is derived for a frozen local field; it is not shown to be a bifurcation set of the full OIM/DIM dynamics, so the barrier scaling may be an artifact of the static conditional landscape.","rationale":"The paper's derivation is internally sound: the parametric solution (19), the astroid elimination (20), and the fold/cusp expansions (32)-(41) check out, and the Stoner-Wohlfarth mapping is an exact identity. The only serious weakness is the transfer from the conditional landscape to statements about real OIM/DIM behavior. The reader's verdict flags the frozen-field assumption; my concern sharpens it: even for a fixed field configuration, the scalar condition used to define the spinodal does not generally coincide with a saddle-node of the coupled phase dynamics. The paper does not attempt to prove this correspondence, and Fig. 3 is not a validation of it. This does not invalidate the mathematical claims about E_i, which are carefully qualified as conditional in the text, so I leave the verdict unchanged; however, any application of the barrier scaling to the dynamics of an actual Ising machine should be treated as unverified until the proposed test is run.","tokens_in":8983,"tokens_out":16517,"duration_ms":155131,"concrete_test":"Integrate the full OIM/DIM phase equations for the same N=10, p=0.5, K=1, K_s=1 network used in Fig. 3. At each time step compute (i) the conditional margin mu_i(t) = 1 - (|h_x|/2K_s)^(2/3) - (|h_y|/2K_s)^(2/3) and (ii) the eigenvalue of the full Jacobian with smallest real part. Test whether mu_i = 0 crossings coincide with zero crossings of that eigenvalue. Separately, initialize one oscillator in the Ising-disfavored binary state, run many realizations, and compare log-mean escape time against mu_i; predicted slopes are 3/2 at smooth spinodal points and 2 at the longitudinal cusp. If either test fails, the conditional astroid is not the dynamical spinodal of the machine.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II defines E_i(theta_i) with h_x,i, h_y,i as instantaneous local fields and derives the astroid from dE_i/dtheta_i = 0 and d^2E_i/dtheta_i^2 = 0. The algebraic steps are correct, and the Stoner-Wohlfarth identity is exact. The load-bearing gap is the step from this frozen-field landscape to the claimed monostable/bistable behavior of the machine. In a coupled network the local field is a function of the other phases; the full flow is not in general a gradient flow on E_i. A zero of the scalar curvature d^2E_i/dtheta_i^2 at a conditional stationary point is not a zero eigenvalue of the full N-oscillator Jacobian, because off-diagonal couplings shift stability boundaries. Thus a node can lie inside the astroid while the corresponding binary state is not a stable fixed point of the coupled system, or outside the astroid while the full system still has a stable state there. Fig. 3 only plots field trajectories against the astroid; it does not verify that astroid crossings coincide with saddle-node events or that escape rates follow the predicted Delta E proportional to mu_i^{3/2} (or mu_i^2 at the cusp). The central quantitative claim therefore rests on an unvalidated quasi-static assumption.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the instantaneous conditional energy E_i(θ_i) = −h_x,i cosθ_i − h_y,i sinθ_i − (K_s/2) cos(2θ_i) of a single oscillator in OIM/DIM networks. It derives the astroid spinodal boundary (Eq. 20), the SHI threshold (Eq. 22), and the barrier scalings ΔE ∝ μ^{3/2} at generic smooth points and ΔE ∝ μ^2 at the longitudinal cusp (Eqs. 35, 41). It also proves an exact correspondence with the Stoner–Wohlfarth energy and shows that OIM and DIM conditional landscapes differ only by a reflection of the transverse field. The derivation is self-contained and algebraically correct within the stated frozen-field framework.","tokens_in":9283,"tokens_out":15319,"duration_ms":136421,"significance":"The result is a clean, parameter-free characterization of local bistability in the conditional landscape of phase-based Ising machines. The Stoner–Wohlfarth correspondence is exact and provides a useful unifying perspective. The barrier scalings are derived from a normal-form expansion rather than fitted, and the astroid boundary is an explicit algebraic result. The main significance is as a local, quasi-static analysis; its direct applicability to the full coupled network dynamics is not established and requires scoping or numerical support.","major_comments":[{"comment":"The astroid is derived from the 1D conditional energy with h_x,i, h_y,i frozen. This is a spinodal of the conditional landscape, but the paper's language ('locally stable', 'trapped', 'bistability in phase-based Ising machines') implies statements about the coupled N-oscillator system. In the full dynamics the local field co-evolves; a zero of ∂²E_i/∂θ_i² at a conditional stationary point is not a zero eigenvalue of the full Jacobian because off-diagonal couplings shift stability. Fig. 3 only plots field trajectories against the astroid; it does not show that crossings coincide with saddle-node events of the full system or that escape rates follow Eqs. (35)/(41). Please either justify a timescale separation or explicitly scope all claims to the frozen-field conditional landscape.","section":"Sec. II, Eq. (2), Eq. (20), Fig. 3"},{"comment":"The paper states that inside the astroid the landscape is bistable and outside monostable. This is plausible—extrema can only appear/disappear at degenerate critical points on the astroid, and h=0 has two minima—but no proof is given. A short argument, or a reference to the Stoner–Wohlfarth analysis, is needed to establish that the number of minima is two for all ℓ<1 and one for ℓ>1. This is load-bearing because the central claim of the paper is that Eq. (20) is the boundary separating monostable from bistable regimes.","section":"Sec. II.B, Eq. (21)"}],"minor_comments":[{"comment":"The notation E_i(q)−E_i(0) with q=θ_i−π conflicts with Eq. (5), where E_i(0) means E_i(θ_i=0). Here E_i(0) is the value at q=0 (i.e., θ_i=π). Please relabel, e.g., E_i(π+q)−E_i(π).","section":"Eq. (39)"},{"comment":"The statement that increasing |h_x,i| 'eventually restores a single-well landscape' is true, but the threshold |h_x,i|=2K_s is not stated. Adding it would make the on-axis discussion more precise.","section":"Sec. II.A, after Eq. (17)"},{"comment":"The caption says 'randomly generated graph with N nodes and edge density p' while the text gives N=10, p=0.5. Please state the graph construction in the caption and describe how the crossing points ('diamonds') are detected.","section":"Fig. 3"},{"comment":"The term 'first-harmonic conditional landscape' in the abstract is not defined. Clarify that it refers to the energy containing cosθ and sinθ terms plus the SHI cos2θ term.","section":"Introduction / abstract"},{"comment":"Minor typographical issue: 'NSF grant #No. 2328961' should read 'NSF grant No. 2328961'.","section":"Acknowledgments"}],"recommendation":"major_revision","confidential_remarks":"The paper is mathematically sound as an analysis of the frozen-field conditional landscape. My main concern is the gap between the conditional-landscape spinodal and the claimed behavior of the full coupled machine. This is fixable either by adding an explicit scope limitation or by a numerical check of the quasi-static assumption. I see no concerns about originality or citation practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the astroid boundary for the conditional one-oscillator energy is derived correctly, the SHI threshold formula is new and useful, and the Stoner–Wohlfarth mapping is a genuinely nice observation. The paper is careful to say all results are for a fixed snapshot of neighboring phases, so the main soft spot is the unvalidated leap from that frozen-field picture to what a coupled network actually does.\n\nThe new material: parametric spinodal (Eq. 19), the astroid boundary (Eq. 20), the generalized SHI threshold (Eq. 22), and the barrier scaling exponents 3/2 and 2 for generic points and the longitudinal cusp. I verified the key algebra myself; it checks out. The paper also cleanly separates longitudinal and transverse fields, shows OIMs and DIMs differ only by transverse-field reversal, and identifies the exact equivalence with the Stoner–Wohlfarth energy. That equivalence is not just a metaphor—the normalized energies are identical up to an additive constant. Good, solid work.\n\nThe soft spots are proportionate. The frozen-field assumption is the big one. The stress-test objection about the full N-oscillator Jacobian is mathematically correct, but it doesn't fully land because the paper never claims the astroid is a bifurcation set of the coupled system. The claims are explicitly about the conditional landscape. The legitimate issue is that the practical significance for choosing SHI strength depends on whether the quasi-static picture is a good approximation in actual dynamics, and the paper doesn't test that. Fig. 3 only plots field trajectories against the astroid; it doesn't show that a crossing coincides with a saddle-node event or that the predicted barrier scaling is observed. The graph-structure discussion is heuristic. There's also no code or data beyond the illustrative figure—\"available upon reasonable request\" is the standard placeholder.\n\nWho this is for: people working on oscillator Ising machines or dynamical Ising machines who want design intuition about SHI strength, and anyone interested in exact analogies between phase-based spin systems and micromagnetic energy landscapes. It's not a breakthrough, but it's a well-made, self-contained contribution.\n\nVerdict: yes, this deserves a serious referee. The math is sound and the SW hook is valuable. A revision should either add a small dynamical validation (e.g., a few-node simulation showing that crossing the astroid correlates with the appearance or disappearance of a metastable state in the full flow) or explicitly temper the framing from \"drives a transition in the machine\" to \"describes the conditional landscape.\" I'd engage with it.","headline":"A clean, honestly framed derivation of the astroid spinodal for the conditional single-oscillator landscape; the gap to full network dynamics is real but not fatal because the paper never claims more than it proves.","tokens_in":9776,"tokens_out":3030,"would_cite":true,"duration_ms":32034,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The conditional energy landscape of a single oscillator in phase-based Ising machines is bistable inside an astroid-shaped boundary, not merely at a critical SHI strength.","keywords":["Oscillator Ising machines","Dynamical Ising machines","spinodal boundary","astroid","bistability","second-harmonic injection","Stoner–Wohlfarth","energy barrier"],"falsifier":"The decisive check is a numerical integration of a small coupled OIM or DIM network where the local field is allowed to evolve dynamically, compared against the frozen-field astroid prediction: if nodes spend time inside the frozen astroid but do not exhibit two stable wells, or vice versa, the quasi-static boundary is violated. Alternatively, escape-time measurements near a smooth crossing versus near the longitudinal cusp could test the predicted 3/2 versus 2 barrier scaling.","tokens_in":8864,"feed_emoji":"🌀","tokens_out":3486,"duration_ms":33253,"temperature":0.7,"pith_summary":"The paper claims that in oscillator Ising machines (OIMs) and dynamical Ising machines (DIMs), the competition between second-harmonic injection (SHI) and the instantaneous local network field reshapes each oscillator's conditional energy landscape, and that this landscape undergoes a monostable-to-bistable transition when the normalized local field crosses an astroid spinodal boundary. The paper derives the astroid boundary from a fold condition and shows the barrier vanishes with exponent 3/2 at generic smooth points and exponent 2 at the longitudinal cusp. Because OIMs and DIMs share the same geometry, and the first-harmonic landscape is mathematically the Stoner–Wohlfarth energy, the result connects phase-based Ising machine design to a classical magnet physics framework. The matter-of-course consequence is that whether a node can get trapped in a wrong phase state is not set by coupling strength alone but by the full orientation of the local field relative to the astroid.","feed_headline":"Astroid spinodal boundary governs bistability in phase-based Ising machines","feed_subtitle":"A single oscillator's landscape turns bistable inside an astroid; barrier grows as the 3/2 power of margin.","key_machinery":"The load-bearing object is the astroid spinodal boundary, obtained by simultaneously imposing stationarity and vanishing curvature on the one-dimensional conditional landscape. Parametrically, the boundary is h_x,i = -2K_s cos^3 θ_sp and h_y,i = 2K_s sin^3 θ_sp, which eliminate to the astroid equation. This object organizes the entire argument: it fixes the SHI threshold K_s,tr,i = (1/2)(|h_x|^{2/3}+|h_y|^{2/3})^{3/2}, determines barrier scaling through a fold expansion, and connects the oscillator problem to the Stoner–Wohlfarth switching astroid.","core_discovery":"The central discovery is that local bistability in phase-based Ising machines is governed by an astroid spinodal boundary in the normalized local-field plane. For a fixed snapshot of neighboring phases, the conditional energy E_i(θ_i) = -h_x,i cos θ_i - h_y,i sin θ_i - (K_s/2) cos 2θ_i has two stable wells exactly when |h_x,i/(2K_s)|^{2/3} + |h_y,i/(2K_s)|^{2/3} < 1. Crossing this boundary creates or destroys a metastable well, with barrier height scaling as K_s μ_i^{3/2} at smooth points and K_s μ_i^2 at the longitudinal cusp, where μ_i is the spinodal margin. The paper also shows OIMs and DIMs differ only by reflection of the transverse field and that the same energy is equivalent (up to a","pith_inferences":["Because the local field is treated as frozen only for the spinodal analysis, the astroid is a quasi-static boundary; in fast networks the effective boundary may be smeared or shifted by field fluctuations, a testable extension being comparison with simulations that include finite relaxation times.","The phasor-sum expression for the local field suggests that graph structure enters through interference, not just degree; wiring loops or motifs that promote partial phase cancellation could systematically delay nodes from entering the bistable region, an effect the paper leaves implicit.","The same first-harmonic landscape appears in other phase-encoded computing schemes, so the astroid criterion could be ported to any system described by a cosine double-well plus a first-harmonic tilt, including some optical or mechanical Ising machines.","The cusp scaling ΔE ∝ μ^2 is slower than the generic 3/2 scaling, implying that longitudinal crossing events create barriers that grow quadratically; one could test this by intentionally biasing nodes along the axis and measuring escape times."],"forward_implications":["The SHI threshold for local bistability is not a single coupling value but a field-orientation-dependent threshold; equal-magnitude fields can require SHI strengths differing by a factor of two.","Since OIMs and DIMs share the astroid geometry, any design rule for avoiding metastable traps applies to both architectures with the transverse field sign reversed.","The Stoner–Wohlfarth equivalence lets phase-based Ising machine bistability be analyzed with known results from single-domain magnet switching, including barrier formulas and switching probabilities.","The barrier scaling (3/2 at smooth points, 2 at cusp) determines how quickly metastable wells become effective traps as the SHI strength is increased beyond threshold.","In a network, nodes repeatedly cross the astroid during dynamics, so local monostable/bistable transitions are a normal part of computation, not just an endpoint effect."],"fun_headline_variants":["Astroid spinodal: bistability frontier in Ising machines","Astroid curve dictates bistability in phase-based Ising machines","Bistability in Ising machines: an astroid spinodal boundary","Astroid spinodal sets bistable regime in Ising machines","Astroid boundary governs bistable switch in Ising machines"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper assumes the local network field (h_x,i, h_y,i) is frozen when computing the spinodal boundary and barrier, whereas in a real coupled network this field evolves together with the oscillator's phase; if it fluctuates quickly, the quasi-static boundary may not describe actual bistability.","fun_headline_variants_meta":{"raw":{"variants":["Astroid spinodal: bistability frontier in Ising machines","Astroid curve dictates bistability in phase-based Ising machines","Bistability in Ising machines: an astroid spinodal boundary","Astroid spinodal sets bistable regime in Ising machines","Astroid boundary governs bistable switch in Ising machines"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001379,"raw_usage":{"total_tokens":5426,"prompt_tokens":748,"completion_tokens":4678,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":492,"completion_tokens_details":{"reasoning_tokens":4591}},"tokens_in":492,"tokens_out":4678,"duration_ms":30466,"temperature":1.0,"reasoning_tokens":4591,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T00:27:48.881379+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The decisive check is a numerical integration of a small coupled OIM or DIM network where the local field is allowed to evolve dynamically, compared against the frozen-field astroid prediction: if nodes spend time inside the frozen astroid but do not exhibit two stable wells, or vice versa, the quasi-static boundary is violated. Alternatively, escape-time measurements near a smooth crossing versus near the longitudinal cusp could test the predicted 3/2 versus 2 barrier scaling.","supporting_citations":[],"review_version":1}