pith. sign in

arxiv: 2507.00689 · v3 · submitted 2025-07-01 · ✦ hep-lat

Polyakov loop model with exact static quark determinant in the 't Hooft-Veneziano limit: U(N) case

S. Voloshyn This is my paper

Pith reviewed 2026-05-19 06:46 UTC · model grok-4.3

classification ✦ hep-lat
keywords Polyakov loop modelunitary matrix model't Hooft-Veneziano limitphase diagramquark condensatedeconfinement transitionstatic quark determinant
0
0 comments X

The pith

In the large N and N_f limit with fixed ratio the Polyakov loop model reduces to an exactly solvable deformed unitary matrix model.

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

The paper studies a d-dimensional U(N) Polyakov loop model that incorporates the exact static determinant for N_f degenerate quark flavors, with explicit dependence on quark mass and chemical potential. In the 't Hooft-Veneziano limit where both N and N_f become large at fixed ratio kappa equal to N_f over N, mean-field theory becomes exact. This reduces the core of the model to a deformed unitary matrix model that admits an exact solution. The solution supplies the free energy, the expectation value of the Polyakov loop, and the quark condensate, and determines how the phase diagram and transition type vary with kappa.

Core claim

In the 't Hooft-Veneziano limit, mean field gives the exact solution for the Polyakov loop model with the exact static quark determinant, reducing its core to a deformed unitary matrix model which is solved exactly to compute the free energy, the expectation value of the Polyakov loop, and the quark condensate. The phase diagram and the type of phase transition depend on the ratio kappa equal to N_f over N.

What carries the argument

Deformed unitary matrix model obtained via exact mean-field reduction of the Polyakov loop model in the large N and N_f limit with fixed kappa.

If this is right

  • The free energy is obtained exactly as a function of temperature, quark mass, and chemical potential.
  • The expectation value of the Polyakov loop is computed precisely and serves as an order parameter for deconfinement.
  • The quark condensate is obtained exactly and tracks chiral symmetry restoration.
  • The nature of the phase transition changes with the value of kappa equal to N_f over N.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The exact solvability offers a benchmark for large-N lattice simulations of related QCD models.
  • The kappa dependence may guide effective descriptions when N and N_f are finite but not asymptotically large.
  • Similar reductions could be tested in SU(N) versions of the model.

Load-bearing premise

The mean-field approximation becomes exact in the 't Hooft-Veneziano limit of large N and N_f with fixed ratio.

What would settle it

Monte Carlo simulations of the full Polyakov loop model at successively larger but finite N and N_f with fixed kappa, checking convergence of the free energy and order parameters to the analytic mean-field predictions.

Figures

Figures reproduced from arXiv: 2507.00689 by S. Voloshyn.

