Pith. sign in

REVIEW 3 major objections 5 minor 17 references

Numerical simulation of fractional topological charge in $SU(N)$ gauge theory coupled with $\mathbb{Z}_N$ 2-form gauge fields

T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Numerical simulations confirm that coupling SU(2) gauge theory to a Z2 2-form gauge field makes the topological charge fractional and suppresses topological freezing.

desk verdict A short proceedings talk that visually confirms half-integer topological charge in SU(2) with a dynamical Z2 2-form gauge field; the result is plausible and reproducible, but the dynamical-B update algorithm that carries the argument is only referenced, not described. read the letter →

arxiv 2501.11438 v1 pith:QOQLVTXK submitted 2025-01-20 hep-lat

Motokazu Abe , Okuto Morikawa This is my paper · ORCID
classification hep-lat PACS 11.15.Ha
keywords fractionaltopologicalchargeZ_N2-formgaugefield1-formsymmetryhybridMonteCarlogradientflowfreezing'tHoofttwistedboundaryconditionSU(2)theory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish, by direct numerical simulation, that fractional topological charge is not only a formal consequence of gauging the $\mathbb{Z}_N$ 1-form symmetry but a real property of the lattice ensembles. The authors simulate SU(2) gauge theory with a dynamical $\mathbb{Z}_2$ 2-form gauge field multiplying each plaquette, smooth the configurations by gradient flow, and measure the topological charge with an improved operator. With the 2-form field present, the charge clusters at half-integers—consistent with $Q = -\frac{1}{N}\int \frac{1}{2} P_2(B_p) \mod 1$—while without it the charge remains integer-valued. The same ensembles show markedly shorter autocorrelation times, suggesting the 2-form field, which acts as a 't Hooft twisted boundary condition, may alleviate the topological freezing that hinders sampling at fine lattice spacings.

What carries the argument

The argument rests on the lattice action (2), in which each plaquette is multiplied by $\exp(-2\pi i B_p / N)$ with $B_p$ a $\mathbb{Z}_N$ 2-form field, and on identity (8), which fixes the fractional part of the topological charge as $Q = -\frac{1}{N} \int \frac{1}{2} P_2(B_p) \mod 1$, with the Pontryagin square $P_2(B_p) = B_p \cup B_p + B_p \cup_1 dB_p$. The charge is measured with a clover-plus-rectangle improved operator ($c_0 = 5/3$, $c_1 = -1/12$) after gradient flow at $t = (0.7L)^2/8$. To simulate the dynamical 2-form field while preserving detailed balance, the authors adopt a 'halfway-updating' HMC algorithm from Ref. [11].

What would settle it

Generate configurations with an independent update scheme—for example, a Metropolis accept/reject step for the $B$-field using the full action—at the same $\beta$ and $L$, and compare the $Q$ distribution and autocorrelation function with the halfway-updating HMC results; a mismatch would indicate the algorithm does not sample the theory intended.

Watch

Extended reading notes

Core claim

The paper claims that fractional topological charge appears in numerical ensembles, not just in formal arguments: after gradient flow, the improved topological charge in $SU(2)$ gauge theory coupled to a $\mathbb{Z}_2$ 2-form field clusters around half-integers, matching the prediction $Q = -\frac{1}{N} \int \frac{1}{2} P_2(B_p) \mod 1$ (Eq. 8). The authors also claim that including the 2-form field dramatically shortens the autocorrelation time of $Q$ relative to the periodic-boundary case, which they interpret as evidence that the 't Hooft twisted boundary condition, realized dynamically through $B$, mitigates topological freezing.

Load-bearing premise

The claim depends on the modified hybrid Monte Carlo update for the 2-form field actually sampling the intended action; the paper gives no description or check of that update, only a citation to another preprint.

Editorial extensions

