Recognition: unknown
Finite-density equation of state of hot QCD using the complex Langevin equation
Pith reviewed 2026-05-10 00:30 UTC · model grok-4.3
The pith
Lattice simulations using the complex Langevin equation determine the QCD equation of state at high baryon densities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The QCD equation of state is obtained by computing the baryon density and the pressure on the lattice as functions of baryon chemical potential and temperature, using the complex Langevin equation for continuum-extrapolated simulations at the physical point and unprecedentedly high densities, with the dynamics shown to converge correctly and to agree with known limits.
What carries the argument
The complex Langevin equation, which evolves the gauge fields in a complexified space to sample the QCD path integral with a complex fermion determinant at finite chemical potential.
If this is right
- Thermodynamic quantities such as energy density and entropy density follow directly from the pressure and baryon density.
- The results extend the range of reliable lattice data for use in hydrodynamic modeling of heavy-ion collisions.
- Consistency with hard-thermal-loop perturbation theory at high temperature supports the method in the deconfined regime.
- The controlled convergence allows systematic studies at densities where the sign problem blocks standard Monte Carlo methods.
Where Pith is reading between the lines
- The same framework could be used to compute higher-order susceptibilities that probe the vicinity of a possible critical endpoint.
- The equation of state tables produced here could be directly imported into simulations of neutron-star mergers.
- Extension of the method to lower temperatures might eventually map the full phase diagram if convergence remains under control.
Load-bearing premise
The complex Langevin dynamics converges to the correct physical distribution rather than a spurious solution.
What would settle it
A direct comparison at overlapping values of temperature and chemical potential showing disagreement with an independent method, such as Taylor expansion results from other groups or resummed perturbation theory at high temperature, would falsify the reliability of the computed equation of state.
Figures
read the original abstract
We present the results of continuum-extrapolated lattice simulations of quantum chromodynamics (QCD) above the crossover temperature and for unprecedentedly high baryon densities at the physical point, employing the complex Langevin equation. In particular, we determine the QCD equation of state by computing the baryon density as well as the pressure as functions of the baryon chemical potential and the temperature. Potential issues with wrong convergence of complex Langevin dynamics are under control and we indeed find agreement with previous lattice studies working at smaller chemical potentials, as well as with perturbative hard-thermal-loop calculations at high temperatures.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports continuum-extrapolated lattice QCD simulations above the crossover temperature at the physical point, using the complex Langevin equation to compute the baryon density and pressure as functions of baryon chemical potential μ_B and temperature T. It reaches unprecedentedly high densities and asserts that wrong-convergence issues of the complex Langevin dynamics are under control, with results agreeing with prior lattice work at smaller μ_B and with hard-thermal-loop perturbation theory at high T.
Significance. If the convergence control is rigorously established, the work would provide a significant extension of lattice QCD into the high-density regime inaccessible to standard Monte Carlo methods due to the sign problem. The direct simulation yields the EOS without fitted parameters or self-referential definitions, offering falsifiable predictions at physical quark masses and continuum limit that can be compared to heavy-ion phenomenology and other non-perturbative approaches.
major comments (2)
- [Abstract] Abstract: the central claim that 'potential issues with wrong convergence of complex Langevin dynamics are under control' is load-bearing for the entire EOS determination at high μ_B, yet the abstract (and by extension the manuscript) provides no explicit quantitative diagnostics, thresholds, or figures showing the distribution of the imaginary part of the action, drift-term norms, or consistency checks with reweighting at the largest accessed densities.
- [Results] Results section: the stated agreement with previous lattice studies at smaller chemical potentials and with HTL perturbation theory lacks reported error budgets, quantitative deviation measures, or details of the continuum-extrapolation procedure (e.g., fit forms, number of lattice spacings, or χ² values), preventing assessment of whether the high-density extrapolation is reliable or merely consistent within large uncertainties.
minor comments (1)
- [Abstract] Abstract: the phrase 'unprecedentedly high baryon densities' should be accompanied by a concrete range or maximum value of μ_B/T to allow immediate gauging of the advance relative to prior work.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each point below and will incorporate revisions to improve the clarity and completeness of the presentation.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claim that 'potential issues with wrong convergence of complex Langevin dynamics are under control' is load-bearing for the entire EOS determination at high μ_B, yet the abstract (and by extension the manuscript) provides no explicit quantitative diagnostics, thresholds, or figures showing the distribution of the imaginary part of the action, drift-term norms, or consistency checks with reweighting at the largest accessed densities.
Authors: We agree that the abstract would benefit from a more explicit reference to the convergence diagnostics. While the main text (particularly the sections on simulation setup and validation) already presents quantitative checks—including histograms of the imaginary part of the action, norms of the drift term, and comparisons to reweighting at moderate densities—we will revise the abstract to briefly mention these controls and their outcomes at the highest μ_B values. This will make the load-bearing claim more self-contained without altering the manuscript's technical content. revision: yes
-
Referee: [Results] Results section: the stated agreement with previous lattice studies at smaller chemical potentials and with HTL perturbation theory lacks reported error budgets, quantitative deviation measures, or details of the continuum-extrapolation procedure (e.g., fit forms, number of lattice spacings, or χ² values), preventing assessment of whether the high-density extrapolation is reliable or merely consistent within large uncertainties.
Authors: We acknowledge that additional quantitative details will aid assessment. The manuscript already performs a continuum extrapolation using multiple lattice spacings and reports statistical errors, but we will expand the results section to include: (i) explicit error budgets separating statistical and systematic contributions for the EOS quantities, (ii) quantitative deviation measures (e.g., relative differences and significance in units of combined uncertainty) for the comparisons to prior lattice data and HTL, and (iii) full details of the extrapolation procedure, including the fit functional form, the number of spacings employed, and the resulting χ² values. These additions will be presented in a revised results section and associated figures/tables. revision: yes
Circularity Check
No circularity: direct numerical lattice computation of the EOS via complex Langevin
full rationale
The paper reports continuum-extrapolated lattice results for the baryon density and pressure obtained by direct stochastic sampling with the complex Langevin equation. Thermodynamic relations are used to obtain the pressure from the density (standard integration dP = n_B dμ_B), but this is an exact identity applied to independently computed observables rather than a self-referential fit or redefinition. Validation against lower-μ lattice data and high-T perturbation theory is presented as consistency check, not as load-bearing input. No derivation step reduces by construction to fitted parameters, self-citations, or ansätze imported from prior work by the same authors; the central output is the numerical evaluation itself under stated convergence diagnostics.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Complex Langevin dynamics converges to the correct complex measure for the QCD action at finite density
- standard math Continuum limit exists and can be reached by extrapolation from finite lattice spacings
Forward citations
Cited by 2 Pith papers
-
Lattice field theories with a sign problem
A review of holomorphic extensions, dual variables, tensor renormalization group, and machine learning approaches for controlling the sign problem in lattice field theories.
-
Lattice field theories with a sign problem
Reviews approaches such as Lefschetz thimbles, complex Langevin dynamics, dual variables, tensor renormalization group, and machine learning to control the sign problem in lattice field theories.
Reference graph
Works this paper leans on
-
[1]
within the region of the phase diagram that is ac- cessible to both studies. In addition, for smallµ B/Twe find agreement with hard-thermal-loop (HTL) perturba- tive results derived in [59–61], as well as with the baryon density of a theory of noninteracting quarks and gluons, only the former of which contribute ton B. It should be noted that we have naiv...
-
[2]
and find no disagreement. Furthermore, as before 4 0 2 4 6 8 10 12 14 16 µB/T 0 10 20 30 40 50 60p/T 4 Free theory T = 600 MeV HTL T = 600 MeV T = 600 MeV T = 360 MeV T = 275 MeV T = 260 MeV Saturation effects 0 1 2 3 4 5 3 4 5 6 7 8 Bors´ anyi et al. T = 275 MeV Bors´ anyi et al. T = 260 MeV FIG. 2: Similar as Fig. 1, but for the pressure. we provide comp...
2000
-
[3]
In [60], HTL results for the pressure difference ∆pas a function ofTwere compared with (older) lattice data for different small chemical potentials and good agree- ment was found
in our entire temperature range using (10) and (11). In [60], HTL results for the pressure difference ∆pas a function ofTwere compared with (older) lattice data for different small chemical potentials and good agree- ment was found. A similar comparison of our results for ∆pwith HTL perturbation theory and a free-theory cal- culation but for larger chemic...
-
[4]
of machine-learning so-called kernels that could sta- bilize simulations in the low-temperature regime of QCD like they have in other systems [66–70]. Further inter- esting research directions include the addition of electro- magnetic fields [71, 72] and/or an isospin chemical po- tential [73]. We would like to thank Szabolcs Bors´ anyi and Jana G¨ unther...
-
[5]
title The QCD phase diagram http://arxiv.org/abs/2512.08843 arXiv:2512.08843
S. Bors´ anyi and P. Parotto, The QCD phase diagram (2025), arXiv:2512.08843 [nucl-th]
-
[6]
Y. Aoki, G. Endr˝ odi, Z. Fodor, S. D. Katz, and K. K. Szab´ o, Nature443, 675 (2006), arXiv:hep-lat/0611014
work page Pith review arXiv 2006
-
[7]
Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
M. Troyer and U.-J. Wiese, Phys. Rev. Lett.94, 170201 (2005), arXiv:cond-mat/0408370
work page Pith review arXiv 2005
- [8]
-
[9]
S. Bors´ anyi, G. Endr˝ odi, Z. Fodor, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szab´ o, JHEP2012(053), arXiv:1204.6710 [hep-lat]. 6
-
[10]
The QCD Equation of State toO(µ 6 B) from Lattice QCD,
A. Bazavov, H.-T. Ding, P. Hegde, O. Kaczmarek, F. Karsch, E. Laermann, Y. Maezawa, S. Mukherjee, H. Ohno, P. Petreczky, H. Sandmeyer, P. Steinbrecher, C. Schmidt, S. Sharma, W. Soeldner, and M. Wagner, Phys. Rev. D95, 054504 (2017), arXiv:1701.04325 [hep- lat]
- [11]
-
[12]
Bazavovet al.(HotQCD), Chiral crossover in QCD at zero and non-zero chemical potentials, Phys
A. Bazavov et al. (HotQCD Collaboration), Phys. Lett. B795, 15 (2019), arXiv:1812.08235 [hep-lat]
-
[13]
M. Giordano, K. Kapas, S. D. Katz, D. Nogradi, and A. Pasztor, JHEP2020(088), arXiv:2004.10800 [hep- lat]
-
[14]
S. Bors´ anyi, Z. Fodor, J. N. Guenther, R. Kara, S. D. Katz, P. Parotto, A. P´ asztor, C. Ratti, and K. K. Szab´ o, Phys. Rev. Lett.126, 232001 (2021), arXiv:2102.06660 [hep-lat]
-
[15]
S. Bors´ anyi, J. N. Guenther, R. Kara, Z. Fodor, P. Parotto, A. P´ asztor, C. Ratti, and K. Szab´ o, Phys. Rev. D105, 114504 (2022), arXiv:2202.05574 [hep-lat]
-
[16]
S. Bors´ anyi, Z. Fodor, M. Giordano, S. D. Katz, D. N´ ogradi, A. P´ asztor, and C. H. Wong, Phys. Rev. D105, L051506 (2022), arXiv:2108.09213 [hep-lat]
-
[17]
D. Bollweg et al. (HotQCD Collaboration), Phys. Rev. D105, 074511 (2022), arXiv:2202.09184 [hep-lat]
-
[18]
S. Bors´ anyi, Z. Fodor, M. Giordano, J. N. Guenther, S. D. Katz, A. P´ asztor, and C. H. Wong, Phys. Rev. D107, L091503 (2023), arXiv:2208.05398 [hep-lat]
-
[19]
D. Bollweg et al. (HotQCD Collaboration), Phys. Rev. D108, 014510 (2023), arXiv:2212.09043 [hep-lat]
-
[20]
S. Bors´ anyi, Z. Fodor, J. N. Guenther, P. Parotto, A. P´ asztor, L. Pirelli, K. K. Szab´ o, and C. H. Wong, Phys. Rev. D110, 114507 (2024), arXiv:2410.06216 [hep- lat]
- [21]
-
[22]
Parisi, Phys
G. Parisi, Phys. Lett. B131, 393 (1983)
1983
-
[23]
J. R. Klauder, J. Phys. A16, L317 (1983)
1983
-
[24]
Parisi and Y.-S
G. Parisi and Y.-S. Wu, Sci. Sin.24, 483 (1981)
1981
- [25]
- [26]
- [27]
-
[28]
Sexty,Simulating full QCD at nonzero density using the complex Langevin equation,Phys
D. Sexty, Phys. Lett B.729, 108 (2014), arXiv:1307.7748 [hep-lat]
- [29]
- [30]
-
[31]
Z. Fodor, S. D. Katz, D. Sexty, and T¨ or¨ ok, Phys. Rev. D 92, 094516 (2015), arXiv:1508.05260 [hep-lat]
- [32]
- [33]
-
[34]
D. Sexty, Phys. Rev. D100, 074503 (2019), arXiv:1907.08712 [hep-lat]
- [35]
- [36]
-
[37]
M. Scherzer, D. Sexty, and I.-O. Stamatescu, Phys. Rev. D102, 014515 (2020), arXiv:2004.05372 [hep-lat]
-
[38]
F. Attanasio, B. J¨ ager, and F. P. G. Ziegler, Eur. Phys. J. A56, 251 (2020), arXiv:2006.00476 [hep-lat]
-
[39]
F. Attanasio, B. J¨ ager, and F. P. G. Ziegler, QCD equa- tion of state via the complex Langevin method (2022), arXiv:2203.13144 [hep-lat]
- [40]
-
[41]
S. Tsutsui, Y. Asano, Y. Ito, H. Matsufuru, Y. Namekawa, J. Nishimura, S. Shimasaki, and A. Tsuchiya, JHEP2025(108), arXiv:2505.06551 [hep- lat]
- [42]
- [43]
-
[44]
G. Aarts, F. A. James, E. Seiler, and I.-O. Stamatescu, Eur. Phys. J. C71, 1756 (2011), arXiv:1101.3270 [hep- lat]
-
[45]
K. Nagata, J. Nishimura, and S. Shimasaki, Phys. Rev. D94, 114515 (2016), arXiv:1606.07627 [hep-lat]
work page Pith review arXiv 2016
-
[46]
M. Scherzer, E. Seiler, D. Sexty, and I.-O. Stamatescu, Phys. Rev. D99, 014512 (2019), arXiv:1808.05187 [hep- lat]
-
[47]
M. Scherzer, E. Seiler, D. Sexty, and I.-O. Stamatescu, Phys. Rev. D101, 014501 (2020), arXiv:1910.09427 [hep- lat]
- [48]
- [49]
-
[50]
A. Joseph and A. Kumar, Thermodynamic Diagnostics for Complex Langevin Simulations: The Role of Configu- rational Temperature (2025), arXiv:2509.08287 [hep-lat]
-
[51]
Correctness criteria for complex Langevin
M. Mandl, Correctness criteria for complex Langevin (2026), arXiv:2604.12388 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[52]
Philipsen, The QCD equation of state from the lattice, Prog
O. Philipsen, Prog. Part. Nucl. Phys.70, 55 (2013), arXiv:1207.5999 [hep-lat]
-
[53]
Mandl, D
M. Mandl, D. Sexty, and D. Unterhuber, in preparation
-
[54]
Ukawa and M
A. Ukawa and M. Fukugita, Phys. Rev. Lett.55, 1854 (1985)
1985
-
[55]
G. Aarts, F. A. James, E. Seiler, and I.-O. Stamatescu, Phys. Lett. B687, 154 (2010), arXiv:0912.0617 [hep-lat]
-
[56]
Analytic Smearing of SU(3) Link Variables in Lattice QCD
C. Morningstar and M. Peardon, Phys. Rev. D69, 054501 (2004), arXiv:hep-lat/0311018
work page Pith review arXiv 2004
-
[57]
R. Bellwied, S. Bors´ anyi, Z. Fodor, S. D. Katz, A. P´ asztor, C. Ratti, and K. K. Szab´ o, Phys. Rev. D 92, 114505¸ (2015), arXiv:1507.04627 [hep-lat]
-
[58]
Properties and uses of the Wilson flow in lattice QCD
M. L¨ uscher, JHEP2010(071), arXiv:1006.4518 [hep-lat]
-
[59]
S. Bors´ anyi, S. D¨ urr, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, T. Kurth, L. Lellouch, T. Lippert, C. McNeile, and K. K. Szab´ o, JHEP2012(010), arXiv:1203.4469 [hep-lat]
- [60]
-
[61]
M. Mandl, D. Sexty, and D. Unterhuber, Zenodo 10.5281/zenodo.19676680 (2026). 7
-
[62]
M. Mandl, D. Sexty, and D. Unterhuber, Zenodo 10.5281/zenodo.19676547 (2026)
- [63]
- [64]
- [65]
-
[66]
S. Borsanyi, Z. Fodor, J. Guenther, K.-H. Kampert, S. D. Katz, T. Kawanai, T. G. Kovacs, S. W. Mages, A. Pasz- tor, F. Pittler, J. Redondo, A. Ringwald, and K. K. Sz- abo, Nature539, 69 (2016), arXiv:1606.07494 [hep-lat]
- [67]
-
[68]
J. Nishimura and S. Shimasaki, Phys. Rev. D92, 011501(R) (2015), arXiv:1504.08359 [hep-lat]
- [69]
- [70]
-
[71]
D. Alvestad, R. Larsen, and A. Rothkopf, JHEP2023 (057), arXiv:2211.15625 [hep-lat]
-
[72]
D. Alvestad, A. Rothkopf, and D. Sexty, Phys. Rev. D 109, L031502 (2024), arXiv:2310.08053 [hep-lat]
-
[73]
K. Boguslavski, P. Hotzy, and D. I. M¨ uller, JHEP2023 (011), arXiv:2212.08602 [hep-lat]
-
[74]
K. Boguslavski, P. Hotzy, and D. I. M¨ uller, Phys. Rev. D109, 094518 (2024), arXiv:2312.03063 [hep-lat]
- [75]
-
[76]
G. Endr˝ odi and G. Mark´ o, Phys. Rev. D109, 034506 (2024), arXiv:2309.07058 [hep-lat]
- [77]
-
[78]
van Rossum and J
G. van Rossum and J. de Boer, CWI Quarterly4, 283 (1991)
1991
-
[79]
van der Walt, S
S. van der Walt, S. C. Colbert, and G. Varoquaux, Com- put. Sci. Eng.13, 22 (2011)
2011
-
[80]
J. D. Hunter, Comput. Sci. Eng.9, 90 (2007)
2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.