Pith. sign in

REVIEW 2 major objections 6 minor 103 references

Topological properties around the Roberge-Weiss transition in $N_f = 2 + 1 + 1$ QCD

T0 review · 2 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read Along the Roberge-Weiss line, the quartic cumulant of QCD's topological charge becomes compatible with the dilute instanton gas value right at the transition, at two lattice spacings.

desk verdict Genuinely new and carefully done measurement of topology along the Roberge-Weiss line; the sharp DIGA onset above T_RW is plausible but rests on two aspect-ratio-2 lattices, so the missing volume scan of b2 is the main thing to push on. read the letter →

arxiv 2608.12883 v1 pith:7E44L3MM submitted 2026-08-13 hep-lat hep-th

classification hep-lathep-th PACS 12.38.Gc11.15.Ha
keywords Roberge-Weisstransitiontopologicalsusceptibilityfourth-ordercumulantb2diluteinstantongaslatticeQCDstaggeredfermionschargeNf=2+1+1
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

The paper argues that in $N_f=2+1+1$ QCD with physical quark masses, the temperature dependence of QCD's $\theta$-angle response changes character sharply at the Roberge-Weiss transition: the fourth-order cumulant $b_2$ of the topological charge distribution becomes compatible with the dilute instanton gas (DIGA) value $-1/12$ as soon as $T \ge T_{\mathrm{RW}}$, rather than approaching it gradually as along the standard thermal crossover. This matters because it would mean that an exact phase transition, not a crossover, is what triggers the onset of cosine-shaped $\theta$-dependence, and that the semiclassical instanton-gas description applies immediately above $T_{\mathrm{RW}}$ in a theory with dynamical quarks. The paper also establishes that the topological susceptibility is roughly constant below $T_{\mathrm{RW}}$ and decays as a power law above it, with exponents less steep than the semiclassical prediction, as expected from remaining lattice artifacts. A sympathetic reader would care because the Roberge-Weiss line then provides a sharply defined laboratory for testing instanton-gas ideas and for relating center symmetry to topological fluctuations.

What carries the argument

The central object is the fourth-order cumulant of the topological charge distribution, $b_2 = -\frac{1}{12}\frac{\langle Q^4\rangle - 3\langle Q^2\rangle^2}{\langle Q^2\rangle}$, which fixes the first deviation of the free-energy density $f(T,\theta)$ from a Gaussian in $\theta$ through $f(T,\theta)-f(T,0)=\frac{1}{2}\chi(T)\theta^2(1+b_2\theta^2+\cdots)$. The dilute instanton gas approximation predicts $f(T,\theta)=\chi(T)(1-\cos\theta)$, which fixes $b_2=-1/12$ and is sensitive only to the non-interacting, integer-charge nature of the topological objects, not to the semiclassical power law for $\chi$. The paper measures $b_2$ from gluonic determinations of the winding number $Q$ on cooled configurations, using a rounding prescription for $Q$ and a linear extrapolation in the number of cooling steps; the simulations use stout-improved staggered fermions with physical quark masses along the Roberge-Weiss line, where fermionic boundary conditions are twisted by the phase $\theta_q=\pi$ and a $\mathbb{Z}_2$ remnant of center symmetry breaks spontaneously.

What would settle it

Measure $b_2$ at fixed temperature just above $T_{\mathrm{RW}}$ on the same $N_t=8$ lattice with three or four spatial volumes; if $b_2$ moves appreciably away from $-1/12$ as the volume grows, or if $\langle Q^2\rangle$ falls below roughly one on any of those volumes at $T$ just above $T_{\mathrm{RW}}$, the claimed sharp DIGA onset would be a finite-volume artifact rather than a physical transition.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that along the Roberge-Weiss line the quartic cumulant $b_2 = -\frac{1}{12}\frac{\langle Q^4\rangle - 3\langle Q^2\rangle^2}{\langle Q^2\rangle}$ is compatible with the DIGA prediction $b_2 = -1/12$ for every temperature at or above $T_{\mathrm{RW}}$ studied on the $16^3 \times 8$ and $20^3 \times 10$ lattices, while below $T_{\mathrm{RW}}$ it remains far from both the chiral perturbation theory and DIGA expectations and suffers sizeable finite-size and taste-breaking effects. The topological susceptibility is approximately constant below $T_{\mathrm{RW}}$ and then drops, with a power-law exponent $c \simeq -4.7$ ($N_t=8$) or $-4.4$ ($N_t=10$) in the fixed-$N_t$ analysis and roughly $-3.9$ in the fixed-scale analysis, softer than the DIGA prediction because of residual discretization effects. This sharp switch of $\theta$-dependence is presented as closer to the behavior of pure gauge theories across deconfinement than to the milder approach observed along the standard thermal line in full QCD, and is interpreted as linked to the spontaneous breaking of the residual $\mathbb{Z}_2$ center symmetry at $T_{\mathrm{RW}}$.

Load-bearing premise

The sharp onset of the dilute-instanton form right above $T_{\mathrm{RW}}$ rests on the assumption that finite-volume suppression of nonzero topological charge is not already faking a DIGA-shaped distribution on the two aspect-ratio-2 lattices used for $b_2$; the paper checks that $\langle Q^2\rangle$ is not much smaller than one there, but it does not carry out a dedicated volume scan of $b_2$ across $T_{\mathrm{RW}}$.

Editorial extensions

If this is right

  • If the claim holds, the free-energy density along the Roberge-Weiss line has the cosine form $\chi(T)(1-\cos\theta)$ for all temperatures above $T_{\mathrm{RW}}$ explored, meaning non-interacting topological objects dominate immediately after the transition.
  • The agreement between $N_t=8$ and $N_t=10$ suggests the sharp onset is not an accidental lattice artifact and should survive the continuum limit.
  • The result brings full QCD with dynamical quarks into line with pure gauge theories, reinforcing the idea that spontaneous breaking of a remnant center symmetry, not merely deconfinement, controls the change in $\theta$-dependence.
  • The measured values $T_{\mathrm{RW}} = 199.0 \pm 1.4$ MeV ($N_t=8$) and $204 \pm 7$ MeV ($N_t=10$), with 3D-Ising scaling, confirm that adding a dynamical charm quark does not alter the Roberge-Weiss critical behavior at the physical point.

