Pith. Sign in

REVIEW

Critical point of QCD at finite T and \mu, lattice results for physical quark masses

High Energy Physics - Lattice · Submitted 8 Feb 2004 · 30 references

Pith's verdict: A real physics update — the endpoint drops from 720 to 360 MeV at physical quark masses — but the reweighting overlap at the endpoint is argued, not measured, and the scale systematic is wider than the quoted error.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-lat/0402006 (version 1)
Pith IDSW4BQUXA, Pith's permanent record of this paper and its review
CitationS.D. Katz, Z. Fodor. "Critical point of QCD at finite T and \mu, lattice results for physical quark masses." arXiv:hep-lat/0402006 (2004). https://pith.science/paper/SW4BQUXA
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The paper locates the QCD critical endpoint at $T_E=162$ MeV and baryon chemical potential $200$ MeV using lattice QCD with physical quark masses.

The reading

The paper asks where, on the plane of temperature and baryon density, QCD stops changing smoothly and begins to jump. That location is the critical endpoint: below it the transition between hadrons and quark-gluon plasma is a crossover, above it a first-order transition, and near it long-wavelength fluctuations could leave experimental signatures in heavy-ion collisions. Using lattice QCD with 2+1 flavours of staggered quarks at approximately physical masses, the authors report the endpoint at $T_E=162$ MeV and baryon chemical potential $200$ MeV, together with a zero-density crossover temperature $T_c=164$ MeV. The result improves their earlier work by lowering the light quark masses by a factor of three and increasing the largest volume by a factor of three; the lower masses moved the endpoint down, while the volume change left it unchanged. The authors state that the continuum extrapolation is still missing, so the numbers are a controlled step rather than a final answer.

What carries the argument

The method is overlap-improving multi-parameter reweighting combined with Lee-Yang zeros. Configurations are generated at $200$ and at the zero-density transition temperature, where importance sampling works; each configuration is then given a Boltzmann weight for nonzero $200$ and for the lower transition temperature at that $200$, with the gauge coupling reweighted by the Ferrenberg-Swendsen procedure. The partition function is continued into the complex $200$ plane, and its Lee-Yang zeros are located. Their finite-volume scaling carries the argument: zeros that stay away from the real axis mean a crossover, zeros that approach it as $1/V$ mean a first-order transition, and the endpoint is the $200$ at which the extrapolated imaginary part vanishes. The physical scale is set from the Sommer scale $R_0$ and the $200$ mass.

Extended reading notes

Core claim

The central claim is that the QCD critical endpoint sits at $T_E=162$ MeV and $200$ MeV for 2+1-flavour QCD with physical quark masses on $L_t=4$ lattices. In lattice units the endpoint is $200$. The evidence comes from the infinite-volume extrapolation of Lee-Yang zeros of the reweighted partition function: at small $200$ the extrapolated imaginary part of the zero, $200$, is inconsistent with zero and the transition is a crossover; as $200$ increases it falls, and at the endpoint the crossover disappears and finite-volume scaling becomes consistent with a first-order transition. The paper also reports a crossover-line curvature $T/T_c=1-C200^2/T_c^2$ with $C=0.0032(1)$. The authors explicitly note that the continuum extrapolation has not been done.

Load-bearing premise

The load-bearing premise is that the reweighting from $200$ ensembles still samples the nonzero-density transition well enough at the endpoint for the Lee-Yang zeros to be reliable; the authors also acknowledge that the continuum extrapolation is still missing.

Editorial extensions

If this is right

  • Heavy-ion searches for critical-point signatures can target the region around $T\approx162$ MeV and baryon chemical potential $\approx360$ MeV rather than the much larger $\mu_E\approx700$ MeV suggested by earlier large-mass estimates.
  • The crossover line's small curvature, $C=0.0032(1)$, gives a concrete input for the QCD equation of state at small baryon density.
  • Because lowering the light quark masses moved $\mu_E$ down from about 720 MeV to 360 MeV, quark-mass effects must be controlled in any comparison between lattice endpoints and experiment.
  • The insensitivity to a threefold volume increase supports the reliability of the finite-volume Lee-Yang analysis, leaving the lattice spacing as the main systematic uncertainty.

