{"id":"a74d52f4-9b1e-41e0-946e-83e049cb527f","arxiv_id":"2512.18917","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a nonreciprocal Blume-Capel magnet, an antisymmetric single-ion anisotropy suppresses oscillatory dynamics and restores a static ordered phase belonging to the 2D Ising universality class.","lead":"This paper studies a two-species spin model with nonreciprocal interactions and shows that an imbalance in chemical potentials between species can restore static magnetic order that would otherwise be destroyed by nonreciprocity. It matters because it identifies vacancy energetics as a practical control knob for stabilizing equilibrium-like order in driven many-body systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermodynamic-limit stability of the restored static order at Δ=1.5 is not fully established; the apparent static phase could still be a slow droplet-induced swap.","rationale":"The reader's verdict was already CONDITIONAL, and the reader's rationale explicitly flagged the unresolved thermodynamic-limit fate of droplet-induced swap phases. My concern sharpens that issue: it applies directly to the parameter point used for the central 2D Ising universality claim, not only to the 3D or off-slice regimes. This is a genuine soft spot because the authors demonstrate at Δ=0 how finite-size data can falsely suggest ordering, and their suppression mechanism at Δ=1.5 is argued physically rather than established by a systematic finite-size/long-time analysis. However, the concern is addressable: a modest set of larger-L, longer-time simulations would likely settle it. The FSS exponents and data collapse in Fig. 11 are otherwise internally consistent, and the Kolmogorov-criterion calculation is clean. I do not see a reason to change the verdict: CONDITIONAL remains appropriate, pending the proposed check. Other issues raised by the reader (single K̃ value, cubic-spline smoothing, single-realization 3D data) are secondary relative to the stability of the static phase itself.","tokens_in":18018,"tokens_out":11284,"duration_ms":123239,"concrete_test":"Perform long-time Monte Carlo simulations at Δ=1.5, K̃=0.3 for L=256 and L=384, at a coupling inside the putative ordered phase (e.g., J̃≈2.5) and near the critical point (J̃≈Jc from Fig. 11), running at least 2×10^7 sweeps. Measure R(t), S(t), and the autocorrelation time of R(t); also check the power spectrum of R(t) for oscillatory peaks. If R extrapolates to a nonzero value as L→∞, S remains zero, and no oscillation period grows with L, the static-order claim survives. If R decays with L or an oscillatory component appears on long timescales, the central claim fails because the phase would be a droplet-induced swap rather than static order.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central 2D claim—that a finite antisymmetric anisotropy Δ restores a robust static ordered phase and that the disorder-to-static transition is 2D Ising—rests on Fig. 7 and the FSS analysis in Fig. 11 at Δ=1.5, K̃=0.3. The evidence is that R converges to a nonzero value for L up to 120 while S tends to zero. But the authors' own Δ=0 results (Fig. 6) show that exactly this kind of finite-size behavior can be misleading: both R and S decrease with L, and the text states that the droplet-induced swap phase occurs only in finite systems and that spiral defects destroy all oscillatory order in the thermodynamic limit. At Δ=1.5 the physical argument in Section IV—that vacancies in species A suppress the nonreciprocal bias on species B—is plausible but not quantitative. The residual nonreciprocal coupling is still present, and the simulations run only 4×10^6 sweeps with L≤160. If the system is actually in a droplet-induced swap phase with a period that grows with L, then R could appear static on accessible timescales while actually decaying in the thermodynamic limit, and the 2D Ising FSS would not describe the true steady state. This is the most load-bearing soft spot because the universality conclusion is built on the assumption that the ordered phase is genuinely static.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a two-species spin-1 Blume–Capel model with fully antisymmetric nonreciprocal onsite coupling (K_AB = -K_BA) and equal-and-opposite single-ion anisotropies (Δ_A/k_BT = -Δ_B/k_BT = Δ). A mean-field theory, Eq. (9), is developed and analyzed by bifurcation theory, yielding disorder, swap (limit-cycle), and static ordered phases separated by Hopf, SNIC, SNLC, pitchfork, and saddle-node bifurcations. The Kolmogorov criterion is shown to be violated with loop affinity 8K̃, Eq. (14). Monte Carlo simulations in 2D show that at Δ=0 spiral defects destroy both swap and static order, while at Δ=1.5, K̃=0.3 a static ordered phase is claimed, with finite-size scaling (Fig. 11) giving exponents close to 2D Ising values. A crossover-to-first-order line within the ordered phase is also reported, ending at a critical point. In 3D, simulations show a stable swap phase for small Δ and a swap→disorder→static sequence as Δ increases.","tokens_in":18307,"tokens_out":5604,"duration_ms":63273,"significance":"If the central 2D claim holds, the paper makes a substantial contribution: it identifies a simple, experimentally plausible local parameter—species-dependent vacancy energetics—that suppresses time-dependent nonreciprocal dynamics and restores equilibrium-like static order and 2D Ising universality in a manifestly non-Hamiltonian model. The mean-field bifurcation analysis is systematic, and the Kolmogorov-criterion calculation (Eq. 14) is clean and explicit. The Monte Carlo evidence for the disorder-to-static transition is strong in the finite-size regime: multiple system sizes, Binder-cumulant crossings, data collapse, and exponents γ=1.779±0.084, ν=1.013±0.047, β=0.120±0.003 are all close to the known 2D Ising values. The authors are also candid about several limitations, including single-realization 3D runs and the open thermodynamic-limit fate of droplet-induced swap regions. The main weakness is that the genuine static nature of the 2D ordered phase at Δ=1.5 is not established as rigorously as the critical scaling itself.","major_comments":[{"comment":"The load-bearing claim is that at Δ=1.5 the 2D ordered phase is genuinely static, not a slow droplet-induced swap. The authors' own Δ=0 analysis (Fig. 6) shows that finite-size behavior of R and S can masquerade as order even when the thermodynamic limit is disordered. At Δ=1.5 the evidence is R converging to a nonzero value while S→0, but S→0 is also observed at Δ=0 in the droplet-induced swap regime. If the swap period grows with L beyond the 4×10^6-sweep observation window, the same diagnostics would look static on accessible timescales. Because the 2D Ising FSS in Fig. 11 is built on the assumption of true static order, please provide a direct test of time-independence: e.g., magnetization autocorrelation functions, the L-dependence of domain-wall or phase-slip times, or runs at least one order of magnitude longer for representative L. A scan in Δ showing an L-independent nonzero R o","section":"IV, Fig. 7 and Fig. 11"},{"comment":"The crossover→first-order→critical-point scenario is under-supported. The FSS slope 1.995±0.013 at J̃=7.5 is consistent with a first-order transition, but no Binder-cumulant dip or order-parameter histogram is shown, and the claimed critical point terminating the first-order line is not located. Since the abstract and conclusion feature the liquid-gas analogy prominently, direct evidence for the critical endpoint—such as histogram shapes or scaling of the discontinuity near the terminus—is needed before this part of the phase diagram is accepted.","section":"Fig. 14 and surrounding text"},{"comment":"The 3D results are averaged over a single realization (stated in Section IV). This is acceptable for a qualitative illustration, but the negative claim that swap→static ordering never occurs directly relies on the absence of a feature in one run. Multiple independent seeds are needed to rule out a rare direct transition and to justify the 'mirrors mean-field expectations' conclusion. Additionally, Fig. 9(b) is described as showing a direct transition from the droplet-induced swap phase to static order, which appears to contradict the earlier statement that swap→static ordering is indirect; please clarify whether the direct path is a property of the droplet-induced swap region only.","section":"IV, 3D simulations"}],"minor_comments":[{"comment":"Captions read 'The provides the value of the critical exponent β'—'This provides' is intended.","section":"Fig. 11(d) and Fig. 12(d)"},{"comment":"The term 'scaling Δ_A/kBt=Δ' appears to be a typo for k_B T.","section":"Fig. 5 caption"},{"comment":"The label 'IVI' in the first row is confusing; it is likely a typographical artifact for 'IV' or a region label and should be cleaned up.","section":"Fig. 1 caption"},{"comment":"The susceptibility maxima used for FSS are obtained from cubic-spline interpolation of already collected data. Please state explicitly how the interpolation error or choice of smoothing affects the reported exponents; for example, show that the results are stable under alternative interpolation or raw-data peak estimates.","section":"IV, FSS methods"},{"comment":"All 2D simulations use K̃=0.3 as a single representative value. The mean-field analysis suggests the phase topology is independent of K̃, but the 2D universality claim is verified only on this slice. A sentence placing K̃=0.3 in the context of the full K̃ dependence would help the reader judge generality.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid candidate for a statistical mechanics journal. The mean-field analysis is sound, the Kolmogorov-criterion violation is clearly derived, and the finite-size scaling at Δ=1.5 is credible within the simulated sizes. My main reservation is the thermodynamic-limit stability of the static phase: the same finite-size diagnostics that distinguish genuine order from droplet-induced swap at Δ=0 are not sufficient at Δ=1.5 without a direct time-scale study. I do not see a circularity problem—the FSS is benchmarked against known Ising exponents—but the universality claim would be greatly strengthened by ruling out a slow-swap alternative. The 3D single-realization issue is secondary but should be addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honestly, this is a solid paper. The new model is a natural Blume-Capel extension of the nonreciprocal Ising model, and the finite-Δ physics is genuinely absent from previous work. The mean-field bifurcation analysis is systematic, and the Kolmogorov-criterion calculation (Eq. 14, loop affinity 8K̃) is clean. The 2D Monte Carlo finite-size scaling is credible: several system sizes up to L=160, Binder crossings, data collapse, and exponents close to Ising values. The contrast between Fig. 6 (Δ=0, R and S both decay with L) and Fig. 7 (Δ=1.5, R converges, S vanishes) makes the qualitative point convincingly.\n\nThe main soft spot is the one the stress-test hangs on: the restored static order at Δ=1.5 may not be truly static in the thermodynamic limit. The authors' own Δ=0 analysis shows that exactly this kind of finite-size behavior can be misleading—the apparent droplet-induced swap phase is finite-size only, destroyed by spiral defects at large L. At Δ=1.5, the physical argument that vacancies in species A suppress the nonreciprocal bias on B is plausible but not quantitative, and the simulations run only 4×10^6 sweeps. If the swap period grows with L, R could look converged while actually decaying. I don't think this is fatal—the FSS and the R-vs-L contrast provide real evidence—but it is the part a referee should push on.\n\nOther soft spots are proportionately minor: the 3D runs are single-realization with no error bars; the critical point at the end of the first-order line is not located (the authors admit this); the spiral-defect mechanism is asserted rather than directly observed; and the universality claim rests on one K̃ value and spline-smoothed susceptibility. None of these sink the paper.\n\nThe authors are honest about the unresolved regime—they explicitly flag the thermodynamic-limit fate of droplet-mediated swap. The paper is worth a serious referee. I'd send it out with a request for more thorough evidence on the Δ=1.5 static phase: longer runs, larger L, or a direct measurement of the swap period scaling.","headline":"The paper's central claim—antisymmetric single-ion anisotropy restores a robust static ordered phase in the 2D nonreciprocal Blume-Capel model—is plausible and mostly well supported, but the thermodynamic-limit stability of that phase is the load-bearing soft spot.","tokens_in":18916,"tokens_out":3336,"would_cite":true,"duration_ms":31506,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A finite antisymmetric single-ion anisotropy restores a stable static ordered phase in the two-dimensional nonreciprocal Blume-Capel model, with the disorder-to-static transition falling in the 2D Ising universality class.","keywords":["nonreciprocal interactions","Blume-Capel model","single-ion anisotropy","swap phase","limit cycle","2D Ising universality","finite-size scaling","vacancy energetics"],"falsifier":"Run 2D Monte Carlo at fixed K̃ = 0.3 and Δ = 1.5 with lattice sizes beyond L = 160 and check whether the order parameter R converges to a size-independent nonzero value; and repeat the 3D swap-to-static sequence at Δ = 0.5, J̃ = 2.2 with many independent realisations instead of one. If R decays with L or the 3D intermediate disorder disappears under averaging, the paper's central stabilization claim fails.","tokens_in":17815,"feed_emoji":"🧲","tokens_out":9569,"duration_ms":85808,"temperature":0.7,"pith_summary":"This paper tries to show that a local, experimentally realistic term — an energy penalty for empty sites that has opposite sign for the two spin species — can reverse the destructive effect of nonreciprocal interactions. In a two-dimensional two-species Blume-Capel model (a lattice spin model with values ±1 and 0, the 0 acting as a vacancy), nonreciprocal couplings alone drive the system to disorder: spiral defects kill both the oscillatory 'swap' phase and static order. When the anisotropy ΔA = −ΔB is switched on, one species preferentially occupies vacancies, the nonreciprocal bias on the other species weakens, and a stable static ordered phase appears; the disorder-to-order transition has exponents matching the 2D Ising universality class. Mean-field theory predicts a richer landscape of phases and bifurcations, and the paper shows how fluctuations and defects reshape that picture. The broader claim is that vacancy energetics is a control knob that can restore equilibrium-like critical behavior in an intrinsically non-equilibrium system.","feed_headline":"Vacancy bias restores static order in a nonreciprocal spin model","feed_subtitle":"Species-asymmetric energy for empty sites suppresses swap oscillations and yields 2D Ising critical exponents.","key_machinery":"Two ingredients carry the argument. First, a compact loop-ratio identity: for the elementary two-species spin cycle ↑↑→↑↓→↓↓→↓↑→↑↑, the product of forward transition rates divided by the backward product equals e^{8K̃}, so when the nonreciprocal coupling K̃ is nonzero the Kolmogorov detailed-balance condition is violated and the system runs out of equilibrium with a cycle affinity 8K̃. Second, the vacancy-bias mechanism: the anisotropy ΔA = −ΔB changes the occupancy statistics of the S=0 state so that species A sits in vacancies more often, reducing the effective nonreciprocal field felt by species B and damping the oscillations; this is what allows static order to survive in 2D. Mean-field","core_discovery":"The central claim is that in the two-dimensional nonreciprocal Blume-Capel model with fully antisymmetric interspecies coupling (K_AB = −K_BA > 0), the equal-and-opposite single-ion anisotropy ΔA/kBT = −ΔB/kBT = Δ acts as a vacancy bias that suppresses nonreciprocal swap dynamics and restores a stable statically ordered phase. At Δ = 0, Monte Carlo simulations show that both the synchronisation order parameter R and the angular-momentum-like order parameter S decay with system size, indicating that in the thermodynamic limit spiral defects destroy global swapping and long-range order. At finite Δ (illustrated for Δ = 1.5, K̃ = 0.3), R converges to a nonzero value with increasing lattice size","pith_inferences":["If the vacancy-bias mechanism is generic, the anisotropy strength needed to restore static order should scale with the nonreciprocal coupling; a systematic scan of the critical line Δc(K̃) would provide a direct test.","The same mechanism may transfer to other nonreciprocal lattice models with an empty state: adding a species-antisymmetric on-site potential could act as a local 'brake' on nonreciprocity and yield equilibrium-like criticality, without requiring fine-tuned interaction asymmetries.","The liquid-gas-like critical point inside the ordered phase may be a generic feature of spin-1 nonreciprocal models; locating it precisely with larger simulations and extracting its exponents would clarify whether it belongs to the Ising or another universality class.","Because the 3D conclusions rest on single-realization averages, an independent multi-realization study would settle whether the intermediate disordered regime between swap and static order is a genuine phase or a finite-sample artifact."],"forward_implications":["If the stabilization claim is correct, an antisymmetric single-ion anisotropy is a sufficient control parameter to suppress the time-dependent swap phase and recover equilibrium-like Ising criticality in a 2D nonreciprocal system.","The measured exponents (γ ≈ 1.78, ν ≈ 1.01, β ≈ 0.12) mean that at the disorder-to-static transition the model behaves like the ordinary 2D Ising model, so standard equilibrium universality applies despite the absence of detailed balance.","The first-order line that emerges inside the static ordered phase at large coupling ends at a critical point; the liquid-gas analogy implies a well-defined endpoint whose location and universality could be probed.","The 3D results indicate that a stable swap phase coexists with static order but connects to it only through an intervening disordered state, so any theory of the swap-to-order transition must account for an intermediate fluctuation-dominated regime.","At zero anisotropy, the thermodynamic-limit disorder in 2D is driven by spiral defects, which means finite-size simulations that appear to show swapping at Δ = 0 are not evidence of a true phase."],"fun_headline_variants":["Antisymmetric anisotropy quells spin swaps, restores order","Vacancy bias suppresses swap phase, yields 2D Ising order","Nonreciprocal model tamed by species-asymmetric vacancy bias","Spin-swap chaos stopped by asymmetric empty-site energy","Defect-driven order in nonreciprocal Blume-Capel model"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the equal-and-opposite anisotropy slice ΔA/kBT = −ΔB/kBT = Δ, plus the particular simulated parameter points, represents the general mechanism; if the vacancy-suppression effect does not persist for arbitrary anisotropies or other couplings, the claimed control-knob interpretation is not established.","fun_headline_variants_meta":{"raw":{"variants":["Antisymmetric anisotropy quells spin swaps, restores order","Vacancy bias suppresses swap phase, yields 2D Ising order","Nonreciprocal model tamed by species-asymmetric vacancy bias","Spin-swap chaos stopped by asymmetric empty-site energy","Defect-driven order in nonreciprocal Blume-Capel model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000228,"raw_usage":{"total_tokens":1362,"prompt_tokens":842,"completion_tokens":520,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":586,"completion_tokens_details":{"reasoning_tokens":442}},"tokens_in":586,"tokens_out":520,"duration_ms":5201,"temperature":1.0,"reasoning_tokens":442,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:49:04.125516+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run 2D Monte Carlo at fixed K̃ = 0.3 and Δ = 1.5 with lattice sizes beyond L = 160 and check whether the order parameter R converges to a size-independent nonzero value; and repeat the 3D swap-to-static sequence at Δ = 0.5, J̃ = 2.2 with many independent realisations instead of one. If R decays with L or the 3D intermediate disorder disappears under averaging, the paper's central stabilization claim fails.","supporting_citations":[],"review_version":1}