REVIEW 2 major objections 4 minor 31 references
Finite-size scaling of Lee-Yang zeros and its application to the 3-state Potts model and heavy-quark QCD
T0 review · 2 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This paper proposes the Lee-Yang-zero ratio method, a finite-size scaling technique that locates critical points from the intersection of ratio curves, and demonstrates it in the 3D Ising model, the 3-state Potts model, and heavy-quark QCD.
desk verdict New crossing observable for finding critical points; the Potts/Ising checks are solid, the QCD central value has an unquantified finite-size correction. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The Lee-Yang-zero ratio (LYZR) is the ratio $R_{nm}(t,L)=h^{(n)}_{LY}(t,L)/h^{(m)}_{LY}(t,L)$ of the imaginary parts of the $n$-th and $m$-th Lee-Yang zeros, the zeros of the partition function in the complex external-field plane. Its finite-size scaling form, $R_{nm}=r_{nm}+c_{nm}L^{y_t}t+\cdots$ with $r_{nm}=X_n/X_m$, is the machinery that carries the argument: it converts the unknown critical point into a volume-independent crossing, exactly the role Binder cumulants play, while the subleading correction $D_{nm}L^{2\bar y}$ is smaller than the Binder counterpart. For systems not written in Ising variables, the partition function is related to the Ising one by a linear coupling transformation $A$ (Eq. (12)); the LYZR then has the same crossing form up to a correction factor that vanishes as $L\to\infty$. The numerical implementation uses reweighting to evaluate the partition function at complex couplings and locate the zeros.
What would settle it
Compute $R_{21}(t,L)$ for the two-dimensional Ising model on lattices of size $L=16,32,64$ using exact transfer-matrix or Monte Carlo reweighting data; if the curves do not intersect at the known $T_c$ (or the intersection drifts with $L$ beyond the predicted $L^{2\bar y}$ correction), the claimed general crossing property is falsified.
Extended reading notes
Core claim
The central claim is that the ratios of Lee-Yang zeros, $R_{nm}=h^{(n)}_{LY}/h^{(m)}_{LY}$, are a critical-point diagnostic with the same crossing structure as Binder cumulants. Starting from the finite-size scaling form $Z(t,h,L^{-1})=\tilde Z(L^{y_t}t, L^{y_h}h)$ and the regularity of the scaled zero functions $\tilde h^{(n)}_{LY}(\tilde t)$, the paper Taylor-expands $\tilde h^{(n)}_{LY}$ at $\tilde t=0$ and obtains $R_{nm}(t,L)=r_{nm}+c_{nm}L^{y_t}t+O(t^2)$, with $r_{nm}=X_n/X_m$ volume-independent. Hence $R_{nm}(t,L)$ curves for different $L$ cross at the critical point $t=0$, and the slope grows as $L^{y_t}$. The paper verifies this in 3D Ising Monte Carlo data (crossing $T\simeq 4.5114$, consistent with the precise $T_c=4.51152322(1)$), then applies the method to the 3-state Potts model, obtaining $\tau_c=0.549375(18)$ from $R_{21}$, and to heavy-quark QCD at $N_t=4$, obtaining $\beta^*_c=5.68595(30)$ from $R_{21}$, both consistent with Binder-cumulant estimates. In general systems the method requires a linear mapping between the physical couplings and the Ising variables $t,h$; this introduces a subleading correction factor $1+D_{nm}L^{2\bar y}$ with $\bar y=y_t-y_h<0$, which the paper notes decays faster than the corresponding Binder correction.
Load-bearing premise
The load-bearing assumption is that the scaled Lee-Yang zero positions $\tilde h^{(n)}_{LY}(\tilde t)$ are analytic functions through the critical point, so the Taylor expansion in Eq. (8) with a nonvanishing leading coefficient $X_n$ is valid; if a zero used in a ratio has $X_n=0$ or the expansion fails, the volume-independent crossing does not occur.
Editorial extensions
If this is right
- In any model in the Ising universality class, a handful of lattice sizes and a reweighting-based zero determination should be enough to read off the critical point from the crossing of $R_{21}$ or $R_{31}$ curves.
- The consistency between $R_{21}$, $R_{31}$, and Binder crossings in the Potts model supports using LYZR as a cross-check or replacement for Binder cumulants, with a smaller leading finite-size correction.
- The heavy-quark QCD result pins the deconfinement critical point at $N_t=4$ to $\beta^*_c \approx 5.686$ with an error comparable to the Binder analysis, so LYZR is usable directly on existing lattice configurations.
- Because the ratios approach 1 for $t>0$ and $(2n-1)/(2m-1)$ for $t<0$, the method also gives a diagnostic of which side of the critical point a simulation is on, not just the location of the crossing.
Reading between the lines
- Beyond the paper, the same crossing logic should apply to Lee-Yang zeros in the complex chemical-potential plane, where reweighting from imaginary $\mu_B$ is sign-problem-free; if so, LYZR could bound the QCD critical point without a real-$\mu_B$ simulation.
- Beyond the paper, the slope $c_{nm}L^{y_t}$ of the ratio curves could be used to extract the thermal exponent $y_t$ from the same fits that locate the crossing, giving a self-contained measurement of critical exponents.
- Beyond the paper, since the method needs only the locations of the leading zeros, it may be more robust to statistical noise in the tail of the reweighted partition function than methods that use the full distribution, a testable advantage on smaller lattices.
- Beyond the paper, applying LYZR to the 2D Ising model, where exact transfer-matrix zeros are available, would provide a clean benchmark of the crossing property and the predicted $L^{2\bar y}$ correction.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript proposes locating critical points via ratios of Lee-Yang zeros (LYZR). Starting from the finite-size scaling form of the partition function, Eq. (2), and a Taylor expansion of the scaled zeros, the paper derives Eq. (9), R_nm(t,L) = r_nm + c_nm L^{y_t} t + O(t^2), so that curves for different L cross at the critical temperature. The method is verified in the 3D Ising model by reweighting, then extended to general Ising-universality systems through the linear mapping Eq. (12), leading to Eq. (16) with a correction factor [1 + D_nm L^{2\bar y}]. The paper applies the method to the 3-state Potts model and to N_f=2 heavy-quark QCD on N_t=4 lattices, reporting tau_c = 0.549375(18)/0.549373(17) and beta*_c = 5.68595(30)/5.68603(31), which are compared with Binder-cumulant results.
Significance. If the LYZR construction is sound, it adds a useful observable to the critical-point-location toolbox: ratios of Lee-Yang zeros have Binder-like crossings, the leading finite-size correction in Eq. (16) decays faster than the corresponding Binder correction, and the method combines naturally with reweighting in theories where complex couplings can be simulated. The paper's strengths are that the scaling argument is concise and testable, and the numerical tests are anchored to independent results: the Ising crossing is compared with the high-precision T_c of Ref. [17], and the Potts and QCD crossings are compared with Binder cumulant analyses. The main weaknesses are that the headline QCD fit values rest on an unquantified neglected correction, and the analyticity assumption behind the Taylor expansion is not justified; both are addressable within the scope of the manuscript.
major comments (2)
- [Sec. 3.3, Eq. (16)] The quoted values beta*_c = 5.68595(30) and 5.68603(31) are obtained by fitting R_n1(beta*,L_T) = r + c(beta* - beta*_c)L_T^{y_t} and explicitly neglecting the second factor [1 + D_nm L_T^{2\bar y} + ...] in Eq. (16). This is not justified at the volumes used: with L_T = 9, 10, 12 and \bar y = y_t - y_h \simeq -0.894, L_T^{2\bar y} is approximately 0.020, 0.016, and 0.012, so the factor varies by about 8 x 10^{-3} across the fitted volumes. Since dR/dbeta* is approximately C L_T^{y_t} with L_T^{y_t} \simeq 30, a coefficient D_nm = O(1) changes the crossing by roughly 2-3 x 10^{-4}, i.e., the size of the quoted statistical errors; the manuscript contains no estimate of D_nm and no QCD check analogous to the 'no clear deviation for L >= 40' statement used in Sec. 3.2. The agreement with the Binder value from Ref. [29] is reassuring but uses the same reweighted configurations and the same action, so it does not by itself bound this finite-size systematic. Please include the correction term in the fit or provide an independent estimate of its size before quoting these CP values.
- [Sec. 2.1, Eq. (8)] The derivation of the crossing property assumes that the scaled zero function \tilde h^{(n)}_{LY}(\tilde t) is analytic at \tilde t = 0 and admits the Taylor expansion of Eq. (8). The text justifies this by the regularity of the finite-L zeros h^{(n)}_{LY}(t,L) in t, but regularity of finite-L zeros does not automatically survive the scaling limit that defines \tilde h^{(n)} in Eq. (3); zeros can collide or the limiting function can develop non-analytic behavior. If \tilde h^{(n)} contains a term such as |\tilde t|^p with 0 < p < 1, Eq. (9) and the claimed volume-independent intersection would fail. The Ising data in Fig. 3 provide empirical support, but the paper should state the analyticity assumption explicitly and add a quantitative check, for example fitting the measured slope dR_n1/dT at the crossing to the predicted L^{y_t} behavior.
minor comments (4)
- [Sec. 3.1, Eqs. (14)-(15)] Please clarify whether the O(L^{2\bar y}) terms in Eqs. (14) and (15) are meant as absolute corrections or as relative corrections. As written, the O(L^{2\bar y}) term in Eq. (15) appears to dominate the displayed leading term L^{-y_h} for large L, which cannot be intended; rewriting these equations with an explicit bracket factor would remove the ambiguity.
- [Sec. 3.3] The fit in Sec. 3.3 is described as a 'three-parameter FSS fit', but four parameters r, c, beta*_c, and y_t are listed; please correct this wording and state whether y_t is fixed to the Ising value or left free.
- [Sec. 2.2] The claim that the crossing point in the Ising model is at T \simeq 4.5114 is made by visual inspection of Fig. 3; please provide a numerical estimate of the crossing point with its statistical error so that the consistency with T_c = 4.51152322(1) can be judged quantitatively.
- [Sec. 2.2, Fig. 2] The data collapse in the right panel of Fig. 2 is shown qualitatively; please specify the values of y_t and y_h used in the rescaling and, if possible, quantify the quality of the collapse.
Circularity Check
No significant circularity: the LYZR crossing is derived from finite-size scaling, and the Potts/QCD critical couplings are fit outputs, not fitted inputs; self-citations are not load-bearing.
full rationale
The paper's derivation chain is self-contained. Eq. (2) is the standard FSS form, Eq. (8) is an explicit analyticity assumption, and Eq. (16) follows by algebra; the crossing property R_nm(tau_c,L)=r_nm is a consequence of the scaling form, not an input. In the Ising test, the crossing T~4.5114 is checked against the independent high-precision Tc=4.51152322(1) from Ref. [17]. For the Potts model and heavy-quark QCD, tau_c and beta*_c are free parameters of least-squares fits to the LYZR data with the derived ansatz; they are not obtained by tuning to the Binder values quoted for comparison. The Binder comparisons in Secs. 3.2 and 3.3 use the same configurations (Refs. [1] and [29]), so they are consistency checks rather than independent external benchmarks, but they are not used to define the LYZR prediction and thus do not make the derivation circular. Refs. [1], [24], [29], and [30] are self-citations, but no load-bearing uniqueness theorem or ansatz is imported from them; Eq. (16) is derived in the text. The neglect of the D_nm L^{2ybar} correction in the finite-L fits is an unquantified finite-size systematic, especially for L_T=9 in QCD, but that is a correctness risk, not circularity, because the extracted critical couplings are not by construction equal to any input parameter. Overall, no circular step is present; the score reflects only minor, non-load-bearing self-citations.
Assumptions & free parameters
free parameters (5)
- r_nm (universal zero ratio at CP) =
not reported
- c (slope coefficient) =
not reported
- y_t (thermal exponent) =
fitted; Ising value 1.588 used for scaling collapse
- tau_c for Potts model =
0.549375(18) from R21, 0.549373(17) from R31
- beta*_c for heavy-quark QCD =
5.68595(30) from R21, 5.68603(31) from R31
assumptions (8)
- domain assumption Finite-size scaling hypothesis: Z(t,h,L^{-1}) = Ztilde(L^{y_t}t,L^{y_h}h) near the critical point (Eq. 2).
- standard math Lee-Yang zeros of the Z2-symmetric Ising model lie on the imaginary h-axis for real t.
- domain assumption For t<0, LYZ are equally spaced as h_n = a(t)(2n-1)/L^3 (Eq. 4).
- domain assumption For t>0, LYZ converge to a finite Lee-Yang edge singularity (Eq. 5).
- standard math Scaled zero functions htilde_n(ttilde) are analytic at ttilde=0 with Taylor expansion Eq. (8).
- domain assumption Pairwise linear mapping A between (tau,xi) and Ising scaling fields (Eqs. 11-12).
- domain assumption Hopping parameter expansion through NLO is accurate near the heavy-quark QCD critical point.
- domain assumption Reweighting overlap is sufficiently controlled in the simulations.
Cite this review
Pith. "Pith review of Finite-size scaling of Lee-Yang zeros and its application to the 3-state Potts model and heavy-quark QCD." pith.science (2026). https://pith.science/paper/FWOHN7ST
@misc{pith2026250118904,
author = {Pith},
title = {Pith review of: Finite-size scaling of Lee-Yang zeros and its application to the 3-state Potts model and heavy-quark QCD},
year = {2026},
howpublished = {\url{https://pith.science/paper/FWOHN7ST}},
note = {Machine review of arXiv:2501.18904}
}
read the original abstract
We propose a new general method to study critical points (CP) using the finite-size scaling of Lee-Yang zeros (LYZ). We first study the LYZ in the three-dimensional Ising model on finite lattices. We show that the ratios of multiple LYZ (Lee-Yang-zero ratios: LYZR) have useful scaling properties similar to the Binder cumulants, providing us with a novel method to study CP. In numerical simulations of the Ising model, we confirm that this method works well. We then apply the method to analyze the CP in the three-dimensional three-state Potts model and finite-temperature QCD in heavy-quark region, which are believed to belong to the same universality class as the Ising model. In these models, the partition function at complex parameters can be evaluated by the reweighting method, which allows us to determine the LYZ by varying coupling parameters continuously around the CP. We demonstrate that the LYZR method is powerful in determining the location of the CP in these models.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[29]
A. Kiyohara, M. Kitazawa, S. Ejiri and K. Kanaya,Finite-size scaling around the critical point in the heavy quark region of QCD,Phys. Rev. D104(2021) 114509 [2108.00118]
arXiv 2021
-
[1]
T. Wada, M. Kitazawa and K. Kanaya,Lee-Yang-zero ratios for locating a critical point, 2410.19345
-
[17]
A.M. Ferrenberg, J. Xu and D.P. Landau,Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model,Phys. Rev. E97(2018) 043301 [1806.03558]
arXiv 2018
-
[2]
P.J. Kortman and R.B. Griffiths,Density of Zeros on the Lee-Yang Circle for Two Ising Ferromagnets,Phys. Rev. Lett.27(1971) 1439
work page 1971
-
[3]
Stephanov,QCD critical point and complex chemical potential singularities, Phys
M.A. Stephanov,QCD critical point and complex chemical potential singularities, Phys. Rev. D73(2006) 094508 [hep-lat/0603014]
arXiv 2006
-
[4]
S.Ejiri,Y.ShinnoandH.Yoneyama, ComplexsingularitiesaroundtheQCDcriticalpointat finite densities,PTEP2014 (2014) 083B02 [1404.6004]
arXiv 2014
-
[5]
X. An, D. Mesterházy and M.A. Stephanov,On spinodal points and Lee-Yang edge singularities, J. Stat. Mech.1803(2018) 033207 [1707.06447]
arXiv 2018
-
[6]
G. Basar, G.V. Dunne and Z. Yin,Uniformizing Lee-Yang singularities,Phys. Rev. D105 (2022) 105002 [2112.14269]
arXiv 2022
Show all 31 references
-
[7]
Rennecke and V.V
F. Rennecke and V.V. Skokov,Universal location of Yang–Lee edge singularity for a one-component field theory in 1≤d≤4, Annals Phys.444 (2022) 169010 [2203.16651]
2022 arXiv
-
[8]
G.Johnson,F.RenneckeandV.V.Skokov, UniversallocationofYang-Leeedgesingularityin classic O(N) universality classes,Phys. Rev. D107(2023) 116013 [2211.00710]
2023 arXiv
-
[9]
Singh, M
S. Singh, M. Cipressi and F. Di Renzo,Exploring Lee-Yang and Fisher zeros in the 2D Ising model through multipoint Padé approximants, Phys. Rev. D109 (2024) 074505 [2312.03178]
2024 arXiv
-
[10]
Karsch, C
F. Karsch, C. Schmidt and S. Singh,Lee-Yang and Langer edge singularities from analytic continuation of scaling functions, Phys. Rev. D109 (2024) 014508 [2311.13530]
2024 arXiv
-
[11]
P.Dimopoulos, L.Dini, F.DiRenzo, J.Goswami,G.Nicotra, C.Schmidtetal., Contribution to understanding the phase structure of strong interaction matter: Lee-Yang edge singularities from lattice QCD, Phys. Rev. D105 (2022) 034513 [2110.15933]
2022 arXiv
-
[12]
Basar,On the QCD critical point, Lee-Yang edge singularities and Pade resummations, 2312.06952
G. Basar,On the QCD critical point, Lee-Yang edge singularities and Pade resummations, 2312.06952
-
[13]
Zambello, D.A
K. Zambello, D.A. Clarke, P. Dimopoulos, F. Di Renzo, J. Goswami, G. Nicotra et al., Determination of Lee-Yang edge singularities in QCD by rational approximations,PoS LATTICE2022(2023) 164 [2301.03952]
2023 arXiv
-
[14]
Clarke, P
D.A. Clarke, P. Dimopoulos, F. Di Renzo, J. Goswami, C. Schmidt, S. Singh et al.,Searching for the QCD critical endpoint using multi-point Padé approximations, 2405.10196
-
[15]
Skokov,Two lectures on Yang-Lee edge singularity and analytic structure of QCD equation of state, 2411.02663
V.V. Skokov,Two lectures on Yang-Lee edge singularity and analytic structure of QCD equation of state, 2411.02663. 9 Finite-size scaling of Lee-Yang zeros – application to Potts model and heavy-quark QCDTatsuya Wada
-
[16]
Binder,Finite size scaling analysis of Ising model block distribution functions,Z
K. Binder,Finite size scaling analysis of Ising model block distribution functions,Z. Phys. B 43 (1981) 119
1981
-
[18]
Lee and C.-N
T.D. Lee and C.-N. Yang,Statistical theory of equations of state and phase transitions. 2. Lattice gas and Ising model,Phys. Rev.87(1952) 410
1952
-
[19]
Ejiri,Lee-Yang zero analysis for the study of QCD phase structure,Phys
S. Ejiri,Lee-Yang zero analysis for the study of QCD phase structure,Phys. Rev. D73(2006) 054502 [hep-lat/0506023]
2006 arXiv
-
[20]
I. Bena, M. Droz and A. LIPOWSKI,Statistical mechanics of equilibrium and nonequilibrium phase transitions: The yang–lee formalism,International Journal of Modern Physics B19 (2005) 4269–4329
2005
-
[21]
Wolff,Collective Monte Carlo Updating for Spin Systems,Phys
U. Wolff,Collective Monte Carlo Updating for Spin Systems,Phys. Rev. Lett.62(1989) 361
1989
-
[22]
Ferrenberg and R.H
A.M. Ferrenberg and R.H. Swendsen,New monte carlo technique for studying phase transitions, Phys. Rev. Lett.61(1988) 2635
1988
-
[23]
F.Kos, D.Poland, D.Simmons-DuffinandA.Vichi, PrecisionIslandsintheIsingand 𝑂(𝑁) Models,JHEP08 (2016) 036 [1603.04436]
2016 arXiv
-
[24]
T. Wada, M. Kitazawa and K. Kanaya,in preparation,
-
[25]
Rehr and N.D
J.J. Rehr and N.D. Mermin,Revised Scaling Equation of State at the Liquid-Vapor Critical Point,Phys. Rev. A8 (1973) 472
1973
-
[26]
Wilding,Simulation studies of fluid critical behaviour,J
N.B. Wilding,Simulation studies of fluid critical behaviour,J. Phys.: Condens. Matter9 (1997) 585
1997
-
[27]
Karsch and S
F. Karsch and S. Stickan,The Three-dimensional, three state Potts model in an external field, Phys. Lett. B488 (2000) 319 [hep-lat/0007019]
2000 arXiv
-
[28]
Cuteri, O
F. Cuteri, O. Philipsen, A. Schön and A. Sciarra,Deconfinement critical point of lattice QCD with𝑁𝑓=2 Wilson fermions,Phys. Rev. D103 (2021) 014513 [2009.14033]
2021 arXiv
-
[30]
Ashikawa, M
R. Ashikawa, M. Kitazawa, S. Ejiri and K. Kanaya,High-precision analysis of the critical point in heavy-quark QCD at Nt=6, Phys. Rev. D110 (2024) 074508 [2407.09156]
2024 arXiv
-
[31]
Wakabayashi, S
N. Wakabayashi, S. Ejiri, K. Kanaya and M. Kitazawa,Scope and convergence of the hopping parameter expansion in finite-temperature quantum chromodynamics with heavy quarks around the critical point,PTEP2022 (2022) 033B05 [2112.06340]. 10
2022 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.