Figure 1
Figure 1. Figure 1: Free energies og PL model as function of [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Quark condensate vs. m κ = 1 (b=0, 0.5, 0.66, 0.85, 1) (left panel) and the average PL in in h-b coordinates at κ = 1 (b=0, 0.36, 0.52, 0.8, 1) (left panel) . Colors counted from top to bottom. 4 Summary and Perspectives U(N) PL model that takes into account exact static quark determinant with Nf degenerate flavors of quarks investigated in the mean-field approximation. First, we solved the deformed unitar… view at source ↗
read the original abstract

I investigate a $d$-dimensional $U(N)$ Polyakov loop model that includes the exact static determinant with $N_f$ degenerate quark flavor and depends explicitly on the quark mass and chemical potential. In the large $N, N_f$ limit mean field gives the exact solution, and the core of the Polyakov loop model is reduced to a deformed unitary matrix model, which I solve exactly. I compute the free energy, the expectation value of the Polyakov loop, and the quark condensate. The phase diagram of the model and the type of phase transition is investigated and shows it depends on the ratio $\kappa =N_f/N$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript investigates a d-dimensional U(N) Polyakov loop model incorporating the exact static quark determinant for N_f degenerate flavors, with explicit dependence on quark mass and chemical potential. In the 't Hooft-Veneziano limit (large N and N_f at fixed ratio κ = N_f/N), the author claims that the mean-field approximation becomes exact, reducing the model to a deformed unitary matrix model that is solved exactly. Analytic expressions are derived for the free energy, the expectation value of the Polyakov loop, and the quark condensate; the phase diagram and the nature of the phase transitions are analyzed as functions of κ.

Significance. If the central reduction and exact solvability hold, the work supplies an analytically tractable effective model for the interplay between deconfinement and chiral symmetry breaking in the presence of static quarks. The κ-dependent phase structure and the explicit mass/chemical-potential dependence constitute a concrete advance over standard Polyakov-loop models, provided the mean-field exactness is placed on a firm footing.

major comments (2)
  1. [Section introducing the 't Hooft-Veneziano limit and mean-field exactness] The assertion that mean field becomes exact in the 't Hooft-Veneziano limit is load-bearing for the entire reduction to a solvable matrix model. The manuscript invokes the standard large-N argument but does not supply an explicit saddle-point stability analysis or 1/N expansion that isolates the fluctuation spectrum induced by the mass- and μ-dependent determinant (see the derivation of the effective action in the section introducing the 't Hooft-Veneziano limit). Without this, it remains unclear whether determinant-induced fluctuations are suppressed at leading order.
  2. [Derivation of the deformed unitary matrix model] The deformation of the unitary matrix model arising from the exact static determinant must be derived explicitly and shown to be consistent with the original Polyakov-loop partition function. The saddle-point equations obtained after the reduction should be cross-checked against the stationarity conditions of the full model to confirm that no additional terms survive at leading order in 1/N.
minor comments (2)
  1. [Notation and model definition] Clarify the precise definition of the deformation term in the matrix-model action and ensure that all parameters (including the ratio κ) are introduced with consistent notation from the outset.
  2. [Introduction] Add a brief comparison with existing large-N treatments of Polyakov-loop models that include dynamical quarks to highlight the novelty of the exact-determinant treatment.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The points raised concerning the justification of mean-field exactness and the explicit derivation of the reduced model are helpful for improving clarity. We address each comment below.

read point-by-point responses

  1. Referee: [Section introducing the 't Hooft-Veneziano limit and mean-field exactness] The assertion that mean field becomes exact in the 't Hooft-Veneziano limit is load-bearing for the entire reduction to a solvable matrix model. The manuscript invokes the standard large-N argument but does not supply an explicit saddle-point stability analysis or 1/N expansion that isolates the fluctuation spectrum induced by the mass- and μ-dependent determinant (see the derivation of the effective action in the section introducing the 't Hooft-Veneziano limit). Without this, it remains unclear whether determinant-induced fluctuations are suppressed at leading order.

    Authors: We agree that an explicit discussion of saddle-point stability strengthens the argument. In the 't Hooft-Veneziano limit with fixed κ = N_f/N, the effective action (including the exact static determinant) is O(N) and the partition function is dominated by the saddle point of the unitary matrix model. Fluctuations around this saddle are suppressed by 1/N in the standard manner for large-N matrix models. The mass- and μ-dependent determinant contributes at leading order but does not alter the suppression of quadratic fluctuations. In the revised manuscript we will add a short appendix or subsection providing the second-variation analysis of the effective action, confirming that the Hessian remains positive definite at leading order with determinant-induced corrections subleading in 1/N. revision: yes

  2. Referee: [Derivation of the deformed unitary matrix model] The deformation of the unitary matrix model arising from the exact static determinant must be derived explicitly and shown to be consistent with the original Polyakov-loop partition function. The saddle-point equations obtained after the reduction should be cross-checked against the stationarity conditions of the full model to confirm that no additional terms survive at leading order in 1/N.

    Authors: We concur that the reduction step benefits from a more explicit presentation. The deformed unitary matrix model is obtained by first writing the full partition function with the exact static quark determinant and then taking the large-N, large-N_f limit at fixed κ; the logarithm of the determinant yields an effective potential for the eigenvalues of the Polyakov-loop matrix. This procedure is consistent with the original model because the determinant is retained exactly before the saddle-point approximation is applied. In the revision we will expand the relevant section to display the intermediate steps and explicitly verify that the saddle-point equations of the reduced model reproduce the stationarity conditions of the full effective action at order N, with no additional terms appearing at this leading order. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard large-N mean-field exactness invoked as external input

full rationale

The paper asserts that in the large N, N_f limit with fixed κ = N_f/N, mean field becomes exact, reducing the Polyakov loop model with exact static quark determinant to a solvable deformed unitary matrix model. This relies on the established large-N mean-field exactness argument common in the literature for matrix models and gauge theories in the 't Hooft-Veneziano limit, rather than any self-referential definition, fitted parameter renamed as prediction, or load-bearing self-citation chain within the paper. The subsequent exact solution for free energy, Polyakov loop vev, and condensate follows directly from this standard approximation without reducing the target quantities to the inputs by construction. No evidence of ansatz smuggling, uniqueness imported from authors, or renaming of known results appears in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that mean-field theory is exact in the large-N, large-N_f limit with fixed ratio; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Mean-field approximation becomes exact in the 't Hooft-Veneziano limit for the Polyakov loop model with exact static quark determinant
    Invoked in the abstract as the step that reduces the model to a solvable deformed unitary matrix model.

pith-pipeline@v0.9.0 · 5636 in / 1356 out tokens · 66089 ms · 2026-05-19T06:46:37.271262+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages · 3 internal anchors

  1. [1]

    Philipsen, PoS LATTICE2019 (2019) 273 [arXiv:1912.04827 [hep-lat]], O

    O. Philipsen, PoS LATTICE2019 (2019) 273 [arXiv:1912.04827 [hep-lat]], O. Philipsen, J. Scheunert, JHEP 11 (2019) 022 [arXiv:1908.03136 [hep-lat]]

    work page arXiv 2019

  2. [2]

    t’ Hooft, Nucl.Phys

    G. t’ Hooft, Nucl.Phys. B 72 (1974) 461

    work page 1974

  3. [3]

    Veneziano, Nucl.Phys

    G. Veneziano, Nucl.Phys. B 117 (1976) 519

    work page 1976

  4. [4]

    D. J. Gross, E. Witten, Phys.Rev. D 21 (1980) 446

    work page 1980

  5. [5]

    S. R. Wadia, Phys.Lett. B 93 (1980) 403

    work page 1980

  6. [6]

    P. H. Damgaard and A. Patk´ os, Phys.Lett. B172 (1986) 369

    work page 1986

  7. [7]

    C. H. Christensen, Phys.Lett. B 714 (2012) 306 [arXiv:1204.2466 [hep-lat]]

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  8. [8]

    Finite-temperature phase transitions of third and higher order in gauge theories at large $N$

    H. Nishimura,, R. D. Pisarski, V. V. Skokov, Phys.Rev. D 97 (2018) 036014 [arXiv:1712.04465 [hep-th]]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  9. [9]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, PoS LATTICE 2021 (2021) 453 [arXiv:2111.07103 [hep-lat]] 15

    work page arXiv 2021

  10. [10]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, Phys.Rev. D 105 (2022) 014501 [arXiv:2111.00474 [hep-lat]]

    work page arXiv 2022

  11. [11]

    Russo, Phases of unitary matrix models and lattice QCD2, arXiv:2010.02950v1

    Jorge G. Russo, Phases of unitary matrix models and lattice QCD2, arXiv:2010.02950v1

    work page arXiv 2010

  12. [12]

    Russo and Miguel Tierz, Multiple phases in a generalized Gross- Witten-Wadia, arXiv:2007.08515v1 [hep-th]

    Jorge G. Russo and Miguel Tierz, Multiple phases in a generalized Gross- Witten-Wadia, arXiv:2007.08515v1 [hep-th]

    work page arXiv 2007

  13. [13]

    Santilli, M

    L. Santilli, M. Tierz, Exact equivalences and phase discrepancies between ran- dom matrix ensembles, arXiv:2003.10475v2 [math-ph]

    work page arXiv 2003

  14. [14]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, P. Yefanov, One-dimensional QCD at finite density and its ’t Hooft-Veneziano limit. JHEP 01 (2025) 008, [arXiv:2410.02328 [hep-lat]]

    work page arXiv 2025

  15. [15]

    Borisenko, S

    O. Borisenko, S. Voloshyn, V. Chelnokov, Rep. Math. Phys. 85 (2020) 129 [arXiv:1812.06069 [hep-lat]]

    work page arXiv 2020

  16. [16]

    Goldschmidt, 1/ N expansion in two-dimensional lattice gauge theory, J

    Y.Y. Goldschmidt, 1/ N expansion in two-dimensional lattice gauge theory, J. Math. Phys. 21 (1980) 1842, DOI: 10.1063/1.524600

    work page doi:10.1063/1.524600 1980

  17. [17]

    Duals of U(N) LGT with staggered fermions

    O. Borisenko, V. Chelnokov, S. Voloshyn, EPJ Web Conf. 175 (2018) 11021 [arXiv:1712.03064 [hep-lat]]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  18. [18]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, Phys.Rev. D 102 (2020) 014502 [arXiv:2005.11073 [hep-lat]]

    work page arXiv 2020

  19. [19]

    Borisenko, V

    O. Borisenko, V. Chelnokov, E. Mendicelli, A. Papa, Nucl.Phys.B 965 (2021) 115332 [arXiv:2011.08285 [hep-lat]]

    work page arXiv 2021

  20. [20]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, P. Yefanov, Phys.Lett. B 827 (2022) 137000 [arXiv:2112.06002 [hep-lat]]

    work page arXiv 2022

  21. [21]

    Borisenko, V

    O. Borisenko, V. Chelnokov, S. Voloshyn, Nucl.Phys. B960 (2020) 115177 [arXiv:2008.00773 [hep-lat]]. 16

    work page arXiv 2020