Pith. Sign in

REVIEW

The density of states method at non-zero chemical potential

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Zoltan Fodor, Sandor D. Katz, Christian Schmidt

High Energy Physics - Lattice · Submitted 25 Jan 2007 · 23 references

Pith's verdict: A real methods advance with a genuinely interesting physics hint, but the triple-point claim rests on a lambda extrapolation performed at the one chemical potential where the claim isn't.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 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 · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-lat/0701022 (version 2)
Pith IDLXIZTXSU, Pith's permanent record of this paper and its review
CitationChristian Schmidt, Sandor D. Katz, Zoltan Fodor. "The density of states method at non-zero chemical potential." arXiv:hep-lat/0701022 (2007). https://pith.science/paper/LXIZTXSU
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

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

By fixing the average plaquette and integrating the complex phase back in afterwards, the authors push lattice QCD to chemical potentials about three times further than reweighting from zero density, and read off three plaquette branches, a

The reading

Lattice QCD stalls at nonzero density because the fermion determinant becomes complex and importance sampling fails; the phase of the determinant fluctuates by more than $O(1)$ and averages almost to zero. This paper asks whether that obstacle can be stepped around far enough to reach the low-temperature, high-density corner of the QCD phase diagram, where colour-superconducting phases are expected. The authors hold the average plaquette fixed at a chosen value, simulate with the magnitude of the determinant only, and then restore the complex phase by numerical integration over the constrained ensemble. On $4^4$ and $6^4$ lattices with four flavours of unimproved staggered fermions this reaches $a^2$-scale statistics at $a^2$; concretely, chemical potentials up to $a^2$, roughly $a^2$ beyond what multi-parameter reweighting from $a^2$ achieved. The payoff is a plaquette expectation value that splits into three branches, which the authors read as three phases meeting at a triple point near $T_{tri}approx 137$ MeV and $a^2$ on the $N_t=6$ lattice, with the new phase carrying a quark number density of 10 to 20 times nuclear density. The authors are explicit that the lattices are coarse, the flavour content unphysical, and neither the continuum nor the thermodynamic limit has been taken.

What carries the argument

The engine is a constrained partition function, the density of states $a^2$, defined by fixing one observable — here the real part of the average plaquette $P$ — and weighting the remaining measure by $g=|det M|e^{-S_G}$ while carrying the complex phase $e^{i heta}$ as the observable to be integrated later. The delta function is replaced by a Gaussian constraint $V(x)=12gamma(x-P)^2$ added to the molecular-dynamics force, and the density itself is recovered from the fluctuation-dissipation relation $a^2$, so no binning in $P$ is needed. A reweighting ratio $R=g(mu,eta)/g(mu_0,eta_0)$ transports each constrained ensemble to nearby chemical potential and coupling, and the final expectation values are one-dimensional integrals such as $a^2$: the density of states supplies the weight, and the phase factor $a^2$, obtained from all eigenvalues of the reduced fermion matrix stored configuration by configuration, supplies the suppression. The three-branch structure is read off from the peak position of the plaquette susceptibility built from these integrals.

Extended reading notes

Core claim

The paper's claim is that the density-of-states (factorization) method, applied to four-flavour unimproved staggered fermions, opens a window on the QCD phase diagram at $a^2$ up to about 3, and that inside that window the plaquette expectation value shows three distinct regimes: a low-temperature hadronic branch near $a^2$, a high-temperature branch near $a^2$, and — for $a^2$ — a new low-temperature branch near $a^2$ that is selected because it is the region least suppressed by the complex phase. The critical coupling falls with increasing $a^2$ below $a^2$ and rises above it, and the quark number density jumps at the same $a^2$ from roughly $(2-3),n_N$ to $(10-20),n_N$. Taken together the authors interpret this as a triple point connecting three phases, located at $a^2$ MeV and $a^2$ MeV on the $N_t=6$ lattice, shifting down from $T_{tri}approx 148$ MeV on the $N_t=4$ lattice. The authors call this the first numerical evidence from lattice QCD for a third phase besides the hadronic phase and the quark-gluon plasma, and nominate the new phase as a candidate colour superconductor, while stating that they make no claim about the order of any of the transition lines.