Reading between the lines

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

  • A dedicated volume scan of $b_2$ at fixed $T$ just above $T_{\mathrm{RW}}$—say aspect ratios 2, 3, and 4 on the same $N_t$—would cleanly separate a genuine DIGA onset from finite-volume suppression of $Q\neq 0$ sectors; the present data check $\langle Q^2\rangle$ but do not include such a scan.
  • If the sharp onset is really tied to the spontaneous breaking of the residual $\mathbb{Z}_2$ symmetry, then other theories with an exact center symmetry at imaginary chemical potential—such as many-flavor QCD near the conformal window, where the Roberge-Weiss transition is used as a probe—should show the same abrupt change in $\theta$-dependence.
  • The slower approach of $b_2$ along the standard thermal line might be at least partly a crossover-smearing effect; a higher-statistics measurement of the standard line, which the authors call for, could reveal that the asymptotic instanton-gas regime begins closer to $T_c$ than previously inferred.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 6 minor

Summary. The manuscript reports a lattice study of topological properties along the Roberge-Weiss (RW) line for N_f = 2 + 1 + 1 QCD at physical quark masses, using stout-improved staggered fermions on lattices with N_t = 8 and N_t = 10. The authors first determine T_RW and the universality class by finite-size scaling of the Polyakov loop and strange quark number observables, finding results consistent with the N_f = 2 + 1 determination of Ref. [21]. They then measure the topological susceptibility and the fourth-order cumulant b_2 using cooling and a rounding prescription. The main results are that chi^{1/4} is roughly constant below T_RW and then drops with an effective power-law exponent c ≈ -4.4 to -4.7, and that b_2 becomes compatible with the DIGA value -1/12 immediately above T_RW on both lattice spacings, in apparent contrast with the milder approach reported along the standard thermal line. The paper concludes that this sharp onset may reflect the spontaneous breaking of the residual Z_2 center symmetry at the RW transition.

Significance. The central observation, if confirmed, is physically interesting: it would indicate that the theta-dependence of the QCD free energy is controlled by the sharp center-symmetry-breaking transition on the RW line, with the dilute-instanton regime setting in as soon as T > T_RW. The paper has clear strengths: it uses high statistics (O(10^6) trajectories for the b_2 measurements), two lattice spacings for the main observable, an independent finite-size-scaling determination of T_RW that does not use topological quantities, and external physical benchmarks (ChiPT b_2 = -0.022, DIGA b_2 = -1/12) that are not fitted to the data. The discussion of zero-smoothing extrapolations and the use of the open-source OpenStaPLE code are also positive. The main caveat is that the sharp-onset claim rests on a single volume per lattice spacing and on a finite-volume diagnostic that is not sufficient to rule out volume-dependent distortions of b_2 just above T_RW; a dedicated volume scan is needed before the claim can be considered established.

major comments (2)
  1. [Sec. 4B, Figs. 11-12] The central claim that b_2 becomes compatible with the DIGA value immediately above T_RW is based on one spatial volume for each N_t (16^3 x 8 and 20^3 x 10), both with aspect ratio two. The only finite-volume control is the simultaneous measurement of <Q^2>, and the assertion that determinations right above T_RW are safe because <Q^2> is not much less than 1 is not a substitute for a volume scan. The mean <Q^2> alone does not identify the distribution: a Skellam distribution with the same mean gives exactly b_2 = -1/12, while a Q in {-1,0,1} truncation yields b_2 = (3<Q^2> - 1)/12, so the measured b_2 must ultimately be checked against the full P(Q) or against results at several spatial volumes. Because the same section reports sizeable finite-size and taste-breaking effects on b_2 below T_RW, the assumption that these effects disappear immediately above T_RW is load-bearing and untested. A dedicated scan, e.g., N_s = 24 and 32 at N_t = 8 for two temperatures just above T_RW, would make the sharp-onset claim robust.
  2. [Sec. 4B] The statement that the result 'seems stable towards the continuum limit' is supported by only two lattice spacings, with no continuum extrapolation. The two spacings also correspond to different physical volumes, and the paper shows that b_2 below T_RW has a visible N_t dependence; it is therefore possible that the same discretization dependence is present above T_RW but masked by statistical errors or by different finite-volume effects. A continuum extrapolation at fixed physical volume, or at least an explicit estimate of O(a^2) corrections, is needed before the sharp onset can be presented as a continuum property rather than as a property of the two simulated lattices.
minor comments (6)
  1. [Abstract; Sec. 3A] There are several typos that should be corrected: 'compactifed' in the abstract, 'spontanous' and 'pesudo-critical' in Sec. 3A, and 'anc' in Sec. 4A.
  2. [Sec. 4A, Fig. 9] The text says simulations were performed on '40^4 x N_t' lattices, while the caption and context indicate the intended notation is 40^3 x N_t; please clarify the volume nomenclature.
  3. [Sec. 4A] The 'rough estimate' of the power-law exponent c approx -3.9 from two fixed-a data points has no quoted error; please either provide an uncertainty or explicitly label it as a qualitative estimate with no statistical significance.
  4. [Sec. 4A, Fig. 9] The comparison with the T = 0 determination from Ref. [86] uses a slightly different discretization; the manuscript notes this, but it would be helpful to state explicitly that the systematic uncertainty from the discretization difference is not estimated.
  5. [Sec. 4B] The statement that finite-size effects become important above about 250 MeV is inferred visually from Figs. 11-12; a quantitative criterion for this threshold, or a reference to a dedicated test, would make the claim more transparent.
  6. [Sec. 4B] Below T_RW the b_2 data do not reach the ChiPT prediction, and the text attributes this to finite-size and taste-breaking effects without numerical support; a brief estimate of the expected size of these effects would improve the discussion.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: b2 is compared against external DIGA/ChiPT benchmarks, and T_RW is fixed by independent finite-size scaling of Polyakov-loop and density observables.

