{"id":"efdf2efa-83c4-442f-8f82-a1c8b4d12d9d","arxiv_id":"2603.17032","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Hierarchical probabilistic counters enable rejection-free Glauber Monte Carlo for the 2D RFIM with O(log N) spin selection and large low-T speedups over Metropolis.","lead":"A new Monte Carlo algorithm simulates the 2D random-field Ising model without rejections, using hierarchical counters to pick spins in logarithmic time. It targets the low-temperature regime where standard Metropolis updates become impractically slow.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"The abstract asserts that hierarchical probabilistic counters yield unbiased Glauber rates and correct continuous-time RFIM dynamics, but supplies no rate-bookkeeping argument once random fields destroy energy-class structure.","rationale":"The reader’s weakest-assumption statement already isolates the precise load-bearing gap: absence of a rate-equivalence argument once random fields break BKL’s energy-class structure. Because only the abstract is available, that gap cannot be closed or refuted; the verdict therefore remains UNVERDICTED. The concrete test above is the minimal analytical/numerical check that would decide whether the concern lands once the full text or code appears. No stronger objection (internal contradiction, numerical artifact, etc.) is visible from the abstract alone, so no verdict shift is warranted.","tokens_in":2048,"tokens_out":476,"duration_ms":10163,"concrete_test":"Extract (or request) the full rate-update and sampling pseudocode; implement a single-spin RFIM (N=1) and a small 4\times4 lattice with known Gaussian fields; compare the empirical flip-time histogram and magnetization autocorrelation against the exact continuous-time master-equation solution (or against a rejection-free Gillespie reference). Any statistically significant deviation in mean waiting time or spectral density falsifies the dynamical-equivalence claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of dynamical fidelity (and therefore of a fair speedup comparison) rests on the hierarchical counters selecting spins with probabilities exactly proportional to the local Glauber rates and advancing continuous time by the correct exponential waiting time. Classical BKL relies on a small number of energy classes whose rates are identical; a Gaussian random field makes every spin’s flip rate unique, so the hierarchy must maintain an exact, updatable sum of heterogeneous rates and sample from it without bias. The abstract states that the construction works and is superior to BKL for RFIM, yet contains neither the recursive update rules, the proof that the sampled rates match the master-equation rates, nor a demonstration that the continuous-time clock remains exact. Without that equivalence the method could be merely a fast approximate sampler, undermining both the “dynamically faithful” claim and the reported >100\times low-T speedups.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript proposes a rejection-free, event-driven Monte Carlo algorithm for the two-dimensional Random Field Ising Model (RFIM). It combines the Bortz–Kalos–Lebowitz (BKL) event-driven framework with Glauber transition probabilities and introduces hierarchical probabilistic counters that select spins in O(log N) operations. The abstract claims that this construction remains dynamically faithful even when a random field destroys the energy-class structure exploited by classical BKL, that Gaussian-RFIM runs reproduce the expected reduction of the pseudo-critical temperature with increasing disorder, and that speedups exceed two orders of magnitude relative to Metropolis in the low-temperature regime.","tokens_in":2233,"tokens_out":825,"duration_ms":17290,"significance":"If the hierarchical-counter construction is rigorously equivalent to continuous-time Glauber dynamics for heterogeneous rates, the method would supply a practically useful tool for both equilibrium and non-equilibrium studies of disordered Ising systems in regimes where Metropolis suffers severe critical slowing down. An O(log N) rejection-free sampler that does not rely on a small number of energy classes would also be of broader algorithmic interest beyond the RFIM. The claimed qualitative consistency with the known disorder-induced drop of the pseudo-critical temperature is a necessary but not sufficient validation target.","major_comments":[{"comment":"The abstract asserts that hierarchical probabilistic counters implement unbiased Glauber rates and preserve exact continuous-time RFIM dynamics once random fields render every flip rate unique. Classical BKL equivalence rests on a finite set of identical energy classes; the abstract supplies neither the recursive rate-update rules, a proof that the sampled probabilities remain proportional to the local Glauber rates, nor a demonstration that the exponential waiting-time clock remains exact. This equivalence is load-bearing for both the “dynamically faithful” claim and the fairness of the reported speedups; without it the method could be a fast approximate sampler.","section":"Abstract"},{"comment":"Speedups “exceeding two orders of magnitude” versus Metropolis are stated without system sizes, temperature and disorder ranges, wall-clock versus Monte-Carlo-step accounting, error bars, or a precise definition of the Metropolis baseline (single-spin vs. sweep, continuous-time vs. discrete-time). These details are required to judge whether the comparison is fair and whether the hierarchical overhead is correctly amortized.","section":"Abstract"},{"comment":"Reproduction of the expected reduction of the pseudo-critical temperature with Gaussian disorder is only a qualitative consistency check. The abstract does not report quantitative comparison to established literature values, finite-size scaling analysis, or any dynamical observable (autocorrelation times, domain-growth exponents) that would test continuous-time fidelity rather than mere equilibrium sampling.","section":"Abstract"}],"minor_comments":[{"comment":"The abstract should briefly distinguish hierarchical probabilistic counters from standard binary trees or Fenwick trees already used in other rejection-free Monte Carlo schemes, so that the claimed novelty is clear.","section":"Abstract"},{"comment":"Clarify whether the continuous-time clock is advanced by the exact sum of all rates or by an approximation, and whether the method is intended for equilibrium sampling only or also for genuine non-equilibrium trajectories.","section":"Abstract"}],"recommendation":"uncertain","confidential_remarks":"Only the abstract was available for this review; the full manuscript could not be examined. The load-bearing rate-equivalence claim and the quantitative benchmarks therefore remain unverified. I recommend obtaining the complete text (including algorithmic pseudocode, complexity analysis, and benchmark tables) before a final editorial decision. If the full paper supplies a correct bookkeeping proof and reproducible speedups, the work is potentially suitable after ordinary revision; if those elements are absent, rejection or major revision would be warranted."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"Punchline: this is a computational methods paper that tries to give you rejection-free, continuous-time Glauber dynamics for the 2D random-field Ising model by replacing classical BKL energy classes with hierarchical probabilistic counters that pick a spin in O(log N). If the bookkeeping is exact, that is a real tool for low-T and non-equilibrium RFIM work; classical BKL does not carry over cleanly once every local field is unique.\n\nWhat looks new is the combination, not the ingredients. BKL, Glauber rates, and tree-based event selection are known. The abstract’s point is that the hierarchy keeps an exact, updatable sum of heterogeneous flip rates so the sampler stays dynamically faithful under Gaussian disorder, where energy-class BKL fails. They report the expected drop of the pseudo-critical temperature with disorder and >100× speedups versus Metropolis at low T. For a methods paper those are the right qualitative checks, and the circularity burden is low: they are not fitting a quantity and then re-predicting it.\n\nThe soft spot is load-bearing and currently invisible. Once random fields destroy energy classes, you need recursive update rules and a short argument that the sampled rates and the exponential waiting time match the master equation. The abstract asserts that and stops. Without the full text we cannot see the rate bookkeeping, complexity analysis, error bars, or code. That is not a manufactured flaw; it is simply what an abstract-only read leaves open. Nothing in the abstract is internally contradictory, and the stress-test concern is exactly the right one to demand from the paper, not a reason to dismiss the idea.\n\nWho this is for: people who actually run disordered Ising dynamics at low T or out of equilibrium. It is not a physics breakthrough; it is subfield infrastructure. I would send it to a serious referee rather than desk-reject. If the proofs and reproducible artifacts are in the manuscript, it is a useful tool paper. If they are hand-waved, the referee will catch it. Worth a look when the full text is available; not something I would cite from the abstract alone.","headline":"Useful methods claim for 2D RFIM dynamics—rejection-free Glauber via hierarchical counters—but the abstract alone cannot verify the rate-equivalence that makes the speedups fair.","tokens_in":2876,"tokens_out":540,"would_cite":false,"duration_ms":10853,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.10.Ln","05.50.+q","75.10.Nr","75.40.Mg"],"model":"grok-4.5","headline":"A rejection-free Glauber Monte Carlo for the 2D random-field Ising model selects spins in O(log N) with hierarchical probabilistic counters, speeding low-temperature dynamics by over 100× while staying faithful to continuous-time evolution.","keywords":["random-field Ising model","Glauber dynamics","rejection-free Monte Carlo","Bortz-Kalos-Lebowitz","hierarchical counters","event-driven simulation","disordered spin systems","critical slowing down"],"falsifier":"Compare the algorithm’s measured flip-rate histogram and long-time magnetization autocorrelation against an exact enumeration (or a carefully equilibrated Metropolis run) on a small lattice with known Gaussian disorder; any systematic deviation in the waiting-time distribution or in the temperature dependence of the susceptibility peak would falsify dynamical fidelity.","tokens_in":2901,"feed_emoji":"⚡","tokens_out":919,"duration_ms":9004,"temperature":0.7,"pith_summary":"This paper introduces a Monte Carlo algorithm for the two-dimensional random-field Ising model that is both rejection-free and event-driven under Glauber rates. Classical Metropolis sampling becomes painfully slow at low temperature because most proposed flips are rejected; the classic Bortz–Kalos–Lebowitz method avoids rejections for the pure Ising model by grouping spins into energy classes, but random fields destroy those classes. The authors restore rejection-free, continuous-time dynamics by replacing the energy-class bookkeeping with hierarchical probabilistic counters that pick the next spin to flip in O(log N) operations. The resulting scheme therefore remains dynamically faithful even when quenched disorder is present. In the low-temperature, low-disorder regime the method yields wall-clock speed-ups of more than two orders of magnitude relative to Metropolis, and it recovers the expected downward shift of the pseudo-critical temperature as the Gaussian random-field strength is increased. The practical consequence is that both equilibrium thermodynamics and non-equilibrium relaxation of disordered Ising systems become accessible on much larger lattices and longer timescales.","feed_headline":"Rejection-free RFIM Monte Carlo flips spins 100× faster","feed_subtitle":"Hierarchical counters restore continuous-time Glauber dynamics for disordered Ising systems at low temperature.","key_machinery":"Hierarchical probabilistic counters: a tree-structured data structure that stores cumulative flip probabilities and permits selection of the next spin (and the waiting time) in O(log N) operations under site-dependent Glauber rates, thereby replacing the energy-class lists of the pure BKL algorithm.","core_discovery":"An event-driven, rejection-free Monte Carlo that employs hierarchical probabilistic counters can implement Glauber dynamics for the two-dimensional random-field Ising model in O(log N) per accepted flip, delivering more than 100-fold acceleration over Metropolis at low temperature while preserving the correct continuous-time master-equation evolution that classical energy-class BKL cannot maintain once random fields are present.","pith_inferences":["The same O(log N) selection structure should extend immediately to three-dimensional RFIM or to models with random bonds, provided the rate tree can be updated after each flip.","Because waiting times are drawn from the exact continuous-time distribution, the method supplies a natural clock for measuring physical aging exponents without arbitrary Monte-Carlo-step rescaling.","A parallel or GPU realization of the hierarchical counters would further enlarge accessible system sizes for finite-size scaling of the RFIM critical line."],"forward_implications":["Low-temperature and low-disorder RFIM dynamics can be followed for wall-clock times two orders of magnitude longer than with Metropolis.","Pseudo-critical temperatures extracted from Gaussian-disorder sweeps can be mapped with higher statistics, confirming the expected depression with increasing field strength.","Both equilibrium sampling and non-equilibrium aging or coarsening protocols become practical for larger two-dimensional lattices.","The same counter hierarchy can be reused for other disordered Ising models whose transition rates lack a simple energy-class partition."],"fun_headline_variants":["Hierarchical counters cut RFIM flips to O(log N)","Rejection-free Glauber MC speeds 2D RFIM 100x at low T","Event-driven counters restore continuous-time RFIM dynamics","Probabilistic counters enable rejection-free RFIM sampling","O(log N) Glauber flips for disordered Ising without rejections"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The hierarchical counters correctly implement the site-dependent Glauber rates without bias and keep the event-driven process equivalent to the underlying continuous-time master equation once random fields destroy the simple energy-class structure.","fun_headline_variants_meta":{"raw":{"variants":["Hierarchical counters cut RFIM flips to O(log N)","Rejection-free Glauber MC speeds 2D RFIM 100x at low T","Event-driven counters restore continuous-time RFIM dynamics","Probabilistic counters enable rejection-free RFIM sampling","O(log N) Glauber flips for disordered Ising without rejections"]},"model":"grok-4.5","effort":"low","cost_usd":0.00557,"raw_usage":{"total_tokens":1487,"prompt_tokens":742,"num_sources_used":0,"completion_tokens":75,"cost_in_usd_ticks":55700000,"prompt_tokens_details":{"text_tokens":742,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":670,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":742,"tokens_out":75,"duration_ms":6860,"temperature":1.0,"reasoning_tokens":670,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T23:23:55.758673+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compare the algorithm’s measured flip-rate histogram and long-time magnetization autocorrelation against an exact enumeration (or a carefully equilibrated Metropolis run) on a small lattice with known Gaussian disorder; any systematic deviation in the waiting-time distribution or in the temperature dependence of the susceptibility peak would falsify dynamical fidelity.","supporting_citations":[],"review_version":1}