{"id":"7fba1c9c-3470-4fef-9004-1ba16a3d3ff8","arxiv_id":"2508.09006","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A diagrammatic criterion based on quark-line disconnected diagrams identifies which multi-hadron excited states dominate contamination in hadron three-point functions in lattice QCD.","lead":"This paper explains why certain excited states, such as a nucleon plus a rho meson, can dominate the signal in lattice QCD calculations of hadron structure, even when those states are heavy. The explanation gives lattice practitioners a recipe for choosing which operators to include in their analyses, potentially reducing a leading source of systematic error.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing step is the partially-quenched mapping from Wick-contraction topology to physical intermediate states; without it, volume-enhanced disconnected diagrams in operator-space three-point functions do not identify the states that contaminate standard three-point functions.","rationale":"Good-faith reading: the paper's mechanism is a heuristic but reasonable way to single out excited-state contributions, supported by ChPT agreements and by the GEVP subtraction in the vector channel. The weakest point is not the volume enhancement itself, which has a transparent non-interacting origin and a cited two-volume check, but the step that converts diagram topology into physical-state selection. That step is imported from a partially quenched framework and is not proven here. The reader identified exactly this assumption, so I agree. The numerical demonstration in App. B is consistent but uses one 400-config ensemble and does not resolve Nrho at zero momentum; it is evidence, not proof. The proposed check—a variational spectral decomposition with and without the predicted current-enhanced operators—would settle the mapping directly. Because the central claim is testable and the existing evidence is suggestive, the conditional verdict stands unchanged.","tokens_in":11419,"tokens_out":17841,"duration_ms":218310,"concrete_test":"On the L=24, m_pi=429 MeV ensemble (or a higher-statistics copy), build variational bases {N, Nsigma, Npi} in the scalar isoscalar channel and {N, Nrho, Npipi} in the isovector vector channel. Solve the GEVP for the two-point matrix, extract the physical excited-state amplitudes A_f = Z_f^* Z'_f <f|J|i> from the standard (sink=N only) three-point function, and separately compute the quark-line disconnected/direct diagram (Eqs. 11-12) at p'_M=q. If the state with the largest |A_f| is always the one selected by a large disconnected diagram, the partially quenched mapping is validated; if any state with no disconnected-diagram support has comparable |A_f|, the criterion fails. Repeating on L=48 tests whether the V-enhancement controls the physical spectral weights.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central criterion—presence/absence of quark-line disconnected diagrams identifies which excited states dominate—rests on the claim that individual Wick contractions can be mapped to specific physical states in the spectral decomposition. The paper invokes a partially quenched framework [15] for this mapping, but in the interacting theory the spectral sum in Eq. (2) is not ordered by contraction topology. The same physical state can receive connected and disconnected contributions, and one contraction (e.g. the direct diagram WD, Eq. (11)) contains all states coupled to the current, not only the meson M; the paper itself notes that rho and sigma couple to their multi-particle decay products and that the dominant state depends on m_pi and L. The volume enhancement (L^3/a^3) of WD is a property of an operator matrix element with a two-hadron sink; it does not directly determine the amplitudes Z_f^* Z'_f <f|J|i> in the standard single-hadron three-point function. The ChPT comparisons in the axial and B->pi channels are consistency checks, not a derivation of the one-to-one mapping. If this mapping fails, the diagrammatic criterion selects the wrong operators and the central claim is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a general diagrammatic criterion for identifying which excited states contaminate lattice QCD three-point functions. The argument is that when a meson operator carries the same momentum as the inserted current, quark-line disconnected Wick contractions, such as the direct diagram in Eq. (11), are enhanced by a factor of order L^3/a^3 in the non-interacting limit, making the corresponding multi-hadron states ('current-enhanced states') dominant at the source-sink separations accessible in present calculations. The authors argue that the presence or absence of such disconnected diagrams therefore identifies which states must be included in multi-state fits or variational analyses, and they support this with examples in the nucleon axial, scalar, and vector channels, in B->pi semileptonic decays, and with a GEVP subtraction using N and N-rho operators in Appendix B. The paper is written as a letter and presents one numerical ensemble test, with additional agreement claimed with ChPT expectations and published lattice results.","tokens_in":11642,"tokens_out":5723,"duration_ms":66582,"significance":"If correct, the criterion would be practically valuable: it gives lattice practitioners a simple rule for choosing operators in variational analyses and for predicting which multi-hadron states dominate excited-state contamination, going beyond ChPT in channels where no chiral prediction exists. The Wick-contraction counting is internally consistent, the non-interacting volume-enhancement argument is transparent, and the GEVP subtraction in Fig. 2 visibly reduces the contamination in the vector channel, corroborating the claimed N-rho dominance. The framework also makes falsifiable predictions, such as the ineffectiveness of N-pi operators in the scalar channel and the importance of N-rho operators in the vector channel, which is a strength. The main weaknesses are that the one-to-one mapping between contraction topology and physical intermediate states is asserted rather than derived, and that the numerical support is limited to one small ensemble subset with several load-bearing checks not shown in the manuscript.","major_comments":[{"comment":"The central criterion that the presence or absence of quark-line disconnected diagrams identifies which excited states contaminate the three-point function rests on a one-to-one mapping between individual Wick contractions and physical intermediate states in the spectral decomposition, introduced via the partially quenched framework [15]. This mapping is not derived, and the spectral decomposition in Eq. (2) is a sum over physical states that is not organized by contraction topology. In particular, the current-meson factor <O_M J> inside W_D in Eq. (11) couples to every state with the meson quantum numbers, not only to the single-meson state, and the paper itself notes after Eq. (14) that resonances can couple to their decay products and that the dominant state depends on m_pi and L. Since the predictions for N-pi, N-sigma, N-rho, and B*-pi all depend on this mapping, the manuscript needs either a derivation of the mapping in the interacting theory or an explicit statement that it is a conjecture whose validation relies on the numerical evidence presented.","section":"Volume enhancement of current-enhanced states, Eq. (11) and following text"},{"comment":"The manuscript states that quark-line disconnected contributions are O(100) larger than connected ones on the L=24a ensemble and that a comparison of L=24a and L=48a volumes 'confirms' the L^3/a^3 scaling, but no numerical results, plots, or tables supporting either statement are shown in the paper. The citation to [14] does not point to a specific figure or table containing this volume-scaling check. Because the volume enhancement is the quantitative mechanism behind the claimed dominance at accessible separations, this evidence should either be presented in the manuscript or the claim should be softened to a heuristic expectation that is not yet numerically verified.","section":"Volume enhancement of current-enhanced states, text after Eq. (12)"},{"comment":"The GEVP-improved three-point function in Eq. (B4) omits the diagonal O_i=O_j=O_{N-rho} contribution, justified only by the statement that it is volume-suppressed because no current-enhanced diagram contributes. No numerical test of this suppression is given, yet the inference that N-rho is the dominant contaminant is drawn from the difference between the standard and this incomplete GEVP-improved correlator. Without demonstrating that the omitted term is negligible for the chosen values of t and tau, the flat plateau in Fig. 2 could be an artefact of an uncontrolled truncation rather than evidence that the physical N-rho state has been removed.","section":"Appendix B, Eqs. (B4)-(B5)"},{"comment":"The central numerical demonstration rests on a single subset of 400 gauge configurations at m_pi=429 MeV and L=24a, with the GEVP analysis shown for one value of t0=3a, no statistical uncertainty bands, no t0 variation, and no systematic check of the fit or operator basis. The claim that excited-state contamination is 'exponentially reduced' is drawn visually from one plot. Since the paper claims a general framework valid across channels and pion masses, either additional ensembles, multiple t0 choices, or a clear statement that the general claim is a prediction supported only by a single proof-of-principle is needed.","section":"Nucleon vector matrix elements and Appendix B"}],"minor_comments":[{"comment":"The notation 'OH and bar-OH' in the text is confusing; the sink interpolator should be written as O_H^dagger(p',t) or an equivalent explicit Hermitian conjugation, especially because Eq. (1) then uses O_H^dagger(p,0) inconsistently.","section":"Introduction, Eq. (1)"},{"comment":"The double-bracket notation <<...>>_F denoting the gauge average of the Wick contraction is not defined in the text; it should be introduced explicitly.","section":"Eq. (17)"},{"comment":"The figure would benefit from error bands on the ratios and from a clear statement of how the statistical uncertainty is estimated, since the main claim of 'exponentially reduced' contamination is currently based on visual inspection of point values.","section":"Figure 2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know about this one: Barca proposes a simple rule for which excited states contaminate hadron three-point functions—look for quark-line disconnected diagrams. When the meson momentum matches the current momentum, those diagrams are enhanced by the volume, and the states they couple to (N rho, N pi, N sigma) dominate the excited-state contamination at accessible source-sink separations. The new piece is the general Wick-contraction criterion, which explains why N pi works in the axial channel but fails in the scalar channel, and why N rho—despite being heavy—dominates the isovector vector channel.\n\nWhat's good: the App. B test is the real evidence. On a 2.4 fm, 429 MeV ensemble, a GEVP with N and N rho operators visibly removes the contamination in both forward and off-forward vector ratios. That is an independent confirmation that the rho's coupling to the vector current makes N rho the dominant contamination, even though N pi and N pi pi are lighter. The ChPT comparisons in axial and scalar channels are consistent with the criterion. The paper is honest about the dependence on m_pi and L.\n\nSoft spots: the general mapping from contraction topology to physical intermediate states is invoked via a partially quenched framework [15], but not derived. In the interacting theory, the spectral sum is not ordered by Wick contraction, so the one-to-one association is a heuristic, not a theorem. That said, the App. B result works without the mapping—it uses an explicit N rho operator—so the central mechanism is not circular. The volume-scaling check is mentioned but not shown; the reader has to trust [14] for that. The numerics are thin: 400 configurations on one ensemble with a heavy pion mass, where the rho is essentially stable. The axial a1 discussion is speculative. And \"prove\" is too strong for a single-ensemble demonstration.\n\nWho it's for: lattice practitioners working on hadron structure and excited-state systematics. It gives a practical operator-selection rule and explains why some standard variational analyses fail to remove contamination. It deserves a serious referee. The referee should push for the volume-scaling data, a lighter-pion check, and softer claims, but the core idea is testable and likely to be useful.\n\nMy verdict: worth engaging. I'd send it to peer review, and I'd cite it if I were working in this area.","headline":"A new diagrammatic criterion for excited-state contamination in lattice three-point functions, with a compelling vector-channel test; the general mapping from contractions to states is plausible but not fully proven.","tokens_in":12130,"tokens_out":8091,"would_cite":true,"duration_ms":90373,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the inserted current and the kinematics select which excited states contaminate hadron three-point functions, through quark-line disconnected diagrams whose contributions grow with the spatial volume.","keywords":["lattice QCD","excited-state contamination","three-point functions","quark-line disconnected diagrams","variational method","current-enhanced states","nucleon matrix elements","vector meson dominance"],"falsifier":"Take a channel predicted to be current-enhanced, e.g. nucleon isovector vector with $N\\rho$, and compute the three-point function at a lighter pion mass where the $\\rho$ lies above the $\\pi\\pi$ P-wave threshold; if the dominant contamination is still $N\\rho$-like rather than $N\\pi\\pi$, or if the quark-line disconnected term does not scale as $L^3$ when two volumes are compared, the argument would be refuted.","tokens_in":11181,"feed_emoji":"⚛️","tokens_out":5836,"duration_ms":60265,"temperature":0.7,"pith_summary":"This paper sets out to explain why certain excited states, not necessarily the lightest ones, dominate excited-state contamination in lattice QCD three-point functions. It argues that the inserted current preferentially creates specific multi-hadron states—$N\\pi$, $N\\sigma$, $N\\rho$, $B^*\\pi$—and that these current-enhanced states enter through quark-line disconnected diagrams whose amplitude is not volume-suppressed. A variational analysis with operators built from the nucleon plus the relevant meson removes the contamination, and numerical evidence shows the predicted states dominate at accessible source–sink separations. If correct, the criterion tells practitioners which operators to include in multi-state fits and which kinematics to avoid, reducing a leading systematic uncertainty in hadron structure calculations.","feed_headline":"Current picks which excited states contaminate lattice QCD fits","feed_subtitle":"A diagrammatic rule identifies the multi-hadron states—Nπ, Nσ, Nρ—to include in three-point fits.","key_machinery":"The mechanism is the current–meson two-point function appearing in quark-line disconnected (direct and disconnected) Wick contractions. A current $J$ inserted at time $\\tau$ can couple directly to a meson with the current's quantum numbers, creating a multi-hadron state with the nucleon; when that meson carries the same momentum as the current transfer, translational invariance gives a spatial-volume factor that is absent in connected diagrams. The variational method (GEVP) with a basis of single- and multi-hadron operators constructs improved operators that overlap with the selected state, so the enhanced contribution can be isolated and subtracted; the paper's numerical tests use this to remove $N\\rho$ contamination from nucleon vector matrix elements.","core_discovery":"The central discovery is a diagrammatic criterion for excited-state contamination: in a hadron three-point function, the excited states that matter are those whose Wick contraction topologies include quark-line disconnected contributions, specifically the direct term where the current creates a meson at time $\\tau$ and the disconnected term where the current is vacuum-averaged. When the meson momentum equals the current momentum transfer, $\\vec p'_M=\\vec q$, these terms are enhanced by the spatial volume $L^3/a^3$ in the non-interacting limit, and numerically are $\\mathcal{O}(100)$ larger than the connected term. Because they decay faster in Euclidean time than ground states, they are exponentially suppressed only at large separations; at the intermediate separations accessible today they dominate. The paper shows, in nucleon axial, scalar and vector channels and in $B\\to\\pi$, that this mechanism identifies $N\\pi$, $N\\sigma$ and $N\\rho$ as the dominant contaminants, even when those states are heavy, and explains why including only $N\\pi$ operators helps in some channels and not others.","pith_inferences":["One testable extension is to use the same diagrammatic criterion to design kinematics with $\\vec p'_M\\neq\\vec q$, suppressing the enhanced contamination before fitting, as a cross-check on the variational subtraction.","The volume-scaling mechanism implies that finite-volume corrections to hadron matrix elements may be larger for current-enhanced channels, so comparisons between volumes or to the infinite-volume limit should watch for enhanced $L^3$ effects in the excited-state sector.","Because the argument only needs a current and a Wick decomposition, it should carry over to four-point functions and to gluonic currents, where the connected contribution is absent and the disconnected one may dominate from even shorter times."],"forward_implications":["The presence or absence of quark-line disconnected diagrams in a three-point function directly indicates which excited states are current-enhanced and must be handled in fits.","In nucleon axial channels, $N\\pi$ states dominate at accessible separations when kinematics allow the pion to carry the current momentum, matching chiral effective theory; in the forward limit $N a_1$-like states may matter.","In the isovector vector channel, $N\\rho$ states dominate even though their energy is well above $N\\pi$ and $N\\pi\\pi$, so variational bases that omit $N\\rho$ leave residual contamination.","In the isoscalar scalar channel $N\\sigma$ states (or $N\\pi\\pi$ in S-wave at lighter pion masses) dominate, while $N\\pi$ operators have no effect because they lack the relevant disconnected diagram.","For $B\\to\\pi\\ell\\nu$, $B^*\\pi$ states are the expected current-enhanced contaminants."],"supporting_citations":[{"why":"Supplies the numerical check that $N\\sigma$ states dominate the isoscalar scalar channel at $m_\\pi=429$ MeV.","marker":"[14]"},{"why":"Provides the partially quenched mapping of individual contractions to intermediate states, the load-bearing connection.","marker":"[15]"},{"why":"Shows $N\\pi$ operators improve nucleon matrix elements only in channels with quark-line disconnected contributions, the comparison the paper explains.","marker":"[17]"},{"why":"The generalized eigenvalue method used to diagonalise the two-point system and construct improved operators.","marker":"[12]"},{"why":"Justifies the exponential suppression of omitted states when a truncated variational basis is used.","marker":"[13]"},{"why":"Reports numerical evidence from variational studies of transition matrix elements consistent with the enhanced contributions.","marker":"[16]"},{"why":"Supplies the $N\\to N\\pi$ matrix element framework and evidence for $N\\pi$ contamination in axial and boosted kinematics.","marker":"[18]"},{"why":"Gives a chiral-effective-theory prediction that $N\\pi$ states contribute large axial-channel contamination, used as a consistency check.","marker":"[23]"},{"why":"Predicts $B\\pi$ excited-state contamination in $B$-meson correlators, the analogue used for $B\\to\\pi\\ell\\nu$.","marker":"[34]"}],"fun_headline_variants":["Current-meson rule exposes Nπ, Nσ, Nρ contamination in QCD fits","Volume boost in current couplings picks leading excited states","Which excited states leak? Current momentum decides in lattice QCD","Disconnected diagrams reveal true excited-state burdens in lattice QCD"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole criterion rests on assuming that each Wick contraction topology can be mapped one-to-one to a physical intermediate state in the interacting theory; if disconnected diagrams mix with other states, the diagrammatic rule may name the wrong contaminant.","fun_headline_variants_meta":{"raw":{"variants":["Current-meson rule exposes Nπ, Nσ, Nρ contamination in QCD fits","Volume boost in current couplings picks leading excited states","Which excited states leak? Current momentum decides in lattice QCD","Disconnected diagrams reveal true excited-state burdens in lattice QCD"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1410,"prompt_tokens":846,"completion_tokens":564,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":490}},"tokens_in":462,"tokens_out":564,"duration_ms":6767,"temperature":1.0,"reasoning_tokens":490,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:31:43.709888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a channel predicted to be current-enhanced, e.g. nucleon isovector vector with $N\\rho$, and compute the three-point function at a lighter pion mass where the $\\rho$ lies above the $\\pi\\pi$ P-wave threshold; if the dominant contamination is still $N\\rho$-like rather than $N\\pi\\pi$, or if the quark-line disconnected term does not scale as $L^3$ when two volumes are compared, the argument would be refuted.","supporting_citations":[{"cited_title":"Nuc- leon sigma terms with a variational analysis from Lattice QCD","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical check that $N\\sigma$ states dominate the isoscalar scalar channel at $m_\\pi=429$ MeV."},{"cited_title":"Parton degrees of freedom in pdfs from the hadronic tensor and large momentum effective theory","cited_arxiv_id":null,"evidence_quote":"Provides the partially quenched mapping of individual contractions to intermediate states, the load-bearing connection."},{"cited_title":"Future methods for studying confined gluons in QCD","cited_arxiv_id":null,"evidence_quote":"Shows $N\\pi$ operators improve nucleon matrix elements only in channels with quark-line disconnected contributions, the comparison the paper explains."},{"cited_title":"On the generalized eigenvalue method for energies and matrix elements in lattice field theory","cited_arxiv_id":null,"evidence_quote":"The generalized eigenvalue method used to diagonalise the two-point system and construct improved operators."},{"cited_title":"On the computation of hadron-to-hadron transition mat- rix elements in lattice QCD.JHEP, 01:140, 2012","cited_arxiv_id":null,"evidence_quote":"Justifies the exponential suppression of omitted states when a truncated variational basis is used."},{"cited_title":"Progress on nucleon transition matrix elements with a lattice QCD variational analysis","cited_arxiv_id":null,"evidence_quote":"Reports numerical evidence from variational studies of transition matrix elements consistent with the enhanced contributions."},{"cited_title":"Investigation of pion- nucleon contributions to nucleon matrix elements.Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the $N\\to N\\pi$ matrix element framework and evidence for $N\\pi$ contamination in axial and boosted kinematics."},{"cited_title":"Axial Vector Form Factors from Lat- tice QCD that Satisfy the PCAC Relation.Phys","cited_arxiv_id":null,"evidence_quote":"Gives a chiral-effective-theory prediction that $N\\pi$ states contribute large axial-channel contamination, used as a consistency check."},{"cited_title":"Pion–Nucleon Sigma Term from Lattice QCD","cited_arxiv_id":null,"evidence_quote":"Predicts $B\\pi$ excited-state contamination in $B$-meson correlators, the analogue used for $B\\to\\pi\\ell\\nu$."}],"review_version":2}