{"id":"5351abdb-0273-4460-9626-3109c864facd","arxiv_id":"2412.00598","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The heavy-quark QCD critical point is measured at N_t = 8 as kappa_c = 0.09024(46), and a first N_t = 4, 6, 8 combination gives m_PS/T_c approximately 17 to 18 with small lattice-spacing dependence.","lead":"Lattice simulations with very large spatial volumes locate the critical endpoint of the heavy-quark deconfinement transition on three lattice spacings. The physical ratio m_PS/T_c at the critical point comes out near 17 to 18 and changes little between lattice spacings.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The N_t=6 central value is fixed by assuming the Z(2) crossing height after the free fit gave b4=1.6297(84); the energy-like contamination model, not the data, sets the critical point, and the quoted systematic does not cover a non-Z2 shift.","rationale":"The paper is careful, and the HPE validation in Figs. 1 and 3 gives real supporting evidence that the effective treatment of high-order terms is sound. However, the central quantitative claim rests on a model assumption that the data themselves disfavor at more than 2 sigma in the unconstrained N_t=6 fit. The N_t=4 result is consistent with Z(2), so the tension appears precisely on the finer lattices that carry the continuum argument. The reader's conditional verdict is appropriate: the quoted errors are conditional on the contamination model being complete, and the proposed large-LT test would either validate that model or expose a missing systematic. No change to the reader's verdict is needed.","tokens_in":8026,"tokens_out":4425,"duration_ms":47814,"concrete_test":"Simulate N_t=6 at LT=21 and 24 (N_s=126 and 144) with the same eff[NLO] method and fit B4 with b4 free, without imposing Z(2). If the crossing height moves toward 1.604 and kappa_c from the free fit agrees with 0.08769 within the quoted +11/-0 systematic, the contamination model is adequate. If b4 remains above about 1.62 or kappa_c shifts by more than the quoted systematic, the fixed-Z(2) analysis is biased and the continuum extrapolation should be reassessed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 4.2 the unconstrained FSS fit at N_t=6 gives b4=1.6297(84), more than 2 sigma above the Z(2) value 1.604. The paper then fixes b4 and nu to Z(2) and adds an energy-like contamination term; the resulting kappa_c=0.08769(7)(+11/-0) is taken as central. The (+11/-0) systematic is only the difference with the fit that disregards this contamination; it does not include uncertainty in the contamination model itself. The same fixed-Z(2) procedure is used for the preliminary N_t=8 point, whose free fit has b4=1.638(28). Because the Binder crossing height is the observable that identifies the transition order, assuming its universal value is the load-bearing step: if the residual Z(3) asymmetry mentioned in Sec. 4.2 produces a non-energy-like correction, kappa_c shifts outside the quoted errors. The continuum claim m_PS/T_c=17-18 then rests on this model-dependent N_t=6/8 pair, with unpublished zero-T mass input. This is not an internal inconsistency but an under-constrained model choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a lattice study of the finite-temperature critical point in the heavy-quark region of two-flavor QCD, extending earlier N_t=4 work to N_t=6 and N_t=8 lattices with spatial aspect ratios up to L_T=18 and L_T=15. To reach these volumes the authors use the hopping parameter expansion combined with an effective shift of the leading-order couplings (eff[LO]/eff[NLO]), with checks of the convergence of the expansion and of the stability of the effective couplings. The Binder cumulant of Re Omega is studied along the transition line, and finite-size scaling fits give kappa_c=0.08769(7)(+11/-0) at N_t=6 and kappa_c=0.09024(46) at N_t=8, corresponding to m_PS^(CP)/T_c=18.07(2)(+11/-0) and 17.2(2), respectively. The paper concludes that the lattice-spacing dependence of this physical ratio is quite small up to N_t=8 and that these values are not far from the continuum limit.","tokens_in":8253,"tokens_out":7467,"duration_ms":73401,"significance":"If the determination is correct, it provides a high-precision anchor for the heavy-quark QCD critical point and demonstrates that the hopping parameter expansion with effective higher-order terms can reach much larger spatial volumes and finer lattices than previously possible. The paper is transparent about the preliminary nature of the N_t=8 result, and the N_t=6 analysis is carefully cross-checked with respect to the definition of the transition line and the convergence of the higher-order terms. The main reservation, which I detail below, is that the central N_t=6 value is obtained by fixing the Binder crossing height b4 to its Z(2) value after the unconstrained fit gives b4=1.6297(84), more than 2 sigma above 1.604. This is not a circularity problem, because b4 is used as an external benchmark rather than a fitted output, but it makes the central value dependent on a contamination model whose systematic uncertainty is not yet quantified.","major_comments":[{"comment":"The unconstrained finite-size scaling fit at N_t=6 gives b4=1.6297(84), more than 2 sigma above the Z(2) value 1.604, and the asymmetry associated with the Z(3) remnant is visible even at L_T=15. The central value is then obtained by fixing nu and b4 to their Z(2) values and adding an energy-like contamination term with a fitted amplitude, because the full six-parameter fit is unstable. The quoted systematic error (+11/-0) is only the difference between fits with and without this contamination term; it does not include uncertainty in the form of the contamination or a possible non-Z(2) correction. Since b4 is the observable that identifies the universality class, this is a load-bearing assumption. I ask the authors to demonstrate quantitatively that the contamination model accounts for the 2-sigma discrepancy, for example by reporting the fitted contamination amplitude and its predicted contribution to b4, by varying the form of the contamination term, or by fitting an independent observable with a different contamination pattern.","section":"Sec. 4.2, Eq. (8), Fig. 4 (right)"},{"comment":"The abstract states that the paper attempts a preliminary continuum extrapolation, but no continuum extrapolation is actually shown. Equations (9)-(11) list three values, m_PS/T_c=16.30(3), 18.07(2), and 17.2(2), which are non-monotonic in N_t, and the text only says that the lattice-spacing dependence is small. The difference between the N_t=6 and N_t=8 values is approximately 0.87, which is about four times the quoted error of the N_t=8 point. To support the claim that these values are not far from the continuum limit, an explicit extrapolation fit (for example linear in 1/N_t^2 or including O(a) corrections) with an assessment of systematic uncertainty is needed; otherwise the wording should be limited to a statement about the available lattice spacings.","section":"Sec. 5, Eqs. (9)-(11), and abstract"},{"comment":"The N_t=8 determination, although labeled preliminary, is used in Sec. 5 as evidence for the small lattice-spacing dependence. The quoted kappa_c=0.09024(46) is obtained after fixing nu (and effectively b4) to the Z(2) values, but the text does not state whether the same energy-like contamination correction is applied at N_t=8 or what systematic error is assigned to the contamination model. Because the N_t=8 point carries a large weight in the comparison underlying the continuum statement, this systematic should be quantified before the comparison is used as evidence.","section":"Sec. 4.3"}],"minor_comments":[{"comment":"The caption refers to 'Eq. (38)' and 'Eq. (40)', but the manuscript contains no equations with these numbers; these references are presumably to equations in ref. [9] and should be replaced with local equations or a clear citation.","section":"Fig. 4 caption"},{"comment":"The phrase 'to be compared with a previous result 1.1135(8)' is unclear: the number 1.1135(8) cannot be kappa_c in the same notation as the rest of the paper, and no conversion or definition is given. Please correct or clarify.","section":"Sec. 4.3"},{"comment":"There are small typographical errors, including 'where M is he Wilson quark kernel' and 'obsereved' in the Fig. 3 caption; these should be corrected.","section":"Sec. 2.1"},{"comment":"The physical conversion uses zero-temperature pseudo-scalar masses from ref. [11], which is listed as 'unpublished'. Please provide the interpolated mass values used in Eqs. (9)-(11) or cite a published source, so that the conversion is reproducible.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proceedings contribution and much of the N_t=6 analysis is already published in ref. [9], while the N_t=8 results depend on ref. [10], listed as 'in preparation'. The editor may wish to ask whether the continuum-related claims should be deferred until the underlying papers are available, and whether the N_t=8 systematic should be completed before the proceedings is finalized."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a straightforward, honest proceedings from WHOT-QCD, and the thing to know is that the new number---the N_t=8 critical point, kappa_c = 0.09024(46)---is real new data but officially preliminary, and the continuum extrapolation built on it is a preliminary three-point exercise, not a settled result. The N_t=6 analysis and the method were published earlier in PRD; the proceedings adds little methodologically but does give the first N_t=8 point and the first N_t=4,6,8 continuum estimate of m_PS/T_c.\n\nWhat the paper does well: the convergence checks for the HPE are concrete and not decorative. They show the relative truncation error as a function of the order, they show the strong linear correlation between high-order terms, and they test the effective resummation against the exact NLO treatment on the N_t=6 phase diagram. The N_t=4 and N_t=6 FSS analyses are already in refereed journals, and the N_t=8 section is explicitly labeled preliminary, with larger errors and fewer volumes. That is decent scientific hygiene for a proceedings.\n\nWhere it is softer, and here the stress-test note mostly lands. The central N_t=6 value comes from a fit in which b4 and nu are fixed to Z(2) after the free fit gave b4 = 1.6297(84), more than two sigma above 1.604. An energy-like contamination term is then introduced to absorb the discrepancy. The quoted systematic (+11 over -0) is simply the difference between fits with and without that term; it does not cover the uncertainty in the contamination model itself. If the Z(3) remnant seen in Fig. 5 induces a correction that is not purely energy-like, kappa_c shifts outside the quoted errors. That is a real model-dependence, and the same fixed-Z(2) procedure is used for the N_t=8 point. The continuum statement at 17-18 therefore rests on an assumption, not on data alone. That said, the paper is fully transparent: it shows the free-fit b4, it shows the Z(3) asymmetry in the complex-plane histogram, and it reports the model-dependence as a systematic. This is an under-constrained choice, not a hidden one. My soundness score matches yours at 6.\n\nWho should read this: anyone who needs the heavy-quark critical point as an anchor for the Columbia plot or as a cross-check for finite-density extrapolations. The N_t=8 number and the continuum estimate will be quoted even in this preliminary form, so the paper deserves to be read carefully. As a full paper, I would send it to referees; as a proceedings, this is the right level of detail and candor.","headline":"Careful WHOT-QCD proceedings: the new N_t=8 critical point is genuinely new but preliminary, and the continuum estimate rests on a Z(2)-fixed Binder analysis whose contamination-model uncertainty is real though transparently disclosed.","tokens_in":8907,"tokens_out":3319,"would_cite":true,"duration_ms":32042,"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":"Using Binder-cumulant scaling on large lattices, this paper locates the finite-temperature critical point of heavy-quark QCD and finds the physical ratio $m_{PS}^{(CP)}/T_c \\simeq 17$–$18$ with only mild lattice-spacing dependence.","keywords":["lattice QCD","heavy-quark region","critical point","Binder cumulant","finite-size scaling","hopping parameter expansion","deconfinement transition","Z(2) universality"],"falsifier":"Measure the same crossing ratio on $N_t=6$ lattices with aspect ratio $L_T=20$ or larger and fit it without fixing the universality class; if the fitted crossing height remains more than two standard deviations above the Ising value, the contamination model used here is not removing the extra contribution and the quoted $\\kappa_c$ would need to shift.","tokens_in":7737,"feed_emoji":"⚛️","tokens_out":10767,"duration_ms":183388,"temperature":0.7,"pith_summary":"On the Columbia plot, the map of the QCD transition as a function of quark masses, the heavy-quark corner is expected to be first-order and to end at a critical point. This paper determines that heavy-quark critical point by finite-size scaling of the Binder cumulant of the Polyakov loop, using the hopping parameter expansion to reach spatial volumes as large as $N_s=120$. On $N_t=6$ and $8$ lattices the extracted critical couplings are $\\kappa_c=0.08769(7)$ and $0.09024(46)$, giving $m_{PS}^{(CP)}/T_c=18.07(2)$ and $17.2(2)$. The paper argues that this physical ratio changes little between $N_t=6$ and $N_t=8$, so the values are probably close to the continuum limit. If correct, this locates the end of the first-order transition line in physical units to about a percent.","feed_headline":"Heavy-quark QCD critical point pinned down near continuum limit","feed_subtitle":"Large-lattice simulations fix the end of QCD's first-order line to a few percent in physical units.","key_machinery":"The load-bearing objects are the Binder cumulant $B_4$ of $\\mathrm{Re}\\,\\hat\\Omega$, treated as the magnetization of a $Z(2)$ spin system, and the hopping parameter expansion of the effective quark action. $B_4$ is a normalized fourth moment whose crossing height is universal at a $Z(2)$ critical point, so the critical coupling is found where $B_4$ becomes independent of the spatial aspect ratio $L_T=N_s/N_t$. The hopping expansion makes large spatial extents affordable, and the central mechanism for incorporating terms beyond next-to-leading order is an effective method that replaces truncated high-order terms by shifted low-order couplings, exploiting the strong linear correlation between the orders; convergence is checked by watching effective couplings stabilize as the truncation order $n_L$ is raised and by comparing with exact next-to-leading-order reweighting. Finite-size-scaling fits that include an energy-like contamination term turn the crossing into the quoted critical couplings.","core_discovery":"The central claim is that in two-flavor heavy-quark QCD the endpoint of the first-order deconfinement line lies at $\\kappa_c=0.08769(7)$ on $N_t=6$ lattices and at $\\kappa_c=0.09024(46)$ on $N_t=8$ lattices, after finite-size-scaling analysis of the Binder cumulant on spatial volumes up to aspect ratios $L_T=18$ and $15$, respectively. Using the zero-temperature pseudo-scalar meson mass, these couplings correspond to $m_{PS}^{(CP)}/T_c=18.07(2)$ and $17.2(2)$. The paper finds this ratio to vary by only about five percent between $N_t=6$ and $N_t=8$, and states that this small lattice-spacing dependence suggests the values are not far from the continuum limit. The analysis treats the Polyakov loop as the magnetization of a $Z(2)$ spin system; an apparent excess of the Binder crossing height at $N_t=6$ is accounted for by an added energy-like contamination term, after which the data are consistent with $Z(2)$ scaling.","pith_inferences":["The near-flatness of $m_{PS}^{(CP)}/T_c$ across $N_t$ suggests a testable prediction: on $N_t=10$ and $12$ lattices the ratio should stay close to $17$–$18$, and a drift beyond a few percent would indicate that the apparent continuum near-independence was accidental.","Because the hopping expansion organizes the fermion determinant as a power series in $\\kappa$, the same large-lattice machinery could be applied at finite baryon density, where the heavy-quark region has a milder sign problem; this could map the critical endpoint in the temperature-density plane.","The $Z(3)$ remnant asymmetry visible in the distribution of $\\hat\\Omega$ implies that finite-volume corrections to $Z(2)$ scaling are dominated by center-symmetry remnants on currently affordable lattices, so future $N_t>8$ studies should check whether even larger $L_T$ are needed before the scaling window opens."],"forward_implications":["If the quoted $\\kappa_c$ values are correct, the first-order deconfinement region in the heavy-quark corner of the Columbia plot ends at a pseudo-scalar meson mass of about $17$–$18$ times $T_c$ in the continuum.","The small $N_t$ dependence of $m_{PS}^{(CP)}/T_c$ between $N_t=6$ and $N_t=8$ means a continuum extrapolation is not yet needed to locate the critical point to about a percent.","The validated effective high-order hopping expansion justifies applying the same method to $N_t\\simeq 10$ lattices, where the hopping parameter is even larger and more orders are needed.","The two-flavor result can be translated to $2+1$ flavor QCD through the analytic flavor dependence of the hopping expansion, which the paper states is straightforward."],"supporting_citations":[{"why":"Supplies the Binder-cumulant finite-size-scaling method used to locate the critical point from the crossing of $B_4$.","marker":"[4]"},{"why":"Establishes the $N_t=4$ critical point with the LO/NLO reweighting method that this paper extends.","marker":"[7]"},{"why":"Provides the convergence analysis and the linear-correlation observation underlying the effective high-order hopping-parameter-expansion method.","marker":"[8]"},{"why":"Gives the $N_t=6$ high-precision analysis that yields the central value $\\kappa_c=0.08769(7)$.","marker":"[9]"},{"why":"Supplies the earlier two-flavor Wilson-fermion determination of the deconfinement critical point that the $N_t=6$ result is compared with and improves upon.","marker":"[5]"},{"why":"Supplies the preliminary $N_t=8$ Binder-cumulant data and the quoted $\\kappa_c=0.09024(46)$.","marker":"[10]"},{"why":"Supplies the zero-temperature pseudo-scalar meson mass used to convert $\\kappa_c$ into the physical ratio $m_{PS}^{(CP)}/T_c$.","marker":"[11]"}],"fun_headline_variants":["Heavy-quark QCD critical point fixed via Binder cumulant scaling","Large-lattice QCD simulation homes in on critical point","Heavy-quark QCD: critical point pinned on Nt=6,8 lattices","Near-continuum critical point in heavy-quark QCD from huge lattices","QCD critical point in heavy-quark region from large-volume scaling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central numbers depend on assuming the transition is of the three-dimensional Ising universality class and that a fitted correction term fully explains why the measured crossing ratio is higher than the Ising value; if that correction is incomplete, the critical coupling would move by more than the quoted uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["Heavy-quark QCD critical point fixed via Binder cumulant scaling","Large-lattice QCD simulation homes in on critical point","Heavy-quark QCD: critical point pinned on Nt=6,8 lattices","Near-continuum critical point in heavy-quark QCD from huge lattices","QCD critical point in heavy-quark region from large-volume scaling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000904,"raw_usage":{"total_tokens":3891,"prompt_tokens":948,"completion_tokens":2943,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":2843}},"tokens_in":564,"tokens_out":2943,"duration_ms":20378,"temperature":1.0,"reasoning_tokens":2843,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:12:25.920581+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same crossing ratio on $N_t=6$ lattices with aspect ratio $L_T=20$ or larger and fit it without fixing the universality class; if the fitted crossing height remains more than two standard deviations above the Ising value, the contamination model used here is not removing the extra contribution and the quoted $\\kappa_c$ would need to shift.","supporting_citations":[{"cited_title":"Binder,Finite size scaling analysis of Ising model block distribution functions, Z","cited_arxiv_id":null,"evidence_quote":"Supplies the Binder-cumulant finite-size-scaling method used to locate the critical point from the crossing of $B_4$."},{"cited_title":"Sugawara et al., for the WHOT-QCD Collaboration, in preparation","cited_arxiv_id":null,"evidence_quote":"Supplies the preliminary $N_t=8$ Binder-cumulant data and the quoted $\\kappa_c=0.09024(46)$."},{"cited_title":"Itagaki et al., for the WHOT-QCD Collaboration, unpublished","cited_arxiv_id":null,"evidence_quote":"Supplies the zero-temperature pseudo-scalar meson mass used to convert $\\kappa_c$ into the physical ratio $m_{PS}^{(CP)}/T_c$."}],"review_version":1}