REVIEW 3 major objections 4 minor 2 cited by
Dense and magnetized QCD from imaginary chemical potential
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read In QCD with background magnetic fields, the leading dense-matter pressure response becomes non-monotonic in temperature once eB reaches about 0.5 GeV², along the strangeness-neutral, isospin-asymmetric line relevant to heavy-ion collisions.
desk verdict A transparent, well-executed proceedings paper whose headline claim outruns the data: the real-axis non-monotonicity rests on an untested linear-in-mu_B^2 fit and a single lattice spacing. 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 observable carrying the argument is the total derivative dP/dμ_B along the constrained trajectory, Eq. (7), which combines the baryon, charge, and strangeness densities with the constrained derivatives of the charge and strangeness chemical potentials. The analytic continuation is done by a multidimensional spline surface in T, B, and μ_B² that is linear in μ_B², so that data taken at imaginary μ_B (j = 0, 3, 4, 5, with μ_B = i π j/8 T) can be extended to real μ_B; node points are stochastically generated and weighted by the Akaike criterion to estimate systematic error.
What would settle it
Compute the O(μ_B⁴) coefficient of the pressure expansion at eB = 0.5 GeV² by adding a simulation at j = 6 or by including a μ_B⁴ term in the spline fit, and check whether the linear-in-μ_B² ansatz is responsible for the turnover; if the curvature term is significant at the real chemical potentials of interest, the non-monotonic peak will shift or vanish. Repeating the measurement at a finer lattice spacing also settles the continuum question, since the current result is at a single lattice spacing and the authors state that a continuum extrapolation is still needed.
Extended reading notes
Core claim
The central claim is that the leading-order dense QCD equation of state, measured by μ_B⁻¹ dP/dμ_B at imaginary chemical potential and analytically continued to real baryon chemical potential, develops a peak and then decreases as temperature rises when the background magnetic field exceeds about 0.5 GeV², so that its temperature dependence becomes non-monotonic. The paper presents this as the main result in Section 6, tied to the crossover region, and notes that a continuum extrapolation is still needed for phenomenological application to heavy-ion collisions.
Load-bearing premise
The extrapolation from imaginary to real chemical potential assumes that the quantity μ_B⁻¹ dP/dμ_B is a purely linear function of μ_B² across the whole range of real chemical potentials of interest; if higher powers of μ_B² matter at those real values, the predicted non-monotonicity is an artifact of the fit rather than a genuine QCD effect.
Editorial extensions
If this is right
- The equation of state of QCD in strong magnetic fields cannot be approximated by the zero-field result in the crossover region; the magnetic field changes the leading density response by more than the statistical errors of this calculation.
- Heavy-ion phenomenology that uses magnetized matter must treat the baryon-density response as non-monotonic in temperature, which will affect hydrodynamic evolution and observables such as directed flow.
- The non-monotonic behavior of μ_B⁻¹ dP/dμ_B at eB ≥ 0.5 GeV² could be a precursor of critical structure expected at larger fields, consistent with previous suggestions of a critical endpoint in the T–B plane.
- Simulations at nonzero magnetic field require careful tuning of μ_Q and μ_S to maintain strangeness neutrality; the paper demonstrates a two-stage procedure (estimate then linear shift) that keeps corrections within errors.
Reading between the lines
- If the non-monotonicity survives the continuum limit, the quark-gluon plasma produced in peripheral heavy-ion collisions may have a density response that changes sign with temperature at fixed beam energy, a feature that could be probed by comparing directed-flow measurements across collision centralities.
- A natural next check is to compute the full μ_B⁴ coefficient of the pressure expansion with the magnetic field turned on; if it grows with B, the linear-in-μ_B² continuation used here systematically underestimates the curvature on the real chemical potential axis.
- The same imaginary-chemical-potential machinery could map where the non-monotonic region begins in the B–T plane, effectively locating a band in which the magnetic field changes the qualitative shape of the dense-matter equation of state before any critical point is reached.
- A hadron resonance gas with magnetic-field-dependent masses could test whether the same non-monotonicity arises from the spin couplings of protons and neutrons, which would distinguish a low-temperature hadronic explanation from a quark-gluon plasma effect.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports lattice QCD results for the leading-order behavior of the dense QCD equation of state in the presence of background magnetic fields. The authors use 2+1+1 flavors of stout-smeared staggered fermions at the physical point on an N_t=8 lattice, with imaginary baryon chemical potentials μ_B = iπj/8 (j = 0,3,4,5), magnetic field strengths eB = 0, 0.3, 0.5, 0.8 GeV^2, and temperatures T = 135–200 MeV. They impose strangeness neutrality and an isospin asymmetry n_Q/n_B = 0.4, determining the required μ_Q and μ_S by a combination of algebraic Taylor expansions and fits, with a linear correction to enforce exact strangeness neutrality. The main observable is μ_B^{-1} dP/dμ_B along this constrained line, which is fitted with a multidimensional spline that is polynomial in T and eB and linear in μ_B^2. The authors report that this quantity develops a non-monotonic temperature dependence for eB ≳ 0.5 GeV^2, and conclude that magnetic fields strongly affect the dense QCD equation of state near the crossover. The paper explicitly identifies the lack of a continuum extrapolation and the need for a full analytic continuation as future work.
Significance. If the reported strong magnetic-field dependence and the non-monotonicity survive a continuum extrapolation and a controlled real-chemical-potential continuation, the result would be an important input for modeling heavy-ion collisions and for mapping the QCD phase diagram in the T-μ_B-B space. The paper is valuable as a first 2+1+1-flavor lattice study along a strangeness-neutral, isospin-asymmetric trajectory at imaginary chemical potential in a magnetic field, and the methods for imposing the experimental constraints are carefully cross-checked (algebraic, fit, and spline-based procedures). The authors are also transparent about the preliminary nature of the analysis, stating that the spline functions are preparatory for a future analytic continuation and that continuum extrapolation is still needed. The main weaknesses are that the abstract and title overstate what has actually been computed, and that the only evidence for the real-axis behavior rests on an untested linear-in-μ_B^2 assumption at a single lattice spacing.
major comments (3)
- [Abstract; Secs. 5–6] The abstract states that "we computed the equation of state of dense QCD" and that the results "suggest a strong change in the equation of state," but the body presents only the ratio μ_B^{-1} dP/dμ_B measured at imaginary μ_B, together with a preliminary spline that the Fig. 5 caption explicitly describes as something "that we will use in a future work to carry out the analytic continuation." No integrated pressure difference ΔP from Eq. (6) and no real-axis equation of state are shown. The claims should be reworded to refer to the leading-order coefficient of the EoS and to a preliminary indication, or the missing integrated quantity should be presented.
- [Sec. 5.2] The linear-in-μ_B^2 continuation is based on fits to j = 0,3,4,5, corresponding to μ_B^2/T^2 = 0, -1.39, -2.47, -3.85, and is then extrapolated to positive μ_B^2, where values of order 4–9 are reached for μ_B/T = 2–3. Since μ_B^{-1} dP/dμ_B = χ2 + (χ4/6)μ_B^2 + ···, the linear fit determines an effective slope that mixes χ4 with higher-order terms, and no estimate of the χ6μ_B^4 contribution is given. Removing j = 5 from the fit only tests stability inside the fitted interval, not the extrapolation domain. If the non-monotonicity claim is meant to apply on the real μ_B axis, this assumption must be tested (or at least quantified); otherwise the claim should be explicitly restricted to imaginary chemical potentials.
- [Sec. 2; Sec. 6] All numerical results are obtained on a single lattice spacing with N_t = 8, and the paper itself states that a continuum extrapolation is still needed. Because magnetic-field effects on the crossover and on thermodynamic quantities can be sensitive to the lattice cutoff, the statement that the EoS changes strongly with B should be labeled as a fixed-N_t result until the continuum limit is available. This is not a request for new simulations in a proceedings article, but the wording should not imply a continuum-physics prediction.
minor comments (4)
- [Abstract; Sec. 5.2] The abstract and the body are inconsistent about the status of the result: the abstract says the equation of state was computed, while Sec. 5.2 calls the spline a "preliminary determination" for future analytic continuation. Please harmonize the wording.
- [Throughout] The manuscript contains numerous typographical and spacing errors, including in the abstract (e.g., "Ourresults") and throughout the introduction and Sec. 4. A careful proofread is needed.
- [Sec. 4, Eq. (8)] The truncation order N in the Taylor expansion (8) is not specified explicitly in the text, and the statement "where a are the parameters we want to determine" should be made more precise by defining the coefficient set and how N is chosen in each step of the iterative tuning procedure.
- [Sec. 5.1, Fig. 4] The text refers to "Figs. 4a and 4b," but the figure panels are labeled (a) and (b); please verify the cross-referencing convention. Also, the captions for Figs. 2 and 4 are very long and would benefit from a short descriptive sentence followed by the details.
Circularity Check
No significant circularity: central lattice results are new direct measurements; self-citations are methodological, not load-bearing.
full rationale
The paper's central object, mu_B^{-1} dP/dmu_B, is a directly measured lattice observable at imaginary chemical potentials (Eq. (7), Sec. 5.2), and the claimed magnetic-field dependence and non-monotonicity are properties of a spline surface fitted to those new data. No fitted parameter is renamed as a prediction: the real-axis continuation is explicitly preliminary and is based on a linear-in-mu_B^2 ansatz, which is a model assumption rather than a definitional tautology. The strangeness-neutrality and isospin-asymmetry constraints are imposed by tuning mu_Q and mu_S using earlier simulation points (Sec. 4), with the final observables corrected by a small linear shift; this is an iterative procedure, not an identity that builds the answer into the input. Self-citations to Refs. [22], [30], and [31] provide the spline-fitting methodology and previous EoS context, but the specific B-dependence and non-monotonicity claims rest on the present simulation data, not on those citations. There is no imported uniqueness theorem, no ansatz smuggled in solely via self-citation, and no renaming of a known empirical pattern. Any concern about the linear extrapolation beyond the fitted imaginary-mu_B range is a correctness/robustness issue, not circularity. Thus the derivation chain is self-contained with respect to the new data, and the paper merits a low circularity score.
Assumptions & free parameters
free parameters (3)
- Taylor coefficients a_{2i+1} for mu_Q and mu_S =
not tabulated
- spline surface coefficients =
not tabulated
- model order choices (N in Eq. 8; polynomial orders in spline; inclusion or exclusion of j=5) =
N and orders chosen by hand; j=5 included or excluded for systematics
assumptions (7)
- domain assumption Rooted-staggered 2+1+1 stout-smeared fermions at the physical point can be used to compute QCD thermodynamics in a magnetic field.
- domain assumption The thermodynamic potential is analytic in mu_B between the imaginary axis and the real axis in the extrapolation region.
- ad hoc to paper mu_B^{-1} dP/dmu_B is linear in mu_B^2 in the range of the data and the extrapolation.
- domain assumption The heavy-ion trajectory is characterized by n_S = 0 and n_Q/n_B = 0.4.
- domain assumption Linear interpolation between integer magnetic flux quanta N_b yields the desired eB values.
- domain assumption A single lattice spacing (N_t = 8) is representative of continuum QCD.
- domain assumption The charm quark chemical potential is set to zero despite dynamical charm quarks.
Cite this review
Pith. "Pith review of Dense and magnetized QCD from imaginary chemical potential." pith.science (2026). https://pith.science/paper/VP5TIDOE
@misc{pith2026250201132,
author = {Pith},
title = {Pith review of: Dense and magnetized QCD from imaginary chemical potential},
year = {2026},
howpublished = {\url{https://pith.science/paper/VP5TIDOE}},
note = {Machine review of arXiv:2502.01132}
}
read the original abstract
In this work, we computed the equation of state of dense QCD in the presence of background magnetic fields using lattice QCD simulations at imaginary baryon chemical potential. Our simulations include 2+1+1 flavors of stout-smeared staggered fermions with masses at the physical point and a tree-level Symanzik-improved gauge action. Using several expansion schemes, we tuned our simulation parameters such that the equation of state satisfies strangeness neutrality and isospin asymmetry constraints, which are relevant to the phenomenology of heavy-ion collisions. Our results suggest a strong change in the equation of state due to the magnetic field, in particular, around the crossover temperature. A continuum extrapolation of our data is still needed for future applications of our equation of state to heavy-ion-collision phenomenology.
Figures
Figures from the paper (2 more)
Forward citations
Cited by 2 Pith papers
-
Leading-Order QCD Equation of State in Strong Magnetic Fields at Nonzero Baryon Chemical Potential
Continuum-estimated leading-order EoS coefficients in magnetized strangeness-neutral QCD at nonzero baryon chemical potential show temperature-band crossings in q1 and P2 and a possible sign change of the trace anomal...
-
QCD Equation of State with Strong Magnetic Fields and Nonzero Baryon Density
Lattice QCD continuum estimates of leading-order baryon-density Taylor coefficients of the magnetized QCD equation of state show deviations from hadron gas and approach to a free gas at strong fields.
Reference graph
Works this paper leans on
-
[22]
S.Borsányi,B.Brandt,G.Endrődi,J.N.Guenther,R.Kara,andA.D.MarquesValois, QCD equation of state in the presence of magnetic fields at low density, PoSLATTICE2023 (2024) 164, [2312.15118]
arXiv 2024
-
[1]
Vachaspati,Magnetic fields from cosmological phase transitions, Phys
T. Vachaspati,Magnetic fields from cosmological phase transitions, Phys. Lett. B265(1991) 258
1991
-
[2]
STAR Collaboration, M. I. Abdulhamidet al.,Observation of the electromagnetic field effect via charge-dependent directed flow in heavy-ion collisions at the Relativistic Heavy Ion Collider, Phys. Rev. X14(2024) 011028, [2304.03430]
arXiv 2024
-
[3]
D. E. Kharzeev, L. D. McLerran, and H. J. Warringa,The Effects of topological charge change in heavy ion collisions: ’Event by event P and CP violation’, Nucl. Phys. A803 (2008) 227, [0711.0950]
arXiv 2008
- [4]
-
[5]
Y. Aoki, G. Endrődi, Z. Fodor, S. D. Katz, and K. K. Szabó,The Order of the quantum chromodynamics transition predicted by the standard model of particle physics, Nature443 (2006) 675, [hep-lat/0611014]
arXiv 2006
-
[6]
W.-j. Fu, J. M. Pawlowski, and F. Rennecke,QCD phase structure at finite temperature and density,Phys. Rev. D101(2020) 054032, [1909.02991]
arXiv 2020
-
[7]
C. S. Fischer,QCD at finite temperature and chemical potential from Dyson–Schwinger equations,Prog. Part. Nucl. Phys.105 (2019) 1, [1810.12938]
arXiv 2019
Show all 33 references
-
[8]
Hippert, J
M. Hippert, J. Grefa, T. A. Manning, J. Noronha, J. Noronha-Hostler, I. Portillo Vazquez, C. Ratti, R. Rougemont, and M. Trujillo,Bayesian location of the QCD critical point from a holographic perspective, Phys. Rev. D110 (2024) 094006, [2309.00579]
2024 arXiv
-
[9]
C. E. Berger, L. Rammelmüller, A. C. Loheac, F. Ehmann, J. Braun, and J. E. Drut,Complex Langevin and other approaches to the sign problem in quantum many-body physics, Phys. Rept.892 (2021) 1, [1907.10183]
2021 arXiv
-
[10]
Roberge and N
A. Roberge and N. Weiss,Gauge Theories With Imaginary Chemical Potential and the Phases of QCD,Nucl. Phys. B275 (1986) 734
1986
-
[11]
Endrődi,Critical point in the QCD phase diagram for extremely strong background magnetic fields,JHEP07 (2015) 173, [1504.08280]
G. Endrődi,Critical point in the QCD phase diagram for extremely strong background magnetic fields,JHEP07 (2015) 173, [1504.08280]
2015 arXiv
-
[12]
T. D. Cohen and N. Yamamoto,New critical point for QCD in a magnetic field,Phys. Rev. D 89 (2014) 054029, [1310.2234]
2014 arXiv
-
[13]
V. V. Braguta, M. N. Chernodub, A. Y. Kotov, A. V. Molochkov, and A. A. Nikolaev, Finite-density QCD transition in a magnetic background field,Phys. Rev. D100 (2019) 114503, [1909.09547]
2019 arXiv
-
[14]
D’Elia, L
M. D’Elia, L. Maio, F. Sanfilippo, and A. Stanzione,Phase diagram of QCD in a magnetic background,Phys. Rev. D105 (2022) 034511, [2111.11237]
2022 arXiv
-
[15]
Zambello, M
K. Zambello, M. D’Elia, L. Maio, and G. Zanichelli,The Roberge-Weiss endpoint in (2+1)-flavor QCD with background magnetic fields,PoSLATTICE2024(2025) 169, [2412.06326]
2025 arXiv
-
[16]
Borsányi, G
S. Borsányi, G. Endrődi, Z. Fodor, A. Jakovác, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szabó, The QCD equation of state with dynamical quarks,JHEP11 (2010) 077, [1007.2580]
2010 arXiv
-
[17]
Bazavovet al., Equation of state in ( 2+1 )-flavor QCD, Phys
HotQCD Collaboration, A. Bazavovet al., Equation of state in ( 2+1 )-flavor QCD, Phys. Rev. D90(2014) 094503, [1407.6387]
2014 arXiv
-
[18]
Bazavovet al.,The QCD Equation of State toO(𝜇6 𝐵) from Lattice QCD,Phys
A. Bazavovet al.,The QCD Equation of State toO(𝜇6 𝐵) from Lattice QCD,Phys. Rev. D95 (2017) 054504, [1701.04325]. 12 Dense and magnetized QCD from imaginary chemical potential M. A. Petri and A. D. M. Valois
2017 arXiv
-
[19]
Borsanyi, J
S. Borsanyi, J. N. Guenther, R. Kara, Z. Fodor, P. Parotto, A. Pasztor, C. Ratti, and K. K. Szabo, Resummed lattice QCD equation of state at finite baryon density: Strangeness neutrality and beyond,Phys. Rev. D105 (2022) 114504, [2202.05574]
2022 arXiv
-
[20]
Astrakhantsev, V
N. Astrakhantsev, V. V. Braguta, A. Y. Kotov, and A. A. Roenko,QCD equation of state at nonzero baryon density in an external magnetic field,Phys. Rev. D109 (2024) 094511, [2403.07783]
2024 arXiv
-
[21]
Endrődi,QCD with background electromagnetic fields on the lattice: A review,Prog
G. Endrődi,QCD with background electromagnetic fields on the lattice: A review,Prog. Part. Nucl. Phys.141 (2025) 104153, [2406.19780]
2025 arXiv
-
[23]
Bellwied, S
R. Bellwied, S. Borsanyi, Z. Fodor, S. D. Katz, A. Pasztor, C. Ratti, and K. K. Szabo, Fluctuations and correlations in high temperature QCD, Phys. Rev. D92 (2015) 114505, [1507.04627]
2015 arXiv
-
[24]
Weisz,Continuum Limit Improved Lattice Action for Pure Yang-Mills Theory
P. Weisz,Continuum Limit Improved Lattice Action for Pure Yang-Mills Theory. 1.,Nucl. Phys. B212(1983) 1
1983
-
[25]
G. S. Bali, F. Bruckmann, G. Endrődi, Z. Fodor, S. D. Katz, S. Krieg, A. Schäfer, and K. K. Szabó, The QCD phase diagram for external magnetic fields,JHEP02 (2012) 044, [1111.4956]
2012 arXiv
-
[26]
Tiesinga, P
E. Tiesinga, P. J. Mohr, D. B. Newell, and B. N. Taylor,CODATA recommended values of the fundamental physical constants: 2018, Rev. Mod. Phys.93 (2021) 025010
2021
-
[27]
Endrődi,QCD equation of state at nonzero magnetic fields in the Hadron Resonance Gas model, JHEP04(2013) 023, [1301.1307]
G. Endrődi,QCD equation of state at nonzero magnetic fields in the Hadron Resonance Gas model, JHEP04(2013) 023, [1301.1307]
2013 arXiv
-
[28]
Vovchenko,Magnetic field effect on hadron yield ratios and fluctuations in a hadron resonance gas,Phys
V. Vovchenko,Magnetic field effect on hadron yield ratios and fluctuations in a hadron resonance gas,Phys. Rev. C110 (2024) 034914, [2405.16306]
2024 arXiv
-
[29]
Marczenko, M
M. Marczenko, M. Szymański, P. M. Lo, B. Karmakar, P. Huovinen, C. Sasaki, and K. Redlich,Magnetic effects in the hadron resonance gas,Phys. Rev. C110 (2024) 065203, [2405.15745]
2024 arXiv
-
[30]
B. B. Brandt and G. Endrődi,QCD phase diagram with isospin chemical potential,PoS LATTICE2016(2016) 039, [1611.06758]
2016 arXiv
-
[31]
B. B. Brandt, F. Cuteri, and G. Endrodi,Equation of state and speed of sound of isospin-asymmetric QCD on the lattice, JHEP07(2023) 055, [2212.14016]
2023 arXiv
-
[32]
Akaike,A new look at the statistical model identification,IEEE Transactions on Automatic Control19(1974) 716
H. Akaike,A new look at the statistical model identification,IEEE Transactions on Automatic Control19(1974) 716
1974
-
[33]
Ding, J.-B
H.-T. Ding, J.-B. Gu, A. Kumar, S.-T. Li, and J.-H. Liu,Baryon Electric Charge Correlation as a Magnetometer of QCD,Phys. Rev. Lett.132 (2024) 201903, [2312.08860]. 13
2024 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.