Load-bearing premise

The load-bearing premise is that the three branches of the plaquette expectation value, together with the coincident jump in quark number density, mark three genuine thermodynamic phases meeting at a triple point, rather than finite-volume or lattice-cutoff artifacts of coarse four-flavour lattices with a plaquette that is not an order parameter.

Editorial extensions

If this is right

  • If the three plaquette branches are genuine phases, the QCD phase diagram acquires a third region at $a^2$ and $a^2$ MeV, bounded above by a transition line from the quark-gluon plasma and to the left by the hadronic phase.
  • The new phase is a candidate colour superconductor, and the paper notes its triple point corresponds to $T_c approx 0.46,E_F$, close to the bound $T_c leq 0.35,E_F$ inferred from trapped cold atoms.
  • The low-temperature phase boundary is steep and nearly temperature independent; extrapolating it to $T=0$ gives $a^2$ MeV, or $a^2$, at the low edge of phenomenological expectations.
  • Because the phase boundary at low temperature is unchanged between $am=0.05$ and $am=0.03$ while the pion mass drops by about $a^2$, the transition is argued to be onset of baryonic matter rather than pion condensation.

Reading between the lines

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

  • The method's real advantage is not just reach in $a^2$ but selectivity: by simulating directly at the plaquette values that survive the phase suppression, it sidesteps the overlap problem that limits reweighting from $a^2$ ensembles, which suggests the same constrained-sampling trick could be paired with newer complex-phase techniques rather than replacing them.
  • The triple point's temperature moves by about 11 MeV between $N_t=4$ and $N_t=6$, which is the size of the cutoff uncertainty on the central result; a single $N_t=8$ determination of the same three branches would be the cheapest decisive test of whether the structure survives.
  • If the third branch is really a dense-matter phase, its density of 10 to 20 times nuclear density at $T approx 100$ MeV should be cross-checkable against canonical-ensemble and strong-coupling calculations on the same lattice geometry, which the paper cites as giving similar densities.
  • The plaquette is not an order parameter, so a sharper signature — an appropriately renormalized Polyakov loop, a baryon-number susceptibility, or a diquark correlator — would be needed before the three-branch reading could be called a phase identification rather than a suggestive coincidence.
Original abstract

We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. For this purpose we employ the density of states method. Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram.

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, 3,472 input tokens, $0.0001 reported.

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

The abstract, sentence by sentence

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