If this is right

  • The numerical confirmation of Eq. (8) means the mixed 't Hooft anomaly for the 1-form symmetry is realized non-perturbatively in a concrete ensemble.
  • Dynamical 2-form fields are equivalent to summing over 't Hooft twisted boundary conditions, so the observed drop in autocorrelation time supports twisted boundary conditions as a practical cure for topological freezing.
  • At the chosen flow time $t=(0.7L)^2/8$, the improved operator yields stable values converging to the nearest half-integer, providing a measurement recipe for fractional $Q$.
  • With additional data at different lattice spacings, the integrated autocorrelation time can be extracted and compared with the open-boundary-condition benchmark of $\sim 1/a^2$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct test of detailed balance for the halfway-updating HMC—comparing, for instance, against a plain Metropolis update of the $B$-field—would confirm whether the observed $Q$ distribution really corresponds to the intended action; the paper does not provide such a check.
  • If the autocorrelation improvement persists at finer lattice spacings, the integrated autocorrelation time may scale as mildly as $1/a^2$, matching the open-boundary-condition result; this is a quantitative prediction worth testing at additional $\beta$ values.
  • The same mechanism should produce fractional $Q$ clustering at multiples of $1/N$ in $SU(N)$ for any $N$, and could be applied to theories with partially broken 1-form symmetries, where the fractional part of $Q$ would diagnose the pattern of symmetry breaking.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This proceedings contribution reports an exploratory lattice simulation of SU(2) gauge theory coupled to a dynamical Z_2 2-form gauge field, i.e., a lattice realization of gauged Z_2 1-form symmetry and of 't Hooft twisted boundary conditions. The authors generate configurations with a Wilson action plus a B-field plaquette coupling (Eq. (2)), using HMC for the link variables and a cited 'halfway-updating' HMC for the dynamical B field. After Wilson flow, they compute an improved topological-charge operator (Eqs. (5)-(7)). Without the B field, Q clusters near integers; with the B field, Q clusters near half-integers, consistent with Eq. (8). They also compare autocorrelation functions of Q with and without the B field and claim that topological freezing is mitigated.

Significance. If the sampling is correct, the paper provides the first numerical observation of fractional topological charge in a non-Abelian lattice gauge theory, and a potentially useful new route to avoiding topological freezing. The numerical check is largely parameter-free: no Q-dependent quantity is fitted, the improved topological-charge operator is defined independently of Eq. (8), and the only genuinely tunable smoothing scale is the flow time t_f. The public Julia code and the direct comparison with the formal prediction are strengths. However, the reliability of the reported distributions hinges on the B-field update, which is neither described nor verified, and the statistical evidence is currently qualitative rather than quantitative.

major comments (3)
  1. [§2.2, halfway-updating HMC] The central load-bearing assumption is the sampling algorithm for the dynamical B field. The paper states that treating B as dynamical 'naively breaks the detailed balance of the HMC method' and that a modified 'halfway-updating' HMC from Ref. [11] is used to restore it, but it gives no description, proof, or numerical verification of this algorithm. Since Figs. 2(b), 3, and 4(b) are only meaningful if the configurations sample the measure of Eq. (2), this is not a peripheral technicality. I ask for a self-contained description of the halfway-updating HMC together with a concrete check of detailed balance (for example, acceptance rates, comparison with a local Metropolis or heat-bath update of B on a small volume, or an explicit demonstration in the cited preprint).
  2. [§3.2, Fig. 2(b) and Fig. 3] The evidence for fractional topological charge is currently presented as a scatter plot of Q versus trajectory number and a single flow-time history. There are no error bars, no histogram or distribution test, and no ensemble-averaged comparison with the prediction Q = -(1/N) ∫ (1/2) P_2(B_p) mod 1 of Eq. (8). The clustering near half-integers is suggestive, but it does not by itself rule out the possibility that the apparent fractional values are an artifact of the B-update or of the choice of flow time, and it does not quantify the statistical significance. Please provide a histogram of Q (or Q mod 1) with errors, an estimate of the central value at the chosen flow time, and an analysis over independent runs or bootstrap resampling.
  3. [§3.3, Fig. 4 and Eq. (12)] The claim that topological freezing is mitigated by the B fields is based only on plots of the normalized autocorrelation function ρ(τ). The paper does not report the integrated autocorrelation time τ_int defined in Eq. (12), nor statistical errors on τ_int, so the comparison between Fig. 4(a) and Fig. 4(b) remains qualitative. Since the freezing problem is specifically about the growth of τ_int with decreasing lattice spacing, please compute τ_int with a stated summation window and error estimate for each β value, with and without the B field.
