Pith. sign in

REVIEW 3 major objections 4 minor 21 references

The confined-deconfined surface tension in SU(N) gauge theories at large N

T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash

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

Pith's one-line read Pure SU(N) gauge theory shows interface tension and latent heat scaling as N^2 at large N.

desk verdict First continuum extrapolation of SU(N) interface tension, with plausible N^2 scaling, but the sigma continuum limit rests on two lattice spacings and should be treated as preliminary. read the letter →

arxiv 2502.01396 v1 pith:VIXUGNXS submitted 2025-02-03 hep-lat

Tobias Rindlisbacher , Kari Rummukainen , Ahmed Salami This is my paper · ORCID
classification hep-lat
keywords SU(N)gaugetheorylarge-Nlimitconfined-deconfinedtransitioninterfacetensionlatentheatlatticePolyakovloopcapillarywaves
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

Pure SU(N) gauge theories with N colors undergo a first-order confinement-deconfinement transition. This paper reports lattice measurements showing that the interface tension between the confined and deconfined phases, and the latent heat released in the transition, both grow as $N^{2}$ at large N. The result matters because the interface tension controls bubble nucleation and the dynamics of the transition, and was previously poorly determined. The paper extracts it from the thermal fluctuation spectrum of the phase interface using a constrained mixed-phase simulation, and takes a continuum limit.

What carries the argument

The central object is the capillary-wave spectrum of the phase interface, Eq. (1): for a nearly flat interface of area $L^2$, the Fourier modes of its height satisfy $\langle |\hat{z}(n_x,n_y)|^2\rangle = T/(4\pi^2 \sigma (n_x^2+n_y^2))$, so the interface tension $\sigma$ is the slope of $1/\langle|\hat{z}|^2\rangle$ versus $k_x^2+k_y^2$. The argument works because the simulation is held in the mixed phase, so the interface can be located by thresholding the smeared Polyakov loop, and the smearing distortion is removed by dividing by the known Fourier transform of the smearing kernel. This turns the measurement of a strongly suppressed mixed-phase probability into a measurement of equilibrium height fluctuations, which are not suppressed and allow large volumes.

What would settle it

Take data at a third lattice spacing, such as $N_t=10$, at one $N$ (for instance $N=8$) and check whether $\sigma/T_c^3$ continues to follow the linear-in-$a^2$ line used here; the current fit predicts the continuum value $1.004(94)$ from $N_t=6,8$, so an $N_t=10$ point that falls off that line, or a continuum value that disagrees with the $N^2$ fit, would falsify the scaling.

Watch

Extended reading notes

Core claim

In the continuum limit, the confined-deconfined interface tension and the latent heat of pure SU(N) gauge theory at large N are described by $\sigma/T_c^3 = 0.0189(11) N^2 - 0.190(19)$ and $L/T_c^4 = 0.354(2) N^2 - 1.65(10)$ for $N \ge 5$. The paper establishes this by simulating $N=4, 5, 8, 10$ and 16 on lattices with inverse temperatures $N_t = 5\ldots 8$ (and $N_t=6$ for $N=16$), constraining the real part of the average Polyakov loop to stay in the mixed-phase region so that two interfaces coexist. The interface tension is read off from the long-wavelength spectrum of interface height fluctuations, after undoing the effect of the smearing kernel analytically. The latent heat comes from the plaquette discontinuity at the critical coupling, using the measured critical couplings to evaluate the $\beta$ function. Both quantities extrapolate linearly in the squared lattice spacing to the continuum, and the continuum values follow the $N^2 + \mathrm{const}$ form.

Load-bearing premise

The result depends on the continuum extrapolation of the interface tension, which is a straight line in the squared lattice spacing drawn through only two lattice spacings ($N_t=6$ and $N_t=8$); if higher-order corrections in the spacing are not negligible or vary with $N$, the quoted $N^2$ coefficient would shift.

Editorial extensions

If this is right

  • The fitted value $\sigma/T_c^3 = 0.0189(11) N^2 - 0.190(19)$ gives a reference interface tension for all $N$, replacing the previously uncertain estimates.
  • The latent heat fit $L/T_c^4 = 0.354(2) N^2 - 1.65(10)$ reduces the error on the large-$N$ transition strength by an order of magnitude compared with earlier work.
  • Both scalings confirm the large-$N$ expectation that the deconfinement transition becomes stronger quadratically with $N$, with finite-$N$ corrections encoded in the small constant terms.
  • The constrained mixed-phase method is shown to work up to $N=16$ and avoids the supercritical slowing down of multicanonical methods, so it can be applied to transitions where the mixed-phase probability is extremely suppressed.

