{"id":"93431af1-0b46-49c4-a7e2-8ee7e37e0bcc","arxiv_id":"2512.02384","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For spiked Wigner inference, Glauber dynamics and AMP reach the same correlation fixed point, with a phase transition at βλ=1 conditional on SK mixing.","lead":"Glauber dynamics on a noisy matrix inference problem is shown to track the same correlation fixed point as Approximate Message Passing, with a sharp transition when the temperature and signal strength cross. The proof works through an auxiliary Gaussian chain and assumes a widely believed mixing condition for the Sherrington–Kirkpatrick model.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 5.2 hinges on unproven SK MLSI for β∈(0.295,1); without it the Glauber-to-RGD transfer in Lemma 4.8 has no polynomial conductance.","rationale":"The reader's weakest_assumption (Condition 1/MLSI) is also the most load-bearing step I see. The paper's unconditional contributions—the RGD reduction to a one-dimensional recursion, the fixed-point characterization, and the magnetization estimates—are substantial and appear internally coherent. But the headline Glauber theorem is conditional on an open MLSI conjecture for SK in 0.295<β<1. I also note two secondary issues that do not change the verdict: (i) Theorem 5.2 states T≥e^{Ω(N^{4-2ε})} whereas the proof's error bound only requires T=Ω(N^{4-2ε}); the exponential statement appears to be a typo but currently contradicts the §1.2 claim of polynomial-time control. (ii) The proof of Lemma 5.6 applies Lemma 5.7 with a Y_t that includes deterministic O(N^{-1/2+ε}) errors; this formally violates Lemma 5.7(ii), though the escape argument seems repairable by a two-phase analysis. Neither issue undermines the central conditional result as much as the open MLSI assumption, which is explicitly acknowledged. Hence the verdict should remain CONDITIONAL.","tokens_in":31514,"tokens_out":40259,"duration_ms":403773,"concrete_test":"Estimate the MLSI constant of Glauber dynamics for the SK model at β=0.5, h=0 and h=1 for N=10^3..10^5 via the Chen-Eldan localization scheme (or spectral gap plus entropy-decay diagnostics). If the constant decays like N^{-c} with c>1 (let alone exponentially), Condition 1 is false and the transfer argument in Theorem 5.2 would need a correspondingly larger polynomial T, so the theorem's stated runtime and the §1.2 polynomial-time claim would have to be revised; if c=1, the stated transfer with T=N^{4-2ε} is consistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 5.2's proof transfers local stationarity from Glauber dynamics to RGD via Lemma 4.8, whose denominator is δ=inf_h ρLS(μβW,h1). This is exactly Condition 1 (MLSI). The paper proves MLSI only for β<0.295 via AKV24; the full high-temperature regime β<1 is a conjecture. If the true MLSI constant is not polynomial (or worse, exponential), then the E_RGD(f,logf)≤O(τ/δ·N/T) bound in the proof of Theorem 5.2 can be made small only with T super-polynomial, and the Corollary 4.7 error term O(logN sqrt(N^3/T)) (or with N^{c+2}/T) cannot be driven to o(N^{-1/2+ε}) while maintaining the paper's polynomial-time narrative. Since Theorem 5.2 is the unique bridge from the RGD analysis to the actual Glauber trajectory, the claim that Glauber recovers the BBP transition for all β>1/λ is only as secure as this conjecture. The paper is transparent about the dependence, but it is still the least-secure load-bearing assumption: the RGD fixed-point analysis may be fully correct while Theorem 5.2 fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the spiked Wigner inference problem and proposes a formal connection between the correlation dynamics of Markov chains and Approximate Message Passing (AMP). The main object is an auxiliary chain, Restricted Gaussian Dynamics (RGD), whose one-step update is shown (informally in Section 2.2) to mirror the AMP state-evolution recursion. The paper derives a one-dimensional recursion for the RGD correlation, proves under certain conditions that RGD reaches the largest fixed point OPT_{β,λ} from a warm start, and then transfers this to Glauber dynamics using locally stationary distributions, conditional on a modified log-Sobolev inequality for the Sherrington-Kirkpatrick model. It also proves high-precision mean-magnetization estimates for the SK model under weak external fields and analyzes the qualitative fixed-point structure of the recursion.","tokens_in":31798,"tokens_out":11586,"duration_ms":121253,"significance":"If the main theorems are correct, this is a substantial step toward a rigorous explanation of the empirical success of MCMC in Bayesian inference: it identifies the exact fixed-point correlation of RGD/Glauber with the AMP state-evolution fixed point and recovers the BBP transition. The RGD analysis is technically rich, particularly the weak-field magnetization estimates (Lemma 6.1 and Lemma 6.3) and the fixed-point analysis of the scalar recursion. The paper is commendably explicit about its conditional dependencies: Theorem 5.2 is stated only under the MLSI conjecture for the SK model, and the RGD-AMP comparison is presented as informal with a deferred numerical condition. These caveats are real and limit the current certainty of the headline Glauber claim, but the underlying RGD fixed-point analysis appears to be a serious contribution.","major_comments":[{"comment":"The proof of Theorem 5.2 transfers local stationarity from Glauber to RGD via Lemma 4.8, whose denominator δ is exactly the MLSI constant of the SK model with external field. The paper uses δ=Ω(1/N), which is Condition 1. As noted in the paper, this is proved only for β<0.295 (via AKV24) and is conjectured for all β<1. Thus the statement that Glauber recovers the BBP transition for all β>1/λ is conditional on an unproved conjecture in exactly the regime where it would extend beyond prior work. If the MLSI constant is not polynomial, the bound E_RGD ≤ O(N^3/T) cannot be made small enough to yield the claimed o(N^{-1/2+ε}) error with the stated exponential T. The paper is transparent about this, but the abstract and introduction should present the Glauber result as explicitly conditional on a conjecture, not as a settled recovery theorem.","section":"§5.2, Theorem 5.2, Condition 1"},{"comment":"The claimed comparison 'one step of RGD ≈ one step of AMP' is never formalized. Theorem 1.7 is marked 'Informal, see Section 2.2', but Section 2.2 contains only algebra and a reference to a numerical condition; there is no formal theorem statement or proof. This is one of the paper's three advertised main contributions and it is the conceptual bridge that justifies the AMP fixed-point identification. Please state a precise formal theorem (with the required numerical condition clearly isolated) or explicitly demote this to a conjecture. As written, the reader cannot verify the central conceptual claim.","section":"§2.2, Theorem 1.7"},{"comment":"The proof of monotonicity states f'_β(h) ≥ βλ( E[S^4] − E[S^2 T]^2 / E[S^2 T^2] ) ≥ 0, with the final step attributed to Cauchy–Schwarz. This inequality is not justified as written: Cauchy–Schwarz gives E[S^2] E[S^2 T^2] ≥ E[S^2 T]^2, not E[S^4] E[S^2 T^2] ≥ E[S^2 T]^2. The displayed term should likely be E[S^2] instead of E[S^4]. Since Lemma 7.5(ii) is used in Lemma 5.3 and hence in the convergence proof (Lemma 5.4) and Proposition 7.1, this is a load-bearing proof step and must be corrected. If the intended argument uses the AT condition in the form β^2 E[S^4]<1, the proof needs to be written out carefully.","section":"§7.1, Lemma 7.5(ii), proof"},{"comment":"The introduction frames Question 1.2 around polynomial-time Glauber dynamics, but Theorem 5.2 requires T ≥ e^{Ω(N^4)} (and e^{Ω(N^3)} for the o(1) version in Remark 1.6). These are exponential runtimes, not polynomial. The paper should not imply that it resolves the polynomial-time question; it proves a characterization at super-polynomial time. If the intended message is 'progress' toward Question 1.2, that is fair, but the abstract and introduction should avoid giving the impression that polynomial-time Glauber success has been established. This is a framing issue, but it affects how the contribution is received.","section":"§1, Question 1.2; §5.2, Theorem 5.2"}],"minor_comments":[{"comment":"The term 'annealed posterior' is used without definition. The distribution μ_{βM} for β≠λ is not the true posterior; clarify the temperature parameterization early.","section":"Abstract / §1"},{"comment":"The statement says 'max_{0≤t≤T} 1/N |⟨x_T,1⟩| ≤ ε' but the subscript of x should be t, not T. Please fix the typo.","section":"§5, Lemma 5.6"},{"comment":"The proof uses the symbol '(∗)' before it is defined. Define the expression explicitly or introduce a numbered equation.","section":"§6, Lemma 6.4 proof"},{"comment":"The four panels are not accompanied by explicit parameter values or the fixed-point locations. Add axis labels and a short caption explaining the 'AT region' boundary.","section":"§2, Figure 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its main conditional assumption, but the introduction somewhat overstates the scope. The missing formal statement of the RGD-AMP comparison and the apparent typo in Lemma 7.5(ii) need to be addressed before I can recommend acceptance. The RGD fixed-point analysis itself is promising and likely correct after the proof repairs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line first: the paper has a genuinely new result — an exact fixed-point characterization of the correlation achieved by restricted Gaussian dynamics (RGD) on the spiked Wigner posterior, shown to match the AMP state evolution fixed points. The one-dimensional recursion in (4) is the right object, and the identification with AMP comes from algebra and standard state evolution, not from fitting; the circularity burden is genuinely low. The escape-from-the-unstable-fixed-point analysis (Lemma 5.6) and the high-precision SK magnetization estimates (Lemma 6.1) are real technical value, and the magnetization bounds are of independent interest. The paper is also honest about its limits: it states plainly that the Glauber-to-RGD transfer breaks for β>1 and that Condition 1 (SK MLSI) is only proved for β<0.295.\n\nThe soft spots, in proportion. The biggest: the Glauber theorem (5.2) is conditional on MLSI for β∈(0.295,1), and that conjecture is load-bearing — δ in Lemma 4.8 is exactly the MLSI constant, so if the constant is not polynomial the transfer is empty. The paper discloses this, so it is not a hidden flaw, but the abstract's 'we recover the phase transition' is stronger than what is conditional. Disclosed-but-oversold.\n\nSecond, the runtime claim and Theorem 5.2 don't agree. The introduction asks whether polynomial-time Glauber matches AMP; Theorem 5.2 requires T ≥ e^{Ω(N^{4−2ε})}, which is exponential. The proof sketch in Section 5 looks like it gives a bound on the order of (C log N)√(N^3/T), which would make polynomial T enough. So something is off — either the statement or the proof — and it needs fixing.\n\nThird, Theorem 1.7 (\"one step of RGD ≈ one step of AMP\") is informal, deferred to a numerical condition. The algebra behind it is solid, but the headline version overstates the formal content.\n\nMinor: RS-AT and AT-fixed-pt for β∈(1,λ) rest on numerical checks; the authors should ship verifiable code. And since the acknowledgments disclose AI-assisted proofs of several lemmas, a referee should check those lemmas especially carefully — the disclosure is fine, but it raises the verification bar, not lowers it.\n\nIf the central RGD analysis holds up — and it looks like it does — this is a serious paper for a serious referee. I'd send it out; I'd expect a revision that reconciles the runtime statement and tightens the informal claims.","headline":"The fixed-point bridge from RGD to AMP is real and the RGD analysis is strong; the Glauber half is honestly conditional on SK mixing, and the paper's own runtime statements don't line up, but this deserves a serious referee.","tokens_in":32262,"tokens_out":10608,"would_cite":true,"duration_ms":103914,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60J10","82B44","60K35","68W20"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that Glauber dynamics on the spiked Wigner posterior recovers the spike with correlation pinned to the largest fixed point of an AMP-like recursion, and that this correlation is positive exactly when β>1/λ, recovering the","keywords":["spiked Wigner","Glauber dynamics","approximate message passing","state evolution","restricted Gaussian dynamics","Sherrington–Kirkpatrick model","phase transition","locally stationary distributions"],"falsifier":"For a fixed λ>1 and β just above 1/λ, numerically compute OPT_{β,λ} from the fixed-point equations, then run RGD from a warm start for T=log²N steps on a planted spiked Wigner instance and measure the mean absolute correlation; if it fails to approach OPT_{β,λ} within N^{-1/2+ε} as N grows, the one-dimensional reduction is wrong.","tokens_in":31398,"feed_emoji":"🎲","tokens_out":4948,"duration_ms":45466,"temperature":0.7,"pith_summary":"The paper establishes a formal equivalence between the short-term behavior of two very different algorithms on the spiked Wigner inference problem: Glauber dynamics, a Markov chain that attempts to sample a high-dimensional posterior, and Approximate Message Passing (AMP), an iterative estimator known to be Bayes-optimal. The bridge is an auxiliary chain, restricted Gaussian dynamics (RGD), whose correlation with the true spike is shown to follow a one-dimensional recursion identical to AMP's state evolution. Using that recursion, the paper proves that after polynomially many steps, Glauber iterates have correlation within N^{-1/2+ε} of the AMP fixed-point value OPT_{β,λ}, with OPT_{β,λ}>0 exactly when β>1/λ. A reader should care because this supplies a sharp, rigorous characterization of what a slow-mixing Markov chain achieves long before it equilibrates, and recovers the BBP transition as a byproduct.","feed_headline":"Glauber dynamics matches AMP's spike recovery at the BBP threshold","feed_subtitle":"Polynomial-time MCMC settles at the same correlation fixed point as message passing, sharp at λ=1.","key_machinery":"The central object is restricted Gaussian dynamics (RGD), an auxiliary Markov chain whose transition from σ samples a Gaussian field z = βλ⟨σ,x⟩/N + sqrt(βλ/N)g and then redraws σ from the SK Gibbs measure with that field. The paper reduces RGD to the one-dimensional map z ↦ E tanh(βg√q_{β,βλz} + βλz) using new high-precision mean-magnetization estimates for the SK model under weak external field (Lemma 6.1). That map has fixed-point equations identical to AMP's state evolution, and its stable fixed point controls the correlation reached by both RGD and—via transfer of local stationarity—Glauber dynamics.","core_discovery":"The central claim is that the correlation between a Glauber iterate and the planted spike is asymptotically pinned to a constant OPT_{β,λ} defined as the largest fixed point of the scalar recursion z = E tanh(βg√q + βλz) and q = E tanh²(βg√q + βλz), where g is a standard Gaussian and q is the SK overlap parameter. The same equations are the state-evolution fixed points of mismatched AMP, so the paper proves that Glauber dynamics—run only for polynomial time, far before mixing—achieves the same recovery performance as AMP. The phase transition is sharp: OPT_{β,λ} is zero for β<1/λ and strictly positive for β>1/λ, recovering the BBP threshold λ=1 when the posterior temperature β=λ is used. The","pith_inferences":["If the MLSI condition is eventually established for all β<1, the Glauber statement becomes unconditional and the technique would likely transfer to a broader class of spiked matrix priors, since the RGD reduction does not use the Boolean prior in an essential way.","The one-dimensional recursion suggests a testable prediction: annealed Glauber dynamics should trace the branch of stable fixed points OPT_{β,λ} as β increases past 1/λ, yielding the same threshold with a slow-heating schedule—going beyond the fixed-temperature result in the paper.","Because OPT_{β,λ} can increase beyond β=λ, samples from the posterior at β>λ may have larger correlation with the spike than the posterior mean, inverting the usual Bayes-optimality intuition; the authors flag this as a subtle consequence.","A direct numerical test of the recursion is feasible: simulate RGD at small N, measure the empirical correlation distribution at T=log²N steps, and compare to the numerically computed fixed point; a mismatch would pinpoint exactly where the one-dimensional reduction breaks."],"forward_implications":["If the MLSI condition holds for all β<1, Glauber dynamics from any initialization has correlation within O(N^{-1/2+ε}) of OPT_{β,λ} after exp(O(N^{4-2ε})) steps.","The threshold for nontrivial recovery by Glauber is β=1/λ, matching the BBP phase transition λ=1 in the Bayes-optimal case β=λ.","RGD from a warm start reaches the stable fixed point in O(log N) steps, implying fast convergence of the correlation process before mixing.","At the posterior temperature β=λ, the RGD fixed point exactly equals the AMP state-evolution fixed point, so a sampling-based estimator can reach AMP-level correlation without waiting for full mixing.","The paper's high-precision magnetization estimates pin the mean magnetization of the SK model under weak field to order N^{-1/2+ε} for all β<1, a spin-glass result with independent uses."],"fun_headline_variants":["Glauber dynamics equals AMP for spiked Wigner recovery","MCMC matches message passing at BBP threshold","Polynomial-time MCMC recovers spike as well as AMP","Glauber dynamics reaches AMP's fixed point in correlation space","Spiked Wigner: Glauber dynamics matches AMP phase transition"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The Glauber-dynamics theorem assumes the Sherrington–Kirkpatrick model at β<1 with any uniform external field satisfies a modified log-Sobolev inequality with constant Ω(1/N), a condition currently proven only for β<0.295 and merely conjectured for all β<1.","fun_headline_variants_meta":{"raw":{"variants":["Glauber dynamics equals AMP for spiked Wigner recovery","MCMC matches message passing at BBP threshold","Polynomial-time MCMC recovers spike as well as AMP","Glauber dynamics reaches AMP's fixed point in correlation space","Spiked Wigner: Glauber dynamics matches AMP phase transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000366,"raw_usage":{"total_tokens":1808,"prompt_tokens":750,"completion_tokens":1058,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":988}},"tokens_in":494,"tokens_out":1058,"duration_ms":12299,"temperature":1.0,"reasoning_tokens":988,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T19:00:21.399143+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a fixed λ>1 and β just above 1/λ, numerically compute OPT_{β,λ} from the fixed-point equations, then run RGD from a warm start for T=log²N steps on a planted spiked Wigner instance and measure the mean absolute correlation; if it fails to approach OPT_{β,λ} within N^{-1/2+ε} as N grows, the one-dimensional reduction is wrong.","supporting_citations":[],"review_version":1}