We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. For this purpose we employ the density of states method. Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram.

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 We study the QCD phase diagram by first principle lattice calculations at so far unreached high densities. contribution 0.98 0.92 0.06 asserted 0.96
s1 For this purpose we employ the density of states method. method 0.99 0.68 0.02 asserted 0.98
s2 Unimproved staggered fermions, which describe four quark flavors in the continuum are used in this analysis. method 1.00 0.57 0.02 asserted 0.98
s3 Though the method is quite expensive, small lattices show an indication for a triple-point connecting three different phases on the phase diagram. result 1.00 0.90 0.74 hedged possibility 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 v2
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. Results for simulations at β = 4 .98, µ = 0 .3, λ = 0 .02, n f = 4, am = 0 .05, and number of lattice points: 4 4 . Shown are the density of states ρ ( x), the phase factor hcos( θ ) i, and their product. 0.0001 0.001 0.01 0.11 10 2.7 2.75 2.8 2.85 2.9 2.95 3 3.05 3.1 3.15 <cos( θ)> P µ=0.30 µ=0.32 µ=0.34 µ=0.36 µ=0.38 µ=0.40
Figure 2. Results for simulations at β = 5 .1, λ = 0 .01, n f = 4, am = 0 .05, and number of lattice points: 6 4 . Shown is the suppression from the complex phase of the fermion determinant hcos ( θ ) i for different chemical potentials. pectation values are strongly altered by the phase factor. figure 2 demonstrates also the advantage of he DOS method over the other approaches of lattice QCD at finite density. Using the DOS method one is able to do simulations directly at those Plaquette values which are relevant at finite density. This results in an even better overlap than that o
Figure 3. Simulation parameters as in figure 1. Shown is the plaquette susceptibility as function of the coupling β . 4.93 4.94 4.95 4.96 4.97 4.98 4.99 5 0 0.1 0.2 0.3 0.4 0.5 β λ/m with phase without phase Fodor, Katz (2002) extrapolated
Figure 4. The λ → 0 extrapolation of the critical couplings on the 4 4 lattice at aµ = 0 .3. 4. The Plaquette expectation value and the phase diagram Let us now discuss results for the plaquette expectation values from the 6 4 lattice as shown in figure 5 and figure 6. At chemical potentials µ . 0 .36, the plaquette signals the QCD transition through a rapid crossover from a low temperature phase of < P > ∼ 2 .9 to a high temperature phase of < P > ∼ 3 .1. For µ & 0 .36 the plaquette expectation value at small temperatures drops to < P > ∼ 2 .85. This new low temperature
Figure 5. Results for simulations at λ = 0 .01, n f = 4, am = 0 .05, and number of lattice points: 6 4 . Shown is the Plaquette expectation value as a function of the coupling β for different chemical potentials 2.83 2.84 2.85 2.86 2.87 2.88 2.89 0.24 0.26 0.28 0.3 0.32 0.34 0.36 µ P β=5.1
Figure 6. Simulation parameters as in figure 5. Shown is the plaquette expectation value at fixed coupling, as a function of the chemical potential. potentials is caused by the fermion determinant. As one can see in figure 2 the region around P ∼ 2 .85 is the region which is less suppressed by the phase factor. Another interesting observation is that the critical coupling, which is decreasing in µ for µ < 0 .36 starts to increase for µ > 0 .36. The plaquette expectation value thus suggests the existence of three different phases in the ( T , µ)-diagram with a triple point, where all
Figure 7. The phase diagram in physical units from Nt = 6 and 8 lattices. the phase diagram in physical units. The phase boundaries were determined by calculating the peaks in the plaquette susceptibility. The points at T = 93 MeV are from calculations on 6 3 ×8 lattices. Note, that we make no statement about the order of the transition lines. To determine the order of the transition one has to perform a finite-size-scaling analysis which is beyond the scope of this article. The triple point is located around µ tri q ≈ 300 MeV, however its temperature ( T tri ) decreases from T tri ≈ 148 MeV…
Figure 8. The quark number density at constant temperature T = 143 MeV (4 4 lattice), T = 124 MeV (6 4 lattice) and T = 93 MeV (6 3 × 8 lattice). The thermodynamic quantity n q are given as usual by nq = 1 a3N3s Nt d ln det M d(aµ )

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

References
How later papers used it

50 later papers cite it

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

Pith read 7 of the 56 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 method or tool 1 of 3 papers with this label · 1 by the same authors

By the same authors

  • Used this paper's density of states method to extend the applicability range, studying n_f=4 staggered QCD on 6^4 and 8·6^3 lattices.

    Recent Result in QCD Thermodynamics from the Lattice (2007) (same authors) · 0.99 · paraphrased: the licence does not allow quoting

How the literature uses it

44 arXiv papers cite this in 50 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 whose main result relates to this one's

Found among the papers closest to this one by title and main claim. Pith's classifier compared each pair's main results and a second model family confirmed the relation; checked against two further families, shown conflicts agreed 118 times in 120 and shown same results 218 times in 239. The quoted sentence is the other paper's main claim.

States the same result 2 papers

Papers citing it · 28 sentences without a confident label

Builds on it 1 paper
Uses its method or tool 3 papers
Cites it as background 18 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.