Reading between the lines

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

  • The exponentially small mixed-phase probability implied by $\sigma \propto N^2$ means that in strongly first-order large-$N$ theories bubble nucleation proceeds through a very suppressed channel; the measured tension is the input needed to quantify supercooling and nucleation rates, which the paper does not compute.
  • A direct test of the large-$N$ prediction would be a continuum extrapolation at $N=16$ from a second lattice spacing; the fit predicts $\sigma/T_c^3 \approx 4.65$ while the single $N_t=6$ value is $5.36(9)$, so the extrapolation must move down by roughly 0.7 if the linear $a^2$ assumption holds.
  • The same fluctuation-spectrum technique could be applied to other order parameters, for example in theories with matter multiplets, where the interface tension is less constrained; nothing in the method ties it to pure glue.
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 / 4 minor

Summary. This proceedings paper reports a high-precision lattice study of the confined-deconfined interface tension and latent heat in pure SU(N) gauge theory at large N. Using the Moore-Turok mixed-phase method, the authors restrict the average Polyakov loop to a narrow window in the coexistence region, determine beta_c by demanding a flat restricted distribution, and extract the interface tension from the capillary-wave spectrum of the interface after applying a semi-analytic smearing-kernel correction. Simulations are performed for N = 4, 5, 8, 10 and 16 on lattices with N_t from 5 to 8. The paper finds sigma/T_c^3 = 0.0189(11) N^2 - 0.190(19) in the continuum limit and L/T_c^4 = 0.354(2) N^2 - 1.65(10) for N >= 5, and concludes that both quantities scale as N^2 at large N.

Significance. If the result holds, this paper provides the first precise lattice determination of the interface tension in the large-N limit and sharpens the existing latent-heat determination by more than an order of magnitude. The methodology has genuine strengths: the mixed-phase restriction avoids supercritical slowing down, the smearing-kernel correction collapses the Fourier spectrum onto a single curve, the beta_c determination is visually and statistically robust, and the latent-heat result is compatible with the older SU(N) determination of Lucini, Teper and Wenger while being much more precise. The use of independent observables for sigma and L avoids a circular N^2 fit. However, the central large-N coefficients currently rest on a two-point continuum extrapolation for the interface tension and on a saturated interpolation ansatz for the beta-function derivative, so the numerical claims are not yet as secure as the paper's wording suggests.

major comments (3)
  1. [Section 3, Table 1, Eq. (9)] The continuum value of sigma/T_c^3 for each N is obtained from exactly two lattice spacings, N_t = 6 and N_t = 8, via a linear extrapolation in a^2. Such a two-point extrapolation has zero internal degrees of freedom, cannot detect O(a^4) or logarithmic corrections, and no systematic error is assigned to it. Since the fitted N^2 coefficient 0.0189(11) is driven by four continuum points, a correlated shift of a few percent in the N = 5 and N = 8 intercepts would move the slope by several quoted errors. Please add at least one more lattice spacing for at least one representative N, or provide a quantitative estimate of the truncation uncertainty and include it in the final error.
  2. [Section 3, Fig. 4] The k -> 0 extrapolation of the kernel-corrected Fourier spectrum is not documented in enough detail. The text does not state the fit ansatz (e.g., constant plus linear or quadratic term in k^2), the momentum range used, the number of Fourier modes, or the chi^2 per degree of freedom. Only a single representative case (SU(16) at N_t = 6) is shown. Because sigma is the inverse of the intercept, a small systematic tilt in this fit propagates directly into every continuum point in Table 1. Please report the fit range, fit form, goodness of fit, and a table showing the stability of the extracted sigma across smearing levels for each N.
  3. [Section 4, Eq. (12)] The derivative d beta / d ln a used in Eq. (11) is obtained from the three-parameter interpolation ansatz of Eq. (12), which the authors themselves state 'saturates the degrees of freedom of the fit.' With only four N_t values per N, this leaves no residual degrees of freedom, and the N_t = 5 point is then excluded from the continuum limit because it deviates. The continuum L/T_c^4 values and hence the N^2 coefficient in Eq. (13) are therefore partly dependent on the chosen functional form. Please add a stability test with a different interpolation ansatz or with the N_t = 5 point included, and propagate any resulting shift as a systematic uncertainty.