Reading between the lines

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

  • A quantitative overlap diagnostic at the endpoint, such as the distribution of reweighting weights or an effective sample size, would test the method independently of the statistical error bars on the zeros.
  • Agreement between this reweighting result and an imaginary-$\mu$ analytic continuation at the same quark masses would be a strong cross-check, because the two methods fail in different ways.
  • Since $T_E$ and $T_c$ differ by only about 2 MeV, within errors, the data do not yet establish whether the endpoint lies exactly on the zero-density transition line or slightly below it.
  • Mapping $T$ and $\mu_B$ to heavy-ion collision energies requires a freeze-out model; the endpoint's location is a thermodynamic statement, not by itself a prediction of which beam energy will see critical fluctuations.
Original abstract

A critical point (E) is expected in QCD on the temperature (T) versus baryonic chemical potential (\mu) plane. Using a recently proposed lattice method for \mu \neq 0 we study dynamical QCD with n_f=2+1 staggered quarks of physical masses on L_t=4 lattices. Our result for the critical point is T_E=162 \pm 2 MeV and \mu_E= 360 \pm 40 MeV. For the critical temperature at \mu=0 we obtained T_c=164 \pm 2 MeV. This work extends our previous study [Z. Fodor and S.D.Katz, JHEP 0203 (2002) 014] by two means. It decreases the light quark masses (m_{u,d}) by a factor of three down to their physical values. Furthermore, in order to approach the thermodynamical limit we increase our largest volume by a factor of three. As expected, decreasing m_{u,d} decreased \mu_E. Note, that the continuum extrapolation is still missing

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,402 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.98 · human participants 0.01 · controlled comparison 0.20 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.93 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

A critical point (E) is expected in QCD on the temperature (T) versus baryonic chemical potential (\mu) plane. Using a recently proposed lattice method for \mu \neq 0 we study dynamical QCD with n_f=2+1 staggered quarks of physical masses on L_t=4 lattices. Our result for the critical point is T_E=162 \pm 2 MeV and \mu_E= 360 \pm 40 MeV. For the critical temperature at \mu=0 we obtained T_c=164 \pm 2 MeV. This work extends our previous study [Z. Fodor and S.D.Katz, JHEP 0203 (2002) 014] by two means. It decreases the light quark masses (m_{u,d}) by a factor of three down to their physical values. Furthermore, in order to approach the thermodynamical limit we increase our largest volume by a factor of three. As expected, decreasing m_{u,d} decreased \mu_E. Note, that the continuum extrapolation is still missing

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 A critical point (E) is expected in QCD on the temperature (T) versus baryonic chemical potential (\mu) plane. background 1.00 0.05 0.07 hedged possibility 0.88
s1 Using a recently proposed lattice method for \mu \neq 0 we study dynamical QCD with n_f=2+1 staggered quarks of physical masses on L_t=4 lattices. method 0.85 0.89 0.04 asserted 0.94
s2 Our result for the critical point is T_E=162 \pm 2 MeV and \mu_E= 360 \pm 40 MeV. result 1.00 0.97 0.86 asserted 1.00
s3 For the critical temperature at \mu=0 we obtained T_c=164 \pm 2 MeV. result 1.00 0.97 0.09 asserted 1.00
s4 This work extends our previous study [Z. contribution 0.58 0.43 0.02 asserted 0.82
s5 Fodor and S.D.Katz, JHEP 0203 (2002) 014] by two means. background 0.51 0.32 0.02 asserted 0.67
s6 It decreases the light quark masses (m_{u,d}) by a factor of three down to their physical values. contribution 0.61 0.89 0.03 asserted 0.99
s7 Furthermore, in order to approach the thermodynamical limit we increase our largest volume by a factor of three. method 0.59 0.93 0.03 asserted 0.99
s8 As expected, decreasing m_{u,d} decreased \mu_E. interpretation 0.67 0.85 0.06 asserted 0.97
s9 Note, that the continuum extrapolation is still missing limitation 1.00 0.37 0.03 asserted 0.98

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. shows Im(β∞ 0 ) as a function of µ enlarged around the endpoint µend. The picture is simple and reflects the physical expectations. For small µ-s the extrapolated Im(β∞ 0 ) is inconsistent with a vanishing value, and the prediction is a crossover. Increas￾ing µ the value of Im(β∞ 0 ) decreases, thus the transition becomes consistent with a first order phase transition. The statistical error was determined by a jackknife analysis using subsamples of the total Ls=6,8,10 and 12 partition functions. Our primary result in lattice units is µend = 0.1825(75)
Figure 2. The phase diagram in physical units. Dotted line illustrates the crossover, solid line the first order phase transition. The small square shows the endpoint. The depicted errors originate from the reweighting procedure. Note, that an overall additional error of 1.3% comes from the error of the scale determination at T=0. Combining the two sources of uncertainties one obtains TE = 162 ± 2 MeV and µE = 360 ± 40 MeV. The above result is a significant improvement on our previous analysis [8] by two means. We increased the physical volume by a factor of three and decreased the light…

