{"id":"4c5f0a56-0a6f-467d-9f1b-00c760b048d6","arxiv_id":"2505.08658","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In SU(3) gauge theory with Nf = 4, 8, and 12 flavors, the normalized eta-prime mass sqrt(8t0) M_eta' takes values near 2.5, 5.0, and 7.5, consistent with n_f^2 scaling in the anti-Veneziano limit.","lead":"This paper uses lattice QCD simulations to measure the flavor-singlet scalar and pseudoscalar masses in SU(3) gauge theories with 4, 8, and 12 light fermions. The central finding is that the eta-prime mass, normalized by the gradient flow scale, follows the predicted 'anti-Veneziano' growth with the number of flavors.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-pole extraction of M_eta' is not validated in the conformal Nf=12 theory, so the claimed n_f^2 scaling may be an effective-mass artifact.","rationale":"The reader's weakest assumption is exactly the right one: the single-pole dominance of the topological-charge density correlator, Eq. (5.8), is the load-bearing step for the central claim. The paper provides substantial lattice evidence, high-statistics ensembles, a data release, and careful plateau checks, and I find no internal inconsistency or obvious numerical error in the presented analysis. However, those plateau checks cannot distinguish a genuine pole from a slowly varying continuum, especially in the Nf=12 conformal theory where the paper itself describes bound states as non-relativistic unparticles. Since the anti-Veneziano scaling is extracted by comparing fitted masses across Nf at a single lattice spacing, the absence of a spectral-function validation is a real limitation. The reader's CONDITIONAL verdict already captures this, so I would not change the verdict; but a concrete refit with a continuum component should be run before the n_f^2 scaling is treated as established.","tokens_in":44404,"tokens_out":3826,"duration_ms":43513,"concrete_test":"Using the released data in Ref. [73], refit the topological-charge correlator for each Nf and mf with a two-component spectral ansatz, e.g. C(r) = A_p e^{-Mr} r^{-3/2} + A_c integral from 4M_pi^2 to infinity of (s - 4M_pi^2)^alpha e^{-sqrt(s) r} r^{-3/2} ds, along with the paper's (1+3/(8r)) correction, leaving A_p, A_c, M, and alpha free. If the continuum component is preferred by the fit or if the pole mass M shifts from the Eq. (5.8) result by more than the quoted systematic error, then the single-pole extraction is not established and the n_f^2 pattern in Eq. (5.10) is not a robust test of anti-Veneziano scaling.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Eq. (5.10)/(6.2), depends entirely on extracting M_eta' by fitting the topological-charge correlator to Eq. (5.8), a one-particle free-propagator form. The paper itself states in Secs. I and II C that in the conformal Nf=12 theory all bound states are non-relativistic 'unparticles' and that no bound states exist in the chiral limit; the spectral function is then not guaranteed to be a single delta function. The plateau checks in Fig. 17 and Appendix D cover only a limited distance window, roughly r in [7,12] lattice units, and a narrow smearing window; a continuum or unparticle component can produce a slowly varying effective mass that looks flat over this range without implying a pole. The paper does not perform a spectral-function reconstruction or an alternative multi-component fit, and it quotes no uncertainty on the values in Eq. (5.10). With only one lattice spacing per Nf and no continuum extrapolation, the observed 1:4:9 pattern could be a property of the fitting window and smearing choices rather than of the theory.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a lattice study of SU(3) gauge theory with Nf = 4, 8, and 12 flavors of HISQ fermions, focusing on the flavor-singlet scalar and pseudoscalar channels. It presents an updated Nf = 8 ensemble at mf = 0.009, new Nf = 4 scalar data, chiral extrapolations of the Nf = 8 spectrum, a comparison of the scalar mass against the WT-identity form M_sigma^2 = d0 + d1 M_pi^2, and the flavor-singlet pseudoscalar mass extracted from the topological-charge correlator. The central new result is Eq. (5.10), restated as Eq. (6.2): M_eta'^2 * 8t0 is reported to be approximately (2.5)^2, (5.0)^2, and (7.5)^2 for Nf = 4, 8, and 12, respectively, independent of the fermion mass, and this is interpreted as evidence for anti-Veneziano n_f^2 scaling. The paper also interprets the Nf = 8 scalar as a walking-technicolor dilaton candidate and compares the Nf = 12 scalar with the conformal-phase expectation M_sigma < M_pi.","tokens_in":44686,"tokens_out":5570,"duration_ms":51164,"significance":"If the n_f^2 scaling of M_eta'^2 * 8t0 is robust, the paper provides a new nonperturbative scaling law connecting the flavor-singlet pseudoscalar mass to the gradient-flow scale across theories with different fermion content, which would be valuable for composite-Higgs model building and for benchmarks of near-conformal gauge theories. The paper's strengths include the use of a common lattice setup for all three theories, high-statistics ensembles, a gluonic operator that avoids pion contamination in the eta-prime channel, public data and workflow release (Ref. [73]), and detailed comparisons with LSD collaboration data. However, the flagship scaling claim currently rests on a single-lattice-spacing determination with an unvalidated single-pole spectral ansatz, and it is quoted without uncertainties; the sigma/dilaton interpretation also relies on a model-dependent relation that is partly circular in the gamma_m comparison.","major_comments":[{"comment":"The extraction of M_eta' is based entirely on fitting the topological-charge correlator to the single free-propagator form C(r) = (A/r^1.5)(1+3/8r)e^{-M r}. In the Nf = 12 theory the paper itself states in Secs. I and II C that all bound states are non-relativistic 'unparticles' and that no bound states exist in the chiral limit, so the spectral function is not guaranteed to be a single delta function. The plateau checks in Fig. 17 and Appendix D cover only the limited distance window r in [7, 12] lattice units and a narrow smearing interval; a continuum or unparticle component can produce a slowly varying effective mass that looks flat over this range without implying a pole. The authors should provide a spectral-function reconstruction or an alternative multi-component/continuum fit, and show explicitly that the fitted exponential corresponds to a pole mass rather than a window-dependent effective mass.","section":"Sec. V, Eq. (5.8) and Appendix D"},{"comment":"The values (2.5)^2, (5.0)^2, and (7.5)^2 in Eq. (5.10) and Eq. (6.2) are quoted without statistical or systematic uncertainties, so the claimed 1:4:9 ratios cannot be quantitatively assessed. In addition, the eta-prime data are taken at a single lattice spacing for each Nf, with no continuum extrapolation; the paper itself states in Sec. VI that continuum and chiral extrapolations are needed. As written, Eq. (5.10) describes a finite-lattice-spacing, fixed-smearing-window result, not a continuum scaling law. The authors should quote uncertainties and explicitly qualify the result as a single-spacing observation pending continuum and chiral extrapolations.","section":"Sec. V, Eq. (5.10), and Sec. VI"},{"comment":"The values of gamma_m inferred from the d1 fits in Table VIII are obtained by inverting Eq. (2.42), which follows from the holographic/linear-sigma-model relation F_sigma^2 = (3-gamma_m)^2 (Nf/2)(F_pi/sqrt2)^2. Comparing those gamma_m values with the hyperscaling values is therefore not an independent check of the model relation; it only tests consistency under that model assumption. For Nf = 12 the agreement between gamma_m = 0.66(31) and the hyperscaling interval 0.4-0.5 is acceptable, but the text should state that this comparison is contingent on Eq. (2.42), not a derivation of it.","section":"Sec. IV C, Eq. (2.42), and Table VIII"},{"comment":"Eq. (5.10) writes M_eta'^2 * 8t0, while Fig. 21 and definition (2.23) plot (M_eta'^2 - M_pi^2) * 8t0. At the simulated pion masses for Nf = 4 and Nf = 8, M_pi^2 * 8t0 is not negligible, so the constants (2.5)^2, (5.0)^2 and (7.5)^2 cannot simultaneously describe both quantities unless M_eta'^2(anomalous) is explicitly identified with the difference. The text and the figure must be made consistent.","section":"Sec. V, Eq. (5.10) vs. Fig. 21"}],"minor_comments":[{"comment":"The sentence 'For Nf = 8, this gives a result consistent with Eq. (2.40), as remarked on in Ref. [16].' is repeated verbatim twice in the same paragraph and should be reduced to a single occurrence.","section":"Sec. II C, before Eq. (2.41)"},{"comment":"The phrase 'loop expansion expansion parameter' contains a duplicated word and should read 'loop expansion parameter'.","section":"Sec. II C, before Eq. (2.58)"},{"comment":"The heading 'Fits of the flavor-singlet psudoscalar' contains a typo; 'psudoscalar' should be 'pseudoscalar'.","section":"Appendix D heading"},{"comment":"The caption quotes sqrt(8t0) = 7.2680(99) without stating the units; if this is in lattice units the text should say so explicitly, as is done for other quantities.","section":"Fig. 17 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a substantial numerical data release from the LatKMI collaboration, and the flavor-singlet pseudoscalar result is interesting. My main concern is that the central scaling claim is currently presented without uncertainties and rests on a spectral ansatz that is not validated in the Nf = 12 conformal phase; the sigma/gamma_m comparison also needs a clearer separation of model input from lattice output. These issues are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe one thing to know: this paper's central claim is that the flavor-singlet pseudoscalar mass, normalized by the gradient flow scale, obeys n_f^2 scaling across Nf=4, 8, and 12 QCD, with M^2_eta' * 8t0 approximately (2.5)^2, (5.0)^2, (7.5)^2. That is a clean, quantitative pattern, and if it survives scrutiny, it offers a useful organizing principle for many-flavor gauge theories.\n\nWhat is actually new: the first full lattice determination of the eta-prime mass in Nf=4, 8, and 12 using a common scale, plus new Nf=4 scalar data and an updated Nf=8 spectrum at a lighter fermion mass. The data and analysis workflows are released, which is good practice and helps reproducibility.\n\nThe paper does several things well. The comparison across theories uses the same lattice action and a gluonic scale, which is sensible. The eta-prime fits are checked for plateaus in both distance and smearing windows, and the systematic error from those choices is estimated. The authors also state plainly that continuum and chiral extrapolations are still needed.\n\nThe soft spots are real but not fatal. First, Eqs. (5.10) and (6.2) quote the normalized masses with no uncertainties attached to the numbers, even though the figures carry error bars. For a headline claim, that is an omission. Second, the eta-prime mass is extracted by fitting the topological charge correlator to a single free-propagator form. In the conformal Nf=12 theory, the paper itself describes bound states as non-relativistic unparticles, with no sharp poles in the chiral limit. With a finite fermion mass and a limited fit window (r between about 7 and 12 lattice units), a continuum contribution could masquerade as a flat effective mass. They do not attempt a spectral-function reconstruction or an alternative multi-component fit, so the 1:4:9 pattern could in principle be a fitting artifact. Third, there is only one lattice spacing per Nf, so discretization effects are uncontrolled.\n\nThe sigma/dilaton analysis is more interpretive. The relation between d1 and gamma_m comes from a holographic/linear-sigma-model assumption, and they then use it to infer gamma_m and compare with hyperscaling. That is a consistency check, not an independent test. The paper does not oversell it, but it is worth flagging.\n\nWho is this for? People working on composite Higgs, walking technicolor, or the conformal window of many-flavor gauge theories. The eta-prime scaling is the kind of result that could shape future model building if it holds up.\n\nMy recommendation: this deserves a serious referee. I would accept it for peer review and ask for (a) uncertainties on the quoted normalized masses, (b) a test or discussion of the single-pole assumption in the conformal phase, and (c) an explicit statement about the one-lattice-spacing limitation. I would cite it for the data and the scaling observation, but with a caveat until the spectral function question is settled.","headline":"The eta-prime n_f^2 scaling is a genuine and useful observation, but the quoted numbers lack uncertainties and the single-pole extraction in the conformal Nf=12 theory needs scrutiny.","tokens_in":45326,"tokens_out":3289,"would_cite":true,"duration_ms":32735,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Ha","12.38.Gc"],"model":"deepseek-v4-flash","headline":"Lattice simulations of SU(3) gauge theory with 4, 8, and 12 light fermions find that the flavor-singlet pseudoscalar mass obeys $M_{\\eta'}^2\\cdot 8t_0\\simeq (2.5)^2, (5.0)^2, (7.5)^2$, matching the anti-Veneziano prediction…","keywords":["lattice QCD","flavor-singlet mesons","eta-prime mass","anti-Veneziano limit","gradient flow","walking technicolor","conformal window","topological charge density"],"falsifier":"Fit the same correlators with a spectral function that includes a continuum or a second state and check whether the quoted $M_{\\eta'}$ values and the $1:4:9$ pattern survive; or extract $M_{\\eta'}$ on the same ensembles from the disconnected fermionic axial-current correlator and require agreement with the gluonic-operator result.","tokens_in":44178,"feed_emoji":"⚛️","tokens_out":17900,"duration_ms":140627,"temperature":0.7,"pith_summary":"The paper uses first-principles lattice simulations to compare SU(3) gauge theories with 4, 8, and 12 light fermion flavors, focusing on the flavor-singlet mesons that are invisible in conventional quark-antiquark analyses. Its central claim is that the mass of the flavor-singlet pseudoscalar, the $\\eta'$, is set almost entirely by the chiral anomaly and obeys a clean counting law: $M_{\\eta'}^2\\cdot 8t_0\\simeq (2.5)^2$, $(5.0)^2$, $(7.5)^2$ for $N_f = 4, 8, 12$, independent of the fermion mass, i.e. $M_{\\eta'}^2/\\Lambda_{\\rm IR}^2\\sim n_f^2$ with $n_f=N_f/N_c$. This 'anti-Veneziano' scaling, the regime of many flavors per color in the large-$N_c$ limit, is the opposite of the familiar suppression of the $\\eta'$ mass at small $N_f/N_c$, and it only becomes visible when masses are normalized by the gradient-flow scale $1/\\sqrt{8t_0}$ rather than by hadronic masses such as $M_\\rho$. The same data also refine the case that the eight-flavor theory walks, with a scalar meson $\\sigma$ as light as the pion that could serve as a composite Higgs, while the twelve-flavor theory sits in the conformal window with $M_\\sigma<M_\\pi$.","feed_headline":"Eta-prime mass scales as flavor count squared","feed_subtitle":"Lattice data for 4, 8, 12 flavors match a square-law prediction once the gradient-flow scale sets the ruler.","key_machinery":"The object that carries the argument is the two-point correlator of the topological charge density, $q(x)=\\frac{1}{32\\pi^2}\\epsilon_{\\mu\\nu\\rho\\sigma}\\mathrm{Tr}\\,G_{\\mu\\nu}G_{\\rho\\sigma}$, averaged over all point pairs at each separation $r=|x-y|$ and smeared by the gradient flow. Because this purely gluonic operator does not couple to pions, the correlator is dominated by the flavor-singlet pseudoscalar, and the mass is extracted by fitting to the single-particle form $C(r)=\\frac{A}{r^{3/2}}\\left(1+\\frac{3}{8r}\\right)e^{-M_{\\eta'}r}$ over a plateau in both the distance window and the smearing scale. The theoretical template is the Ward–Takahashi identity for the flavor-singlet axial current: in the anti-Veneziano limit, the regime $N_c\\to\\infty$ with $N_c\\alpha$ and $n_f=N_f/N_c$ fixed and large, the fermion-loop diagram dominates the anomaly correlator and yields $M_{\\eta'}^2/\\Lambda_{\\rm IR}^2\\sim n_f^2$, the law the data confirm. For the scalar channel, the corresponding machinery is the scale-symmetry Ward–Takahashi identity condensed into the fit form $M_\\sigma^2 = d_0 + d_1 M_\\pi^2$, whose slope $d_1$ is linked to the mass anomalous dimension through $d_1=(1+\\gamma_m)/(3-\\gamma_m)$ and to the dilaton decay constant through $F_\\sigma/(F_\\pi/\\sqrt2)=C_{\\gamma_m}/\\sqrt{d_1}$.","core_discovery":"The authors' central claim, in their own terms, is that the anomalous part of the $\\eta'$ mass in large-$N_f$ SU(3) gauge theory obeys the anti-Veneziano scaling law $M_{\\eta'}^2/\\Lambda_{\\rm IR}^2\\sim n_f^2$, and that the lattice data realize this law when the infrared scale is chosen to be the gradient-flow scale, $\\Lambda_{\\rm IR}=1/\\sqrt{8t_0}$. Concretely, they find $M_{\\eta'}^2\\cdot 8t_0\\simeq (2.5)^2$ for $N_f=4$, $(5.0)^2$ for $N_f=8$, and $(7.5)^2$ for $N_f=12$, with no significant fermion-mass dependence within each theory, so the ratios of the three values are $1:4:9$, matching $(N_f/N_c)^2$. The same normalized quantity is flat and ordered only in this scale: using $M_\\rho$ as the reference scale destroys the pattern, which the authors read as evidence that the $\\eta'$ mass is anchored to a gluonic infrared scale that tracks the number of flavors, rather than to the chiral-symmetry-breaking scale.","pith_inferences":["If the $n_f^2$ pattern survives the continuum limit, it gives model builders a sharp diagnostic: near-conformal composite-Higgs candidates with $N_f/N_c\\simeq 3$ should carry an anomalously heavy flavor-singlet pseudoscalar, about three times heavier (nine times in squared mass) than the four-flavor case in units of $1/\\sqrt{8t_0}$.","The single-pole fit implies a definite value of the topological susceptibility in each theory; measuring $\\chi_{\\rm top}$ directly from the flowed gauge configurations and checking $\\chi_{\\rm top}\\sim n_f^2$ would independently test the assumption that underlies the mass extraction.","The contrast between $M_\\rho$-normalization, which fails, and $1/\\sqrt{8t_0}$, which succeeds, suggests the $\\eta'$ mass is tied to a glueball-like scale that responds to flavor count differently from chiral-breaking scales; testing the same scaling against the string tension or glueball masses would clarify which gluonic scale the $\\eta'$ actually tracks."],"forward_implications":["The $\\eta'$ mass in multi-flavor SU(3) theories is fixed by the flavor content: $M_{\\eta'}^2\\cdot 8t_0\\simeq 2.5^2\\,(N_f/4)^2$, a parameter-free target for any theory with a given $N_f/N_c$.","The scaling is invisible when hadronic reference scales such as $M_\\rho$ are used, so comparisons of many-flavor spectra must be anchored to a gluonic scale like $1/\\sqrt{8t_0}$ rather than to hadron masses.","The eight-flavor walking scenario survives the new lightest-mass ensemble: $M_\\sigma\\lesssim M_\\pi$, $d_1\\simeq 1$, and $F_\\sigma/(F_\\pi/\\sqrt2)\\simeq 4$, keeping $\\sigma$ viable as a composite-Higgs dilaton.","The twelve-flavor theory shows $M_\\sigma<M_\\pi$ with $d_1\\simeq 0.71$, consistent with a pseudo-dilaton picture inside the conformal window, where $\\sigma$ is a pseudo-Nambu–Goldstone boson of scale symmetry but $\\pi$ is not."],"supporting_citations":[{"why":"supplies the anti-Veneziano-limit calculation predicting $M_{\\eta'}^2/\\Lambda_{\\rm IR}^2\\sim n_f^2$, the scaling law the lattice data are compared against.","marker":"[12]"},{"why":"defines the gradient-flow scale $t_0$ used both as the smearing device for the topological charge operator and as the infrared scale $1/\\sqrt{8t_0}$ that makes the scaling manifest.","marker":"[41]"},{"why":"provides the established method of extracting the $\\eta'$ mass from the topological charge density correlator, which this paper adopts for many-flavor theories.","marker":"[70]"},{"why":"derives the dilaton mass formula $M_\\sigma^2=d_0+d_1M_\\pi^2$ that organizes the scalar spectrum analysis across all three theories.","marker":"[26]"},{"why":"contains the authors' preliminary flavor-singlet results that this complete lattice study extends and finalizes.","marker":"[37, 38]"},{"why":"provides the previous eight-flavor spectra, chiral extrapolations, and analysis procedures that the new lightest-mass ensemble updates.","marker":"[17]"},{"why":"establishes the anomalous mass decomposition $M_{\\eta'}^2=M_{\\eta'}^{2(\\mathrm{anom})}+M_\\pi^2$ that separates the anomaly-driven part under study.","marker":"[39, 40]"}],"fun_headline_variants":["Eta-prime mass scales as Nf squared in lattice QCD","Anti-Veneziano scaling realized for eta-prime on the lattice","Square-law for eta-prime mass across 4, 8, 12 flavors","Multi-flavor lattice QCD sees eta-prime mass follow Nf^2","Gradient-flow scale uncovers eta-prime mass square-law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the topological charge density correlator is dominated by a single $\\eta'$ pole over the fitted distance and smearing windows; if a continuum or unparticle component contributes, the extracted exponential mass is a window-dependent effective quantity and the $n_f^2$ scaling could be an artifact of the fitting procedure rather than a property of the theory.","fun_headline_variants_meta":{"raw":{"variants":["Eta-prime mass scales as Nf squared in lattice QCD","Anti-Veneziano scaling realized for eta-prime on the lattice","Square-law for eta-prime mass across 4, 8, 12 flavors","Multi-flavor lattice QCD sees eta-prime mass follow Nf^2","Gradient-flow scale uncovers eta-prime mass square-law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":3019,"prompt_tokens":1110,"completion_tokens":1909,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":726,"completion_tokens_details":{"reasoning_tokens":1811}},"tokens_in":726,"tokens_out":1909,"duration_ms":13127,"temperature":1.0,"reasoning_tokens":1811,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:49:24.996024+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same correlators with a spectral function that includes a continuum or a second state and check whether the quoted $M_{\\eta'}$ values and the $1:4:9$ pattern survive; or extract $M_{\\eta'}$ on the same ensembles from the disconnected fermionic axial-current correlator and require agreement with the gluonic-operator result.","supporting_citations":[{"cited_title":"Walking on the Ladder: 125 GeV Technidilaton, or Conformal Higgs -Dedicated to the late Professor Yoichiro Nambu-","cited_arxiv_id":"1508.07688","evidence_quote":"supplies the anti-Veneziano-limit calculation predicting $M_{\\eta'}^2/\\Lambda_{\\rm IR}^2\\sim n_f^2$, the scaling law the lattice data are compared against."},{"cited_title":"Bando, K.-i","cited_arxiv_id":null,"evidence_quote":"derives the dilaton mass formula $M_\\sigma^2=d_0+d_1M_\\pi^2$ that organizes the scalar spectrum analysis across all three theories."}],"review_version":1}