minor comments (4)
  1. [Throughout] The rendered text contains numerous missing spaces (e.g., 'Wepresentresultsfrom...', 'Theresultsoftheextrapolationgivestheinverseofsigma'), which makes the manuscript difficult to read. These should be fixed in the final version.
  2. [Abstract and Conclusion] The phrase 'we observe unambiguously' is stronger than the current two-point continuum extrapolation for sigma supports; a wording such as 'consistent with N^2 scaling' would better match the evidence presented.
  3. [Fig. 5, Eq. (13)] The large-N latent-heat fit is stated to apply for N >= 5, but the text does not explicitly list which continuum points enter the fit (N = 5, 8, 10) and whether N = 4 is excluded from both the sigma and latent-heat fits. Please state this explicitly in a caption or in the text.
  4. [Section 3, Fig. 5] The SU(3) comparison points from refs. [20] and [21] are shown in the figures but not discussed in the text. Since they provide an important external consistency check, a brief sentence describing the comparison would be useful.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the N^2 scalings are fitted outcomes from independent lattice measurements, not constructed inputs.

full rationale

The paper's central claim is a fit, not a construction. The interface tension values are extracted from the capillary-wave spectrum of the measured interface height fluctuations via Eq. (1), and the latent-heat values are extracted from the plaquette discontinuity combined with d beta / d ln a obtained from the measured beta_c(N_t) series via Eq. (12). Neither observable is defined in terms of the other or in terms of the final N^2 fits. The N^2 + const forms in Eqs. (9) and (13) are fitted after the continuum limits are taken and are not used as constraints in the measurements. The beta-function ansatz in Eq. (12) uses the known 1-loop N^2 leading term and fitted c_i parameters, but it is fitted to beta_c data, not to the latent-heat values, so the N^2 coefficient of L is not forced by the fit to L itself. The self-citations (refs. [13], [19], [20]) are methodological or comparison-only and do not supply the large-N scaling. The main weaknesses, such as the two-point linear a^2 continuum extrapolation for sigma and the saturating three-parameter beta_c interpolation, are statistical or model-selection risks rather than circular reductions. The exclusion of the N_t = 5 latent-heat point is likewise a data-selection concern, not evidence that the quoted scaling was put in by hand.

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

The central results are measurements, not derivations. They rest on the capillary-wave model for interface fluctuations, a specific threshold and smearing definition of the interface position, linear a^2 continuum extrapolation for sigma, and an interpolation ansatz for the beta function used in the latent heat. The N^2 scaling is an empirical fit, not a first-principles consequence derived in the paper.

free parameters (3)
  • Interface threshold fraction = 0.5
    Interface position in each column is defined by the crossing of half the deconfined Polyakov line expectation value; no variation of this threshold or systematic error is reported (Section 3).
  • Smearing parameter rho in Eq (2) = not stated
    Controls the UV filtering of the Polyakov line before interface detection; the paper does not give its value, though the stability over the number of smearing steps is shown in Fig 4.
  • c_i (i=1,2,3) in the beta-function interpolation ansatz (Eq 12) = not reported
    Three parameters per SU(N), fitted to beta_c(N_t) at N_t=5,6,7,8 to compute d beta / d ln a for the latent heat; the fit saturates the degrees of freedom, making the derivative model-dependent.
assumptions (4)
  • domain assumption Capillary wave spectrum Eq (1) relates the interface Fourier modes to the inverse interface tension, assuming a thin sheet with small amplitude fluctuations at long wavelengths.
    This is the central formula used to extract sigma from the measured interface fluctuations; it is standard but assumes no interface stiffness corrections or nonlinear effects in the fitted momentum range.
  • ad hoc to paper The interface position is determined by the threshold crossing of the smeared Polyakov line at half the deconfined expectation value.
    This definition is chosen for practical detection; the paper does not demonstrate independence of the result from the threshold value or from the smearing level beyond the shown data collapse.
  • domain assumption The continuum limit for the interface tension is linear in a^2 using only N_t=6 and N_t=8 lattice spacings.
    Standard O(a^2) scaling is assumed; with only two lattice spacings, O(a^4) terms or logarithmic corrections cannot be excluded and no systematic error is assigned (Section 3, Table 1).
  • ad hoc to paper The interpolation ansatz Eq (12) for aTc as a function of beta_c, containing the 1-loop beta function plus three polynomial terms, is used to compute the derivative d beta / d ln a.
    The leading term is motivated by asymptotic freedom, but the three polynomial terms are ad hoc and saturate the degrees of freedom of the fit; the latent heat result depends on this interpolation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The confined-deconfined surface tension in SU(N) gauge theories at large N." pith.science (2026). https://pith.science/paper/VIXUGNXS

