{"id":"4cdba477-617f-40ea-823a-7c5258d11b96","arxiv_id":"2506.03002","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A high-statistics lattice-QCD calculation reports the kaon gluon momentum fraction as 0.557(18)(24)(56) at 2 GeV in the MS scheme, about twice the pion value from the same ensemble.","lead":"Physicists used supercomputer simulations of quarks and gluons to measure how gluons share momentum inside a kaon, a meson containing a strange quark. They report the kaon gluon momentum fraction as about 56 percent at 2 GeV, roughly twice the pion value from the same lattice ensemble.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline moment and the 1.92(18) kaon/pion ratio both rest on an uncomputed quark-gluon mixing term; the 10% mixing error may not cover the 20% upper end cited in the paper.","rationale":"The reader's weakest_assumption is exactly the missing Z_gq term, and I agree that it is the most load-bearing assumption. The bare measurement and the renormalization factor are transparently reported; the smearing study is extensive; and the Pz = 1 exclusion is openly discussed with a physical rationale. The one place where the quoted uncertainty is not derived from the calculation is the mixing term. If that term is at the upper end of the range the paper itself cites (20%), the central number shifts by roughly twice the assigned error, and the kaon/pion ratio—quoted with statistical errors only—would move by enough to change the physical conclusion. This is not an ad hominem or a dispute with the lattice method; it is a missing internal computation. The continuum extrapolation and smearing systematics are also real limitations, but they are explicitly caveated as future work and do not undercut the stated ensemble-level claim as directly. The proposed test—computing the disconnected quark contribution and the full 2×2 renormalization on the same ensemble—would settle the question. Until then, I leave the reader's conditional verdict unchanged: the analysis is careful and reproducible in structure, but the central number's uncertainty depends on an external estimate.","tokens_in":16652,"tokens_out":10527,"duration_ms":130968,"concrete_test":"Compute the disconnected quark-loop contribution to the O_g^OPE matrix element on a subset of the 1013 a12m310 configurations, renormalize with the full 2×2 RI/MOM mixing matrix including Z_gq from Ref. [20], and evaluate the omitted term Z_gq⟨x⟩_q^bare in Eq. (14) at μ = 2 GeV. If that term is larger in magnitude than 0.056 (10% of 0.557), the quoted central value and the ratio 1.92(18) move beyond the assigned mixing uncertainty, and the conditional verdict should be tightened. If it is below 0.056, the 10% systematic is adequate and the central claim stands as reported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (14) is the load-bearing renormalization relation: ⟨x⟩_g^MS = Z_gg^MS ⟨x⟩_g^bare + Z_gq^MS ⟨x⟩_q^bare. The paper sets Z_gq⟨x⟩_q^bare to zero and assigns a 10% systematic to the final moment, citing previous work reporting mixing between 2% and 20% (Refs. [18,20,31]). A 20% mixing contribution would shift 0.557 by about ±0.11, nearly twice the quoted 0.056 mixing error, so the stated uncertainty depends on choosing the 10% rather than the upper end of the cited range. The same omission enters the PDF matching in Eq. (15), where the R_gq kernel is dropped. Moreover, the headline ratio ⟨x⟩_g^K/⟨x⟩_g^π = 1.92(18) is quoted with only the statistical error, so the unknown mixing term is not propagated into the ratio at all. The central value, the ratio, and the PDF comparison with DSE therefore all rest on an external estimate rather than a computed quantity on this ensemble. This is a correctness risk, not a disagreement with consensus: the quoted number may be right, but the uncertainty assigned to the dominant neglected term is not derived from the calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a high-statistics lattice-QCD calculation of the kaon gluon momentum fraction and gluon PDF on a single MILC 2+1+1 HISQ ensemble at a≈0.12 fm and Mπ≈310 MeV. Using 1,296,640 kaon-correlator measurements, the authors extract bare gluon matrix elements via two-state fits to three-point correlators, study the dependence on HYP, Stout, and Wilson smearing, and quote a nonperturbatively renormalized kaon gluon momentum fraction ⟨x⟩_g^{MS,K}=0.557(18)_stat(24)_NPR(56)_mixing at μ=2 GeV in the MS scheme. They also extract the kaon gluon PDF from reduced pseudo-Ioffe-time distributions (RpITDs) using pseudo-PDF matching, compare it with a DSE prediction and with a pion gluon PDF from the same ensemble, and report a kaon-to-pion gluon momentum fraction ratio of 1.92(18)_stat.","tokens_in":16944,"tokens_out":9992,"duration_ms":108484,"significance":"If the result withstands scrutiny, the kaon-to-pion ratio ⟨x⟩_g^K/⟨x⟩_g^π≈1.92 would be a notable finding, contradicting the near-equality suggested by ETMC and DSE and with implications for emergent-mass phenomenology. The paper has clear strengths: the statistics are very large for gluonic observables; the smearing study is systematic and informative; the ground-state matrix-element extractions are checked against variations in fit windows; the OPE-based moment is independent of the PDF model; and the error budget is reported transparently. The main caveat is that the central renormalized value and the PDF both rest on an uncomputed quark-gluon mixing term that is replaced by an ad hoc 10% systematic, and the ratio is quoted without propagating that uncertainty.","major_comments":[{"comment":"The headline moment 0.557(18)_stat(24)_NPR(56)_mixing is obtained by setting Z_gq⟨x⟩_q^bare to zero and assigning a 10% systematic based on Refs. [18,20,31]. The paper itself cites those works as reporting mixing between 2% and 20%, so the quoted 56_mixing uncertainty is not an envelope of the cited range but a choice near its lower end. A 20% mixing contribution would shift the central value by roughly ±0.11, nearly twice the quoted mixing error. The authors should either compute Z_gq⟨x⟩_q^bare, or quote the central value with a conservative uncertainty derived from the full cited range, and should state explicitly that the current 10% is an assumption rather than a derived systematic.","section":"III.A, Eq. (14)"},{"comment":"The ratio ⟨x⟩_g^K/⟨x⟩_g^π = 1.92(18)_stat is quoted with only the statistical error. The kaon moment has additional 24_NPR and 56_mixing uncertainties, and the pion moment from Ref. [21] carries its own 13_NPR uncertainty; none of these are propagated into the ratio. Adding just the kaon mixing uncertainty in quadrature changes the ratio uncertainty to about 0.19, and including the NPR terms brings the total to roughly 0.29. The conclusion that the kaon gluon carries significantly more momentum than the pion is therefore not supported by the quoted error. The ratio must be accompanied by a full uncertainty budget, including the quadrature combination or, preferably, an estimate of the correlation between the two moments.","section":"III.A, ratio paragraph"},{"comment":"The pseudo-PDF matching used to extract the kaon gluon PDF drops the quark-gluon kernel R_gq, and the paper justifies this by citing previous nucleon and pion gluon PDF studies [23,24,29]. Those studies do not automatically cover the kaon, whose quark disconnected contributions can differ. Because the same neglected mixing term affects both the renormalized moment and the RpITD entering the fit, the comparisons with DSE for x≳0.4 and the higher-moment ratios are not independent of the mixing assumption. The authors should provide a quantitative estimate of the R_gq contribution to M(ν,z²), or propagate a corresponding systematic into the extracted PDF band.","section":"III.B, Eq. (15)"},{"comment":"The fit excludes Pz=1 because those data cannot be properly described by the matching, yet no systematic uncertainty is assigned to this exclusion. The paper itself notes that small-Pz contamination has become visible only at the present statistics, which suggests that the remaining Pz∈[2,5] data may also carry a residual, smaller version of the same effect. At minimum, the authors should report a fit that includes Pz=1 with a modified treatment, or vary the minimum Pz in the fit and include the spread in the quoted PDF uncertainties.","section":"III.B, fit-range paragraph"}],"minor_comments":[{"comment":"There are typos that should be corrected, including 'in in the MS scheme' in the abstract and 'koan' for 'kaon' in Sec. I.","section":"Abstract and Sec. I"},{"comment":"The text refers to 'WILSON2' although the labels used elsewhere are WILSON1 and WILSON3; also, the Wilson-flow labels should be explicitly defined (flow time in units of a²) since the text states Nsteps=100 for both chosen flow times.","section":"Sec. II, Fig. 5 discussion"},{"comment":"Equation (20) as written displays ⟨x^n⟩_g on the left-hand side, but the expression B(α+n,β+1)/B(α+1,β+1) actually gives the normalized ratio ⟨x^n⟩_g/⟨x⟩_g. The notation should be corrected to avoid implying that the quoted numbers are absolute moments.","section":"III.B, Eq. (20)"},{"comment":"The y-axis label appears as 'M( , z2)' with the Ioffe-time variable ν missing; this should be fixed for clarity.","section":"Fig. 8"},{"comment":"The fitted values of α and β, along with the χ²/dof of the RpITD fits, are not reported. Reporting these values would greatly improve reproducibility and allow readers to judge the influence of the Bayesian priors.","section":"III.B, PDF fit"},{"comment":"The bands in Fig. 10 exclude the uncertainties in ⟨x⟩_g; this is stated in the caption but should be emphasized in the main text, and ideally a version with those uncertainties propagated should be provided so that the visual separation between pion and kaon curves is not overinterpreted.","section":"Fig. 10 and surrounding text"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper is transparent and the bare-moment extraction is careful, but the load-bearing systematic is externally assigned rather than derived. The stress-test concern about the 10% mixing error is valid and is compounded by the ratio being quoted with only statistical error. I recommend major revision rather than rejection because these issues are fixable within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll cut to it: this is a solid, transparent lattice-QCD analysis that delivers a new high-statistics kaon gluon moment and PDF. If the central value 0.557(18)(24)(56) holds, it's a useful data point, but the error budget has a hole: the quark-gluon mixing term is not computed and is replaced by a 10% systematic. That is the thing to look at carefully.\n\nThe paper does quite well on the lattice side. 1.3M measurements, a clean smearing comparison across HYP/Stout/Wilson, and stability checks on the ground-state extraction. The bare moment is fit over several momenta and is robust. The NPR uncertainty is separated from statistical, which is good practice. The comparison of the PDF with DSE and with the pion PDF from the same ensemble is a nice way to frame the result.\n\nThe soft spots are real. Eq. 14 includes the mixing term Z_gq, and they drop it. They cite prior work reporting 2% to 20% mixing, then use 10% as the error. If the actual mixing is closer to 20%, the shift is about ±0.11 on 0.557, which is nearly twice their 0.056 mixing error. The ratio 1.92(18) is quoted with statistical error only, so the mixing uncertainty is not propagated there either; it might partially cancel in the ratio, but they haven't shown that. Also, this is one lattice spacing, no continuum extrapolation, and the pion mass is 310 MeV; the paper acknowledges these limits. The higher moments in Eq. 20 are just outputs of the fitted alpha and beta, not independent predictions.\n\nNone of these kill the paper, but the mixing issue is load-bearing. The central value may be right, but the uncertainty on the dominant neglected term is not derived from this calculation. A referee should ask for either a computation of the disconnected diagrams on this ensemble, or at least a demonstrably conservative prior that covers the cited range.\n\nAt the end of the day, this is the kind of paper I'd send out. It's a careful, honest calculation with a new result, and the critique is addressable. I'd read it again after the mixing treatment improves. Worth a serious referee; conditional acceptance at best.","headline":"A careful, high-statistics kaon gluon moment and PDF, but the headline ratio rests on an uncomputed quark-gluon mixing term.","tokens_in":17514,"tokens_out":3361,"would_cite":true,"duration_ms":35866,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t"],"model":"deepseek-v4-flash","headline":"Kaon gluons carry about twice the pion's momentum fraction: $\\langle x\\rangle_g^{K}=0.557(18)(24)(56)$ at $\\mu=2$ GeV, with ratio $1.92(18)$.","keywords":["kaon","gluon parton distribution function","gluon momentum fraction","lattice QCD","pseudo-PDF method","reduced pseudo-Ioffe-time distribution","nonperturbative renormalization","meson structure"],"falsifier":"Compute the quark disconnected diagrams and the $Z_{gq}\\langle x\\rangle_q$ mixing contribution on the same ensemble and renormalize with the full mixing matrix: if the resulting $\\langle x\\rangle_g^{\\overline{\\text{MS}},K}$ moves by more than 0.056, or if the kaon-to-pion ratio moves outside $1.92(18)$, the central claim is falsified. A physical-continuum extrapolation with full mixing that brings the kaon moment down to the earlier value near 0.42 would also falsify the factor-of-two conclusion.","tokens_in":16425,"feed_emoji":"⚛️","tokens_out":15814,"duration_ms":130993,"temperature":0.7,"pith_summary":"This paper asks how much of the kaon's momentum is carried by gluons, using lattice QCD with 1,296,640 kaon-correlator measurements on one ensemble with lattice spacing about 0.12 fm and a 310 MeV pion mass. It reports a nonperturbatively renormalized gluon momentum fraction $\\langle x\\rangle_g^{\\overline{\\text{MS}},K}=0.557(18)_{\\text{stat}}(24)_{\\text{NPR}}(56)_{\\text{mixing}}$ at $\\mu=2$ GeV in the $\\overline{\\text{MS}}$ scheme. The central result is the kaon-to-pion ratio $\\langle x\\rangle_g^K/\\langle x\\rangle_g^\\pi = 1.92(18)$: on the same ensemble, gluons carry about twice as much momentum in the kaon as in the pion. The paper also extracts the kaon gluon PDF through pseudo-PDF matching and finds its normalized shape agrees with a Dyson-Schwinger calculation for $x\\gtrsim 0.4$. If the ratio holds up, it would mean the gluon momentum share is not the same for kaon and pion as earlier calculations suggested, and future electron-ion collider measurements have a concrete prediction to test.","feed_headline":"Kaon gluons carry nearly twice the pion's share","feed_subtitle":"Kaon gluon momentum fraction comes out 0.557(18)(24)(56) at 2 GeV, about 1.9 times the pion's on the same ensemble.","key_machinery":"The argument is carried by the reduced pseudo–Ioffe-time distribution (RpITD), a double ratio of fitted ground-state gluon matrix elements $\\mathcal{M}(\\nu,z^2)=\\frac{M(z\\cdot P_z,z^2)/M(0\\cdot P_z,0)}{M(z\\cdot0,z^2)/M(0\\cdot0,0)}$ that cancels renormalization factors and ultraviolet divergences. The gluon operator $F^{ti}(z)W(z,0)F^t_i(0)-F^{ij}(z)W(z,0)F_{ij}(0)$, with a gauge-invariant Wilson line $W$, produces the matrix elements; for the first moment a separate OPE gluon operator is used and related to $\\langle x\\rangle_g$ by a kinematic factor. The RpITD is matched to the light-cone gluon PDF through a convolution with the gluon-gluon kernel $R_{gg}$, using the model $x^\\alpha(1-x)^\\beta/B(\\alpha+1,\\beta+1)$ with Bayesian priors. A smearing study over hypercubic, Stout, and Wilson-flow steps leads to the choice of five steps of HYP smearing as the conservative balance between noise reduction and preservation of the physics. Quark-gluon mixing in the renormalization is not computed; it is replaced by a 10 percent systematic error.","core_discovery":"The paper's central claim is that, on a 2+1+1-flavor HISQ ensemble with $a\\approx0.12$ fm and $M_\\pi\\approx310$ MeV, the nonperturbatively renormalized kaon gluon momentum fraction is $\\langle x\\rangle_g^{\\overline{\\text{MS}},K}=0.557(18)_{\\text{stat}}(24)_{\\text{NPR}}(56)_{\\text{mixing}}$ at $\\mu=2$ GeV. Taking the pion's value from the same ensemble, the paper finds $\\langle x\\rangle_g^K/\\langle x\\rangle_g^\\pi=1.92(18)$, which is in tension with earlier lattice and Dyson-Schwinger results that put the two moments nearly equal. Using the pseudo-PDF method, the paper extracts the kaon gluon PDF and reports that its shape, divided by the first moment, agrees with the Dyson-Schwinger prediction for $x\\gtrsim0.4$ and with the pion's shape within one standard deviation; the absolute distributions differ because of the factor-of-two momentum fraction. Higher normalized moments are $\\langle x^2\\rangle_g^K/\\langle x\\rangle_g^K=0.123(16)$ and $\\langle x^3\\rangle_g^K/\\langle x\\rangle_g^K=0.0277(46)$.","pith_inferences":["One consequence the paper does not pursue: if the factor-of-two ratio survives the physical-continuum limit, then SU(3) flavor breaking in the gluonic sector is much larger than in the quark sector, which would be a new constraint on models of emergent hadronic mass.","A direct test is to compute the quark disconnected diagrams and the full quark-gluon mixing term on the same ensemble; if the correction exceeds the assumed 10 percent, both the central moment and the kaon-to-pion ratio would move, possibly reconciling the result with earlier near-equality predictions.","The observation that the normalized kaon and pion gluon PDF shapes agree within 1$\\sigma$, while the absolute fractions differ by a factor of two, suggests the $x$-dependence of gluon distributions in light pseudoscalar mesons may be approximately universal even when the total gluon momentum share is not; a physical-pion-mass comparison would test this.","The small-$P_z$ contamination identified in the RpITD data implies that previous lower-statistics pseudo-PDF results may have underestimated similar systematics; re-analyzing older ensembles with the same fit-window and momentum cuts would quantify that."],"forward_implications":["If the kaon-to-pion ratio $1.92(18)$ is confirmed, the kaon's gluon momentum fraction is about twice the pion's, so by the momentum sum rule the quark momentum fraction in the kaon is correspondingly smaller.","The extracted gluon PDF, matched at $\\mu=2$ GeV, gives a lattice prediction for kaon gluon distributions that future electron-ion collider measurements can test directly for $x\\gtrsim0.4$.","The high-statistics data reveal low-momentum contamination in the RpITD that was not visible at lower statistics, so future pseudo-PDF analyses should exclude or model the smallest $P_z$ points.","The result provides a single-ensemble data point toward a physical-continuum extrapolation of the kaon gluon moment; combining it with other ensembles would sharpen the comparison with earlier continuum-extrapolated lattice results.","The smearing study identifies five steps of HYP smearing as a conservative default for gluonic matrix elements, which can inform the design of future gluon-PDF calculations."],"supporting_citations":[{"why":"It defines the gluon pseudo-distribution operator and the pseudo-PDF framework used for the extraction.","marker":"[5]"},{"why":"It is the earlier lattice determination of pion and kaon gluon momentum fractions with continuum extrapolation that this work compares against.","marker":"[13]"},{"why":"It is the previous kaon gluon PDF calculation that this work extends with higher statistics and a moment determination.","marker":"[14]"},{"why":"It reports quark-gluon mixing magnitudes used to set the 10 percent mixing systematic.","marker":"[18]"},{"why":"It supplies the nonperturbative renormalization constant $Z_{gg}^{-1}=1.512(65)$ used for the bare moment.","marker":"[19]"},{"why":"It is the clover-fermion study finding mixing no larger than 10 percent, the basis for the quoted mixing uncertainty.","marker":"[20]"},{"why":"It gives the pion gluon moment and PDF on the same ensemble, used for the kaon-to-pion ratio and the direct PDF comparison.","marker":"[21]"},{"why":"It is one of the previous lattice studies reporting the quark-gluon mixing magnitude used to bound the systematic error.","marker":"[31]"},{"why":"It provides the Dyson-Schwinger kaon and pion gluon PDFs compared against the extracted PDF and used for higher moments.","marker":"[32]"},{"why":"It provides the gluon matching relation connecting the RpITD to the light-cone gluon PDF.","marker":"[33]"}],"fun_headline_variants":["Kaon gluons double pion's share in lattice QCD","Kaon gluon momentum fraction nearly twice pion's","Lattice QCD: kaon gluons carry 1.9x pion's share","High-stat kaon lattice QCD finds gluon boost over pion","Kaon's gluon momentum beats pion's by nearly factor 2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that neglecting the quark-gluon mixing term $Z_{gq}$ in the renormalization and absorbing it into a 10 percent systematic error is adequate; if the true mixing for the kaon exceeds 10 percent, the central value 0.557 and the kaon-to-pion ratio shift by more than the quoted mixing uncertainty.","fun_headline_variants_meta":{"raw":{"variants":["Kaon gluons double pion's share in lattice QCD","Kaon gluon momentum fraction nearly twice pion's","Lattice QCD: kaon gluons carry 1.9x pion's share","High-stat kaon lattice QCD finds gluon boost over pion","Kaon's gluon momentum beats pion's by nearly factor 2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000273,"raw_usage":{"total_tokens":1697,"prompt_tokens":1068,"completion_tokens":629,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":684,"completion_tokens_details":{"reasoning_tokens":550}},"tokens_in":684,"tokens_out":629,"duration_ms":5569,"temperature":1.0,"reasoning_tokens":550,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:11:21.516443+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the quark disconnected diagrams and the $Z_{gq}\\langle x\\rangle_q$ mixing contribution on the same ensemble and renormalize with the full mixing matrix: if the resulting $\\langle x\\rangle_g^{\\overline{\\text{MS}},K}$ moves by more than 0.056, or if the kaon-to-pion ratio moves outside $1.92(18)$, the central claim is falsified. A physical-continuum extrapolation with full mixing that brings the kaon moment down to the earlier value near 0.42 would also falsify the factor-of-two conclusion.","supporting_citations":[{"cited_title":"Gluon Pseudo-Distributions at Short Distances: Forward Case","cited_arxiv_id":null,"evidence_quote":"It defines the gluon pseudo-distribution operator and the pseudo-PDF framework used for the extraction."},{"cited_title":"Quark and Gluon Momen- tum Fractions in the Pion and in the Kaon","cited_arxiv_id":null,"evidence_quote":"It is the earlier lattice determination of pion and kaon gluon momentum fractions with continuum extrapolation that this work compares against."},{"cited_title":"Quark and Gluon Momen- tum Fractions in the Pion from Nf=2+1+1 Lattice QCD","cited_arxiv_id":null,"evidence_quote":"It reports quark-gluon mixing magnitudes used to set the 10 percent mixing systematic."},{"cited_title":"Nonperturbatively renormalized nucleon gluon momen- tum fraction in the continuum limit of Nf=2+1+1 lattice QCD","cited_arxiv_id":null,"evidence_quote":"It supplies the nonperturbative renormalization constant $Z_{gg}^{-1}=1.512(65)$ used for the bare moment."},{"cited_title":"Hackett, Patrick R","cited_arxiv_id":null,"evidence_quote":"It is the clover-fermion study finding mixing no larger than 10 percent, the basis for the quoted mixing uncertainty."},{"cited_title":"Gluon moment and parton distribution func- tion of the pion from Nf=2+1+1 lattice QCD","cited_arxiv_id":null,"evidence_quote":"It gives the pion gluon moment and PDF on the same ensemble, used for the kaon-to-pion ratio and the direct PDF comparison."},{"cited_title":"Gluon momentum frac- tion of the nucleon from lattice QCD","cited_arxiv_id":null,"evidence_quote":"It is one of the previous lattice studies reporting the quark-gluon mixing magnitude used to bound the systematic error."},{"cited_title":"Kaon and pion parton distributions","cited_arxiv_id":null,"evidence_quote":"It provides the Dyson-Schwinger kaon and pion gluon PDFs compared against the extracted PDF and used for higher moments."},{"cited_title":"Investigating Gluonic Operators in Coordinate Space","cited_arxiv_id":null,"evidence_quote":"It provides the gluon matching relation connecting the RpITD to the light-cone gluon PDF."}],"review_version":1}