minor comments (5)
  1. [Eq. (6)] The displayed definitions of C^P_{μν} and C^R_{μν} are garbled in the manuscript; the reader cannot see the actual expressions inside the Im{...} braces. Please replace them with explicit, correctly typeset formulas for the clover and rectangular field-strength clovers.
  2. [§2.2, Eq. (4)] When the gradient flow is applied in the presence of the B field, please state explicitly whether B_p is held fixed during the flow and how the derivative ∂_{n,μ} S_W[V_t, B] is defined for the coupled action; the current sentence 'we substitute S_W[V_t,B] instead of S_W[V_t]' leaves this ambiguous.
  3. [Table 1 and Fig. 4] The 'number of configurations' in Table 1 should specify whether these are independent configurations or HMC trajectories (Fig. 2 says configurations are separated by 10 trajectories), how many trajectories were discarded for thermalization, and the total HMC length. In addition, Fig. 4 varies both β and L simultaneously, so the observed autocorrelation differences are a mixture of volume and lattice-spacing effects.
  4. [Fig. 4 caption] The phrase 'From the top of the label' in the caption of Fig. 4 is unclear; a conventional legend or an explicit table entry for (L, a√σ) would be easier to read.
  5. [§2.2, text near Eq. (8)] The sentence 'we can ensure that the construction in this talk is consistent with Ref. [6, 12]' should read 'Refs. [6, 12]'; more importantly, please spell out one sentence on how the equivalence to 't Hooft twisted boundary conditions is obtained, since the current statement is only asserted.

Circularity Check

1 steps flagged · score 4.0 of 10

The fractional-charge measurement itself is an independent check, but the dynamical Z_N 2-form sampling rests entirely on a same-author 'halfway-updating' HMC algorithm that is neither described nor verified.

  1. self citation load bearing [Section 2.2, paragraph beginning 'Next, we generate gauge configurations using an HMC method with the action coupled with B-fields']
    "Here, we treat the B field as a dynamical flux rather than a background field, which naively breaks the detailed balance of the HMC method. To restore the detailed balance, we employ a modified HMC algorithm, the “halfway-updating” HMC given in Ref. [11]."

    The numerical distributions of Q (Fig. 2(b)) and the autocorrelation comparison (Fig. 4(b)) are only meaningful if the update of the B field samples the action of Eq. (2) with correct detailed balance. The paper's only justification for that sampling is a citation to Ref. [11], a same-author preprint; no description, proof, pseudocode, or numerical test of the halfway-updating HMC is given here. Thus the central numerical claim reduces to the credibility of the self-cited algorithm rather than to an independently verifiable argument in this paper. The Q measurement itself is not circular — it uses a standard improved operator and is compared with, not fitted to, Eq. (8) — but the ensemble on which it is applied is asserted by self-citation.

full rationale

The derivation chain for fractional Q is: define action (2) with Z_N 2-form field; generate configurations by HMC; smooth by gradient flow; measure Q with clover-plus-rectangular improved operator. The measured Q values cluster near 0.5 and other half-integers for N=2, matching the formula Q = -1/N ∫ 1/2 P2(B_p) mod 1 from Eq. (8). This is a genuine consistency check rather than a fitted result: the improved operator's coefficients (c0=5/3, c1=-1/12) come from standard Lüscher-Weisz improvement, the flow time is a fixed choice, and no parameter is tuned to force fractional Q. The main circularity concern is the 'halfway-updating' HMC algorithm, which is the load-bearing mechanism for sampling the dynamical B field and is cited only to a same-author preprint (Ref. [11]) without an explanation or check in this paper. If that algorithm does not satisfy detailed balance for action (2), the Q distribution and autocorrelation improvement are artifacts. Since the fractional-charge measurement has independent content but rests on this unverified self-citation, the circularity score is 4.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central numerical claim rests on imported formal results (Eq. (8), twisted-boundary equivalence) and on an unverified algorithmic assumption (halfway-updating HMC). No new free parameters were fitted to make Q fractional; the only hand-chosen scale is the flow time. No new particles or fields are invented.

