{"id":"b2004e08-56d4-414b-9489-5a6ba567ae4c","arxiv_id":"2411.13782","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New next-to-leading-order mass-dependent results for the flowed quark condensate, with small- and large-mass expansions and a full numerical evaluation, supply the perturbative input for a proposed gauge-invariant lattice quark mass determination.","lead":"This paper proposes a new way to measure quark masses by comparing lattice measurements of a smoothed quark condensate with new perturbative calculations. It computes the needed mass-dependent corrections at next-to-leading order in the strong coupling and develops an expansion method valid for both small and large values of the flow-time mass parameter.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (8)-(9) rest on an unproven residue expansion (Eqs. (4)-(7)); meromorphy, arc-vanishing, and pole completeness are not demonstrated, and the Fig. 1 check is limited-range, so a missed singularity or residue error would invalidate the new O(alpha_s) coefficients.","rationale":"Good-faith reading: the paper proposes a lattice method for heavy-quark masses and supplies one perturbative ingredient, the O(alpha_s) mass-dependent flowed condensate. The claim is scoped honestly in the body; the lattice feasibility and higher orders are explicitly left open. The reader's CONDITIONAL verdict is appropriate, and my analysis lands on the same load-bearing assumption: the Laplace-residue expansion (Eqs. (4)-(7)) is the foundation of every new analytic number in Eqs. (8)-(9), and it is asserted with details deferred. I identify no additional fatal flaw; the internal consistency of the analytic expansions with the ftint numerics (Fig. 1) and the availability of ftint as published, reproducible code (Ref. [20]) count as genuine supporting evidence, as does the known closed-form free-theory result from Ref. [17] that anchors the leading terms. The residual risk is a missed singularity or a residue-arithmetic error that a limited-range plot-based comparison could conceal, particularly in the intermediate m^2t region where only the numerics apply. Because the concern is absence of a checkable derivation rather than any identified error, the verdict should not change: CONDITIONAL (UNCHANGED). Secondary issues—the abstract's 'two-loop level' versus the body's 'O(alpha_s)', and the lack of precision statements on the constants in Eqs. (8)-(9)—are worth fixing but do not affect the central claim. My requested test, an independent expansion-by-regions derivation of the small-m^2t coefficient plus high-precision numerical checks of the large-mass coefficients, would settle the concern if it lands.","tokens_in":9906,"tokens_out":23451,"duration_ms":221970,"concrete_test":"Independently recompute the full O(alpha_s) coefficient of the m^2t term in Eq. (8) using expansion-by-regions, which the paper itself states is valid in the small-m^2t limit (only the large-m^2t limit is problematic): the coefficient including the log(m^2t)*log(m^2/mu^2) terms should be compared digit-by-digit with 2.3334, -2.16804, 2.22817, -0.95493. Separately, verify the three O(alpha_s) large-mass coefficients in Eq. (9) by high-precision (relative error < 1e-4) numerical evaluation with ftint or an independent integrator at 8m^2t = 50, 100, 200, where the 1/(m^2t) and 1/(m^2t)^2 series converges quickly; any quoted digit that shifts beyond rounding demonstrates that the residue expansion has missed a singularity or arc contribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2's central result is that the O(alpha_s) mass-dependent flowed quark condensate is known for a wide range of m^2t, via the new analytic expansions (8)-(9) and the ftint numerics. All analytic content follows from the Laplace-residue method of Eqs. (4)-(7): after the Laplace transform, the inverse contour (5) is closed and the integral is asserted to equal the residue sums (6)-(7). The derivations are deferred ('Further details will be explained in a future publication'). For Eqs. (6)-(7) to be correct for each of the eight scalar integrals, three conditions must hold: (i) the Mellin transform tilde-I(s,t) is meromorphic in the relevant half-plane with only pole singularities; (ii) the contribution of the arc at |s| -> infinity vanishes, so the integrand decays faster than (m^2t)^s grows; and (iii) the catalogue of singularities contributing in Eqs. (6)-(7) is complete. None of these is demonstrated. If a branch cut in the s-plane (e.g., from a threshold-like branch point in m^2t) crosses the closure contour, or if an arc contribution fails to vanish, then the quantities claimed as new—the O(alpha_s) numbers 0.759581, 2.3334, 2.16804, 2.22817, 0.95493 in Eq. (8) and the entire O(alpha_s) content of Eq. (9)—would be incomplete, and the paper's central claim would fail. The only check offered is the agreement of the expansions with ftint shown in Fig. 1. That agreement is meaningful but bounded: it tests the final ratio over a limited flow-time window (especially narrow in the charm case, where only the small-m^2t expansion applies for a modest range), it carries no stated numerical precision, and it is a joint test of the two computations, so a common-mode error in the residue arithmetic could be masked by plot resolution. The paper also contains an internal wording conflict (abstract 'two-loop level' versus body 'O(alpha_s)') and the constants in Eqs. (8)-(9) have no precision statement; these are secondary.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper proposes a new method to determine heavy-quark MS-bar masses from lattice QCD by matching the ratio of flowed quark condensates, Eq. (1), to perturbation theory. The perturbative input is the O(alpha_s) mass-dependent flowed quark condensate, which the authors compute in two complementary ways: analytic expansions in m^2t and 1/(m^2t), Eqs. (8)-(9), and a full numerical evaluation using ftint. The paper claims that the next-to-leading small-m^2t O(alpha_s) term in Eq. (8) and all O(alpha_s) terms in Eq. (9) are new, and it shows results for charm and bottom quarks in Fig. 1.","tokens_in":10150,"tokens_out":3151,"duration_ms":33204,"significance":"If the results are correct, they provide a genuinely useful perturbative ingredient for a novel, gauge-invariant, one-point-function-based lattice determination of heavy-quark masses. The proposed Laplace-residue expansion technique for massive gradient-flow integrals is also of methodological interest, and the explicit numerical cross-check with ftint is a strength. The novelty claims for Eqs. (8)-(9) are concrete and falsifiable, and the paper is refreshingly direct about what remains to be done. However, the central technical tool is only asserted, not demonstrated, so the significance is conditional on the missing derivation.","major_comments":[{"comment":"The inverse-Laplace residue expansion is the load-bearing element of the paper, but its validity is not established. After Eq. (5), the text simply says the contour is closed and the result is given by residue sums, with 'Further details will be explained in a future publication.' For Eqs. (6)-(7) to hold, one must show that the Laplace-transformed integrals are meromorphic in the relevant half-plane, that the arc contributions at |s| -> infinity vanish, and that the catalog of singularities is complete for each of the eight scalar integrals. None of these conditions is demonstrated. Since every new O(alpha_s) coefficient in Eqs. (8)-(9) depends on this method, the central claim is not yet fully supported. Please include the derivation, or at least a rigorous statement of the analyticity assumptions and a proof sketch, in the manuscript.","section":"Section 2, Eqs. (4)-(7)"},{"comment":"The agreement between the expansions and the ftint evaluation shown in Fig. 1 is a useful internal cross-check, but it is not an independent validation of the residue method: ftint evaluates the same scalar integrals that the expansions are supposed to represent. The plotted flow-time range is also limited, and the large-m^2t expansion is not shown for the charm case. The text states that Fig. 1 provides a 'convincing verification' of the expansion method; this overstates what a self-consistency check in a finite window can establish. Please quantify the plotted range, the order of neglected terms, the numerical uncertainties of ftint, and, if possible, test the expansions against an independent method (for example, expansion by regions or direct high-precision integration) at selected values of m^2t.","section":"Section 2, Fig. 1 and surrounding text"},{"comment":"The novelty claims--'the next-to-leading term in m^2t at O(alpha_s) of Eq. (8) is new' and 'in Eq. (9), all the O(alpha_s) terms are new'--cannot be verified from the manuscript because no intermediate steps are shown. A reader cannot see which residues produce which coefficients, how the MS-bar renormalization of the flowed operator and mass is implemented beyond the quoted factors R_chi and Z_m^MS, or how the numerical values such as 0.759581 and 2.3334 are obtained. For a proceedings contribution this level of detail may be acceptable, but for a refereed journal article the derivation must be present or the claims must be explicitly labeled as relying on a forthcoming publication.","section":"Section 2, Eqs. (8)-(9) and the paragraph following them"}],"minor_comments":[{"comment":"The momentum-space integration measure is garbled in the displayed integrals; please typeset the d-dimensional measures ∫_p ∫_k clearly so that the reader can distinguish the loop integrations from the flow-time integrations.","section":"Eq. (3)"},{"comment":"The inequalities '0.1 ≪ 8m_c^2 t ≪ 20' and '1.0 ≪ 8m_b^2 t ≪ 200' are written without explicitly stating the unit convention; adding a sentence that t has mass dimension -2 and m^2t is dimensionless would help the reader.","section":"Section 1, flow-time window"},{"comment":"The figure axes and legends are difficult to read in the current rendering, and it is not clear over which t-interval each expansion is expected to be accurate. Please improve the figure and add a statement of the approximate validity ranges of the small- and large-m^2t expansions.","section":"Fig. 1"},{"comment":"The footnotes marked 'one.sup' and 'two.sup' in the text are not rendered as ordinary footnotes, and the reference to the latest arXiv version of Ref. [17] is made only in a footnote; please incorporate these corrections into the reference list and remove the placeholder superscripts.","section":"References and footnotes"}],"recommendation":"major_revision","confidential_remarks":"This is clearly a proceedings-style contribution, and the authors are transparent that full derivations will appear elsewhere. For a journal version, the deferred proof of the Laplace-residue method is the main obstacle; if the journal does not require self-contained derivations, the paper could be acceptable after the presentation issues are fixed. I also note that the numerical check in Fig. 1 is internal to the same computation, so the referee should not be asked to treat it as independent verification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Name],\n\nPunchline: this proceedings paper contains a genuinely new perturbative result—the O(alpha_s) mass-dependent term for the flowed quark condensate, in both small- and large-m^2t expansions, plus a full numerical evaluation. If the coefficients are right, this is a useful building block for a proposed gauge-invariant, noise-suppressed heavy-quark mass scheme. The paper is honestly scoped: it calls the computation a first step and explicitly lists lattice feasibility and higher orders as open issues.\n\nWhat's good: the authors derive the expansions and cross-check them against an independent numerical evaluation (ftint) in Fig. 1, with agreement in the charm and bottom windows. That is real internal evidence, not fitting. They also clearly mark which terms are new (Eq. 8's NLO small-mass term, all of Eq. 9's NLO terms) and which come from earlier work. The citation pattern is clean; the self-citations are to published LO results.\n\nNow the soft spots, in proportion. The load-bearing expansion method—Laplace transform, contour closure, residue sums in Eqs. (4)–(7)—is asserted, with details deferred to 'a future publication.' The stress-test note is correct that meromorphy, arc vanishing, and pole completeness are not demonstrated here. That means the new numerical constants rest on an unproven method. The Fig. 1 check is meaningful but bounded: it is a joint test of the analytic and numerical pipelines, covers a limited flow-time window (especially narrow for charm), and carries no stated precision. A common-mode error could hide, though I think that is unlikely given the method's analogy to prior work and the standard structure of the integrals.\n\nMinor: the abstract says 'two-loop level' while the body consistently says O(alpha_s); that wording should be fixed. The constants in Eqs. (8)–(9) lack precision statements, a common proceedings omission.\n\nVerdict: solid, worth a serious referee. The referee's main request should be the expansion-method derivation, either in an appendix or by pointing to the promised follow-up. The paper is not a complete proof, but the result is new, useful, and probably correct. I'd take it for peer review and ask for the missing details.","headline":"New O(alpha_s) mass-dependent condensate terms with a promising but unproven expansion method; solid proceedings paper worth refereeing.","tokens_in":10940,"tokens_out":3660,"would_cite":true,"duration_ms":32889,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T25","81T13","81V05"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes that the flowed quark condensate, computed at O(alpha_s) with full $m^2t$ dependence, supplies the perturbative input needed to extract heavy-quark MS-bar masses by matching a gauge-invariant ratio of flowed quark…","keywords":["gradient flow","quark mass determination","quark condensate","heavy quarks","perturbative expansion","lattice QCD","Laplace transform","MS-bar scheme"],"falsifier":"Evaluate one of the eight integrals in Eq. (3) at an intermediate value such as $m^2t=10$ with an independent high-precision numerical integrator and compare with the large-$m^2t$ expansion of Eq. (9) truncated at the $1/(m^2t)^2$ term; a disagreement larger than the $O(1/(m^2t)^3)$ remainder would falsify the residue-expansion method.","tokens_in":9522,"feed_emoji":"⚛️","tokens_out":16686,"duration_ms":109856,"temperature":0.7,"pith_summary":"The paper proposes determining heavy-quark masses from the gradient flow by comparing the lattice value of a ratio of flowed quark condensates to perturbation theory. The ratio is gauge invariant, involves only one-point functions, and is finite after renormalization, which should give clean control of perturbative errors and suppressed lattice noise. The missing ingredient was the quark-mass dependence of the flowed quark condensate; the paper computes it at $O(\\alpha_s)$ in both the small-$m^2t$ and large-$m^2t$ limits, and also numerically over the full mass range. The next-to-leading small-mass term and all the $O(\\alpha_s)$ terms of the large-mass expansion are new. If the result is right, lattice data on this ratio can be matched to determine MS-bar charm and bottom masses at next-to-leading order.","feed_headline":"Flowed quark condensate expansions open a new path to quark masses","feed_subtitle":"One-point functions now have the next-to-leading-order input for MS-bar quark masses from lattice data.","key_machinery":"The central mechanism is the Laplace-transform residue expansion for massive gradient-flow integrals. Given an integral $I(m^2,t)$, one forms $\\tilde I(s,t)=\\int_0^\\infty d(m^2)\\,(m^2)^{-s-1}I(m^2,t)$, inverts with a contour integral, and obtains the $m^2t\\ll 1$ series as a residue sum over positive singularities and the $m^2t\\gg 1$ series as a residue sum over negative singularities of $\\tilde I(s,t)(m^2)^s$. This bypasses the failure of expansion by regions in the large-mass limit, where the flow-time integral covers hard and soft momentum regions simultaneously. The observable it is applied to is the ratio $R(t,m_1,m_2)=\\langle\\bar{\\chi}_1\\chi_1\\rangle/\\langle\\bar{\\chi}_2{\\overleftrightarrow{D}}\\chi_2\\rangle$, which is UV finite because the only divergent part of the flowed fermion fields is the wave-function renormalization.","core_discovery":"At $O(\\alpha_s)$, the flowed quark condensate $\\langle[\\bar{\\chi}(t,x)\\chi(t,x)]_R\\rangle$ is given by Eq. (8) for $m^2t\\ll 1$ and by Eq. (9) for $m^2t\\gg 1$. In the small-mass expansion the coefficient of the $m^2t$ term at $O(\\alpha_s)$ is new; in the large-mass expansion all $O(\\alpha_s)$ terms are new. These analytic results are supplemented by a numerical evaluation with full $m^2t$ dependence, so the perturbative prediction for the ratio $R(t,m_1,m_2)$ can be evaluated for charm and bottom quarks across their physical flow-time windows. The calculation uses a Laplace transform in $m^2$ followed by a contour inversion in the transform variable; the small- and large-$m^2t$ series are obtained from residues at positive and negative singularities, respectively, avoiding the failure of expansion by regions when the flow-time integral spans hard and soft momenta.","pith_inferences":["The Laplace-transform residue technique should transfer directly to the denominator observable $\\langle\\bar{\\chi}\\overleftrightarrow{D}\\chi\\rangle$ and to $O(\\alpha_s^2)$ calculations, which would raise the matching precision to next-to-next-to-leading order.","Because expansion by regions is known to fail in the large-$m^2t$ flow-time case, the same method may also unlock analytic expansions of other massive gradient-flow quantities, such as flowed observables involving massive quarks.","A stronger check than numerical self-consistency would be to apply the residue expansion to a simpler massive gradient-flow integral with a known closed form; passing that test would harden the evidence before the detailed proof appears.","If lattice data reach the required precision, this gauge-invariant one-point observable could serve as an independent cross-check of existing lattice averages for charm and bottom MS-bar masses, which rest on different renormalization schemes."],"forward_implications":["Lattice groups can match $R(t,m,0)$ or $R(t,m,m)$ to the new $O(\\alpha_s)$ expressions and extract MS-bar charm and bottom masses at next-to-leading order, once sufficiently precise flowed-condensate data exist.","The large-$m^2t$ expansion is valid in the physical bottom-quark window $1.0 \\ll 8m_b^2t \\ll 200$, where the previously known small-mass results do not apply.","The numerical full-mass-dependence result covers the intermediate $m^2t$ region, so a mass determination does not depend on which asymptotic expansion happens to converge.","Because the observable is gauge invariant, the leading nonperturbative corrections start at dimension-four condensates, giving a controlled perturbative uncertainty that gauge-dependent schemes such as RI-MOM lack.","The same framework should extend to light-quark masses, as the paper notes, since the ratio does not require the quark to be heavy."],"supporting_citations":[{"why":"Supplies the wave-function renormalization $Z_\\chi$ and the argument that only this factor is UV divergent, making the ratio in Eq. (1) finite.","marker":"[8]"},{"why":"Provides the $O(\\alpha_s)$ linear-mass term for the flowed quark condensate that Eq. (8) extends to the next order in $m^2t$.","marker":"[12]"},{"why":"Supplies the previous small-mass $O(\\alpha_s^2)$ results and the gradient-flow higher-order techniques that frame the calculation.","marker":"[13]"},{"why":"Defines expansion by regions, the standard technique whose failure in the large-$m^2t$ flow-time case motivates the new Laplace-transform approach.","marker":"[16]"},{"why":"Provides the $O(\\alpha_s^0)$ terms of the expansions in Eqs. (8)-(9) and describes the earlier application of expansion by regions to gradient-flow integrals.","marker":"[17]"},{"why":"Offers the scale-setting and momentum-flow idea analogous to the Laplace-transform residue expansion used here.","marker":"[18]"},{"why":"Together with [18], supplies the analogous residue and Laplace idea underlying the new expansion method.","marker":"[19]"},{"why":"The numerical evaluation tool described there is used to compute the full-mass-dependence result and to verify the small- and large-$m^2t$ expansions.","marker":"[20]"}],"fun_headline_variants":["Gradient flow yields new quark mass route","Two-loop flow condensate paves way to quark masses","New expansions make flowed quark condensate precise","Flow-time expansions unlock quark mass determination"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that closing the inverse-Laplace contour and summing the residues at the singularities of $\\tilde I(s,t)$ reproduces the full massive gradient-flow integral for small and large $m^2t$, a claim whose proof is deferred to a future publication and which is only checked by numerical self-consistency.","fun_headline_variants_meta":{"raw":{"variants":["Gradient flow yields new quark mass route","Two-loop flow condensate paves way to quark masses","New expansions make flowed quark condensate precise","Flow-time expansions unlock quark mass determination"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000179,"raw_usage":{"total_tokens":1284,"prompt_tokens":910,"completion_tokens":374,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":317}},"tokens_in":526,"tokens_out":374,"duration_ms":854825,"temperature":1.0,"reasoning_tokens":317,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:54:20.911643+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate one of the eight integrals in Eq. (3) at an intermediate value such as $m^2t=10$ with an independent high-precision numerical integrator and compare with the large-$m^2t$ expansion of Eq. (9) truncated at the $1/(m^2t)^2$ term; a disagreement larger than the $O(1/(m^2t)^3)$ remainder would falsify the residue-expansion method.","supporting_citations":[{"cited_title":"Chiral symmetry and the Y ang–Mills gradien t flow,","cited_arxiv_id":null,"evidence_quote":"Supplies the wave-function renormalization $Z_\\chi$ and the argument that only this factor is UV divergent, making the ratio in Eq. (1) finite."},{"cited_title":"Lattice energy–momentum tens or from the Y ang–Mills gradient flow—inclusion of fermion fields,","cited_arxiv_id":null,"evidence_quote":"Provides the $O(\\alpha_s)$ linear-mass term for the flowed quark condensate that Eq. (8) extends to the next order in $m^2t$."},{"cited_title":"Results and techniques for higher order calculations within the gradient-flow formali sm,","cited_arxiv_id":null,"evidence_quote":"Supplies the previous small-mass $O(\\alpha_s^2)$ results and the gradient-flow higher-order techniques that frame the calculation."},{"cited_title":"Quantum electrodynamics on t he lattice and numerical perturbative computation of g− 2,","cited_arxiv_id":null,"evidence_quote":"Together with [18], supplies the analogous residue and Laplace idea underlying the new expansion method."},{"cited_title":"ftint: Calculating gradient-flow integrals with pySecDec,","cited_arxiv_id":null,"evidence_quote":"The numerical evaluation tool described there is used to compute the full-mass-dependence result and to verify the small- and large-$m^2t$ expansions."}],"review_version":1}