{"id":"1adb1306-bf0f-4bde-85f0-dbf6ff4c26a0","arxiv_id":"2501.19336","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"Scaling fits give a chiral transition temperature around 144 MeV at one lattice spacing and a not yet continuum extrapolated QCD critical endpoint near T=105 MeV, mu_B=422 MeV.","lead":"By analyzing lattice QCD simulations near the chiral limit and at imaginary chemical potential, this report extracts the transition temperature, its curvature, and a candidate location for the QCD critical point. The reason to read it is that a reliable map of the QCD phase diagram is central to heavy-ion collisions, neutron stars, and cosmology.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Linear scaling-field ansatz (Eq. 7) is the key unchecked assumption; the wide fit range in Fig. 3 may bias the extracted (T_cep, mu_cep).","rationale":"The reader correctly identified the linear mixing ansatz (Eq. 7) as the weakest assumption in the endpoint extraction. My analysis confirms this is the most load-bearing concern: every other step in the Lee-Yang edge analysis (Pade construction, pole identification, universal scaling fit) either has been benchmarked in known models or is explicitly acknowledged as preliminary, but the mapping from the universal z_LY location to the (T, mu_B) trajectory is entirely determined by the assumed linear forms. The paper provides no test of the ansatz's validity in the fit region, and the plotted temperature range suggests that data far from the critical point are included, where linear dominance is not guaranteed. The concrete test I propose directly probes whether nonlinear terms or a narrower fit window alter the endpoint coordinates by more than the quoted uncertainties; if so, the result is underdetermined and the CONDITIONAL verdict is appropriate, but no stronger conclusion should be drawn until such a test is performed. I therefore agree with the reader's assessment and recommend no change to the conditional verdict.","tokens_in":13999,"tokens_out":9033,"duration_ms":89953,"concrete_test":"Using the data behind Fig. 3, refit the Lee-Yang edge trajectory with an extended scaling-field ansatz that adds quadratic terms, e.g., t = alpha_t (T - T_cep) + beta_t (mu_B - mu_cep) + gamma_t (T - T_cep)^2 and h = alpha_h (T - T_cep) + beta_h (mu_B - mu_cep) + gamma_h (mu_B - mu_cep)^2, and also repeat the fit with the temperature window restricted to |T - T_cep| < 15 MeV and < 25 MeV. If the extracted (T_cep, mu_cep) moves by more than the quoted errors, or the chi^2 improves significantly with the added quadratic terms, the linear ansatz is not controlling the systematics and the endpoint claim is not reliable as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central endpoint claim rests on the linear mixing ansatz Eq. (7), which assumes that the Z(2) scaling fields t and h are linear in (T - T_cep) and (mu_B - mu_cep). This is the leading-order Taylor expansion, but the fit in Fig. 3 uses Lee-Yang edge data over a wide temperature range (the axis spans roughly 80-170 MeV), while T_cep is determined to be 105 MeV with a lower uncertainty of -18 MeV. For |T - T_cep| up to about 50 MeV, regular and nonlinear scaling-field contributions are not necessarily negligible; the universal power law Im[mu_LY] ~ (T - T_cep)^{beta delta} and the linear relation Re[mu_LY] = mu_cep - (beta_h/alpha_h)(T - T_cep) hold only asymptotically. The paper does not state the temperature interval used in the fit nor test the stability of the result under a restricted fit window, and it relies solely on this assumed form to convert the measured Pade poles into a critical-point coordinate. If the linear ansatz is inadequate, both T_cep and mu_cep shift, and the quoted asymmetric errors do not reflect this systematic. This is the most load-bearing concern because it is the only path from the Pade poles to the claimed endpoint; the Pade pole identification itself has been benchmarked in the 2d Ising and Roberge-Weiss cases, but the linear scaling-field mapping has no such cross-check in QCD.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings article reports two quantitative results from lattice QCD with (2+1)-flavor HISQ fermions: (i) a determination of the chiral transition temperature and curvature coefficients from O(2) scaling fits to the improved order parameter and mixed susceptibilities at N_tau=8, and (ii) a determination of the QCD critical endpoint (T_cep, mu_cep) = (105(+8,-18), 422(+80,-35)) MeV at N_tau=6 from the universal scaling of Lee-Yang edge singularities obtained via multi-point Pade approximants to the baryon number density at imaginary chemical potential. The Lee-Yang edge analysis uses a linear mixing ansatz for the Z(2) scaling fields (Eq. (7)) and identifies the closest Pade pole with the Lee-Yang edge, a strategy benchmarked in earlier work on the Roberge-Weiss transition and the 2d Ising model. The paper also quotes curvature coefficients for the chiral transition line and compares its endpoint estimate with a HotQCD [4,4] Pade analysis.","tokens_in":14453,"tokens_out":7021,"duration_ms":64458,"significance":"If the endpoint result holds, it is a noteworthy step toward locating the QCD critical point from first-principles lattice calculations, and the agreement with the HotQCD [4,4] Pade data in Fig. 3 is encouraging. The chiral Tc determination is consistent with previous continuum results and uses universal O(2) scaling functions imported from Ref. [11], which gives the chiral analysis an external, non-perturbative anchor. The multi-point Pade method has been validated in controlled settings, so the methodological pipeline is credible. However, the endpoint value is presented with only a companion-paper reference for the detailed analysis, and the present manuscript does not provide the fit diagnostics needed to judge the robustness of the linear scaling-field ansatz. The paper is a proceedings summary, and its main value is in drawing attention to a specific, falsifiable prediction for the critical endpoint.","major_comments":[{"comment":"The endpoint extraction rests entirely on the linear mixing ansatz t = alpha_t (T - T_cep) + beta_t (mu_B - mu_cep) and h = alpha_h (T - T_cep) + beta_h (mu_B - mu_cep), with the four mixing parameters fitted to data. The paper states that 'the scaling directions are not known a priori,' so the linear form is an assumption rather than a derived relation. The fit shown in Fig. 3 spans a wide temperature range (the horizontal axis runs from 80 to 170 MeV), while T_cep is about 105 MeV; for |T - T_cep| up to roughly 50 MeV, regular and nonlinear contributions to the scaling fields are not guaranteed to be negligible. The manuscript does not state the temperature interval used in the fit, the number of data points, or the chi^2/dof, and no stability check under a restricted fit window is reported. Because the quoted asymmetric errors (+8/-18 MeV in T, +80/-35 MeV in mu_B) do not include this systematic, the central claim is underevidenced. Please report the fit range and diagnostics and test the stability of (T_cep, mu_cep) when the fit window is narrowed around T_cep and when quadratic terms are added to Eq. (7).","section":"Sec. 3, Eq. (7)"},{"comment":"The chiral transition temperature is reported inconsistently. The Fig. 2 caption gives (T_c, z_0, h_0^{-1/delta}) = (143.8(2) MeV, 1.45(3), 39.0(3)), while the text states 'The fit yields a critical temperature of Tc = 143.9(5) MeV for Ntau=8' and later says the joint pseudo-critical fit gives 'yet another independent determination ... Tc = 143.9(5) MeV.' The figure itself contains 'z0=1.42 Tc=143.7' in the right panel. The reader cannot tell which of these numbers corresponds to which fit (order-parameter vs pseudo-critical, with or without H=1/160). Since the abstract lists the critical temperature as a reported result, these values must be reconciled and clearly attributed to the specific fits.","section":"Sec. 2 and Fig. 2"},{"comment":"The order-parameter fit excludes the H=1/160 data point because of finite-size concerns, and the text later says differences are included as a systematic error, but the magnitude of that systematic is not given in the manuscript. Similarly, the curvature coefficients kappa_2^l, kappa_2^s, kappa_11^{ls} are said to be determined from ratios of mixed susceptibilities, but their numerical values are not tabulated here (only the combination kappa_B^{mu_S=0}=0.015(1) appears). Since the abstract lists the critical temperature and curvature as reported results, the manuscript should state the final values and the systematic error attached to the H=1/160 exclusion, or clearly point to the specific table in Ref. [8] that contains them.","section":"Sec. 2, H=1/160 exclusion"}],"minor_comments":[{"comment":"There are several typographical errors: 'Genvea' should be 'Geneva', 'weather' should be 'whether', 'easely' should be 'easily', and 'chrial' in the Fig. 1 caption should be 'chiral'.","section":"General"},{"comment":"Ref. [20] is cited only as arXiv:2405.10196 without a journal or version; since the endpoint analysis relies on it, please provide the preprint version or a DOI if available.","section":"References"},{"comment":"The labels 'HotQCD [4,4]' and 'BiPar Multi' are not defined in the caption; 'BiPar' should be spelled out as the Bielefeld-Parma collaboration and the meaning of '[4,4]' stated.","section":"Fig. 3"},{"comment":"The notation kappa_11^{ls} is not defined explicitly; the subscript 11 convention (presumably kappa times mu_l mu_s normalized by T^2) should be stated for readers not familiar with the companion papers.","section":"Eq. (2)"},{"comment":"The sentence 'We performed fits in which Tc and t1,x are kept as free fit parameters as well as using for Tc and t1,x the values determined from the scaling fits to the order parameter M given in Eqs. 36, 37' lacks the actual numbers from these two procedures; please provide the resulting Tc values or refer to a table.","section":"Sec. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings-style summary that leans heavily on two companion papers (Refs. [8] and [20]). My recommendation is based on the manuscript as submitted: the chiral Tc inconsistency and the absence of fit diagnostics for the central endpoint claim are fixable in revision, but they are load-bearing enough that the present version is not yet acceptable. I would not reject the paper; the methodological benchmarking and external scaling functions are genuine strengths."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here’s my take: this is a competent proceedings summary, not a new research paper. All the quantitative results come from Ref. [8] and Ref. [20], and the paper says so. Its value is a compact, readable exposition of the scaling framework and the Lee-Yang analysis, useful for someone wanting the Bielefeld-Parma picture in one place.\n\nThe chiral scaling section is solid. The order-parameter scaling collapse in Fig. 2 is convincing, and the T_c values are consistent with the earlier continuum estimate. The treatment of the H=1/160 data—excluded or shifted—is honest and gives a sense of the systematic uncertainty. The Lee-Yang endpoint analysis is genuinely strengthened by the benchmarks in 2d Ising and Roberge-Weiss; those checks are evidence the multipoint Padé method finds the right singularities.\n\nThe main soft spot is the linear mixing ansatz in Eq. (7). The paper acknowledges the scaling directions are not known a priori and then assumes linear fields. The stress-test concern is fair: Fig. 3 spans a wide T range, and the paper never states the fit window or checks stability under a restricted range. Since T_cep=105 MeV and the data extend well above that, the universal power laws are being used in a regime where regular terms might matter. That makes the endpoint preliminary, which the authors say, but a fit-window stability test would turn a plausible claim into a credible one. Second, there’s a minor internal inconsistency: Fig. 2’s caption gives T_c=143.8(2) while the text reports 143.9(5); consistent within errors but it should be fixed. Third, the endpoint is at a single lattice spacing, N_tau=6, with no continuum extrapolation, so the comparison to N_tau=8 is suggestive rather than conclusive.\n\nNone of this is fatal for what the paper is: a proceedings report that clearly states the methods, the assumptions, and the uncertainties. The central argument holds up at that level. I would not cite this for the numbers—I’d cite the primary papers—but I might cite it as an overview. Who is it for: a newcomer to the field, or a colleague who wants the scaling analysis summarized without wading through three long papers.\n\nRecommendation: yes, it deserves a referee. The topic is important, the exposition is competent, and a referee could help the author fix the T_c discrepancy and add a fit-range stability check. As a proceedings, that is the right level of engagement.","headline":"A competent proceedings summary of the Bielefeld-Parma chiral scaling and Lee-Yang endpoint work; no new results, but the endpoint claim rests on an untested linear mixing ansatz.","tokens_in":14901,"tokens_out":4693,"would_cite":false,"duration_ms":44544,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V05","82B26","82B27"],"pacs":["12.38.Gc","64.60.F-","11.15.Ha"],"model":"deepseek-v4-flash","headline":"This paper claims that the QCD critical endpoint lies at $(T_{\\mathrm{cep}}, \\mu_{\\mathrm{cep}}) = (105^{+8}_{-18}, 422^{+80}_{-35})$ MeV, extracted from universal Lee-Yang scaling of multipoint Pad\\'e approximants to…","keywords":["QCD phase diagram","Lee-Yang zeros","critical endpoint","chiral phase transition","universal scaling","multipoint Pade approximant","imaginary chemical potential","lattice QCD"],"falsifier":"A multipoint Pad\\'e analysis on $N_\\tau = 8$ lattices that gives an endpoint outside the quoted $N_\\tau = 6$ uncertainty would show the result is a cutoff artefact; a direct check would also test whether the closest pole obeys $\\Im\\mu_{LY} \\sim (T-T_{\\mathrm{cep}})^{\\beta\\delta}$ with the O(N) value of $\\beta\\delta$ as $T$ approaches $T_{\\mathrm{cep}}$ from above.","tokens_in":13775,"feed_emoji":"⚛️","tokens_out":8262,"duration_ms":70784,"temperature":0.7,"pith_summary":"Using universal scaling relations, this paper reports two load-bearing results for the QCD phase diagram: a chiral-limit transition temperature $T_c = 143.9(5)$ MeV from O(2) scaling fits at $N_\\tau = 8$, and a first determination of the QCD critical endpoint at $(T_{\\mathrm{cep}}, \\mu_{\\mathrm{cep}}) = (105^{+8}_{-18}, 422^{+80}_{-35})$ MeV from $N_\\tau = 6$ lattices. The endpoint is obtained by fitting the universal scaling of Lee-Yang edge singularities, whose locations are read off as the closest poles of multipoint Pad\\'e approximants to the baryon number density computed at imaginary chemical potentials. The paper's most consequential claim is the endpoint coordinate: it is not yet continuum extrapolated, but it is consistent between the multipoint Pad\\'e and a $[4,4]$ Pad\\'e analysis at $N_\\tau = 8$. If confirmed, this would anchor the location of the QCD critical point for heavy-ion collision searches and for models of dense matter.","feed_headline":"QCD critical point pinned near T=105 MeV, mu=422 MeV","feed_subtitle":"Lee-Yang zero scaling from imaginary chemical potential lattice data locates the endpoint.","key_machinery":"The load-bearing machinery is the multi-point Pad\\'e approximant $R_n^n(\\hat\\mu_B)$ to the first cumulant of the baryon density, constructed by matching the function and its derivative at imaginary simulation points in the range $\\hat\\mu_B \\in [0, i\\pi]$; its closest pole in the complex chemical-potential plane is identified with the Lee-Yang edge singularity. The temperature flow of that pole is then mapped to the universal scaling of the Lee-Yang edge, $z_{LY} = t/h^{1/(\\beta\\delta)}$, through the linear mixing ansatz for the scaling fields $t$ and $h$ (Eq. 7). In the chiral sector, the key object is the improved order parameter $M = M_l - H\\chi_l$, which removes additive UV divergences and the leading regular $H$ contribution; its scaling is described by O(2) scaling functions, and the pseudo-critical temperatures are obtained from maxima of five different susceptibilities.","core_discovery":"The central quantitative claim of the paper is that the QCD critical endpoint is located at $(T_{\\mathrm{cep}}, \\mu_{\\mathrm{cep}}) = (105^{+8}_{-18}, 422^{+80}_{-35})$ MeV, determined on $N_\\tau = 6$ lattices by fitting the universal scaling of Lee-Yang edge singularities. The edges are obtained as the closest poles of multi-point Pad\\'e approximants to the baryon number density computed at imaginary chemical potentials $\\mu_B = i\\theta_B$; the scaling fit uses the linear mixing ansatz $t = \\alpha_t(T-T_{\\mathrm{cep}}) + \\beta_t(\\mu_B-\\mu_{\\mathrm{cep}})$ and $h = \\alpha_h(T-T_{\\mathrm{cep}}) + \\beta_h(\\mu_B-\\mu_{\\mathrm{cep}})$, together with the universal location of the Lee-Yang edge in the scaling variable $z = t/h^{1/(\\beta\\delta)}$. The paper also reports a separate chiral-limit determination $T_c = 143.9(5)$ MeV on $N_\\tau = 8$ lattices from O(2) scaling fits to the improved order parameter and pseudo-critical temperatures, and curvature coefficients $\\kappa_2^{B,\\mu_S=0} = 0.015(1)$. The endpoint result is not continuum extrapolated, but is consistent between the multi-point Pad\\'e analysis and a $[4,4]$ Pad\\'e analysis of eight-order Taylor coefficients at $N_\\tau = 8$.","pith_inferences":["A natural next step is a finer-lattice ($N_\\tau = 8$) multipoint Pad\\'e analysis of the same baryon-density data; agreement with the $N_\\tau = 6$ endpoint would be the first evidence that the location survives continuum extrapolation.","The linear mixing ansatz for the scaling fields is the softest point of the endpoint analysis; replacing it with a nonlinear ansatz or extracting the mixing angles from auxiliary observables would provide a direct test of the coordinate's stability.","Because only the closest Pad\\'e pole is used, spurious poles due to numerical noise are a risk; requiring the pole trajectory to follow the predicted $\\Im\\mu_{LY} \\sim (T - T_{\\mathrm{cep}})^{\\beta\\delta}$ power law over several temperatures would sharpen the identification.","If the endpoint is real, freeze-out curves from heavy-ion collisions at moderate energies would cross the critical region at $\\mu_B$ near 422 MeV, giving a concrete target for fluctuation measurements in low-energy collision programs."],"forward_implications":["Multi-point Pad\\'e approximants to imaginary-chemical-potential data can locate Lee-Yang edges in QCD, as validated against the Roberge-Weiss transition and the 2D Ising model.","The Lee-Yang edge scaling gives a direct handle on the endpoint; systematic errors from varying the Pad\\'e order are included in the quoted coordinate uncertainties.","The chiral scaling analysis, using the improved order parameter and pseudo-critical temperatures from five susceptibilities, yields a consistent $T_c = 143.9(5)$ MeV that agrees with the earlier continuum estimate.","The curvature coefficients $\\kappa_2^{B} = 0.015(1)$ quantify how mildly the transition temperature bends with baryon chemical potential at zero strangeness chemical potential and at zero strangeness density."],"supporting_citations":[{"why":"Supplies the lattice data and O(2) scaling fits for the chiral transition temperature and curvature.","marker":"[8]"},{"why":"Provides the earlier continuum estimate of $T_c$ that the $N_\\tau = 8$ result is consistent with.","marker":"[7]"},{"why":"Supplies the O(2) (and Z(2), O(4)) scaling functions used in the chiral fits.","marker":"[11]"},{"why":"Introduces the revised scaling equation of state underlying the linear mixing ansatz.","marker":"[12]"},{"why":"Applies the linear mixing ansatz to the QCD critical endpoint, the form used in Eq. (7).","marker":"[13]"},{"why":"Provides the universal location of the Lee-Yang edge in the scaling variable $z$ used in the endpoint fit.","marker":"[16]"},{"why":"Derives the temperature behavior of the Lee-Yang edge singularities near the QCD critical point.","marker":"[18]"},{"why":"Provides the $N_\\tau = 8$ $[4,4]$ Pad\\'e results to which the multipoint Pad\\'e endpoint is compared.","marker":"[19]"},{"why":"Companion paper with details of the multipoint Pad\\'e analysis and AIC estimates.","marker":"[20]"},{"why":"Validated the pole identification method on the Roberge-Weiss transition.","marker":"[14]"}],"fun_headline_variants":["QCD critical point at T≈105 MeV, μ≈422 MeV","Lee-Yang zero scaling places QCD endpoint at T=105, μ=422","Universal scaling fixes QCD critical point near T=105, μ=422","Lattice QCD uses Lee-Yang zeros to find critical endpoint"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The endpoint extraction assumes that, in the sampled region around the critical point, the scaling fields $t$ and $h$ are strictly linear functions of $T - T_{\\mathrm{cep}}$ and $\\mu_B - \\mu_{\\mathrm{cep}}$ with four fitted coefficients, and that the closest Pad\\'e pole is the Lee-Yang edge; nonlinear scaling fields or a misidentified pole would move the fitted endpoint coordinates.","fun_headline_variants_meta":{"raw":{"variants":["QCD critical point at T≈105 MeV, μ≈422 MeV","Lee-Yang zero scaling places QCD endpoint at T=105, μ=422","Universal scaling fixes QCD critical point near T=105, μ=422","Lattice QCD uses Lee-Yang zeros to find critical endpoint"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000462,"raw_usage":{"total_tokens":2307,"prompt_tokens":940,"completion_tokens":1367,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":1285}},"tokens_in":556,"tokens_out":1367,"duration_ms":11208,"temperature":1.0,"reasoning_tokens":1285,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T20:30:39.019452+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A multipoint Pad\\'e analysis on $N_\\tau = 8$ lattices that gives an endpoint outside the quoted $N_\\tau = 6$ uncertainty would show the result is a cutoff artefact; a direct check would also test whether the closest pole obeys $\\Im\\mu_{LY} \\sim (T-T_{\\mathrm{cep}})^{\\beta\\delta}$ with the O(N) value of $\\beta\\delta$ as $T$ approaches $T_{\\mathrm{cep}}$ from above.","supporting_citations":[],"review_version":1}