Figures reproduced from arXiv:hep-lat/0402006 by the authors.

References
How later papers used it

160 later papers cite it

Evidence: substantial. 141 of the citing papers share no author with it, and 37 of those carry a confident label from Pith's citation classifier.

Pith read 27 of the 217 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its data or values 2 of 20 papers with this label
How the literature uses it

130 arXiv papers cite this in 184 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 107 sentences without a confident label

Contests it 1 paper
  • Our results about the critical surface appear to be in qualitative contradiction with those of [ref] , [ref] , which both conclude for the existence of a critical point $(μ_E,T_E)$ at small chemical potential $μ_E/T_E 1$.

    Exploring the QCD phase diagram · Discussion · contest 0.72

Builds on it 1 paper
Uses its method or tool 8 papers
Uses its data or values 17 papers
  • There is as yet no complete agreement as to the location of the quark matter critical endpoint, as obtained from lattice studies: in [ref] one finds a discussion of various proposed points in the phase diagram, such as one (see [ref] ) with coordinates $(μ_B^CEP/T^CEP, T^CEP/T^C_0) ≈(1.68, 0.94)$, where $T^C_0$ is the crossover temperature at $μ_B = 0$, for which we follow [ref] and take the value $≈$ 150 MeV; but also another with coordinates $(μ_B^CEP/T^CEP, T^CEP/T^C_0) ≈(2.2, 0.99)$ (see [ref] ).

    The Trajectory of the Cosmic Plasma Through the Quark Matter Phase Diagram · Section · 0.80

  • The lattice data includes a first order phase transition at $T=270$ MeV and a pressure functional P(T) similar to Refs. [ref] at $μ=0$ for pure gauge, a crossover at vanishing chemical potential with a transition temperature of $171$ MeV (determined as the peak of the change of the chiral condensate and $Φ$) and the location of the critical end-point (at $μ_c=354$ MeV , $T_c=167$ MeV for symmetric matter in accordance with one of the existent calculations [ref] ).

    A Novel Approach to Model Hybrid Stars · Introduction · 0.97

  • The case of physical quark masses, after conversion to continuum units, is shown in [ref] (right) [ref] .

    Status of Lattice Studies of the QCD Phase Diagram · Calculations at finite density · 0.92

Compares with it 9 papers
Cites it as support 1 paper
  • In particular, it includes: $∙$ lattice computations of topological susceptibility in the vicinity of the phase transition [ref] ; $∙$ lattice computations of the critical temperature $T_c(μ)$ as a function of $μ$ at small $μ$ [ref] ; $∙$ analysis of the holographic models of QCD in vicinity of the phase transition at nonzero temperature [ref] .

    Confinement- Deconfinement Phase Transition in Hot and Dense QCD at Large N · Confinement- Deconfinement Phase Transition in hot QCD at large $N$. · 1.00

Cites it as background 39 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.