REVIEW 3 major objections 4 minor 19 references
Equation of state of isospin asymmetric QCD with small baryon chemical potentials
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read This paper reports the first lattice QCD equation of state along the electric charge chemical potential axis, obtained by expanding the pressure in a leading-order Taylor series around nonzero isospin chemical potential.
desk verdict First lattice pressure on the mu_Q axis via Taylor expansion around nonzero isospin, but the BEC-phase numbers rest on an uncontrolled spline derivative. 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 identity is $c_{LL} = c_{II}$ between the connected parts of the second-order Taylor coefficients in the $\mu_L$ and $\mu_I$ directions, which follows from the trace representation of the coefficients. Using this identity, the coefficient $\chi_2^L$ is computed as $\chi_2^I$ (obtained as $\partial n_I / \partial \mu_I$ at zero pion source, $\lambda = 0$), plus the difference of the disconnected contributions of $\chi_2^I$ and $\chi_2^L$. Only the latter needs a $\lambda$-extrapolation, which greatly reduces uncertainties inside the BEC phase. The expansion itself is the leading-order Taylor series in $\mu_L$ and $\mu_s$ around simulation points on the isospin axis, with coefficients $\chi_2^L$, $\chi_2^s$ and $\chi_{11}^{Ls}$ interpolated in $T$ and $\mu_I$ by a two-dimensional spline and a Silver-Blaze boundary condition at $T=0$.
What would settle it
Compute the pressure on the charge axis to next-to-leading order in the same Taylor expansion (order $\mu^4$ terms) or by an independent method such as reweighting or imaginary chemical potential; if the difference from the leading-order result exceeds the statistical errors at $\mu_Q/m_\pi = 1.5$, the leading-order equation of state shown here is not valid up to that point.
Extended reading notes
Core claim
The central claim is that the QCD equation of state at pure charge chemical potential can be obtained from the leading-order Taylor expansion around the isospin axis, despite the pion-condensed BEC phase where standard expansions around zero chemical potential fail. To make this possible, the authors show that the coefficient $\chi_2^L$ can be computed reliably inside the BEC phase by exploiting the identity $c_{LL} = c_{II}$, which lets them obtain the connected contribution from $\chi_2^I = \partial n_I / \partial \mu_I$ evaluated directly at vanishing pion source via spline interpolation, instead of extrapolating the noisy fully connected trace. They present the resulting pressure $p/T^4$ on the $\mu_Q$ axis for $T \approx 123$–$165$ MeV and $\mu_Q/m_\pi$ up to $\sim 1.5$, with the largest deviation from the isospin-axis pressure deep inside the BEC phase. The expansion is valid only as long as leading order is sufficient, and the paper explicitly notes it will break down when expanding through a phase boundary.
Load-bearing premise
The leading-order Taylor expansion in the baryon and strangeness directions around the isospin axis is accurate enough up to $\mu_Q/m_\pi \sim 1.5$ that omitting fourth-order and higher terms does not change the pressure meaningfully; the expansion breaks down at a phase boundary.
Editorial extensions
If this is right
- The charge-axis equation of state can serve as input for early-Universe models with large lepton flavour asymmetries, where the trajectory runs near the $\mu_Q$ axis.
- The Taylor coefficients at nonzero isospin chemical potential open a route to the full three-dimensional parameter space of light-quark chemical potentials in the vicinity of the isospin axis.
- Inside the pion-condensed BEC phase, the density-improvement method yields significant results where the standard improved observable is too noisy; the same technique can be applied to other Taylor coefficients.
- The equation of state along the charge axis differs most strongly from the isospin axis deep in the BEC phase, indicating that charge chemical potential effects are not negligible there.
- The Silver-Blaze boundary condition at $T=0$ provides a useful constraint for future interpolations of the Taylor coefficients.
Reading between the lines
- The identity $c_{LL}=c_{II}$ likely generalizes to higher-order Taylor coefficients, so future work could compute light-quark coefficients at $\lambda=0$ from improved density derivatives, cutting the dominant systematic of the $\lambda$ extrapolation.
- A natural test is to compare the leading-order charge-axis pressure with a next-to-leading-order ($\mu^4$) calculation; if the difference is within errors up to $\mu_Q/m_\pi\sim1.5$, the expansion window is established empirically, and the breakdown near the phase boundary can be mapped.
- Because the leading-order expansion is expected to fail through a phase boundary, the same framework with imaginary chemical potentials could probe the radius of convergence and locate the transition on the charge axis.
- The charge-axis equation of state could be plugged into cosmic-QCD transition codes to compute gravitational-wave signatures of pion condensation; the quantitative impact on those signatures is testable once the equation-of-state table is released.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript extends the lattice QCD equation of state for isospin-asymmetric matter to small baryon and strangeness chemical potentials by Taylor expanding around simulation points on the pure isospin axis. To control the lambda->0 extrapolation for the leading baryon susceptibility chi_L^2, the authors use the exact identity c_LL=c_II and compute chi_I^2 as a numerical derivative of the spline-interpolated isospin density from Ref. [17], which they call the density-improved estimator. They then combine two-dimensional spline interpolations of chi_L^2, chi_s^2, and chi_Ls^11 with the isospin-axis EoS to present first results for the pressure on the pure charge chemical potential axis, p/T^4 at T~123-165 MeV and mu_Q/m_pi up to about 1.5, including inside the pion BEC phase.
Significance. If the reported results are correct, this is the first lattice QCD equation of state on the mu_Q axis in a regime relevant for early-Universe scenarios with large lepton flavour asymmetries, and it demonstrates a useful technique for expanding around non-zero isospin chemical potentials despite the complex-action problem. The algebraic reduction of the leading baryon susceptibility to the isospin susceptibility, Eq. (16), is clean and the improvement plots in Figs. 1 and 2 show a genuine reduction of uncertainties in the BEC phase. However, the central numerical result inherits uncontrolled systematics from a spline derivative, from the leading-order truncation, and from the absence of a continuum extrapolation; these need to be quantified before the result can be used as a quantitative prediction.
major comments (3)
- [Sec. 2.3, Eq. (16)] The density-improved chi_L^2 is obtained by substituting chi_I^2 = d n_I / d mu_I, evaluated as a numerical derivative of a spline interpolation of the isospin density from Ref. [17]. No systematic error for the spline fit or for the derivative is quoted, no comparison with a direct lambda=0 estimate of chi_I^2 or with the connected part c_II is shown, and the spline is most strained exactly near the BEC boundary mu_I = m_pi/2, where n_I has its strongest mu_I dependence and where the new results in Fig. 4 rely on it most. Because this derivative enters the O(mu^2) coefficient itself, the statement in the Conclusions that the Taylor expansion is only leading order does not cover a possible bias in the spline derivative. Please add a systematic error estimate or an explicit cross-check for this derivative.
- [Sec. 3, Fig. 4] The charge-axis pressure is presented at fixed lattice spacing without a continuum extrapolation or a comparison between the available lattice spacings (e.g., 24^3 x 6 and 24^3 x 8 used elsewhere in the paper). Since the result is advertised as the first lattice QCD EoS on the mu_Q axis and is aimed at cosmological applications, an estimate of discretization effects is necessary before this can be considered a quantitative result.
- [Sec. 3, Eq. (3)] The expansion is truncated at leading order, O(mu^2), but no higher-order coefficients or radius-of-convergence estimate are given for the mu_Q axis. At the largest mu_Q/m_pi=1.5 shown in Fig. 4, the offsets from the isospin axis are mu_L/m_pi=0.25 and mu_s/m_pi=-0.5; the paper does not demonstrate that the neglected O(mu^4) terms are small in this region. The qualitative caveat in the Conclusions is appropriate but is not a substitute for a quantitative convergence check if the Fig. 4 results are to be used as an EoS.
minor comments (4)
- [Conclusions] The sentence 'valid as long the expansion to this order is sufficient' should read 'valid as long as the expansion to this order is sufficient'.
- [Fig. 2 caption] The labels in the caption ('standard impr.' and 'improved') do not match the labels in the text and main body ('standard impr.' and 'density impr.'); please harmonize them so the figure is self-contained.
- [Fig. 1, left panel] The y-axis label uses a ratio of expectation values, but the text describes the improvement term as a difference; please define the normalization explicitly and state that the plotted quantity is negative, as shown in the figure.
- [Sec. 2.2] The notation chi_L^2 and chi_s^2 in Eq. (4) is inconsistent with the more common chi_2^L and chi_2^s notation; a brief note defining the index convention would help the reader.
Circularity Check
No significant circularity: the charge-axis EoS is a new combination of newly computed Taylor coefficients and the prior isospin-axis EoS; Eq. (16) is an exact rearrangement, not a fit to the target.
full rationale
The derivation chain is self-contained. The new pressure on the mu_Q axis in Fig. 4 is obtained by substituting Eq. (17) into the Taylor expansion Eq. (3), using newly computed coefficients chi_L^2, chi_s^2, and chi_Ls^11 together with the prior isospin-axis EoS from Ref. [17]. The density-improved extraction of chi_L^2 via Eq. (16) uses the exact identity c_LL = c_II from Eq. (15); it replaces one connected part with another, not with the target pressure. The spline derivative of n_I is a numerical estimator for chi_I^2, which is an input coefficient, not the predicted quantity; any spline bias is a systematic-error concern, not circularity. Self-citations to Refs. [7], [16], and [17] supply prior data and methodology, but none of these prior works presupposes the charge-axis result. Thus no prediction reduces to its own inputs by construction.
Assumptions & free parameters
free parameters (2)
- Spline fit for isospin density derivative (chi_I2 = dn_I/dmu_I) =
not reported (spline knot and smoothing parameters not given)
- Two-dimensional spline interpolation of Taylor coefficients for pressure integration =
not reported (spline coefficients not tabulated)
assumptions (4)
- domain assumption Rooted improved staggered quarks at physical quark masses and finite lattice spacing reproduce continuum QCD in the simulated temperature and chemical potential range.
- domain assumption Silver-Blaze property: at T=0 the Taylor coefficients vanish for all mu_I because the lightest excitations (neutron, kaon) have nonzero mass.
- domain assumption The leading-order Taylor expansion in mu_L and mu_s around the isospin axis is sufficient on the mu_Q axis up to the largest mu_Q considered.
- domain assumption The isospin density data from Ref. [17] and its spline interpolation are sufficiently precise and smooth to yield a reliable numerical derivative chi_I2 = dn_I/dmu_I.
Cite this review
Pith. "Pith review of Equation of state of isospin asymmetric QCD with small baryon chemical potentials." pith.science (2026). https://pith.science/paper/HRH5WVUX
@misc{pith2026241112918,
author = {Pith},
title = {Pith review of: Equation of state of isospin asymmetric QCD with small baryon chemical potentials},
year = {2026},
howpublished = {\url{https://pith.science/paper/HRH5WVUX}},
note = {Machine review of arXiv:2411.12918}
}
read the original abstract
We extend our measurement of the equation of state of isospin asymmetric QCD to small baryon and strangeness chemical potentials, using the leading order Taylor expansion coefficients computed directly at non-zero isospin chemical potentials. Extrapolating the fully connected contributions to vanishing pion sources is particularly challenging, which we overcome by using information from isospin chemical potential derivatives evaluated numerically. Using the Taylor coefficients, we present, amongst others, first results for the equation of state along the electric charge chemical potential axis, which is potentially of relevance for the evolution of the early Universe at large lepton flavour asymmetries.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[17]
B. B. Brandt, F. Cuteri and G. Endrődi,Equation of state and speed of sound of isospin-asymmetric QCD on the lattice, JHEP07(2023) 055 [2212.14016]
arXiv 2023
-
[1]
M. M. Wygas, I. M. Oldengott, D. Bödeker and D. J. Schwarz,Cosmic QCD Epoch at Nonvanishing Lepton Asymmetry, Phys. Rev. Lett.121 (2018) 201302 [1807.10815]
arXiv 2018
-
[2]
M. M. Middeldorf-Wygas, I. M. Oldengott, D. Bödeker and D. J. Schwarz,The cosmic QCD transition for large lepton flavour asymmetries, 2009.00036
arXiv 2009
-
[3]
V. Vovchenko, B. B. Brandt, F. Cuteri, G. Endrődi, F. Hajkarim and J. Schaffner-Bielich, Pion Condensation in the Early Universe at Nonvanishing Lepton Flavor Asymmetry and Its Gravitational Wave Signatures, Phys. Rev. Lett.126 (2021) 012701 [2009.02309]
arXiv 2021
-
[4]
D. T. Son and M. A. Stephanov,QCD at finite isospin density, Phys. Rev. Lett.86(2001) 592 [hep-ph/0005225]
arXiv 2001
-
[5]
J. B. Kogut and D. K. Sinclair,Quenched lattice QCD at finite isospin density and related theories, Phys. Rev. D66(2002) 014508 [hep-lat/0201017]
arXiv 2002
-
[6]
J. B. Kogut and D. K. Sinclair,Lattice QCD at finite isospin density at zero and finite temperature, Phys. Rev. D66(2002) 034505 [hep-lat/0202028]
arXiv 2002
-
[7]
B. B. Brandt, G. Endrődi and S. Schmalzbauer,QCD phase diagram for nonzero isospin-asymmetry, Phys. Rev. D97(2018) 054514 [1712.08190]
arXiv 2018
Show all 19 references
-
[8]
B. B. Brandt and G. Endrődi,Reliability of Taylor expansions in QCD,Phys. Rev. D99 (2019) 014518 [1810.11045]
2019 arXiv
-
[9]
B. B. Brandt, F. Cuteri, G. Endrődi and S. Schmalzbauer,The Dirac spectrum and the BEC-BCS crossover in QCD at nonzero isospin asymmetry,Particles3 (2020) 80 [1912.07451]
2020 arXiv
-
[10]
Cuteri, B
F. Cuteri, B. B. Brandt and G. Endrődi,Searching for the BCS phase at nonzero isospin asymmetry, PoSLATTICE2021(2022) 232 [2112.11113]
2022 arXiv
-
[11]
B. B. Brandt, V. Chelnokov, F. Cuteri and G. Endrődi,Pion condensation at lower than physical quark masses, PoSLATTICE2022(2023) 146 [2301.08607]
2023 arXiv
-
[12]
B. B. Brandt, G. Endrődi and S. Schmalzbauer,QCD at finite isospin chemical potential,EPJ Web Conf.175 (2018) 07020 [1709.10487]
2018 arXiv
-
[13]
B. B. Brandt, G. Endrődi, E. S. Fraga, M. Hippert, J. Schaffner-Bielich and S. Schmalzbauer, New class of compact stars: Pion stars,Phys. Rev. D98(2018) 094510 [1802.06685]
2018 arXiv
-
[14]
B. B. Brandt, G. Endrődi and S. Schmalzbauer,QCD at nonzero isospin asymmetry, PoS Confinement2018 (2018) 260 [1811.06004]
2018 arXiv
-
[15]
B. B. Brandt, F. Cuteri and G. Endrődi,QCD thermodynamics at non-zero isospin asymmetry, PoSLATTICE2021(2022) 132 [2110.14750]. 9 EoS of isospin asymmetric QCD with small baryon chemical potentials Bastian B. Brandt
2022 arXiv
-
[16]
B. B. Brandt, F. Cuteri and G. Endrődi,Equation of state and Taylor expansions at nonzero isospin chemical potential,PoSLATTICE2022(2023) 144 [2212.01431]
2023 arXiv
-
[18]
S. A. Gottlieb, W. Liu, D. Toussaint, R. L. Renken and R. L. Sugar,Fermion Number Susceptibility in Lattice Gauge Theory,Phys. Rev. D38 (1988) 2888
1988
-
[19]
Adhikari and J
P. Adhikari and J. O. Andersen,Pion and kaon condensation at zero temperature in three-flavor𝜒PPT at nonzero isospin and strange chemical potentials at next-to-leading order,JHEP 06(2020) 170 [1909.10575]. 10
2020 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.