full rationale

The paper's central claim is that b2 becomes compatible with the DIGA prediction just above T_RW. This is a comparison of lattice-measured b2, defined in Eq. (12) from the measured topological charge distribution, against the external analytic prediction b2 = -1/12 that follows from the non-interacting cos(theta) form of Eq. (13). No parameter in the DIGA or ChiPT benchmarks is fitted to the lattice data, and b2 is not used to determine T_RW. The transition temperature is instead located independently via finite-size scaling of the imaginary Polyakov loop and the strange-quark number density susceptibilities, Eqs. (7)-(8), using externally tabulated 3D-Ising critical exponents. The same-group references, such as Ref. [21] for the N_f=2+1 T_RW value and Ref. [86] for the T=0 topological susceptibility baseline, provide context and cross-checks, but the paper independently reproduces T_RW and the b2 comparison does not reduce to those citations. The finite-volume discussion, including the <Q^2> check, is a systematic uncertainty rather than a circular step, since it neither defines the DIGA value nor tunes the comparison to force agreement. No equation or procedure in the paper exhibits a reduction of the claimed result to its own inputs, so the appropriate finding is no significant circularity.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

No invented entities are introduced. The central claim depends on standard lattice-QCD assumptions (rooting, line of constant physics, 3D-Ising exponents, DIGA/ChiPT benchmarks) and on the paper-specific assumption that finite-volume effects are negligible just above T_RW. Fitted parameters are the pseudo-critical couplings used to locate T_RW, the per-ensemble rounding factor alpha, and the chi power-law exponent; the latter is a fit, not an input to the b2 claim.

free parameters (3)
  • Pseudo-critical coupling beta_RW(N_t, N_s) from finite-size fits = T_RW = 199.0(1.4) MeV (N_t=8), 204(7) MeV (N_t=10)
    Sec. 3: the transition location used to separate below/above behavior is obtained from tanh and Lorentzian peak fits plus FSS extrapolation; it is an output of fits, but the b2 comparison depends on it.
  • Rounding factor alpha in topological charge assignment = per-ensemble, not tabulated, chosen to minimize cluster offset
    Sec. 2, Eq. (4): alpha calibrates Q to integers before computing cumulants; different choices could bias b2 and chi.
  • Power-law exponent c in chi^{1/4} = A T^{c/4} = c ~ -4.7 (N_t=8), ~ -4.4 (N_t=10), ~ -3.9 (fixed-a)
    Sec. 4A: fitted to the high-T drop of chi; used for comparison with DIGA, not as an input to the b2 claim.
assumptions (6)
  • domain assumption Rooted staggered fermion determinant correctly reproduces the light-quark and charm physics of QCD after rooting.
    Invoked in Sec. 2, Eq. (1); the paper notes charm effects are small and the continuum limit is only marginally approached at a ~ 0.1 fm.
  • domain assumption The line of constant physics from the scale-setting points of Ref. [30] applies to this stout-staggered 2+1+1 discretization and fixes physical quark masses.
    Sec. 2: bare masses and beta are tuned by spline interpolation of Ref. [30]; a wrong LCP would shift T_RW and all temperature values.
  • standard math 3D Ising critical exponents from Ref. [71] are the correct universality class for a second-order Z2 RW transition.
    Used in Sec. 3 FSS collapse plots; the paper uses published exponents rather than deriving them.
  • domain assumption DIGA prediction b2=-1/12 and ChiPT prediction b2=-0.022 are valid external benchmarks in the relevant temperature regimes.
    Sec. 4 uses these as comparison standards; DIGA is expected valid for T much larger than Lambda_QCD, and the paper tests it at T just above T_RW ~ 200 MeV.
  • domain assumption Cooling plus linear extrapolation in cooling steps removes UV renormalization effects and does not bias the integer topological charge assignment.
    Sec. 4, Eqs. (3)-(4) and Fig. 8; if the extrapolation is invalid, chi and b2 values carry an unquantified systematic.
  • ad hoc to paper Finite-volume effects on b2 are negligible above T_RW on 16^3 x 8 and 20^3 x 10 lattices.
    Sec. 4B: the paper argues safety because Q^2 is not small, but no dedicated volume scan of b2 is provided; this is the weakest assumption.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Topological properties around the Roberge-Weiss transition in $N_f = 2 + 1 + 1$ QCD." pith.science (2026). https://pith.science/paper/7E44L3MM

@misc{pith2026260812883,
  author       = {Pith},
  title        = {Pith review of: Topological properties around the Roberge-Weiss transition in $N_f = 2 + 1 + 1$ QCD},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7E44L3MM}},
  note         = {Machine review of arXiv:2608.12883}
}
abstract

We investigate the topological properties of QCD across the finite temperature Roberge-Weiss transition, which is found for particular values of the imaginary baryon chemical potential. Our study is conducted for $N_f = 2+1+1$ QCD with physical quark masses, discretized via stout improved staggered fermions and considering mostly two different values of the compactifed dimension, $N_t = 8$ and $N_t = 10$. Results for $T_{RW}$ and for the associated universality class are consistent with those found in the $N_f = 2 + 1$ case with a slightly different discretization. The topological susceptibility appears to be practically constant for $T \lesssim T_{RW}$, then rapidly decaying for higher temperatures. The analysis of the fourth order cumulant of the topological charge distribution, $b_2$, reveals that it is compatible with the prediction of the Dilute Instanton Gas Approximation right after $T_{RW}$, showing thus a sharp transition, which is more similar to what observed in pure gauge theories rather than to full QCD along the standard thermal line, where instead a slower transition was observed in previous studies.

Figures

