REVIEW 3 major objections 5 minor 1 cited by
The QCD phase diagram, universal scaling, and Lee-Yang zeros
T0 review · 3 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This paper claims that the QCD critical endpoint lies at $(T_{\mathrm{cep}}, \mu_{\mathrm{cep}}) = (105^{+8}_{-18}, 422^{+80}_{-35})$ MeV, extracted from universal Lee-Yang scaling of multipoint Pad\'e approximants to…
desk verdict A competent proceedings summary of the Bielefeld-Parma chiral scaling and Lee-Yang endpoint work; no new results, but the endpoint claim rests on an untested linear mixing ansatz. 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 load-bearing machinery is the multi-point Pad\'e approximant $R_n^n(\hat\mu_B)$ to the first cumulant of the baryon density, constructed by matching the function and its derivative at imaginary simulation points in the range $\hat\mu_B \in [0, i\pi]$; its closest pole in the complex chemical-potential plane is identified with the Lee-Yang edge singularity. The temperature flow of that pole is then mapped to the universal scaling of the Lee-Yang edge, $z_{LY} = t/h^{1/(\beta\delta)}$, through the linear mixing ansatz for the scaling fields $t$ and $h$ (Eq. 7). In the chiral sector, the key object is the improved order parameter $M = M_l - H\chi_l$, which removes additive UV divergences and the leading regular $H$ contribution; its scaling is described by O(2) scaling functions, and the pseudo-critical temperatures are obtained from maxima of five different susceptibilities.
What would settle it
A multipoint Pad\'e analysis on $N_\tau = 8$ lattices that gives an endpoint outside the quoted $N_\tau = 6$ uncertainty would show the result is a cutoff artefact; a direct check would also test whether the closest pole obeys $\Im\mu_{LY} \sim (T-T_{\mathrm{cep}})^{\beta\delta}$ with the O(N) value of $\beta\delta$ as $T$ approaches $T_{\mathrm{cep}}$ from above.
Extended reading notes
Core claim
The central quantitative claim of the paper is that the QCD critical endpoint is located at $(T_{\mathrm{cep}}, \mu_{\mathrm{cep}}) = (105^{+8}_{-18}, 422^{+80}_{-35})$ MeV, determined on $N_\tau = 6$ lattices by fitting the universal scaling of Lee-Yang edge singularities. The edges are obtained as the closest poles of multi-point Pad\'e approximants to the baryon number density computed at imaginary chemical potentials $\mu_B = i\theta_B$; the scaling fit uses the linear mixing ansatz $t = \alpha_t(T-T_{\mathrm{cep}}) + \beta_t(\mu_B-\mu_{\mathrm{cep}})$ and $h = \alpha_h(T-T_{\mathrm{cep}}) + \beta_h(\mu_B-\mu_{\mathrm{cep}})$, together with the universal location of the Lee-Yang edge in the scaling variable $z = t/h^{1/(\beta\delta)}$. The paper also reports a separate chiral-limit determination $T_c = 143.9(5)$ MeV on $N_\tau = 8$ lattices from O(2) scaling fits to the improved order parameter and pseudo-critical temperatures, and curvature coefficients $\kappa_2^{B,\mu_S=0} = 0.015(1)$. The endpoint result is not continuum extrapolated, but is consistent between the multi-point Pad\'e analysis and a $[4,4]$ Pad\'e analysis of eight-order Taylor coefficients at $N_\tau = 8$.
Load-bearing premise
The endpoint extraction assumes that, in the sampled region around the critical point, the scaling fields $t$ and $h$ are strictly linear functions of $T - T_{\mathrm{cep}}$ and $\mu_B - \mu_{\mathrm{cep}}$ with four fitted coefficients, and that the closest Pad\'e pole is the Lee-Yang edge; nonlinear scaling fields or a misidentified pole would move the fitted endpoint coordinates.
Editorial extensions
If this is right
- Multi-point Pad\'e approximants to imaginary-chemical-potential data can locate Lee-Yang edges in QCD, as validated against the Roberge-Weiss transition and the 2D Ising model.
- The Lee-Yang edge scaling gives a direct handle on the endpoint; systematic errors from varying the Pad\'e order are included in the quoted coordinate uncertainties.
- The chiral scaling analysis, using the improved order parameter and pseudo-critical temperatures from five susceptibilities, yields a consistent $T_c = 143.9(5)$ MeV that agrees with the earlier continuum estimate.
- The curvature coefficients $\kappa_2^{B} = 0.015(1)$ quantify how mildly the transition temperature bends with baryon chemical potential at zero strangeness chemical potential and at zero strangeness density.
Reading between the lines
- A natural next step is a finer-lattice ($N_\tau = 8$) multipoint Pad\'e analysis of the same baryon-density data; agreement with the $N_\tau = 6$ endpoint would be the first evidence that the location survives continuum extrapolation.
- The linear mixing ansatz for the scaling fields is the softest point of the endpoint analysis; replacing it with a nonlinear ansatz or extracting the mixing angles from auxiliary observables would provide a direct test of the coordinate's stability.
- Because only the closest Pad\'e pole is used, spurious poles due to numerical noise are a risk; requiring the pole trajectory to follow the predicted $\Im\mu_{LY} \sim (T - T_{\mathrm{cep}})^{\beta\delta}$ power law over several temperatures would sharpen the identification.
- If the endpoint is real, freeze-out curves from heavy-ion collisions at moderate energies would cross the critical region at $\mu_B$ near 422 MeV, giving a concrete target for fluctuation measurements in low-energy collision programs.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings article reports two quantitative results from lattice QCD with (2+1)-flavor HISQ fermions: (i) a determination of the chiral transition temperature and curvature coefficients from O(2) scaling fits to the improved order parameter and mixed susceptibilities at N_tau=8, and (ii) a determination of the QCD critical endpoint (T_cep, mu_cep) = (105(+8,-18), 422(+80,-35)) MeV at N_tau=6 from the universal scaling of Lee-Yang edge singularities obtained via multi-point Pade approximants to the baryon number density at imaginary chemical potential. The Lee-Yang edge analysis uses a linear mixing ansatz for the Z(2) scaling fields (Eq. (7)) and identifies the closest Pade pole with the Lee-Yang edge, a strategy benchmarked in earlier work on the Roberge-Weiss transition and the 2d Ising model. The paper also quotes curvature coefficients for the chiral transition line and compares its endpoint estimate with a HotQCD [4,4] Pade analysis.
Significance. If the endpoint result holds, it is a noteworthy step toward locating the QCD critical point from first-principles lattice calculations, and the agreement with the HotQCD [4,4] Pade data in Fig. 3 is encouraging. The chiral Tc determination is consistent with previous continuum results and uses universal O(2) scaling functions imported from Ref. [11], which gives the chiral analysis an external, non-perturbative anchor. The multi-point Pade method has been validated in controlled settings, so the methodological pipeline is credible. However, the endpoint value is presented with only a companion-paper reference for the detailed analysis, and the present manuscript does not provide the fit diagnostics needed to judge the robustness of the linear scaling-field ansatz. The paper is a proceedings summary, and its main value is in drawing attention to a specific, falsifiable prediction for the critical endpoint.
major comments (3)
- [Sec. 3, Eq. (7)] The endpoint extraction rests entirely on the linear mixing ansatz t = alpha_t (T - T_cep) + beta_t (mu_B - mu_cep) and h = alpha_h (T - T_cep) + beta_h (mu_B - mu_cep), with the four mixing parameters fitted to data. The paper states that 'the scaling directions are not known a priori,' so the linear form is an assumption rather than a derived relation. The fit shown in Fig. 3 spans a wide temperature range (the horizontal axis runs from 80 to 170 MeV), while T_cep is about 105 MeV; for |T - T_cep| up to roughly 50 MeV, regular and nonlinear contributions to the scaling fields are not guaranteed to be negligible. The manuscript does not state the temperature interval used in the fit, the number of data points, or the chi^2/dof, and no stability check under a restricted fit window is reported. Because the quoted asymmetric errors (+8/-18 MeV in T, +80/-35 MeV in mu_B) do not include this systematic, the central claim is underevidenced. Please report the fit range and diagnostics and test the stability of (T_cep, mu_cep) when the fit window is narrowed around T_cep and when quadratic terms are added to Eq. (7).
- [Sec. 2 and Fig. 2] The chiral transition temperature is reported inconsistently. The Fig. 2 caption gives (T_c, z_0, h_0^{-1/delta}) = (143.8(2) MeV, 1.45(3), 39.0(3)), while the text states 'The fit yields a critical temperature of Tc = 143.9(5) MeV for Ntau=8' and later says the joint pseudo-critical fit gives 'yet another independent determination ... Tc = 143.9(5) MeV.' The figure itself contains 'z0=1.42 Tc=143.7' in the right panel. The reader cannot tell which of these numbers corresponds to which fit (order-parameter vs pseudo-critical, with or without H=1/160). Since the abstract lists the critical temperature as a reported result, these values must be reconciled and clearly attributed to the specific fits.
- [Sec. 2, H=1/160 exclusion] The order-parameter fit excludes the H=1/160 data point because of finite-size concerns, and the text later says differences are included as a systematic error, but the magnitude of that systematic is not given in the manuscript. Similarly, the curvature coefficients kappa_2^l, kappa_2^s, kappa_11^{ls} are said to be determined from ratios of mixed susceptibilities, but their numerical values are not tabulated here (only the combination kappa_B^{mu_S=0}=0.015(1) appears). Since the abstract lists the critical temperature and curvature as reported results, the manuscript should state the final values and the systematic error attached to the H=1/160 exclusion, or clearly point to the specific table in Ref. [8] that contains them.
minor comments (5)
- [General] There are several typographical errors: 'Genvea' should be 'Geneva', 'weather' should be 'whether', 'easely' should be 'easily', and 'chrial' in the Fig. 1 caption should be 'chiral'.
- [References] Ref. [20] is cited only as arXiv:2405.10196 without a journal or version; since the endpoint analysis relies on it, please provide the preprint version or a DOI if available.
- [Fig. 3] The labels 'HotQCD [4,4]' and 'BiPar Multi' are not defined in the caption; 'BiPar' should be spelled out as the Bielefeld-Parma collaboration and the meaning of '[4,4]' stated.
- [Eq. (2)] The notation kappa_11^{ls} is not defined explicitly; the subscript 11 convention (presumably kappa times mu_l mu_s normalized by T^2) should be stated for readers not familiar with the companion papers.
- [Sec. 2] The sentence 'We performed fits in which Tc and t1,x are kept as free fit parameters as well as using for Tc and t1,x the values determined from the scaling fits to the order parameter M given in Eqs. 36, 37' lacks the actual numbers from these two procedures; please provide the resulting Tc values or refer to a table.
Circularity Check
No significant circularity: the endpoint is a direct fit of Pade-derived Lee-Yang edge locations to an explicitly quoted linear scaling ansatz, with universal exponents and edge constant imported from external references; the chiral T_c is cross-checked against an earlier continuum result.
full rationale
The paper's Sec. 2 chiral-scale analysis is a conventional fit: Eq. (6) uses O(2)/O(4) scaling functions determined externally in Ref. [11], and the non-universal constants (T_c, z0, h0^{-1/beta}) are fit parameters. The resulting T_c=143.9(5) MeV for N_tau=8 is compared with the previous continuum estimate in Ref. [7], providing an independent benchmark. No quantity is defined in terms of the target result. Sec. 3's endpoint extraction is also a fit rather than a disguised tautology: Eq. (7) states the linear scaling-field ansatz explicitly, with T_cep and mu_cep as parameters; the Lee-Yang edge positions come from multipoint Pade approximants to lattice data, and the universal scaling form t/h^{1/beta delta}=z_LY is taken from external Ref. [16] with the asymptotic behaviors from Refs. [17,18]. The quoted endpoint (105(+8,-18), 422(+80,-35)) MeV is the output of fitting that ansatz to those pole locations, not a separately fitted input renamed as a prediction. The self-citations [14,15,20] support the method via earlier benchmarks (Roberge-Weiss transition, 2d Ising model) and provide details of the same analysis; such internal consistency checks are not load-bearing in a circular sense. The main caveats -- the linear mixing assumption, the temperature range of the Lee-Yang fit, the pole identification, and the lack of continuum extrapolation -- are substantive systematic/correctness concerns, and the paper itself notes the continuum extrapolation is still pending, but they are not instances of circular reasoning. No equation reduces to its own input by construction.
Assumptions & free parameters
free parameters (7)
- Chiral critical temperature T_c (N_tau=8, staggered) =
143.8(2) MeV in Fig. 2 caption; 143.9(5) MeV in text
- z0 =
1.45(3)
- h0^(-1/delta) =
39.0(3)
- Chiral curvature coefficients kappa2_l, kappa2_s, kappa2_ls =
kappa_B^{muS=0}=0.015(1); kappa_B^{nS=0}=0.893(35) kappa_B^{muS=0}
- Pseudo-critical correction-to-scaling coefficients t1_x, tc_x, t2_x, t3_x =
not quoted individually
- Finite-volume shift for H=1/160 T_pc values =
0.25 MeV
- CEP scaling-field mixing parameters (alpha_t, beta_t, alpha_h, beta_h) =
T_cep=105(+8,-18) MeV; mu_cep=422(+80,-35) MeV
assumptions (6)
- domain assumption O(4) universality class for the two-flavor chiral transition in the continuum, and O(2) on staggered lattices at finite lattice spacing
- ad hoc to paper The scaling-field ansatz in Eqs. (1)-(2): t is a linear combination of temperature and squared chemical potentials, and h is proportional to H
- ad hoc to paper The linear mixing ansatz for CEP scaling fields, Eq. (7), with unknown scaling directions treated as linear in T and mu_B
- domain assumption The closest pole of the multipoint Pade approximant is identified with the Lee-Yang edge and obeys universal scaling
- domain assumption Finite-size effects are controlled by adjusting the spatial extent and are small for H>=1/40; for H=1/160 they are corrected by a global shift
- domain assumption Universal scaling functions of the 3d O(2) model apply to HISQ staggered lattice QCD at N_tau=8 in the scaling window
Cite this review
Pith. "Pith review of The QCD phase diagram, universal scaling, and Lee-Yang zeros." pith.science (2026). https://pith.science/paper/MK6PD4QS
@misc{pith2026250119336,
author = {Pith},
title = {Pith review of: The QCD phase diagram, universal scaling, and Lee-Yang zeros},
year = {2026},
howpublished = {\url{https://pith.science/paper/MK6PD4QS}},
note = {Machine review of arXiv:2501.19336}
}
read the original abstract
We will report on current progress in the understanding of the QCD phase diagram, including universal scaling in the chiral limit and the vicinity of the QCD critical point. In the latter case we will discuss the universal scaling of Lee-Yang zeros and their determination from multi-point Pad\'e approximations to the baryon number density at imaginary chemical potentials. In particular, reported results include the critical phase transition temperature, the curvature of the critical and pseudo-critical transition temperature with respect to the chemical potential and the location of the QCD critical point.
Figures
Forward citations
Cited by 1 Pith paper
-
Lee--Yang edge singularities in Nonlocal Nambu--Jona-Lasinio Model
In a nonlocal NJL model, Lee-Yang edge singularities follow trajectories that end at the QCD critical point, with critical exponent 1.494(1) matching mean-field scaling.
Reference graph
Works this paper leans on
-
[8]
H. T. Ding, O. Kaczmarek, F. Karsch, P. Petreczky, M. Sarkar, C. Schmidt, S. Sharma, Curvature of the chiral phase transition line from the magnetic equation of state of (2+1)-flavor QCD, Phys. Rev. D 109 (11) (2024) 114516. arXiv:2403.09390, doi: 10.1103/PhysRevD.109.114516
arXiv 2024
-
[20]
D. A. Clarke, P. Dimopoulos, F. Di Renzo, J. Goswami, C. Schmidt, S. Singh, K. Zambello, Searching for the QCD crit- ical endpoint using multi-point Pad´ e approximations (5 2024). arXiv:2405.10196. 5
arXiv 2024
- [11]
- [1]
-
[2]
A. Pasztor, The QCD phase diagram at finite temperature and density - a lattice perspective, PoS LATTICE2023 (2024) 108. doi:10.22323/1.453.0108
-
[3]
R. D. Pisarski, F. Wilczek, Remarks on the Chiral Phase Tran- sition in Chromodynamics, Phys. Rev. D 29 (1984) 338–341. doi:10.1103/PhysRevD.29.338
-
[4]
R. D. Pisarski, F. Rennecke, Conjectures about the Chiral Phase Transition in QCD from Anomalous Multi-Instanton In- teractions, Phys. Rev. Lett. 132 (25) (2024) 251903. arXiv: 2401.06130, doi:10.1103/PhysRevLett.132.251903
arXiv 2024
-
[5]
A. Pelissetto, E. Vicari, Relevance of the axial anomaly at the finite-temperature chiral transition in QCD, Phys. Rev. D 88 (10) (2013) 105018. arXiv:1309.5446, doi:10.1103/ PhysRevD.88.105018
arXiv 2013
Show all 20 references
-
[6]
Giacosa, G
F. Giacosa, G. Kov´ acs, P. Kov´ acs, R. D. Pisarski, F. Ren- necke, Anomalous U(1)A couplings and the Columbia plot, Phys. Rev. D 111 (1) (2025) 016014. arXiv:2410.08185, doi: 10.1103/PhysRevD.111.016014
2025 arXiv
-
[7]
H. T. Ding, et al., Chiral Phase Transition Temperature in ( 2+1 )-Flavor QCD, Phys. Rev. Lett. 123 (6) (2019) 062002. arXiv:1903.04801, doi:10.1103/PhysRevLett.123.062002
2019 arXiv
-
[9]
L. Dini, P. Hegde, F. Karsch, A. Lahiri, C. Schmidt, S. Sharma, Chiral phase transition in three-flavor QCD from lattice QCD, Phys. Rev. D 105 (3) (2022) 034510. arXiv:2111.12599, doi: 10.1103/PhysRevD.105.034510
2022 arXiv
-
[10]
A. Y. Kotov, M. P. Lombardo, A. Trunin, QCD transition at the physical point, and its scaling window from twisted mass Wilson fermions, Phys. Lett. B 823 (2021) 136749. arXiv:2105.09842, doi:10.1016/j.physletb.2021.136749
2021 arXiv
-
[12]
J. J. Rehr, N. D. Mermin, Revised Scaling Equation of State at the Liquid-Vapor Critical Point, Phys. Rev. A 8 (1973) 472–480. doi:10.1103/PhysRevA.8.472
1973 doi
-
[13]
Nonaka, M
C. Nonaka, M. Asakawa, Hydrodynamical evolution near the QCD critical end point, Phys. Rev. C 71 (2005) 044904. arXiv: nucl-th/0410078, doi:10.1103/PhysRevC.71.044904
2005 arXiv
-
[14]
Dimopoulos, L
P. Dimopoulos, L. Dini, F. Di Renzo, J. Goswami, G. Nico- tra, C. Schmidt, S. Singh, K. Zambello, F. Ziesch´ e, Contri- bution to understanding the phase structure of strong inter- action matter: Lee-Yang edge singularities from lattice QCD, Phys. Rev. D 105 (3) (2022) 034513....
2022 arXiv
-
[15]
Singh, M
S. Singh, M. Cipressi, F. Di Renzo, Exploring Lee-Yang and Fisher zeros in the 2D Ising model through multipoint Pad´ e approximants, Phys. Rev. D 109 (7) (2024) 074505. arXiv: 2312.03178, doi:10.1103/PhysRevD.109.074505
2024 arXiv
-
[16]
Connelly, G
A. Connelly, G. Johnson, F. Rennecke, V. Skokov, Universal Location of the Yang-Lee Edge Singularity in O(N ) Theories, Phys. Rev. Lett. 125 (19) (2020) 191602. arXiv:2006.12541, doi:10.1103/PhysRevLett.125.191602
2020 arXiv
-
[17]
C.-N. Yang, T. D. Lee, Statistical theory of equations of state and phase transitions. 1. Theory of condensation, Phys. Rev. 87 (1952) 404–409. doi:10.1103/PhysRev.87.404
1952 doi
-
[18]
M. A. Stephanov, QCD critical point and complex chemical potential singularities, Phys. Rev. D 73 (2006) 094508. arXiv: hep-lat/0603014, doi:10.1103/PhysRevD.73.094508
2006 arXiv
-
[19]
Bollweg, J
D. Bollweg, J. Goswami, O. Kaczmarek, F. Karsch, S. Mukher- jee, P. Petreczky, C. Schmidt, P. Scior, Taylor expansions and Pad´ e approximants for cumulants of conserved charge fluctua- tions at nonvanishing chemical potentials, Phys. Rev. D 105 (7) (2022) 074511. arXiv:2202.0...
2022 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.