{"id":"faeb2787-0bd3-49d7-89ce-20c38963100b","arxiv_id":"2505.08236","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"The purely relaxational Metropolis dynamics of the 3D inverted XY universality class has dynamic critical exponent z=2.59(3), determined from quench simulations of Z6 and Z8 gauge models.","lead":"The authors measure the critical slowing down of 3D lattice gauge models with Z6 and Z8 symmetry and obtain the dynamic critical exponent z=2.59(3) for the inverted XY universality class. This shows that these topological transitions relax more slowly than the ordinary 3D XY model even though both share the same static exponent nu about 0.6717.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (15) is a conjectured energy-density scaling form; the energy flow may measure a faster mode than the slowest topological mode, so z=2.59(3) is not yet pinned to the slowest critical dynamics.","rationale":"The paper is careful, and the consistency between Z6 and Z8, the stability of the fits, and the plausible physical picture are real evidence. However, the central claim is a universal dynamic exponent of the slowest critical modes, and the out-of-equilibrium FSS ansatz of Eq. (15) is neither derived nor independently verified for the IXY class in this work. The energy-density observable may not project onto the slowest modes, and the authors explicitly report a smaller effective exponent from energy-density autocorrelations. The Polyakov-loop check is the natural control, but its uncertainty (0.12) is too large to validate the headline precision of 0.03. This does not make the paper wrong; it makes the acceptance conditional on a direct slowest-mode measurement. I therefore recommend CONDITIONAL rather than ACCEPT, with the condition being a Polyakov-loop exponential-autocorrelation or out-of-equilibrium Polyakov-loop FSS analysis at larger sizes.","tokens_in":16799,"tokens_out":5213,"duration_ms":59930,"concrete_test":"Measure the exponential autocorrelation time tau_exp of the zero-momentum Polyakov loop P_x (Eq. 24) at Kc for both Z6 and Z8, at L = 32, 40, 48, 56, 64, using long equilibrium runs and fitting tau_exp = a L^z from the exponential tail of C_P(t), excluding L < 40. Alternatively, apply the same out-of-equilibrium FSS analysis to the Polyakov-loop susceptibility rather than the energy density. If the resulting z differs from 2.59 by more than 2 sigma, the energy-flow estimate in Eq. (19) is biased; if z = 2.59(5) is reproduced, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that z=2.59(3) is the universal slowest-mode dynamic exponent of the IXY universality class. The extraction relies on Eq. (15), Omega = L^{3-y_r} E_s ~ A(Theta, Upsilon), which the paper itself labels 'conjectured' and refers to Refs. [55,61] for numerical verification; no independent derivation is given for the IXY class. The same energy-density dataset is used both to fit z and to demonstrate collapse in Fig. 2, so the collapse is not an independent test of the ansatz. The risk is concrete and acknowledged in Sec. IV: the integrated autocorrelation time of the energy density gives z_i ~ 2.3, smaller than the claimed z, which the authors attribute to weak coupling of the energy density to the slowest modes. If that is the case, the out-of-equilibrium energy flow can saturate to an effective exponent below the true asymptotic z. The only independent slowest-mode probe, the Polyakov-loop integrated autocorrelation, yields z = 2.52(12), a much larger error than the headline and compatible with a wide range. Thus the slowest-mode exponent, which is the physically claimed quantity, is not directly constrained at the claimed precision.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the critical relaxational (Model A) dynamics of three-dimensional Z_N lattice gauge models with N=6 and N=8, whose topological transitions belong to the inverted XY (IXY) universality class. The authors perform large-scale Monte Carlo simulations of a local Metropolis dynamics and analyze out-of-equilibrium relaxational flows after instantaneous quenches to the critical point. Using an out-of-equilibrium finite-size scaling framework adapted from Refs. [55,61], they fit the L-dependence of the time at fixed rescaled energy density, and of time differences, to aL^z, obtaining z=2.595(35) for N=6 and z=2.590(35) for N=8. They combine these into a final estimate z=2.59(3). They also report an equilibrium study of integrated autocorrelation times of Polyakov loops, giving a consistent but less precise value z=2.52(12). The central claim is that z=2.59(3) characterizes the relaxational critical slowing down of the entire IXY universality class, including topological transitions in 3D Abelian Higgs models.","tokens_in":17042,"tokens_out":10188,"duration_ms":105392,"significance":"If the result holds, this is the first determination of the Model A dynamic exponent for the IXY universality class, filling a clear gap in the literature. The consistency between N=6 and N=8 provides strong evidence for universality, and the comparison with the standard XY value z≈2.02 highlights a striking difference in dynamics despite shared static exponents. The paper is careful in several ways: fits are tested for stability as L_min is increased, scaling-correction fits with L^{-ω} are included, and an independent equilibrium probe based on Polyakov loops is used as a cross-check. The prediction is falsifiable in other IXY-class models, such as lattice Abelian Higgs theories, making the paper a useful reference for future work.","major_comments":[{"comment":"The final estimate z=2.59(3) rests on the conjectured out-of-equilibrium scaling form (15) for the subtracted energy density. The collapse shown in Fig. 2 is not an independent test, since the same energy-density data determine z. The paper itself reports in Sec. IV that the integrated autocorrelation time of the energy density gives z_i≈2.3, smaller than the claimed z, which the authors attribute to weak coupling of the energy density to the slowest modes. If that is the case, the out-of-equilibrium energy flow could similarly measure an effective exponent rather than the slowest-mode exponent. The only independent slowest-mode probe, the Polyakov-loop integrated autocorrelation time, yields z=2.52(12), whose uncertainty is considerably larger than the claimed 0.03. The error quoted in Eq. (19) therefore does not include a possible systematic uncertainty from observable choice. I ask the authors to either provide an out-of-equilibrium analysis of a nonlocal topological observable (e.g., Polyakov-loop related quantities) or to widen the final error to a conservative value and explicitly state this limitation.","section":"Sec. III C and Eq. (15); Sec. IV and Eq. (19)"}],"minor_comments":[{"comment":"In the sentence preceding Eq. (19), 'we consider' should be 'we obtain', and the capitalization 'We' after the comma should be lower-case.","section":"Sec. III D"},{"comment":"The acceptance ratio of about 2% at the critical point is quite low; it would be helpful to state the corresponding autocorrelation times or the number of Metropolis sweeps per trajectory to help the reader judge the computational efficiency of the proposed update.","section":"Sec. III A"},{"comment":"Reference [63] appears to have a formatting error in the author list: 'E. V. Ivanovae. M. V. Kompaniets' should have a space and period, e.g., 'E. V. Ivanova, E. M. V. Kompaniets' or similar.","section":"References"},{"comment":"The text says 'in Figs. 2 and we show plots' in the third paragraph of Sec. III D; the figure number for the second panel is missing and should be 'Figs. 2 and 3' or the sentence should be rewritten.","section":"Fig. 2 caption"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid numerical study from an experienced group. The requested revision is primarily about conservative treatment of systematic errors rather than any detected technical error in the analysis. I do not have concerns about scope or citation practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe short summary: this is a careful numerical paper that fills a genuine gap. For the purely relaxational (Model A) critical dynamics of the 3D inverted XY universality class, no dynamic exponent had been measured. The authors extract z=2.59(3) from out-of-equilibrium FSS of the energy density after quenches to Kc, using two different models (Z6 and Z8 gauge), and the result is stable across fits of time differences with Lmin≥24. The equilibrium autocorrelation of the Polyakov loop gives z=2.52(12), consistent but less precise. That consistency is the key supporting evidence, because the energy density alone is not a clean probe of the slowest modes.\n\nThe paper does several things well. It uses the difference observable Δ(Ω, Υ1, Υ2, L) to suppress scaling corrections, and the fits look honest with acceptable χ². The data collapse in Fig. 2 for both models is reassuring, though it uses the fitted z so it is not an independent test. The authors also openly state that the integrated autocorrelation of the energy density gives zi≈2.3, smaller than their z, and explain why this does not invalidate the estimate; they then bring in the Polyakov-loop result as the slowest-mode check. That is exactly the right way to handle the issue.\n\nThe soft spot is the out-of-equilibrium FSS ansatz of Eq. (15). It is labelled conjectured and relies on the authors' own earlier numerical work (Refs. [55,61]), not on a derivation or an independent confirmation. For this paper, the distribution of Θ and the collapse show the ansatz works self-consistently, but the framework is not independently verified for the IXY class. The stress-test concern that the energy flow might track a faster mode than the slowest topological mode is real in principle; the paper's own Sec. IV shows it should be taken seriously. However, the Polyakov-loop autocorrelation data, though limited to L≤44 for Z6, do point to the same z. So I don't think the caveat sinks the main claim; it mainly sets the precision and the reliance on the conjectured ansatz.\n\nMinor issues: no data/code repository, and the acceptance rate of the Metropolis updates is very low (~2%) but that's a detail. The citation pattern is fine; the self-citations are to the framework they previously built, and the Z2 gauge model results are relevant.\n\nWho is this for: anyone working on critical dynamics of lattice gauge theories, topological deconfinement transitions, or dynamic universality classes. It's a solid numerical contribution, not a breakthrough, but it deserves a serious referee. I would recommend acceptance with the usual request for a bit more discussion of the ansatz's status and perhaps a clearer statement of the Polyakov-loop check's limitations.\n\nRecommendation: send it to peer review; I'd accept it if the above caveats are addressed in the revision.","headline":"Solid first estimate of the dynamic exponent for the 3D IXY universality class, with an honest treatment of the main caveat; the conjectured out-of-equilibrium scaling ansatz is the only real soft spot.","tokens_in":17610,"tokens_out":2857,"would_cite":true,"duration_ms":26906,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.70.Jk","64.60.Ht","11.15.Ha"],"model":"deepseek-v4-flash","headline":"The paper establishes that the purely relaxational critical dynamics of the three-dimensional inverted XY universality class has dynamic exponent z = 2.59(3), obtained from out-of-equilibrium finite-size scaling of Z6 and Z8 gauge models.","keywords":["dynamic critical exponent","critical slowing down","inverted XY universality class","topological phase transition","lattice gauge theory","Z_N gauge model","out-of-equilibrium finite-size scaling","Metropolis dynamics"],"falsifier":"Run the same local Metropolis quench protocol on the 3D IXY gauge model (or on a 3D Abelian Higgs model) at its critical point and extract $z$ from the out-of-equilibrium scaling of $t(\\Omega, \\Upsilon, L)$; if the result lies outside $2.59(3)$ — for instance near $2.3$, as the energy-density integrated autocorrelation naively suggests — the claimed universality of the dynamic exponent across the whole IXY class would be disproved.","tokens_in":16574,"feed_emoji":"⏳","tokens_out":8903,"duration_ms":78845,"temperature":0.7,"pith_summary":"This paper determines the dynamic critical exponent that controls critical slowing down under local, purely relaxational Metropolis dynamics at topological transitions in the three-dimensional inverted XY (IXY) universality class. Simulating Z6 and Z8 lattice gauge models quenched instantaneously to their critical points, it analyzes the relaxational flow of the gauge-invariant energy density within an out-of-equilibrium finite-size scaling framework and obtains $z = 2.59(3)$. This matters because the IXY class includes the continuous transitions of 3D lattice Abelian Higgs models, effective theories of superconductors, and because it shows that gauge-theory topological transitions slow down far more severely than the standard 3D XY model, whose relaxational exponent is $z \\approx 2.02$ despite sharing the same correlation-length exponent $\\nu \\approx 0.6717$.","feed_headline":"Critical slowing down in 3D gauge models: z = 2.59","feed_subtitle":"Local relaxational dynamics in the inverted-XY class is slower than standard XY, despite sharing its static exponents.","key_machinery":"The central object is the subtracted post-quench energy density $E_s(t,r,L) = E(t,r,L) - E_{c,\\infty}$ and its rescaled form $\\Omega = L^{3-y_r} E_s$, with $y_r = 1/\\nu$. The argument runs on the conjectured out-of-equilibrium finite-size scaling ansatz $\\Omega(t,r,L) \\approx A(t L^{-z}, r L^{y_r})$, which removes the analytic background that dominates the equilibrium energy density and lets the time variable $\\Theta = t L^{-z}$ expose the dynamic exponent. In practice the paper fits $t(\\Omega, \\Upsilon, L) \\approx L^z F(\\Omega, \\Upsilon)$ and, more stably, the difference $\\Delta(\\Omega, \\Upsilon_1, \\Upsilon_2, L) \\approx L^z F_\\Delta(\\Omega, \\Upsilon_1, \\Upsilon_2)$ at fixed $\\Omega$, while the Polyakov-loop integrated autocorrelation time supplies the equilibrium cross-check.","core_discovery":"The central claim is that the dynamic exponent $z$ of the IXY universality class under Model A (locally reversible Metropolis) dynamics is $z = 2.59(3)$, a value the paper treats as universal for all topological transitions in this class. The evidence comes from consistent out-of-equilibrium finite-size scaling analyses of the Z6 and Z8 gauge models, yielding $z = 2.595(35)$ and $z = 2.590(35)$ respectively, with equilibrium autocorrelation data of the Polyakov loop giving the less precise but consistent $z = 2.52(12)$. Because the duality relating the IXY free energy to the XY model is nonlocal, local dynamics in the gauge model maps to nonlocal dynamics in the spin model, so the dynamic universality class differs even though the static thermal sector is shared.","pith_inferences":["If $z = 2.59(3)$ carries over to Abelian Higgs models, numerical studies of superconducting transitions will need roughly $L^{0.57}$ more sweeps than XY spin-model equivalents to decorrelate at criticality.","A direct quench experiment on the 3D IXY gauge model itself, rather than on the dual ZN clock or XY systems, would test the universality claim; agreement at the level of a few percent would close the loop on the class.","The observable dependence of integrated autocorrelation exponents ($z_i \\approx 2.3$ for energy density versus $z = 2.52(12)$ for Polyakov loops) is a caution that future dynamic studies in gauge theories should report both the observable and the estimator.","Applying the same out-of-equilibrium protocol to noncompact U(1) Higgs or other ZN models with $N > 8$ could show whether $z$ stays pinned near $2.59$ or drifts with $N$, sharpening the meaning of 'universal' for this class."],"forward_implications":["The value $z = 2.59(3)$ should govern critical slowing down for any local relaxational dynamics at a 3D IXY transition, including the lattice IXY gauge model and the 3D Abelian Higgs models.","Relaxational dynamics in the IXY class is significantly slower than in the standard 3D XY class: $z = 2.59(3)$ versus $z \\approx 2.02$, despite identical static exponents.","The nonlocal duality between the IXY and XY partition functions does not imply equal dynamic exponents, because it does not preserve locality of the dynamics.","The integrated autocorrelation time of the energy density grows with an effective exponent near $2.3$, smaller than the true $z$; the Polyakov loop, not the energy density, couples to the slowest critical mode.","Equilibrium autocorrelation analyses confirm but do not improve on the out-of-equilibrium estimate, so the quench-based method is the more powerful route to $z$ in gauge systems."],"supporting_citations":[{"why":"Review that frames the inverted XY universality class of topological deconfinement transitions, placing the ZN and Abelian-Higgs models in the class.","marker":"[8]"},{"why":"Duality relations showing the IXY free energy maps to the XY model with Villain action, giving the shared static critical exponents.","marker":"[13, 14]"},{"why":"Earlier out-of-equilibrium finite-size scaling study of the Z2 gauge model that establishes the quench protocol and z-extraction method now applied to Z6 and Z8.","marker":"[55]"},{"why":"Source of the conjectured out-of-equilibrium scaling ansatz for the subtracted energy density, Eq. (15), which underpins the z determination.","marker":"[61]"},{"why":"Provides the critical couplings Kc for the Z6 and Z8 models that define the quench targets.","marker":"[65]"},{"why":"Precise value of the 3D Ising dynamic exponent z = 2.0245(15) used as the spin-model baseline for comparison.","marker":"[59]"},{"why":"High-order perturbative calculation giving the 3D XY relaxational dynamic exponent z = 2.0246(10) that the paper compares against.","marker":"[62]"},{"why":"Sokal's remark that integrated autocorrelation exponents can be smaller than the true dynamic exponent, used to set the energy-density zi ≈ 2.3 aside in favor of Polyakov loops.","marker":"[73]"}],"fun_headline_variants":["Inverted XY dynamics: z=2.59, slower than standard XY","3D gauge models exhibit slower critical dynamics: z=2.59","Same static, slower dynamic: IXY z=2.59 vs XY z=2.02","Topological transitions: critical slowing down z=2.59","Gauge models: dynamic exponent 2.59, not 2.02, for IXY"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The determination of $z$ assumes that the out-of-equilibrium finite-size scaling form of Eq. (15) for the subtracted energy density is exact, a conjecture whose numerical support comes from earlier work; if that form is contaminated by analytic background or corrections stronger than $O(L^{-\\omega})$, the quoted exponent inherits a bias.","fun_headline_variants_meta":{"raw":{"variants":["Inverted XY dynamics: z=2.59, slower than standard XY","3D gauge models exhibit slower critical dynamics: z=2.59","Same static, slower dynamic: IXY z=2.59 vs XY z=2.02","Topological transitions: critical slowing down z=2.59","Gauge models: dynamic exponent 2.59, not 2.02, for IXY"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000788,"raw_usage":{"total_tokens":3509,"prompt_tokens":1015,"completion_tokens":2494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":2387}},"tokens_in":631,"tokens_out":2494,"duration_ms":18174,"temperature":1.0,"reasoning_tokens":2387,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:59:53.322819+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same local Metropolis quench protocol on the 3D IXY gauge model (or on a 3D Abelian Higgs model) at its critical point and extract $z$ from the out-of-equilibrium scaling of $t(\\Omega, \\Upsilon, L)$; if the result lies outside $2.59(3)$ — for instance near $2.3$, as the energy-density integrated autocorrelation naively suggests — the claimed universality of the dynamic exponent across the whole IXY class would be disproved.","supporting_citations":[{"cited_title":"Bonati, H","cited_arxiv_id":null,"evidence_quote":"Earlier out-of-equilibrium finite-size scaling study of the Z2 gauge model that establishes the quench protocol and z-extraction method now applied to Z6 and Z8."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"High-order perturbative calculation giving the 3D XY relaxational dynamic exponent z = 2.0246(10) that the paper compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sokal's remark that integrated autocorrelation exponents can be smaller than the true dynamic exponent, used to set the energy-density zi ≈ 2.3 aside in favor of Polyakov loops."}],"review_version":1}