Figures reproduced from arXiv: 2608.12883 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figures from the paper (7 more)
Figure 6
Figure 6. Figure 6: ), where systematics related to the determination of the physical scale have been added; taking into ac￾count the systematics from the different observables, we find TRW = 199.0±1.4 MeV. In the same picture we also report the critical temperature obtained for Nf = 2 + …
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 9
Figure 9. Figure 9: FIG. 9 [PITH_FULL_IMAGE:figures/full_fig_p008_9.png]
Figure 8
Figure 8. Figure 8: FIG. 8 [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p009_10.png]
Figure 9
Figure 9. Figure 9: the agreement is reasonable, given also that, [PITH_FULL_IMAGE:figures/full_fig_p009_9.png]
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p010_11.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

103 extracted references · 5 canonical work pages

  1. [21]

    Borsanyi, Z

    S. Borsanyi, Z. Fodor, D. A. Godzieba, R. Kara, P. Parotto, D. Sexty, and R. Vig, Topological features of the deconfinement transition, Phys. Rev. D107, 054514 (2023), arXiv:2212.08684 [hep-lat]

  2. [1]

    INTRODUCTION Many of the non-perturbative features which charac- terize QCD and QCD-like theories are related to the existence of non-trivial topological properties. Most of these properties are encoded in the dependence of the free (vacuum) energy density on the topological parame- terθ, i.e., the coupling to the winding numberQthat can be added to the Q...

  3. [2]

    NUMERICAL SET-UP We consider the finite temperature partition function ofN f = 2 + 1 + 1 QCD discretized via stout-improved staggered fermions as follows: Z= Z [DU]e−SYM Y f=u,d,s,c det (Mst,f) 1 4,(1) where SYM =− β 3 X i µ̸=ν 5 6W 1×1 i,µν− 1 12W 1×2 i,µν (Mst,f)ij =am fδij + 4X ν=1 ηi;ν 2 eaµfδν,4U (4) i;νδij−ˆν− e−aµfδν,4U (4)† i−ˆν;νδij+ˆν (2) are re...

  4. [3]

    THE ROBERGE-WEISS LINE The Roberge-Weiss (RW) symmetry is a remnant of center symmetry in full QCD. For aSU(Nc) lattice gauge theory with periodic boundary conditions in the Eu- clidean time, a center transformation consists of mul- tiplying the temporal links on a given time slice by a given element of the gauge group center, namely, by exp(2πik/Nc),k= 0...

  5. [5]

    Non-perturbative aspects of fundamental inter- actions, in the Standard Model and beyond

    CONCLUSIONS In this study, we have investigated some properties of Nf = 2 + 1 + 1 QCD and physical quark masses along the RW line, which is a modification of the standard ther- mal line along which an exact residual center symmetry exists and gets spontaneously broken at a critical temper- atureT RW . We have adopted a stout-smeared staggered discretizati...

  6. [6]

    Vicari and H

    E. Vicari and H. Panagopoulos, Theta dependence of SU(N) gauge theories in the presence of a topological term, Phys. Rept.470, 93 (2009), arXiv:0803.1593 [hep- th]

  7. [7]

    Bonanno, C

    C. Bonanno, C. Bonati, and M. D’Elia, Strong CP problem, theta term and QCD topological properties, (2025), arXiv:2510.03059 [hep-lat]

  8. [8]

    Bonanno, C

    C. Bonanno, C. Bonati, and M. D’Elia, Topological susceptibility and qcd at finite theta angle (2026), arXiv:2604.28035 [hep-lat]

