{"id":"09ca9f22-7b2d-4ecc-9c5e-dfaf14163f1c","arxiv_id":"2411.12918","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A Taylor expansion around nonzero isospin chemical potentials yields the first lattice QCD pressure estimates along the electric charge chemical potential axis.","lead":"The paper computes the equation of state of strongly interacting matter along the electric charge chemical potential axis using lattice QCD, for the first time from simulations carried out at nonzero isospin density. The result matters for early Universe scenarios with large lepton flavour asymmetries, where charge density can dominate the dynamics.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The charge-axis EoS inherits an uncontrolled systematic from the spline derivative used for chi_L^2 via Eq. (16); this is not covered by the paper's LO-truncation caveat.","rationale":"The Reader's stated weakest assumption is the sufficiency of the leading-order Taylor expansion, which the paper explicitly acknowledges in the Conclusions. That is a real limitation, but it is a range-of-validity caveat rather than a threat to the leading-order term itself. The more load-bearing point is the unvalidated numerical derivative in Eq. (16): chi_L^2 is the dominant new input for the charge-axis pressure, and inside the BEC phase it is obtained by differentiating a spline fit to n_I rather than by a direct, controlled lambda=0 measurement. The algebraic identity behind the density improvement is exact, and the use of the well-established Ref. [17] isospin EoS is a strength, but the spline derivative carries an unquantified systematic that propagates directly into Fig. 4. A single direct comparison on one BEC ensemble would settle whether this systematic is real. Since the Reader already calls for validating the spline derivative and quantifying truncation errors, the conditional verdict stands without adjustment.","tokens_in":8340,"tokens_out":13454,"duration_ms":146128,"concrete_test":"On a representative 24^3 x 8 ensemble inside the BEC phase (e.g., T = 123 MeV, mu_I/m_pi ~ 0.77), compute chi_I^2 directly at lambda = 0 via (T/V)[<c_II> + <n_I^2> - <n_I>^2] using the improved operators of Sec. 2.2, and compare it with the spline derivative d n_I / d mu_I used in Eq. (16). If the two agree within combined errors, the density-improved chi_L^2 and hence the Fig. 4 pressure are unaffected; if they disagree beyond errors, the charge-axis EoS needs a systematic error or a different estimator.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 2.3 the leading baryon-susceptibility coefficient chi_L^2 is obtained from Eq. (16), which replaces the connected part c_LL by c_II and evaluates chi_I^2 as the numerical derivative d n_I / d mu_I of a spline interpolation of the isospin density from Ref. [17]. The algebraic identity c_LL = c_II is exact, but the numerical derivative is not a controlled estimator: no spline-systematic error is quoted, no comparison with a direct lambda=0 estimate of chi_I^2 is shown, and no cross-check is given that the spline resolves the rapid mu_I-dependence of n_I near the BEC boundary at mu_I = m_pi/2. This derivative enters directly into the leading pressure correction on the charge axis in Fig. 4, and it is most relied upon exactly inside the BEC phase where the density has its strongest variation. A bias here is not mitigated by the Conclusions statement that the Taylor expansion is only leading order; the LO-truncation caveat concerns neglected O(mu^4) terms, while a spline-derivative bias would corrupt the O(mu^2) term itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":8609,"tokens_out":6037,"duration_ms":59722,"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":[{"comment":"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.","section":"Sec. 2.3, Eq. (16)"},{"comment":"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.","section":"Sec. 3, Fig. 4"},{"comment":"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.","section":"Sec. 3, Eq. (3)"}],"minor_comments":[{"comment":"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'.","section":"Conclusions"},{"comment":"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.","section":"Fig. 2 caption"},{"comment":"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.","section":"Fig. 1, left panel"},{"comment":"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.","section":"Sec. 2.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a proceedings contribution and the scope fits the venue. The main new result, the charge-axis EoS, is genuinely interesting and the algebraic construction is sound. My concerns are about the numerical systematics of the spline derivative and the lack of a continuum estimate; these are fixable in a revision and do not require new conceptual work. I do not see a novelty or citation-pattern issue."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe thing to know: this proceedings paper contains a genuinely new lattice result—the pressure on the pure electric charge chemical potential axis—obtained by Taylor expanding around nonzero isospin chemical potential rather than around zero. That axis is relevant for the early Universe with large lepton flavor asymmetries, and I don't know of a prior lattice calculation. The method is also notable: they use the exact relation c_LL = c_II to replace the hard connected part of chi_L^2 with a numerical derivative of the isospin density, computed directly at zero pion source. That reduces the statistical errors in the BEC phase, where the standard improvement was failing.\n\nWhat the paper does well: the identity is exact, the implementation is described in enough detail to reproduce, and the result in Fig. 4 is a real first. The authors are candid that this is leading order only and will break at a phase boundary.\n\nThe soft spots are real, though not disqualifying for a proceedings. The density-improved chi_L^2 values in the BEC phase depend on a spline interpolation of n_I from a previous paper, differentiated numerically. No uncertainty from the spline is quoted, and the paper shows no cross-check against a direct computation of chi_I^2 at vanishing pion source. The stress-test note makes a fair point: this is a systematic in the O(mu^2) coefficient itself, so the conclusions' caveat about LO truncation doesn't cover it. I'd want to see a spline-sensitivity study or a direct comparison before trusting the points deep in the BEC phase. The mu_Q-axis pressure also inherits this, plus the unknown size of O(mu^4) terms. No continuum extrapolation is shown—fine for a proceedings, but worth stating explicitly.\n\nIn short: the central new result is plausible and the method is a genuine step forward, but the BEC-phase numbers need more validation than the paper provides. It is a good conference paper that should grow into a stronger full publication.\n\nTake it to peer review—it is a legitimate first result, not a desk reject. I would send it to a referee who knows the isospin-density literature. I wouldn't cite the numbers yet; I'd wait for the journal version.","headline":"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.","tokens_in":9116,"tokens_out":3801,"would_cite":false,"duration_ms":35487,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["lattice QCD","equation of state","isospin chemical potential","Taylor expansion","charge chemical potential","pion condensation","lepton flavour asymmetry","BEC phase"],"falsifier":"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.","tokens_in":8159,"feed_emoji":"⚛️","tokens_out":7326,"duration_ms":65005,"temperature":0.7,"pith_summary":"This paper extends lattice QCD equation-of-state results from the pure isospin chemical potential axis into the directions of small baryon and strangeness chemical potentials. It does so by computing the leading-order Taylor expansion coefficients directly at nonzero isospin chemical potential, treating the isospin axis as the expansion point. The central technical obstacle, extrapolating the fully connected two-point contributions to zero pion source, is met with a singular-value based valence improvement combined with a density-improved method that reduces fluctuations. Using these coefficients, the paper obtains, for the first time, the QCD pressure on the pure electric charge chemical potential axis, a regime relevant for early-Universe evolution with large lepton flavour asymmetries. The calculation covers temperatures around 123–165 MeV and charge chemical potentials up to about 1.5 pion masses, including the pion-condensed phase.","feed_headline":"First lattice QCD pressure at nonzero charge chemical potential","feed_subtitle":"Taylor expansion from the isospin axis reaches pion-condensed matter, a regime for early Universe lepton asymmetries.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the isospin-axis equation of state and the spline-interpolated isospin density used to build the charge-axis pressure.","marker":"[17]"},{"why":"Provides the lattice action, simulation setup, and the valence quark improvement program that this paper adapts to Taylor coefficients.","marker":"[7]"},{"why":"Introduces the Taylor expansion method on the lattice that this paper re-centres on nonzero isospin expansion points.","marker":"[18]"},{"why":"Establishes the pion-condensed BEC phase and the zero-temperature Silver-Blaze structure used as boundary conditions.","marker":"[4]"},{"why":"Introduces the pionic source regulator used for simulations in the BEC phase.","marker":"[5]"},{"why":"Provides the finite-temperature finite-isospin lattice framework with the source regulator.","marker":"[6]"},{"why":"First account of Taylor expansions at nonzero isospin chemical potential, extended here to leading order and to the charge axis.","marker":"[16]"}],"fun_headline_variants":["Isospin-axis Taylor expansion yields QCD pressure at charge","Lattice QCD pressure at charge chemical potential via isospin","Pion-condensed phase reached for charge chemical potential in QCD","First QCD equation of state along charge chemical potential axis","QCD EoS at charge from isospin expansion reaches pion-condensed regime"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Isospin-axis Taylor expansion yields QCD pressure at charge","Lattice QCD pressure at charge chemical potential via isospin","Pion-condensed phase reached for charge chemical potential in QCD","First QCD equation of state along charge chemical potential axis","QCD EoS at charge from isospin expansion reaches pion-condensed regime"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000636,"raw_usage":{"total_tokens":2894,"prompt_tokens":871,"completion_tokens":2023,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":1940}},"tokens_in":487,"tokens_out":2023,"duration_ms":14126,"temperature":1.0,"reasoning_tokens":1940,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:04:19.509002+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Taylor expansion method on the lattice that this paper re-centres on nonzero isospin expansion points."},{"cited_title":"Equation of state and Taylor expansions at nonzero isospin chemical potential","cited_arxiv_id":"2212.01431","evidence_quote":"First account of Taylor expansions at nonzero isospin chemical potential, extended here to leading order and to the charge axis."}],"review_version":1}