@misc{pith2026250201396,
  author       = {Pith},
  title        = {Pith review of: The confined-deconfined surface tension in SU(N) gauge theories at large N},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/VIXUGNXS}},
  note         = {Machine review of arXiv:2502.01396}
}
abstract

We present results from an investigation of the $N$-dependency of the confined-deconfined interface tension and latent heat in pure SU($N$) gauge theory at large $N$. The interface tension is determined by measuring the transverse fluctuations of the phase interface on large lattices with coexisting confined and deconfined phases. We observe unambiguously that both the interface tension and latent heat scale as $N^2$ at large $N$.

Figures

Figures reproduced from arXiv: 2502.01396 by the authors.

Figure 1
Figure 1. Left: schematic Polyakov line 𝑃𝐿 distribution at 𝛽𝑐 (dashed curve), with restriction to the central region (red lines). Right: two-phase configuration with two planar interfaces corresponding to the central region of 𝑃𝐿 distribution. 0,181 0,182 0,183 0,184 <PL> 0,6 0,8 1 1,2 1,4 1,6 Probability SU(8) 502 * 160 * 6 β = 44.5605 0,181 0,182 0,183 0,184 <PL> 0,6 0,8 1 1,2 1,4 1,6 Probability SU(8) 502 * 160 * 6 β = 44.… view at source ↗
Figure 2
Figure 2. Restricted average Polyakov line 𝑃𝐿 probability distributions of SU(8) gauge theory on a 6 × 502 × 160 lattice slightly below, at, and above the critical 𝛽. In the simulation 𝑃𝐿 was restricted to range 0.1805 ≤ 𝑃𝐿 ≤ 0.1845. range in the center of the distribution, see [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Capturing the confined-deconfined interface from a coexisting-phase-restricted simulation of SU(10) theory on a 602 × 240 × 6 lattice at 𝛽 = 69.9225, after 10 (left) and 70 (right) smearing steps. Dislocations still present at 10 smearing steps disappear at sufficiently large number of steps. describe the shape of the volume occupied by the deconfined phase, neglecting bulk fluctuations, so that 𝐴0 = |⟨𝑃𝐿⟩deconf.| i… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Left: extrapolation of the interface Fourier modes to 𝑘 → 0 limit at different number of smearing steps for SU(16) at 𝛽 = 179.855 on a 402 × 160 × 6 lattice. Right: same after kernel correction. Data collapses on a single curve, except for the lowest smearing level. Th…
Figure 5
Figure 5. Figure 5: Left: interface tension 𝜎/𝑇 3 𝑐 as a function of 𝑁 shown on a log-log plot. Dashed lines are large-𝑁 fits to 𝑁𝑡 = 6 and continuum limit results. The SU(3) 𝑁𝑡 = 6 point is from ref. [20] and not used in the fit. Right: latent heat 𝐿/𝑇 4 𝑐 , with the continuum large-𝑀 fi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 5 canonical work pages

  1. [6]

    Lucini, M

    B. Lucini, M. Teper and U. Wenger,Properties of the deconfining phase transition in SU(N) gauge theories, JHEP 02(2005) 033 [hep-lat/0502003]

    arXiv 2005

  2. [20]

    Iwasaki, K

    Y. Iwasaki, K. Kanaya, L. Karkkainen, K. Rummukainen and T. Yoshie,Interface tension in quenched QCD,Phys. Rev. D49(1994) 3540 [hep-lat/9309003]

    arXiv 1994

  3. [1]

    Classical Field Dynamics of the Electroweak Phase Transition

    G.D. Moore and N. Turok,Classical field dynamics of the electroweak phase transition, Phys. Rev. D55(1997) 6538 [hep-ph/9608350]

    work page Pith review arXiv 1997

  4. [2]

    ’t Hooft,A Planar Diagram Theory for Strong Interactions, Nucl

    G. ’t Hooft,A Planar Diagram Theory for Strong Interactions, Nucl. Phys. B72(1974) 461

    work page 1974

  5. [3]

    Witten,Anti-de Sitter space, thermal phase transition, and confinement in gauge theories, Adv

    E. Witten,Anti-de Sitter space, thermal phase transition, and confinement in gauge theories, Adv. Theor. Math. Phys.2(1998) 505 [hep-th/9803131]

    arXiv 1998

  6. [4]

    F.R. Ares, O. Henriksson, M. Hindmarsh, C. Hoyos and N. Jokela,Effective actions and bubble nucleation from holography, Phys. Rev. D105 (2022) 066020 [2109.13784]

    arXiv 2022

  7. [5]

    Lucini, A

    B. Lucini, A. Rago and E. Rinaldi,SU(𝑁𝑐) gauge theories at deconfinement, Phys. Lett. B 712 (2012) 279 [1202.6684]

    arXiv 2012

  8. [7]

    Wilson,Confinement of Quarks,Phys

    K.G. Wilson,Confinement of Quarks,Phys. Rev. D10(1974) 2445

    1974

Show all 21 references
  1. [8]

    Cabibbo and E

    N. Cabibbo and E. Marinari,A New Method for Updating SU(N) Matrices in Computer Simulations of Gauge Theories, Phys. Lett. B119 (1982) 387

    1982

  2. [9]

    Kennedy and B.J

    A.D. Kennedy and B.J. Pendleton,Improved Heat Bath Method for Monte Carlo Calculations in Lattice Gauge Theories,Phys. Lett. B156 (1985) 393

    1985

  3. [10]

    Brown and T.J

    F.R. Brown and T.J. Woch,Overrelaxed Heat Bath and Metropolis Algorithms for Accelerating Pure Gauge Monte Carlo Calculations, Phys. Rev. Lett.58(1987) 2394

    1987

  4. [11]

    de Forcrand and O

    P. de Forcrand and O. Jahn,Monte Carlo overrelaxation for SU(N) gauge theories, inQCD and Numerical Analysis III, (Berlin, Heidelberg), pp. 67–73, Springer Berlin Heidelberg, 3, 2005, DOI [hep-lat/0503041]

    2005 arXiv

  5. [12]

    Hila lattice simulation framework

    “Hila lattice simulation framework.”https://github.com/CFT-HY/HILA

  6. [13]

    Rummukainen, R

    K. Rummukainen, R. Seppä and D.J. Weir,Resolving the critical bubble in SU(8) deconfinement transition, in41st International Symposium on Lattice Field Theory, 1, 2025 [2501.17593]

    2025 arXiv

  7. [14]

    B.A.BergandT.Neuhaus, Multicanonicalensemble: ANewapproachtosimulatefirstorder phase transitions, Phys. Rev. Lett.68(1992) 9 [hep-lat/9202004]

    1992 arXiv

  8. [15]

    Wang and D.P

    F. Wang and D.P. Landau,Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of States,Phys. Rev. Lett.86 (2001) 2050 [cond-mat/0011174]

    2001 arXiv

  9. [16]

    Bennett, B

    E. Bennett, B. Lucini, D. Mason, M. Piai, E. Rinaldi and D. Vadacchino,The density of states method for symplectic gauge theories at finite temperature, 2409.19426. 9 The confined-deconfined surface tension in SU(N) gauge theories at large N Ahmed Salami

    arXiv

  10. [17]

    Caselle, G

    M. Caselle, G. Costagliola, A. Nada, M. Panero and A. Toniato,Jarzynski’s theorem for lattice gauge theory,Phys. Rev. D94(2016) 034503 [1604.05544]

    2016 arXiv

  11. [18]

    Kajantie, M

    K. Kajantie, M. Laine, K. Rummukainen and M.E. Shaposhnikov,The Electroweak phase transition: A Nonperturbative analysis, Nucl. Phys. B466 (1996) 189 [hep-lat/9510020]

    1996 arXiv

  12. [19]

    Moore and K

    G.D. Moore and K. Rummukainen,Electroweak bubble nucleation, nonperturbatively, Phys. Rev. D63(2001) 045002 [hep-ph/0009132]

    2001 arXiv

  13. [21]

    Giusti, M

    L. Giusti, M. Hirasawa, M. Pepe and L. Virzì,A precise study of the SU(3) Yang-Mills theory across the deconfinement transition, 2501.10284. 10

Pith tools

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