Show all 103 references
  1. [9]

    D. J. Gross, R. D. Pisarski, and L. G. Yaffe, QCD and Instantons at Finite Temperature, Rev. Mod. Phys.53, 43 (1981)

  2. [10]

    R. D. Pisarski and L. G. Yaffe, The density of instantons at finite temperature, Phys. Lett. B97, 110 (1980)

  3. [11]

    Teper, Axial Anomaly Suppression (and Axial U(1) Symmetry Restoration) at High Temperatures: A Lat- tice Monte Carlo S tudy, Phys

    M. Teper, Axial Anomaly Suppression (and Axial U(1) Symmetry Restoration) at High Temperatures: A Lat- tice Monte Carlo S tudy, Phys. Lett. B171, 81 (1986)

  4. [12]

    Teper, The Topological Susceptibility in SU(2) Lat- tice Gauge Theory: An Exploratory Study, Phys

    M. Teper, The Topological Susceptibility in SU(2) Lat- tice Gauge Theory: An Exploratory Study, Phys. Lett. B171, 86 (1986)

  5. [13]

    All´ es, M

    B. All´ es, M. D’Elia, and A. Di Giacomo, Topology at zero and finiteTin SU(2) Yang-Mills theory, Phys. Lett. B412, 119 (1997), arXiv:hep-lat/9706016

  6. [14]

    All´ es, M

    B. All´ es, M. D’Elia, and A. Di Giacomo, Topological susceptibility at zero and finiteTin SU(3) Yang-Mills theory, Nucl. Phys. B494, 281 (1997), [Erratum: Nucl. Phys. B 679, 397–399 (2004)], arXiv:hep-lat/9605013

  7. [15]

    Del Debbio, H

    L. Del Debbio, H. Panagopoulos, and E. Vicari, Topo- logical susceptibility of SU(N) gauge theories at finite temperature, JHEP09, 028, arXiv:hep-th/0407068

  8. [16]

    Berkowitz, M

    E. Berkowitz, M. I. Buchoff, and E. Rinaldi, Lattice QCD input for axion cosmology, Phys. Rev. D92, 034507 (2015), arXiv:1505.07455 [hep-ph]

  9. [17]

    Borsanyi, M

    S. Borsanyi, M. Dierigl, Z. Fodor, S. D. Katz, S. W. Mages, D. Nogradi, A. Redondo, J. a nd Ringwald, and K. K. Szabo, Axion cosmology, lattice QCD and the dilute instanton gas, Phys. Lett. B752, 175 (2016), arXiv:1508.06917 [hep-lat]

  10. [18]

    Frison, R

    J. Frison, R. Kitano, H. Matsufuru, S. Mori, and N. Ya- mada, Topological susceptibility at high temperature on the lattice, JHEP09, 021, arXiv:1606.07175 [hep-lat]. 12

  11. [19]

    Kitano and N

    R. Kitano and N. Yamada, Topology in QCD and the axion abundance, JHEP10, 136, arXiv:1506.00370 [hep-ph]

  12. [20]

    Borsanyi and D

    S. Borsanyi and D. Sexty, Topological susceptibility of pure gauge theory using Density of States, Phys. Lett. B815, 136148 (2021), arXiv:2101.03383 [hep-lat]

  13. [22]

    Bonanno, M

    C. Bonanno, M. D’Elia, and L. Verzichelli, Theθ- dependence of the SU(N) critical temperature at large N, JHEP02, 156, arXiv:2312.12202 [hep-lat]

  14. [23]

    Bonati, M

    C. Bonati, M. Cardinali, and M. D’Elia,θdepen- dence in trace deformedSU(3) Yang-Mills theory: a lattice study, Phys. Rev. D98, 054508 (2018), arXiv:1807.06558 [hep-lat]

  15. [24]

    Bonati, M

    C. Bonati, M. Cardinali, M. D’Elia, and F. Mazziotti, θ-dependence and center symmetry in Yang-Mills theo- ries, Phys. Rev. D101, 034508 (2020), arXiv:1912.02662 [hep-lat]

  16. [25]

    Roberge and N

    A. Roberge and N. Weiss, Gauge Theories With Imagi- nary Chemical Potential and the Phases of QCD, Nucl. Phys. B275, 734 (1986)

  17. [26]

    Bonati, M

    C. Bonati, M. D’Elia, M. Mariti, M. Mesiti, F. Negro, and F. Sanfilippo, Roberge-Weiss endpoint at the phys- ical point ofN f = 2 + 1 QCD, Phys. Rev. D93, 074504 (2016), arXiv:1602.01426 [hep-lat]

  18. [27]

    Bonati, E

    C. Bonati, E. Calore, M. D’Elia, M. Mesiti, F. Negro, F. Sanfilippo, S. F. Schifano, G. Silvi, and R. Tripic- cione, Roberge-Weiss endpoint and chiral symmetry restoration inNf = 2+1 QCD, Phys. Rev. D99, 014502 (2019), arXiv:1807.02106 [hep-lat]

  19. [28]

    Cuteri, J

    F. Cuteri, J. Goswami, F. Karsch, A. Lahiri, M. Neu- mann, O. Philipsen, C. Schmidt, and A. Sciarra, Toward the chiral phase transition in the Roberge-Weiss plane, Phys. Rev. D106, 014510 (2022), arXiv:2205.12707 [hep-lat]

  20. [29]

    Cardinali, M

    M. Cardinali, M. D’Elia, F. Garosi, and M. Giordano, Localization properties of Dirac modes at the Roberge- Weiss phase transition, Phys. Rev. D105, 014506 (2022), arXiv:2110.10029 [hep-lat]

  21. [30]

    D’Elia, M

    M. D’Elia, M. Nacci, and K. Zambello, The Roberge- Weiss transition as a probe for conformality in many- flavor QCD, (2026), arXiv:2604.07284 [hep-lat]

  22. [31]

    Weisz, Continuum Limit Improved Lattice Action for Pure Yang-Mills Theory

    P. Weisz, Continuum Limit Improved Lattice Action for Pure Yang-Mills Theory. 1., Nucl. Phys. B212, 1 (1983)

  23. [32]

    Curci, P

    G. Curci, P. Menotti, and G. Paffuti, Symanzik’s Im- proved Lagrangian for Lattice Gauge Theory, Phys. Lett. B130, 205 (1983), [Erratum: Phys.Lett.B 135, 516 (1984)]

  24. [33]

    J. B. Kogut and L. Susskind, Hamiltonian Formulation of Wilson’s Lattice Gauge Theories, Phys. Rev. D11, 395 (1975)

  25. [34]

    Morningstar and M

    C. Morningstar and M. J. Peardon, Analytic smearing of SU(3) link variables in lattice QCD, Phys. Rev. D 69, 054501 (2004), arXiv:hep-lat/0311018

  26. [35]

    Borsanyiet al., Leading hadronic contribution to the muon magnetic moment from lattice QCD, Nature593, 51 (2021), arXiv:2002.12347 [hep-lat]

    S. Borsanyiet al., Leading hadronic contribution to the muon magnetic moment from lattice QCD, Nature593, 51 (2021), arXiv:2002.12347 [hep-lat]

  27. [36]

    Berg, Dislocations and topological background in the lattice O(3)σ-model, Physics Letters B104, 475 (1981)

    B. Berg, Dislocations and topological background in the lattice O(3)σ-model, Physics Letters B104, 475 (1981)

  28. [37]

    Iwasaki and T

    Y. Iwasaki and T. Yoshie, Instantons and Topological Charge in Lattice Gauge Theory, Phys. Lett. B131, 159 (1983)

  29. [38]

    S. Itoh, Y. Iwasaki, and T. Yoshie, Stability of Instan- tons on the Lattice and the Renormalized Trajectory, Phys. Lett. B147, 141 (1984)

  30. [39]

    Teper, Instantons in the Quantized SU(2) Vacuum: A Lattice Monte Carlo Investigation, Phys

    M. Teper, Instantons in the Quantized SU(2) Vacuum: A Lattice Monte Carlo Investigation, Phys. Lett. B162, 357 (1985)

  31. [40]

    Ilgenfritz, M

    E.-M. Ilgenfritz, M. laursen, M. M¨ uller-Preuβker, G. Schierholz, and H. Schiller, First evidence for the existence of instantons in the quantized SU(2) lattice vacuum, Nuclear Physics B268, 693 (1986)

  32. [41]

    Campostrini, A

    M. Campostrini, A. Di Giacomo, H. Panagopoulos, and E. Vicari, Topological Charge, Renormalization and Cooling on the Lattice, Nucl. Phys. B329, 683 (1990)

  33. [42]

    Alles, L

    B. Alles, L. Cosmai, M. D’Elia, and A. Papa, Topol- ogy in 2-D CP**(N-1) models on the lattice: A Critical comparison of different cooling techniques, Phys. Rev. D62, 094507 (2000), arXiv:hep-lat/0001027

  34. [43]

    Bonati and M

    C. Bonati and M. D’Elia, Comparison of the gradient flow with cooling inSU(3) pure gauge theory, Phys. Rev. D89, 105005 (2014), arXiv:1401.2441 [hep-lat]

  35. [44]

    Alexandrou, A

    C. Alexandrou, A. Athenodorou, K. Cichy, A. Dro- mard, E. Garcia-Ramos, K. Jansen, U. Wenger, and F. Zimmermann, Comparison of topological charge def- initions in Lattice QCD, Eur. Phys. J. C80, 424 (2020), arXiv:1708.00696 [hep-lat]

  36. [45]

    Alexandrou, A

    C. Alexandrou, A. Athenodorou, and K. Jansen, Topo- logical charge using cooling and the gradient flow, Phys. Rev. D92, 125014 (2015), arXiv:1509.04259 [hep-lat]

  37. [46]

    Del Debbio, H

    L. Del Debbio, H. Panagopoulos, and E. Vicari, theta dependence of SU(N) gauge theories, JHEP08, 044, arXiv:hep-th/0204125

  38. [47]

    Bonati, M

    C. Bonati, M. D’Elia, and A. Scapellato,θdependence in SU(3) Yang-Mills theory from analytic continuation, Phys. Rev. D93, 025028 (2016), arXiv:1512.01544 [hep- lat]

  39. [48]

    Bonanno, C

    C. Bonanno, C. Bonati, M. Papace, and D. Vadacchino, Theθ-dependence of the Yang-Mills spectrum from an- alytic continuation, JHEP05, 163, arXiv:2402.03096 [hep-lat]

  40. [49]

    M. G. Alford, A. Kapustin, and F. Wilczek, Imagi- nary chemical potential and finite fermion density on the lattice, Phys. Rev. D59, 054502 (1999), arXiv:hep- lat/9807039

  41. [50]

    Lombardo, Finite density (might well be easier) at finite temperature, Nucl

    M.-P. Lombardo, Finite density (might well be easier) at finite temperature, Nucl. Phys. B Proc. Suppl.83, 375 (2000), arXiv:hep-lat/9908006

  42. [51]

    de Forcrand and O

    P. de Forcrand and O. Philipsen, The QCD phase diagram for small densities from imaginary chemical potential, Nucl. Phys. B642, 290 (2002), arXiv:hep- lat/0205016

  43. [52]

    D’Elia and M.-P

    M. D’Elia and M.-P. Lombardo, Finite density QCD via imaginary chemical potential, Phys. Rev. D67, 014505 (2003), arXiv:hep-lat/0209146

  44. [53]

    de Forcrand and O

    P. de Forcrand and O. Philipsen, The QCD phase di- agram for three degenerate flavors and small baryon density, Nucl. Phys. B673, 170 (2003), arXiv:hep- lat/0307020

  45. [54]

    D’Elia and M

    M. D’Elia and M. P. Lombardo, QCD thermodynamics from an imaginary mu(B): Results on the four flavor lat- tice model, Phys. Rev. D70, 074509 (2004), arXiv:hep- lat/0406012. 13

  46. [55]

    Giudice and A

    P. Giudice and A. Papa, Real and imaginary chemical potential in two color QCD, Phys. Rev. D69, 094509 (2004), arXiv:hep-lat/0401024

  47. [56]

    Chen and X.-Q

    H.-S. Chen and X.-Q. Luo, Phase diagram of QCD at finite temperature and chemical potential from lattice simulations with dynamical Wilson quarks, Phys. Rev. D72, 034504 (2005), arXiv:hep-lat/0411023

  48. [57]

    Azcoiti, G

    V. Azcoiti, G. Di Carlo, A. Galante, and V. Laliena, Phase diagram of QCD with four quark flavors at finite temperature and baryon density, Nucl. Phys. B723, 77 (2005), arXiv:hep-lat/0503010

  49. [58]

    D’Elia, F

    M. D’Elia, F. Di Renzo, and M. P. Lombardo, The Strongly interacting quark gluon plasma, and the crit- ical behaviour of QCD at imaginary mu, Phys. Rev. D 76, 114509 (2007), arXiv:0705.3814 [hep-lat]

  50. [59]

    D’Elia and F

    M. D’Elia and F. Sanfilippo, The Order of the Roberge- Weiss endpoint (finite size transition) in QCD, Phys. Rev. D80, 111501 (2009), arXiv:0909.0254 [hep-lat]

  51. [60]

    P. Cea, L. Cosmai, M. D’Elia, C. Manneschi, and A. Papa, Analytic continuation of the critical line: Sug- gestions for QCD, Phys. Rev. D80, 034501 (2009), arXiv:0905.1292 [hep-lat]

  52. [61]

    de Forcrand and O

    P. de Forcrand and O. Philipsen, Constraining the QCD phase diagram by tricritical lines at imaginary chem- ical potential, Phys. Rev. Lett.105, 152001 (2010), arXiv:1004.3144 [hep-lat]

  53. [62]

    Bonati, G

    C. Bonati, G. Cossu, M. D’Elia, and F. Sanfilippo, The Roberge-Weiss endpoint inN f = 2 QCD, Phys. Rev. D 83, 054505 (2011), arXiv:1011.4515 [hep-lat]

  54. [63]

    P. Cea, L. Cosmai, M. D’Elia, A. Papa, and F. Sanfil- ippo, The critical line of two-flavor QCD at finite isospin or baryon densities from imaginary chemical potentials, Phys. Rev. D85, 094512 (2012), arXiv:1202.5700 [hep- lat]

  55. [64]

    Wu and X.-F

    L.-K. Wu and X.-F. Meng, Nature of the Roberge- Weiss transition end points in two-flavor lattice QCD with Wilson quarks, Phys. Rev. D87, 094508 (2013), arXiv:1303.0336 [hep-lat]

  56. [65]

    Philipsen and C

    O. Philipsen and C. Pinke, Nature of the Roberge-Weiss transition inN f = 2 QCD with Wilson fermions, Phys. Rev. D89, 094504 (2014), arXiv:1402.0838 [hep-lat]

  57. [66]

    Bonati, P

    C. Bonati, P. de Forcrand, M. D’Elia, O. Philipsen, and F. Sanfilippo, Chiral phase transition in two-flavor QCD from an imaginary chemical potential, Phys. Rev. D90, 074030 (2014), arXiv:1408.5086 [hep-lat]

  58. [67]

    Wu and X.-F

    L.-K. Wu and X.-F. Meng, Nature of Roberge-Weiss transition endpoints for heavy quarks inN f = 2 lattice QCD with Wilson fermions, Phys. Rev. D90, 094506 (2014), arXiv:1405.2425 [hep-lat]

  59. [68]

    Czaban, F

    C. Czaban, F. Cuteri, O. Philipsen, C. Pinke, and A. Sciarra, Roberge-Weiss transition inN f = 2 QCD with Wilson fermions andN τ = 6, Phys. Rev. D93, 054507 (2016), arXiv:1512.07180 [hep-lat]

  60. [69]

    Philipsen and A

    O. Philipsen and A. Sciarra, Finite Size and Cut-Off Effects on the Roberge-Weiss Transition inNf = 2 QCD with Staggered Fermions, Phys. Rev. D101, 014502 (2020), arXiv:1909.12253 [hep-lat]

  61. [70]

    B. B. Brandt, A. Chabane, V. Chelnokov, F. Cuteri, G. Endr˝ odi, and C. Winterowd, Light Roberge-Weiss tricritical endpoint at imaginary isospin and baryon chemical potential, Phys. Rev. D109, 034515 (2024), arXiv:2207.10117 [hep-lat]

  62. [71]

    B. B. Brandt, V. Chelnokov, F. Cuteri, G. Endrodi, and C. Winterowd, The light Roberge-Weiss tricritical endpoint at imaginary isospin chemical potential, PoS LA TTICE2021, 164 (2022), arXiv:2207.10586 [hep- lat]

  63. [72]

    Endrodi, G

    G. Endrodi, G. D. Moore, A. Pieczynski, and A. Scia- rra, Exact center symmetry and first-order phase tran- sition in QCD with three degenerate dynamical quarks, (2026), arXiv:2603.22046 [hep-lat]

  64. [73]

    Endr˝ odi, G

    G. Endr˝ odi, G. D. Moore, and A. Sciarra, First- order phase transition in dynamical 3-flavor QCD at imaginary isospin, PoSLA TTICE2024, 159 (2025), arXiv:2501.00792 [hep-lat]

  65. [74]

    D’Elia, L

    M. D’Elia, L. Maio, K. Zambello, and G. Zanichelli, Roberge-Weiss transition for QCD in a magnetic background, Phys. Rev. D111, 094509 (2025), arXiv:2502.19294 [hep-lat]

  66. [75]

    Zambello, M

    K. Zambello, M. D’Elia, L. Maio, and G. Zanichelli, The Roberge-Weiss endpoint in (2+1)-flavor QCD with background magnetic fields, PoSLA TTICE2024, 169 (2025), arXiv:2412.06326 [hep-lat]

  67. [76]

    H. W. J. Blote, E. Luijten, and J. R. Heringa, Ising universality in three dimensions: a monte carlo study, Journal of Physics A: Mathematical and General28, 6289–6313 (1995)

  68. [77]

    A. D. Kennedy, I. Horvath, and S. Sint, A New exact method for dynamical fermion computations with non- local actions, Nucl. Phys. B Proc. Suppl.73, 834 (1999), arXiv:hep-lat/9809092

  69. [78]

    M. A. Clark and A. D. Kennedy, Accelerating dynamical fermion computations using the rational hybrid Monte Carlo (RHMC) algorithm with multiple pseudofermion fields, Phys. Rev. Lett.98, 051601 (2007), arXiv:hep- lat/0608015. [74]https://github.com/GPU-LQCD/OpenStaPLE

  70. [79]

    Bonati, E

    C. Bonati, E. Calore, S. Coscetti, M. D’Elia, M. Mesiti, F. Negro, S. F. Schifano, G. Silvi, and R. Tripiccione, Design and optimization of a portable LQCD Monte Carlo code using OpenACC, Int. J. Mod. Phys. C28, 1750063 (2017), arXiv:1701.00426 [hep-lat]

  71. [80]

    Bonati, E

    C. Bonati, E. Calore, M. D’Elia, M. Mesiti, F. Negro, F. Sanfilippo, S. F. Schifano, G. Silvi, and R. Tripic- cione, Portable multi-node LQCD Monte Carlo simula- tions using OpenACC, Int. J. Mod. Phys. C29, 1850010 (2018), arXiv:1801.01473 [hep-lat]

  72. [81]

    DeGrand and R

    T. DeGrand and R. Hoffmann, QCD with one com- pact spatial dimension, JHEP02, 022, arXiv:hep- lat/0612012

  73. [82]

    DeGrand, R

    T. DeGrand, R. Hoffmann, and J. Najjar, More about QCD on compact spaces, JHEP01, 032, arXiv:0711.4290 [hep-lat]

  74. [83]

    Lucini, A

    B. Lucini, A. Patella, and C. Pica, Baryon currents in QCD with compact dimensions, Phys. Rev. D75, 121701 (2007), arXiv:hep-th/0702167

  75. [84]

    Lucini and A

    B. Lucini and A. Patella, Discrete symmetry breaking and baryon currents in U(N) and SU(N) gauge theories, Phys. Rev. D79, 125030 (2009), arXiv:0904.3479 [hep- th]

  76. [85]

    Dong and K.-F

    S.-J. Dong and K.-F. Liu, Stochastic estimation with Z(2) noise, Phys. Lett. B328, 130 (1994), arXiv:hep- lat/9308015

  77. [86]

    Vafa and E

    C. Vafa and E. Witten, Parity Conservation in QCD, Phys. Rev. Lett.53, 535 (1984)

  78. [87]

    Di Vecchia and G

    P. Di Vecchia and G. Veneziano, Chiral Dynamics in the Large N Limit, Nucl. Phys. B171, 253 (1980). 14

  79. [88]

    Grilli di Cortona, E

    G. Grilli di Cortona, E. Hardy, J. Pardo Vega, and G. Villadoro, The QCD axion, precisely, JHEP01, 034, arXiv:1511.02867 [hep-ph]

  80. [89]

    Boccaletti and D

    A. Boccaletti and D. Nogradi, The semi-classical ap- proximation at high temperature revisited, JHEP03, 045, arXiv:2001.03383 [hep-ph]

  81. [90]

    Bonati, M

    C. Bonati, M. D’Elia, M. Mariti, G. Martinelli, M. Mesiti, F. Negro, F. Sanfilippo, and G. Villadoro, Axion phenomenology andθ-dependence fromN f = 2 + 1 lattice QCD, JHEP03, 155, arXiv:1512.06746 [hep-lat]

  82. [91]

    Borsanyiet al., Calculation of the axion mass based on high-temperature lattice quantum chromodynamics, Nature539, 69 (2016), arXiv:1606.07494 [hep-lat]

    S. Borsanyiet al., Calculation of the axion mass based on high-temperature lattice quantum chromodynamics, Nature539, 69 (2016), arXiv:1606.07494 [hep-lat]

  83. [92]

    Alexandrou, A

    C. Alexandrou, A. Athenodorou, K. Cichy, M. Con- stantinou, D. e. P. Horkel, K. Jansen, G. Koutsou, and C. Larkin, Topological susceptibility from twisted mass fermions using spectral projectors and the gradient flow, Phys. Rev. D97, 074503 (2018), arXiv:1709.06596 [hep- lat]

  84. [93]

    Athenodorou, C

    A. Athenodorou, C. Bonanno, C. Bonati, G. Clemente, F. . D’Angelo, M. D’Elia, L. Maio, G. Martinelli, F. Sanfilippo, and A. Todaro, Topological susceptibil- ity ofN f = 2 + 1 QCD from staggered fermions spec- tral projectors at high tem peratures, JHEP10, 197, arXiv:2208.08921...

  85. [94]

    Alles, M

    B. Alles, M. D’Elia, and A. Di Giacomo, Topologi- cal susceptibility in full QCD at zero and finite tem- perature, Phys. Lett. B483, 139 (2000), arXiv:hep- lat/0004020

  86. [95]

    Petreczky, H.-P

    P. Petreczky, H.-P. Schadler, and S. Sharma, The topological susceptibility in finite temperature QCD and axion cosmology, Phys. Lett. B762, 498 (2016), arXiv:1606.03145 [hep-lat]

  87. [96]

    Bonati, M

    C. Bonati, M. D’Elia, G. Martinelli, F. Negro, F. San- filippo, and A. ˜Todaro, Topology in full QCD at high temperature: a multicanonical approach, JHEP11, 170, arXiv:1807.07954 [hep-lat]

  88. [97]

    Burger, E.-M

    F. Burger, E.-M. Ilgenfritz, M. P. Lombardo, M. M¨ uller- Preussker, and A. T˜ runin, Topology (and axion’s prop- erties) from lattice QCD with a dynamical charm, Nucl. Phys. A967, 880 (2017), arXiv:1705.01847 [hep-lat]

  89. [98]

    Chen, T.-W

    Y.-C. Chen, T.-W. Chiu, and T.-H. Hsieh (TWQCD), Topological susceptibility in finite temperature QCD with physical (u/d,s,c) domain-wall quarks, Phys. Rev. D106, 074501 (2022), arXiv:2204.01556 [hep-lat]

  90. [99]

    A. Y. Kotov, M. P. Lombardo, and A. Trunin, Topo- logical observables andθdependence in high temper- ature QCD from lattice simulations, JHEP09, 045, arXiv:2502.15407 [hep-lat]

  91. [100]

    Bonati, M

    C. Bonati, M. D’Elia, H. Panagopoulos, and E. Vi- cari, Change ofθdependence in 4d su(n) gauge theories across the deconfinement transition, Phys. Rev. Lett. 110, 252003 (2013), arXiv:1301.7640 [hep-lat]

  92. [101]

    Bonati, M

    C. Bonati, M. D’Elia, P. Rossi, and E. Vicari,θdepen- dence of 4DSU(N) gauge theories in the large-Nlimit, Phys. Rev. D94, 085017 (2016), arXiv:1607.06360 [hep- lat]

  93. [102]

    D’Elia, Field theoretical approach to the study ofθ dependence in Yang-Mills theories on the latti ce, Nucl

    M. D’Elia, Field theoretical approach to the study ofθ dependence in Yang-Mills theories on the latti ce, Nucl. Phys. B661, 139 (2003), arXiv:hep-lat/0302007 [hep- lat]

  94. [103]

    Giusti, S

    L. Giusti, S. Petrarca, and B. Taglienti, Theta depen- dence of the vacuum energy in the SU(3) gauge the- ory from the lattice, Phys. Rev. D76, 094510 (2007), arXiv:0705.2352 [hep-th]

  95. [104]

    Panagopoulos and E

    H. Panagopoulos and E. Vicari, The 4D SU(3) gauge theory with an imaginaryθterm, JHEP11, 119, arXiv:1109.6815 [hep-lat]

Pith tools

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