free parameters (1)
  • Gradient-flow time t_f = (0.7L)^2/8 (e.g., 7.84 for L=8)
    Selected by hand as the time before oversmearing; used for all Q measurements. Fig. 3 shows a plateau near 0.5, so the choice is not tightly constrained, but it is still a chosen scale.
assumptions (5)
  • domain assumption The lattice construction of Refs. [3,4] and Eq. (8) correctly give the fractional topological charge for SU(N) theory coupled to Z_N 2-form fields.
    The paper uses Eq. (8) as the expected value to compare against; it is imported from prior same-author work and not re-derived here.
  • ad hoc to paper The halfway-updating HMC of Ref. [11] restores detailed balance when B is dynamical.
    This is the key algorithmic assumption enabling the simulations; it is only cited, not described or verified.
  • domain assumption Gradient flow at the chosen flow time removes UV fluctuations without changing the topological-charge sector.
    Standard Wilson flow, but the specific time is hand-picked and no systematic check is shown for all configuration samples.
  • domain assumption Coupling B is equivalent to imposing 't Hooft twisted boundary conditions, per Refs. [6,12].
    This equivalence is used to interpret the autocorrelation improvement as a twisted-boundary effect.
  • standard math The improved topological charge operator with c0=5/3, c1=-1/12, from Ref. [10], approximates the continuum topological charge after flow.
    Coefficients are taken from the literature; the paper provides no new derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Numerical simulation of fractional topological charge in $SU(N)$ gauge theory coupled with $\mathbb{Z}_N$ 2-form gauge fields." pith.science (2026). https://pith.science/paper/QOQLVTXK

@misc{pith2026250111438,
  author       = {Pith},
  title        = {Pith review of: Numerical simulation of fractional topological charge in $SU(N)$ gauge theory coupled with $\mathbbZ_N$ 2-form gauge fields},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QOQLVTXK}},
  note         = {Machine review of arXiv:2501.11438}
}
abstract

The pure $SU(N)$ gauge theory with a $\theta$ term has the $\mathbb{Z}_N$ $1$-form global symmetry. When this symmetry is gauged, it is formally established that the topological charge becomes fractional. In this talk, we generate gauge configurations using the HMC method with coupling to the gauged $\mathbb{Z}_N$ $2$-form gauge field. After smoothing these configurations via the gradient flow method, we numerically confirm that the topological charge has a fractional value. We also anticipate that these higher-form fields can solve the topological freezing problem.

Figures

Figures reproduced from arXiv: 2501.11438 by the authors.

