{"id":"e12547b4-280f-45d3-be21-4d47bf6b7763","arxiv_id":"2501.19251","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Lattice QCD simulations with up to seven flavors indicate that the first-order chiral transition seen on coarse lattices is a cutoff artifact, so the continuum chiral transition is second-order for all flavors up to the conformal window.","lead":"Using lattice simulations with many quark flavors, this paper studies whether the QCD phase transition becomes first-order as more massless quarks are added. The results suggest the first-order behavior seen on coarse lattices is an artifact, and the transition stays second-order in the continuum, possibly down to zero temperature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim 'second order for all Nf<=7' is not supported by the presented fits: Nf=3,4 have no extrapolation, and the nonzero intercepts for Nf=5-7 rest on an unvalidated three-point tricritical fit.","rationale":"The reader's weakest-assumption analysis already identifies Eq. (4) and the missing extrapolation for Nf=3,4 as the load-bearing point, and my stress-test converges on the same issue. The physical reasoning behind Eq. (4) is plausible: near a tricritical point, the wing-line scaling exponent 2/5 is expected from tricritical mean-field theory, and the Nf=6 data showing leading-order scaling across three lattice spacings is genuine supporting evidence. The trend aT^tric -> 0 as Nf increases is also internally consistent with the physical picture of a tricritical point at Nf>7. Nevertheless, the paper overstates what the fits establish. For Nf=3 and 4 the paper explicitly says extrapolation is impossible, so including them in the 'all Nf<=7' conclusion is not supported by the presented analysis. For Nf=5,6,7 the entire conclusion rests on a three-point fit in which the coarsest lattice spacing strongly influences the intercept; without fit-range stability or an independent check of the assumed functional form, a spurious nonzero intercept cannot be excluded. The final sentence of Section 4 concedes exactly this residual possibility, yet the abstract and Section 3.1 present the conclusion without that caveat. A concrete refit excluding the coarsest point or adding higher-order terms would settle whether the intercepts are robust. If they survive, the conditional acceptance is appropriate; if not, the claim should be restricted to Nf=5-7 or downgraded further. My recommendation is therefore to keep the reader's CONDITIONAL verdict, with the condition that the fit stability and the Nf=3,4 gap be addressed.","tokens_in":10090,"tokens_out":8588,"duration_ms":87835,"concrete_test":"Refit the Nf=5,6,7 Z2-boundary points using Eq. (4) in three variants: (a) omitting the coarsest Ntau=4 point, (b) adding a (am)^{6/5} term, and (c) restricting to Ntau=6,8,10 only. If aT^tric(Nf) becomes compatible with zero in any variant, the continuum conclusion does not follow. In addition, require the authors either to supply a fit for Nf=3 and 4 or to remove these flavours from the headline claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central conclusion in Section 3.1 states that extrapolations terminate for all Nf<=7 at nonzero aT^tric(Nf), implying that the first-order region is a cutoff effect. However, the same section says that for Nf=3 and 4 'we only have two data points making an extrapolation impossible'; those flavours are asserted rather than fitted. For Nf=5,6,7 the conclusion depends entirely on Eq. (4), aT_c(am,Nf)=aT^tric(Nf)+A(Nf)(am)^{2/5}+B(Nf)(am)^{4/5}+..., with at most three data points per flavour at Ntau=4,6,8,10. The nonzero intercept is the key evidence that the Z2-boundary does not reach the continuum. If the asymptotic 2/5 behaviour is not yet dominant at these lattice spacings, or if O(a^2) discretisation effects mimic the 2/5 curvature over the fitted mass range, the fitted intercepts are unreliable. No fit parameters, chi^2 values, or fit-range stability checks are reported, so this cannot be verified from the paper. The conclusion is also weakened by the paper's own final caveat that a first-order region could still emerge 'at very small lattice spacings', which is logically inconsistent with the earlier statement that such a transition 'can be excluded for all numbers of flavours'.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the order of the chiral phase transition in massless many-flavour QCD using unimproved staggered fermions. For each N_f in [2,8] and lattice spacings N_tau in {4,6,8,10}, the authors map the Z2 boundary separating the first-order region from the crossover in the lattice parameter space, extract critical masses via finite-size scaling of the kurtosis with 3D Ising critical exponents, and then fit the tricritical scaling form Eq. (4) to obtain a tricritical lattice spacing aT^tric(N_f). The fitted nonzero intercepts for N_f <= 7 lead the authors to conclude that the first-order region seen on coarse lattices is a cutoff effect and that the continuum chiral-limit transition is second order for all N_f <= 7. For N_f = 8 the apparent critical points at N_tau = 8 and 10 are attributed to a lattice bulk transition. Scale setting via the Sommer parameter is used to convert the results to physical temperatures, yielding tricritical temperatures T^tric(N_f) that decrease toward zero, consistent with a physical tricritical point at T=0 between N_f=7 and 8.","tokens_in":54,"tokens_out":3893,"duration_ms":97990,"significance":"If the central claim is correct, it is an important result: the first-order chiral transition observed on coarse lattices at several flavour numbers would not survive the continuum limit, and the chiral phase transition in massless QCD would be second order for all N_f below the conformal window. The paper also gives a falsifiable prediction of a physical tricritical point between 7 and 8 flavours at T=0. The approach is sound in conception: the finite-size scaling analysis uses standard 3D Ising critical values, the tricritical scaling form Eq. (1) is externally motivated by tricritical theory, the N_f=6 data show leading-order scaling over three lattice spacings in Fig. 3b, and the codebase CL2QCD is openly available. The main weakness is that the central extrapolation is not documented quantitatively: no fit parameters, chi^2 values, fit ranges, or stability checks are reported, and for N_f=3 and 4 no extrapolation is possible with only two data points.","major_comments":[{"comment":"The conclusion that the first-order region is absent in the continuum for 'all flavours N_f <= 7' is not supported for N_f=3 and 4. The same section states that for these flavours 'we only have two data points making an extrapolation impossible.' No aT^tric values are obtained by fitting for N_f=3 and 4, yet these flavours are included in the summary claim. Either additional data or an explicitly restricted conclusion (e.g., N_f=5,6,7) is needed.","section":"Section 3.1"},{"comment":"The paper does not report the fit parameters, chi^2/dof, fit ranges, or fit-range stability for the aT^tric(N_f) extrapolations. The NLO form in Eq. (4) has three free parameters (aT^tric, A, B); with only two or three data points per flavour, the fit is underdetermined or has zero degrees of freedom, so the plotted lines in Fig. 3b cannot be quantitatively assessed. The authors should provide a table of fit results and demonstrate that the nonzero intercepts are stable under variations of the fit range.","section":"Section 3.1, Eq. (4)"},{"comment":"The same fit-documentation problem applies to the T^tric(N_f) values used to support the statement that T^tric(N_f) decreases to zero. In addition, the scale-setting procedure relies on the improved Sommer parameter, but no lattice spacing values, statistical or systematic uncertainties, or checks of the N_f-independence of r_1 are presented. The claim T^tric(N_f) -> 0 is therefore not quantitatively supported by the reported data.","section":"Section 3.2, Eq. (5)"},{"comment":"The final caveat that a first-order region could still emerge 'at very small lattice spacings' is logically in tension with the earlier claim that a first-order transition in the chiral limit 'can be excluded for all numbers of flavours.' If the caveat stands, the data constrain only the simulated lattice spacings, and the conclusion should be rephrased as a statement about the absence of evidence in the investigated range, not an exclusion.","section":"Section 4"},{"comment":"For N_f=8, the conclusion that the observed first-order region is not connected to the continuum rests on identifying the N_tau=8 and 10 points as a bulk transition and on the general expectation that bulk transitions are lattice artifacts. This is plausible, but it does not demonstrate that no thermal Z2 line exists at finer lattice spacings. The statement that the first-order region is 'not connected to the continuum for any N_f' therefore exceeds what the N_f=8 data alone can show.","section":"Section 3.3"}],"minor_comments":[{"comment":"The word 'Szenario' should be 'Scenario'.","section":"Figure 2 caption"},{"comment":"The x-axis label aT = N_\\tau^{-1} is introduced only in the text; the figure caption should state that aT is in lattice units and that the continuum limit is at the origin.","section":"Section 3.1, Figure 3b"},{"comment":"The caption mentions 'dotted lines are LO-fits, while dashed lines are NLO-fits,' but the fit ranges and the number of points included in each fit are not specified; this information should be given in the caption or in the text.","section":"Section 3.1, Figure 3b caption"},{"comment":"All numerical results are presented only through figures; a table of critical masses, critical couplings, lattice spacings, and fit results would substantially improve reproducibility and allow independent verification.","section":"General"},{"comment":"The sentence stating that critical beta-values for N_tau=8 and 10 are identical but 'not shown' should either be supported by a table or figure or be removed.","section":"Section 3.3"},{"comment":"The axis label 'm^{2/5} [r1]^{2/5}' should be explained; the mass variable and its units are not defined in the figure caption.","section":"Figure 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a proceedings-style contribution with a strong central idea, but the published version currently overreaches the reported evidence. The main fixable issues are the unsupported inclusion of N_f=3 and 4 in the central claim, the lack of fit documentation for Eqs. (4) and (5), and the internally inconsistent final caveat. These can be addressed by adding a fit-results table, providing stability checks, and carefully rewording the conclusion. If the authors are unable to supply additional data for N_f=3 and 4, the claim should be restricted to the flavour numbers for which extrapolations are actually performed. With those changes the paper could be acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new content is the temperature determination at the tricritical points for Nf=5,6,7 and the suggestion that the physical tricritical point lies at T=0 between seven and eight flavours. That would tie the chiral transition to the conformal window, so if it holds it matters. The paper is honest about the Nf=8 bulk obstruction and gives a plausible argument that the first-order region on coarse lattices shrinks with lattice spacing. The Nf=6 case does show leading-order tricritical scaling across three lattice spacings, which is the strongest piece of evidence.\n\nThe soft spots are real and need to be named. The headline conclusion 'second order for all Nf<=7' is not supported by the fits shown. In Section 3.1 the paper itself says Nf=3 and 4 have only two data points, making extrapolation impossible, yet those flavours are included in the all-Nf claim. For Nf=5,6,7 the fits have at most three data points per flavour, with no error bars, no chi-square, and no fit-range stability checks. With three parameters in Eq. (4), a three-point fit is a tautology: the intercept is whatever the curve wants it to be. The nonzero aT_tric values could be real, but the paper does not demonstrate that the asymptotic 2/5 behaviour is already dominant at these lattice spacings, and O(a^2) discretisation effects could mimic the curvature. The final caveat about a possible first-order region emerging at very small lattice spacings is also in tension with the earlier statement that such a transition 'can be excluded for all numbers of flavours.' The T=0 extrapolation is likewise based on only a few T_tric values, so it is a hint, not a determination.\n\nThe paper is not circular: the critical masses are measured, then fitted with externally cited tricritical exponents, and the T=0 statement uses separately scale-set temperatures. Self-citations to the group's earlier work are reasonable here because that work set up the method.\n\nWho is this for? Lattice practitioners working on the chiral transition and conformal window. They should read it, but with the caveat that the central claim is overstated. A serious referee should not desk-reject it; the strategy is sound and the data are real. The manuscript needs revision to restrict the conclusion to what the fits actually cover, add numerical tables and fit-range stability checks, and reconcile the final caveat with the exclusion claim.","headline":"A useful step in mapping the many-flavour chiral transition, but the all-Nf<=7 conclusion reaches beyond the fits shown.","tokens_in":10963,"tokens_out":2465,"would_cite":false,"duration_ms":23453,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T25","81V05"],"pacs":["11.15.Ha","12.38.Gc","11.30.Rd"],"model":"deepseek-v4-flash","headline":"The paper claims that the first-order chiral transition observed on coarse lattices is a cutoff effect for up to seven quark flavours, leaving a second-order transition in the continuum chiral limit.","keywords":["chiral phase transition","many-flavour QCD","tricritical point","chiral limit","lattice QCD","staggered fermions","Z2 boundary","conformal window"],"falsifier":"Run the same Z2-boundary determination on finer lattices (N_tau = 12 or 16) for Nf = 6, where three lattice spacings already support leading-order scaling. If the critical masses on those finer lattices fall on the fitted curve and keep the intercept at the same nonzero aT^tric, the paper's conclusion holds; if the intercept moves toward zero as N_tau grows, the first-order region survives in the continuum.","tokens_in":9887,"feed_emoji":"⚛️","tokens_out":6821,"duration_ms":53529,"temperature":0.7,"pith_summary":"This paper asks whether the chiral phase transition in massless QCD turns first order as the number of quark flavours grows, and answers no for every flavour number up to seven. On coarse lattices a first-order region does appear, but the authors show it is a cutoff effect: for each Nf <= 7, the boundary line between the first-order and crossover regions terminates at a nonzero tricritical lattice spacing, so the continuum extrapolation reaches a second-order transition. The critical temperature at the tricritical point falls as Nf increases, and the data are consistent with a physical tricritical point at vanishing temperature, which would coincide with the onset of the conformal window. If correct, the result rules out a first-order chiral transition in the continuum and ties the order of the transition to conformal dynamics.","feed_headline":"First-order chiral region vanishes in continuum QCD","feed_subtitle":"Data with up to seven flavours put the tricritical point at zero temperature, pointing to second order everywhere.","key_machinery":"The load-bearing object is the Z2 boundary: the line of second-order transitions that separates the first-order region from the crossover in the extended parameter space of gauge coupling, quark mass, flavour number, and lattice spacing. The boundary is located by finite-size scaling of the Binder cumulant of the chiral condensate, using the 3D Ising critical value B4 = 1.6044(10) and the Ising exponents y_t and y_h. Its extrapolation to the chiral limit is governed by tricritical scaling, aT_c(am,Nf) = aT^tric(Nf) + A(Nf)(am)^{2/5} + B(Nf)(am)^{4/5} + ..., whose intercept aT^tric is the tricritical lattice spacing below which the transition becomes second order. That intercept is what decides whether the first-order region reaches the continuum.","core_discovery":"On its own terms, the paper establishes that with standard staggered fermions and the Wilson gauge action, the first-order chiral transition seen for many flavours on coarse lattices does not survive the continuum limit. Fitting the Z2-boundary critical masses with the tricritical scaling form aT_c(am,Nf) = aT^tric(Nf) + A(Nf)(am)^{2/5} + B(Nf)(am)^{4/5} + ... gives nonzero intercepts aT^tric(Nf) for each Nf from 3 to 7, meaning the first-order region pinches off at a finite lattice spacing; the continuum chiral limit is therefore second order for all Nf <= 7. For Nf = 8 the same analysis is blocked by a bulk transition, but the observed first-order region is again disconnected from the continuum. Temperatures at the tricritical points decrease with rising Nf, and the authors find the data compatible with T = 0 at the physical tricritical point, located between seven and eight flavours.","pith_inferences":["If the T = 0 tricritical point is real, the chiral transition order becomes a clean observable for locating the conformal-window onset in many-flavour gauge theories.","The two-point fits for Nf = 3 and 4 leave the extrapolation unconstrained; finer lattices or additional masses would provide a direct test of the assumed scaling.","The same technique of extending the Columbia plot into the lattice-spacing direction could be used to separate other thermal transitions from bulk artifacts.","A confrontation with improved lattice actions at Nf = 8 would help distinguish the bulk transition from any physical signal, sharpening the boundary of the claim."],"forward_implications":["For every Nf <= 7, the first-order region is a lattice artifact, and the continuum chiral limit transition is second order.","The physical tricritical point, where the tricritical line meets the continuum, has T = 0 and lies between Nf = 7 and 8, possibly coinciding with the onset of the conformal window.","No first-order chiral transition occurs for any number of flavours, unless an unknown first-order region appears at very small lattice spacings.","The tricritical temperature decreases with flavour number, confirming the expected trend and connecting the transition to conformal scaling.","For Nf = 8, the thermal Z2 line ends at the bulk transition, so no tricritical scaling applies, but the first-order region is still not connected to the continuum."],"supporting_citations":[{"why":"Earlier mapping of the Z2 line showed the first-order region shrinks with lattice spacing, providing the starting point this work extends.","marker":"[6]"},{"why":"Established the method of analysing the chiral transition with non-integer numbers of flavours, used here for extrapolations.","marker":"[18]"},{"why":"Supplies the tricritical scaling form aT_c = aT^tric + A(am)^{2/5} + B(am)^{4/5} + ... used in all fits.","marker":"[32]"},{"why":"Provides the 3D Ising critical exponent y_t = 1/v entering the finite-size scaling fit.","marker":"[34]"},{"why":"Provides the 3D Ising exponent y_h and the critical Binder ratio used to extract the Z2 critical mass.","marker":"[35]"},{"why":"Introduces the improved Sommer scale r1 used to set physical units and convert lattice spacings to fm.","marker":"[36]"},{"why":"Summarises scale-setting practice in lattice QCD, supporting the temperature determination.","marker":"[37]"},{"why":"Documents the bulk transition in many-flavour QCD, the lattice artifact that blocks tricritical scaling for Nf = 8.","marker":"[40]"}],"fun_headline_variants":["First-order chiral region is a lattice artifact","Tricritical point at zero temperature, so second order","Many-flavour QCD: continuum transition is second order","Chiral limit is second order for up to seven flavours"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes that the tricritical scaling formula aT_c = aT^tric + A(am)^{2/5} + B(am)^{4/5} accurately describes the Z2 boundary over the simulated mass range; if it does not, the nonzero intercepts do not prove the first-order region vanishes in the continuum.","fun_headline_variants_meta":{"raw":{"variants":["First-order chiral region is a lattice artifact","Tricritical point at zero temperature, so second order","Many-flavour QCD: continuum transition is second order","Chiral limit is second order for up to seven flavours"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001586,"raw_usage":{"total_tokens":6370,"prompt_tokens":1034,"completion_tokens":5336,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":650,"completion_tokens_details":{"reasoning_tokens":5272}},"tokens_in":650,"tokens_out":5336,"duration_ms":35467,"temperature":1.0,"reasoning_tokens":5272,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T20:48:29.014950+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Z2-boundary determination on finer lattices (N_tau = 12 or 16) for Nf = 6, where three lattice spacings already support leading-order scaling. If the critical masses on those finer lattices fall on the fitted curve and keep the intercept at the same nonzero aT^tric, the paper's conclusion holds; if the intercept moves toward zero as N_tau grows, the first-order region survives in the continuum.","supporting_citations":[{"cited_title":"Lawrie and S","cited_arxiv_id":null,"evidence_quote":"Supplies the tricritical scaling form aT_c = aT^tric + A(am)^{2/5} + B(am)^{4/5} + ... used in all fits."},{"cited_title":"Ising Universality in Three Dimensions: A Monte Carlo Study","cited_arxiv_id":"cond-mat/9509016","evidence_quote":"Provides the 3D Ising critical exponent y_t = 1/v entering the finite-size scaling fit."},{"cited_title":"Scale setting in lattice QCD","cited_arxiv_id":"1401.3270","evidence_quote":"Summarises scale-setting practice in lattice QCD, supporting the temperature determination."},{"cited_title":"The bulk transition of many-flavour QCD and the search for a UVFP at strong coupling","cited_arxiv_id":"1012.5971","evidence_quote":"Documents the bulk transition in many-flavour QCD, the lattice artifact that blocks tricritical scaling for Nf = 8."}],"review_version":1}