Figure 1
Figure 1. Z𝑁 1-form symmetry transformation on the lattice. The left panel 1(a) illustrates the global symmetry, while the right panel 1(b) depicts the gauged symmetry. The red line represents the symmetry operator 𝑈(Σ), the purple line indicates the Wilson loop 𝑊 (𝐶), and the orange-shaded region highlights the plaquette coupled with the 𝐵-fields. 2.2 Topological charge At first, we generate gauge configurations using the HM… view at source ↗
Figure 2
Figure 2. Topological charge 𝑄 vs. the trajectory number 𝜏 (molecular dynamics time). Each configuration is separated by 10 trajectories. The left panel shows that 𝑄 is almost integral, while the right panel indicates that 𝑄 is fractional. 𝛽 = 2.4 and 𝐿 = 8. how 𝑄 evolves with the flow time for a configuration. In this example, 𝑄 approaches the fractional value 0.5 because we set 𝑁 = 2 in Eq. (8). We also adapted some definit… view at source ↗
Figure 3
Figure 3. Topological charge 𝑄 vs. the flow time 𝑡. The configuration is selected from the region where thermalization has occurred and sufficient time has passed, specifically at 𝜏 = 1000. From this, it becomes evident that the rectangular improved method we chose exhibits the best convergence. Here, the purple line represents a specific flow time, 𝑡 = (0.7𝐿) 2 /8, chosen as the stage before oversmearing in the flow equation… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Normalized autocorrelation function 𝜌(𝜏) of the topological charge 𝑄 without ( [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 11 canonical work pages

  1. [11]

    M. Abe, O. Morikawa and H. Suzuki,Monte Carlo simulation of the𝑆𝑈(2)/Z2 Yang–Mills theory, 2501.00286

    arXiv

  2. [1]

    Gaiotto, A

    D. Gaiotto, A. Kapustin, N. Seiberg and B. Willett,Generalized Global Symmetries,JHEP 02 (2015) 172. 7 Numericalsimulationoffractionaltopologicalchargein 𝑆𝑈(𝑁)gaugetheorycoupledwith Z𝑁 2-formgauge fields Motokazu Abe

    work page 2015

  3. [2]

    Gaiotto, A

    D. Gaiotto, A. Kapustin, Z. Komargodski and N. Seiberg,Theta, Time Reversal, and Temperature,JHEP 05 (2017) 091

    work page 2017

  4. [3]

    M. Abe, O. Morikawa and H. Suzuki,Fractional topological charge in lattice Abelian gauge theory, PTEP 2023(2023) 023B03 [2210.12967]

    work page Pith review arXiv 2023

  5. [4]

    M. Abe, O. Morikawa, S. Onoda, H. Suzuki and Y. Tanizaki,Topology of SU(N) lattice gauge theories coupled withZ𝑁 2-form gauge fields,JHEP 08 (2023) 118 [2303.10977]

    arXiv 2023

  6. [5]

    ’t Hooft,Confinement and Topology in Nonabelian Gauge Theories,Acta Phys

    G. ’t Hooft,Confinement and Topology in Nonabelian Gauge Theories,Acta Phys. Austriaca Suppl. 22(1980) 531

    work page 1980

  7. [6]

    M. Abe, O. Morikawa and S. Onoda,Note on lattice description of generalized symmetries in 𝑆𝑈(𝑁)/Z𝑁 gauge theories,Phys. Rev. D108 (2023) 014506 [2304.11813]

    work page Pith review arXiv 2023

  8. [7]

    Lüscher,Topology of Lattice Gauge Fields, Commun

    M. Lüscher,Topology of Lattice Gauge Fields, Commun. Math. Phys.85(1982) 39

    work page 1982

Show all 17 references
  1. [8]

    Lüscher,Properties and uses of the Wilson flow in lattice QCD, JHEP 08(2010) 071 [1006.4518]

    M. Lüscher,Properties and uses of the Wilson flow in lattice QCD, JHEP 08(2010) 071 [1006.4518]

    2010 arXiv

  2. [9]

    CP-PACScollaboration, Topological susceptibility in lattice QCD with two flavors of dynamical quarks, Phys. Rev. D64 (2001) 114501 [hep-lat/0106010]

    2001 arXiv

  3. [10]

    Lüscherand P

    M. Lüscherand P. Weisz,On-shell improvedlattice gauge theories,Commun. Math.Phys. 98 (1985) 433

    1985

  4. [12]

    van Baal,Some Results for SU(N) Gauge Fields on the Hypertorus,Commun

    P. van Baal,Some Results for SU(N) Gauge Fields on the Hypertorus,Commun. Math. Phys. 85 (1982) 529

    1982

  5. [13]

    Lüscher and S

    M. Lüscher and S. Schaefer,Lattice QCD without topology barriers,JHEP 07 (2011) 036 [1105.4749]

    2011 arXiv

  6. [14]

    Teper,Physics from the lattice: glueballs in qcd; topology; su(n) for all n, 1997

    M. Teper,Physics from the lattice: glueballs in qcd; topology; su(n) for all n, 1997

    1997

  7. [15]

    Nguyen, Y

    M. Nguyen, Y. Tanizaki and M. Ünsal,Study of gapped phases of 4d gauge theories using temporal gauging of theZ𝑁 1-form symmetry, JHEP 08(2023) 013 [2306.02485]

    2023 arXiv

  8. [16]

    Morikawa and H

    O. Morikawa and H. Suzuki,Direct Monte Carlo computation of the ’t Hooft partition function, 2501.07042

    arXiv

  9. [17]

    Nagai and A

    Y. Nagai and A. Tomiya,JuliaQCD: Portable lattice QCD package in Julia language, 2409.03030. 8

    arXiv

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.