{"id":"26d6b1ad-9f6a-4b22-b5b7-226ab80955e4","arxiv_id":"2608.12467","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A new dimension-seven operator is generated by the gluon equation of motion in the presence of the axion-gluon coupling, and its spectator-scattering contribution to B to K a is roughly 25 percent of the soft-overlap contribution.","lead":"This paper builds a specialized effective-theory framework for producing a heavy axion in rare B-meson decays, focusing on B to K plus axion. It shows that a previously missed hard-scattering contribution is about a quarter of the usual soft contribution when the axion couples only to gluons at the bottom-quark scale.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 25–30% spectator/soft ratio inherits an uncontrolled E0(q²) error: the soft matching uses the constant E0≈5 at hard virtuality ~m_b², where the paper's own footnote says the local approximation fails, while the spectator term is evaluated at E0(0).","rationale":"The reader identified the E0 locality assumption as the weakest point, and I agree. This is genuinely load-bearing because the central claim is the ratio in Eq. (67), not merely the existence of the new operator. The paper deserves credit for a clear footnote acknowledging the limitation and for constructing the SCET factorization in a parameter-light way; the derivation of O_∂ag and the spectator factorization are internally consistent. But the numerical hierarchy between soft and spectator is exactly where the approximation bites. The soft amplitude is a one-loop matching calculation at the hard scale, where E0 is not a constant, and this cannot be treated as a small refinement: no parametric suppression prevents E0(m_b²) from differing substantially from E0(0). Since the ratio is a ratio of two terms using the same (incorrect) constant, the error does not cancel. The suggested numerical test is straightforward and should settle whether the 25–30% claim survives. If it does, the CONDITIONAL verdict remains; if it does not, the paper's central quantitative message needs revision even though the structural SCET construction may still be valid.","tokens_in":25266,"tokens_out":12403,"duration_ms":116971,"concrete_test":"Recompute Eq. (45) with the momentum-dependent b→s g vertex: replace the constant C_Dg inside the O_∂ag insertion by the full charm-loop form factor C_Dg((p−ℓ)²) (e.g., from Ref. [35] or a direct one-loop evaluation with m_c ≠ 0), then perform the ℓ integral and the SCET matching to obtain M_soft, while keeping Eq. (44) unchanged because the hard-collinear spectator gluon has ℓ² ≪ 4m_c². Re-evaluate Eq. (67) with this corrected M_soft and compare with 0.25–0.30. A simpler diagnostic is to compute the ratio E0(m_b²)/E0(0) from the known Inami-Lim or charm-loop functions; if this ratio deviates from 1 by more than about 30%, the quoted R_spec interval is not settled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim is Eq. (67), R_spec ≈ 0.25–0.30. Both M_spec and M_soft are proportional to C_∂ag = 4 C_gg C_Dg, and C_Dg contains E0. The spectator diagram (Sec. III A) has an internal n-hard-collinear gluon with ℓ² ~ -2ω ū E_K ~ -Λ_QCD m_b, far below 4m_c², so using the q²→0 local value |E0|≈5 is justified. The soft-overlap amplitude (Sec. III B, Eqs. (45)–(49)) instead integrates a hard gluon with virtuality up to m_b². There the constant-E0 treatment is not valid; E0 is a q²-dependent charm-loop form factor and should be evaluated or integrated at q² ~ m_b². Because the same |E0|=5 is used in both amplitudes, R_spec is effectively E0(0)/E0(m_b²) times the true ratio. If the charm-loop form factor at hard virtuality differs from its low-energy value by an O(1) factor, the headline 25–30% range is not robust. The effect can go in either direction: if E0(m_b²) < E0(0), the spectator term becomes relatively larger or even dominant; if E0(m_b²) > E0(0), it drops toward the 6–7% level of Case 2. The paper's own footnote (page 3) flags exactly this: the local approximation holds only for gluon virtualities well below 4m_c² and 'for hard gluons it is only an approximation.'","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs a soft-collinear effective theory (SCET) description of B to K a decays for a heavy QCD axion coupled only through a G tilde-G at low scales. The authors show that eliminating the redundant operator O_Dg via the gluon equation of motion, modified by the axion-gluon interaction, generates a dimension-seven operator O_{partial a g}. They match this operator onto SCET and calculate the leading-power spectator-scattering amplitude (Eq. (44)) and the one-loop soft-overlap amplitude (Eq. (51)). The central numerical result is that the spectator term amounts to about 25-30 percent of the soft term (Eq. (67)) when aG tilde-G is present only below the electroweak scale. In the alternative UV realization where aG tilde-G is present above the electroweak scale, the induced direct operator O_bsa dominates and the gluonic spectator term is only about 6-7 percent (Eq. (73)). The paper also derives f_a sensitivity estimates for both cases and discusses ALP-meson mixing and RG effects in the appendices.","tokens_in":25692,"tokens_out":21828,"duration_ms":200781,"significance":"If the central result holds, the paper identifies a genuinely new leading-power spectator-scattering contribution to B to K a that is absent in the conventional form-factor description and is large enough to matter for precision studies of the low-scale gluonic axion scenario. The operator derivation is transparent, the SCET matching is presented in detail, and the factorization formula in Eq. (44) is explicit and reproducible from the stated inputs. The paper also gives a clean separation of two UV realizations and makes falsifiable numerical predictions. The main quantitative claim, however, depends on treating the charm-loop coefficient E0 as a local constant in the soft-overlap amplitude, an approximation that the paper itself flags as invalid for hard gluons; this is the main obstacle to accepting the 25-30 percent ratio.","major_comments":[{"comment":"The central numerical claim in Eq. (67) is not robust under the momentum dependence of E0, which the footnote after Eq. (6) identifies as a q^2-dependent charm-loop form factor. In the spectator amplitude (Eq. (44)) the internal hard-collinear gluon has virtuality l^2 = -2 omega u-bar E_K ~ -Lambda_QCD m_b, below 4 m_c^2, so using the low-energy value |E0| = 5 is justified there. In the soft-overlap matching (Eqs. (45)-(49)) the loop integral samples hard gluon virtualities up to O(m_b^2), exactly the region where the paper states that the local approximation fails. Because both M_spec and M_soft are proportional to C_{partial a g} = 4 C_gg C_Dg, the ratio R_spec in Eq. (67) implicitly contains E0(q_low^2)/E0(q_hard^2), and an O(1) variation of this form factor at q^2 ~ m_b^2 would move R_spec well outside the quoted 0.25-0.30 range. Please recompute the soft-overlap amplitude with the actual q^2-dependent E0 (for instance, using the charm-loop dispersion treatment of Ref. [35]) or provide a quantitative bound on the induced error.","section":"Sec. III B, Eq. (67)"},{"comment":"The discussion of renormalization-group effects in Appendix C is qualitative: the statement that resummation changes the individual amplitudes by 10-20 percent and partially cancels in R_spec is not supported by a calculation. The quoted range for R_spec in Eq. (67) is itself of order 25 percent, so an unquantified 10-20 percent effect is comparable in size to the central claim. The paper should either present a resummed result for the ratio, or explicitly label Eq. (67) as a fixed-order leading-order estimate and assign a conservative uncertainty from the missing resummation.","section":"Appendix C"}],"minor_comments":[{"comment":"The heading 'F uture Directions' contains a spacing typo and should read 'Future Directions'.","section":"Sec. VI"},{"comment":"The relative sign between M_spec and M_soft is not computed, so the enhancement factor in Eq. (69) and the statement that the spectator contribution strengthens the f_a limit by 25-30 percent assume constructive interference. This assumption should be stated in the main text together with the branching-fraction predictions, since a negative relative sign would instead reduce the total rate.","section":"Sec. VI, Eqs. (67)-(69)"},{"comment":"The text says x = 2E_K/m_B with x <= 1, but for small axion masses x is very close to 1 (e.g., x = 0.9996 for m_a = 0.5 GeV). The wording 'x < 1' should be clarified so that the large-recoil window and the x to 1 limit are not confused.","section":"Eq. (26)"},{"comment":"The derivation leading to Eq. (44) drops the term 1/(2 u-bar E_K) relative to 1/omega in Eq. (42). This is justified at leading power but only after an endpoint-region check; a brief statement that endpoint contributions are power-suppressed for LCDAs vanishing linearly at the endpoint would make the step self-contained.","section":"Eq. (44)"}],"recommendation":"major_revision","confidential_remarks":"This is a serious SCET calculation with a clear and interesting central claim. The main obstacle is the E0(q^2) issue, which is fixable within the scope of the paper by evaluating or bounding the charm-loop form factor at hard virtuality. The self-citation to Ref. [6] is used appropriately for the Case-2 comparison, and I see no novelty or scope problems."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth reading. It identifies a new dimension-seven operator O_∂ag, generated when the redundant flavor-changing gluon operator O_Dg is eliminated using the axion-modified gluon equation of motion, and matches it onto SCET for B→K a. That is a real extension: earlier heavy-QCD-axion analyses in rare B decays used the total B-to-K form factor and the electroweak-generated b→s a operator, missing this gluonic spectator mechanism. The SCET work is careful—spinor traces, LCDA convolutions, color factors all check out, and the factorization of the spectator term is explicit. The second scenario (gluon coupling above the electroweak scale) sensibly reproduces the known dominance of the form-factor treatment and yields a spectator contribution of only 6–7%. The paper is also honestly written: it flags its own approximations, including the one that matters.\n\nThe soft spot is exactly where the stress-test lands. The central ratio in Eq. (67), R_spec ≈ 0.25–0.30, uses |E0|=5 in both the spectator and soft-overlap amplitudes. The spectator gluon is hard-collinear, with ℓ² ~ −Λ_QCD m_b, so the local value E0(0) is fine there. But the soft-overlap matching in Sec. III B integrates a hard gluon with virtuality up to m_b², and the footnote to Eq. (6) concedes that the local-E0 treatment is only an approximation for hard gluons. Since both amplitudes are proportional to C_∂ag ∝ E0, the ratio is effectively E0(0)/E0(m_b²) times the true value. If the charm-loop form factor at q² ~ m_b² differs from its low-energy value by an O(1) factor—not unlikely, since 4m_c² ≈ 3.2 GeV² is far below m_b²—the headline 25–30% could drop toward the 6–7% level or grow so that the spectator dominates. The direction is genuinely unknown without computing the charm loop at hard virtuality.\n\nMinor issues: the RG resummation is not performed and is estimated at 10–20%, and the ALP-meson mixing interference discussed in Appendix B can be as large as ~50% of the direct term in Case 1, but is set aside. Both are acknowledged in the text.\n\nWho is this for? People working on ALP phenomenology in rare B decays and on SCET factorization. The structural result—the existence of a leading-power gluonic spectator contribution generated by O_∂ag—is robust and will need to be incorporated in future precision analyses. The numerical 25–30% estimate, however, should be treated as preliminary until the E0(q²) dependence is handled properly.\n\nRecommendation: send it to referees, with instructions to check whether the soft-overlap matching can be reformulated with a q²-dependent E0 and whether the hard-virtual charm-loop form factor changes the ratio. The paper deserves referee time; the headline number does not yet deserve to be quoted.","headline":"Genuinely new SCET operator and spectator contribution for heavy QCD axions in B→K a, but the headline 25–30% ratio rests on a local-E0 approximation the paper itself flags as doubtful at the relevant hard virtuality; the structural claim holds, the number does not yet.","tokens_in":26221,"tokens_out":2311,"would_cite":true,"duration_ms":21249,"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":"A new gluonic operator adds a 25–30% spectator correction to B→K a in the low-scale axion scenario.","keywords":["heavy QCD axion","axion-like particles","soft-collinear effective theory","B to K a decay","spectator scattering","soft-overlap form factor","gluon equation of motion","flavor-changing neutral currents"],"falsifier":"Recompute the soft-overlap matching of $O_{\\partial ag}$ using the momentum-dependent charm-loop form factor $E_0(q^2)$ evaluated at $q^2\\sim m_b^2$ instead of the constant $|E_0|=5$; if the resulting soft amplitude changes by more than the quoted uncertainties, the central $25$–$30\\%$ spectator-to-soft ratio and the rate enhancement $1.5$–$1.7$ shift accordingly. Alternatively, measure the $m_a$ dependence of the $B\\to K a$ rate: the low-scale gluonic scenario predicts a total SCET amplitude that decreases with $m_a$, whereas the conventional form-factor prediction increases.","tokens_in":25092,"feed_emoji":"⚛️","tokens_out":13651,"duration_ms":111689,"temperature":0.7,"pith_summary":"This paper develops a soft-collinear effective theory (SCET) treatment of rare $B\\to K a$ decays induced by a heavy QCD axion coupled to gluons through $aG\\tilde G$. Its central claim is that when that gluonic coupling is the only axion interaction at the weak scale, eliminating the redundant flavor-changing gluon operator in the low-energy theory generates a new dimension-seven operator, and this operator produces a hard-spectator-scattering contribution to $B\\to K a$ that sits at the same power as the soft-overlap term and equals roughly 25–30% of it. In the alternative, more standard realization in which the gluon coupling is imposed above the electroweak scale, renormalization-group evolution generates a direct $b\\to s a$ operator whose soft form-factor contribution dominates and reduces the spectator piece to 6–7%. The paper thereby identifies when the usual form-factor description of $B\\to K a$ is reliable and when it must be replaced by a factorized SCET treatment that includes spectator scattering.","feed_headline":"Rare B-meson axion decay gets a 25% spectator boost","feed_subtitle":"A new gluonic operator makes hard spectator scattering a 25% correction in the low-scale gluonic axion scenario.","key_machinery":"The load-bearing mechanism is the modified gluon equation of motion in the presence of $aG\\tilde G$: $(D^\\mu G_{\\alpha\\mu})^A = g_s J_\\alpha^A + (4C_{gg}/f_a)(\\partial^\\mu a)\\tilde G_{\\alpha\\mu}^A$. Substituting this into the redundant operator $O_{Dg}$ is what generates $O_{\\partial ag}$ and, after tree-level matching onto SCET (the effective theory that separates soft, collinear, and hard-collinear momentum modes in heavy-to-light decays), the leading SCET operator $O^{\\rm SCET}_{\\partial ag}\\propto [\\bar\\chi_n\\gamma^\\mu_\\perp P_L T^A Y_s^\\dagger h_v](\\partial_\\nu a)\\epsilon^{\\mu\\nu\\rho\\sigma} n_\\rho G^A_{\\bar n\\perp\\sigma}$. The spectator amplitude then factorizes as $M^{\\partial ag}_{\\rm spec}= g_s (C_{\\partial ag}/f_a)(C_F E_K/N_c)(f_K f_B m_B^2)/(8\\lambda_B)$, with the inverse moment $1/\\lambda_B$ encoding the $B$-meson light-cone distribution amplitude, while the soft amplitude is the matrix element of the radiatively generated SCET current proportional to $(n\\cdot q)E_K\\zeta(E_K)$. The same machinery, with different Wilson coefficients, controls the chromomagnetic operator $O_{8g}$ and the direct $b\\to s a$ operator.","core_discovery":"The central discovery is that the axion–gluon operator $aG\\tilde G$ modifies the gluon equation of motion, so the standard elimination of the redundant operator $O_{Dg}$ from the weak effective theory no longer produces only the familiar QCD-penguin four-quark operators. It also produces a new dimension-seven operator $O_{\\partial ag} = \\bar{s} T^A \\gamma^\\mu P_L b\\,(\\partial_\\nu a)\\tilde G^{A\\mu\\nu}$ with Wilson coefficient $C_{\\partial ag}/f_a = 4 C_{gg} C_{Dg}/f_a$. Matching $O_{\\partial ag}$ onto SCET at leading power yields a soft-overlap amplitude proportional to $(n\\cdot q) E_K \\zeta(E_K)$ and a spectator-scattering amplitude that factorizes into a perturbative hard kernel and the $B$- and $K$-meson light-cone distribution amplitudes. In the low-scale gluonic realization the spectator amplitude is $25$–$30\\%$ of the soft one (Eq. 67), so under constructive interference the total rate is enhanced by a factor $1.5$–$1.7$ and the conventional form-factor treatment is incomplete; in the above-electroweak realization the direct $b\\to s a$ operator makes the soft term dominant and the spectator correction only $6$–$7\\%$ (Eq. 73), so the conventional treatment remains valid. The same framework yields the chromomagnetic-operator contributions and updated bounds on the axion decay constant $f_a$ for both realizations.","pith_inferences":["A testable extension is to compute the charm-loop form factor $E_0(q^2)$ at hard virtuality $q^2\\sim m_b^2$ and redo the one-loop soft matching; if it departs from $|E_0|=5$, the 25–30% ratio is the first quantity to shift.","The same modified-gluon-equation mechanism should generate analogous spectator operators in other $b\\to s$ modes with pseudoscalar final states, such as $B_s\\to\\phi a$ or $\\Lambda_b\\to\\Lambda a$, where the soft-versus-spectator balance could be probed in differential distributions.","Because the SCET total amplitude in the low-scale scenario decreases with $m_a$ while the LCSR form-factor result increases, a shape measurement of the $B\\to K a$ rate as a function of $m_a$ could distinguish the spectator mechanism from the conventional treatment."],"forward_implications":["In the low-scale gluonic scenario, precision analyses of $B\\to K a$ must include the $O_{\\partial ag}$ spectator term: omitting it changes the amplitude by roughly 25–30% and strengthens the inferred bound on $f_a$ by the same amount under constructive interference.","When the gluon coupling is imposed above the electroweak scale, the RG-induced direct $b\\to s a$ operator makes the soft-overlap form factor dominant and the spectator correction only 6–7%, so the conventional $f_0(q^2)$ form-factor treatment remains valid within its uncertainties.","The chromomagnetic operator $O_{8g}$ has a larger spectator-to-soft ratio (about 0.5–0.9), but its small Wilson coefficient makes its absolute contribution only about 15% of the $O_{\\partial ag}$ spectator term.","The two realizations have very different experimental reach: the low-scale gluonic scenario probes only $f_a\\sim O(\\text{GeV})$, while the above-electroweak scenario with the direct $b\\to s a$ operator gives branching fractions of order $10^{-6}$ for $f_a\\sim 200$ GeV and correspondingly stronger limits."],"supporting_citations":[{"why":"It supplies the equation-of-motion elimination of the redundant operator $O_{Dg}$ that produces the new operator $O_{\\partial ag}$.","marker":"[30]"},{"why":"It defines the redundant operator $O_{Dg}$ and the chromomagnetic operator $O_{8g}$ whose coefficients drive the leading contributions.","marker":"[29]"},{"why":"It provides the soft-overlap form factor $\\zeta(E_K)$ and its $m_a$ dependence used in the numerical ratios.","marker":"[12]"},{"why":"It establishes the large-recoil decomposition of the $B\\to K$ form factor into soft overlap plus spectator correction that the paper refines.","marker":"[9]"},{"why":"It supplies the SCET factorization of heavy-to-light form factors and the power counting that places soft and spectator contributions at the same order.","marker":"[10]"},{"why":"It provides the matching coefficients for the $B$-type SCET currents used in the spectator amplitude of the direct $b\\to s a$ operator.","marker":"[24]"},{"why":"It defines the heavy QCD axion scenario in $b\\to s$ transitions and supplies the two-loop $b\\to s a$ Wilson coefficient for the above-electroweak realization.","marker":"[3]"},{"why":"It supplies the conventional light-cone-sum-rule form-factor treatment and the branching-fraction estimates against which the SCET results are compared.","marker":"[6]"}],"fun_headline_variants":["Axion B decay gets 25% spectator boost from new gluonic operator","New gluonic operator shifts B→K axion rate by 25%","SCET uncovers 25% spectator term in heavy-axion B decay","Low-scale gluonic axion: B decay boosted 25% by spectator"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the coefficient $E_0\\simeq5$ of the redundant operator $O_{Dg}$ can be treated as a local constant even though it is really a $q^2$-dependent charm-quark loop form factor; the approximation is trusted below the charm threshold, while the soft-overlap amplitude receives hard gluons at the $b$-quark scale.","fun_headline_variants_meta":{"raw":{"variants":["Axion B decay gets 25% spectator boost from new gluonic operator","New gluonic operator shifts B→K axion rate by 25%","SCET uncovers 25% spectator term in heavy-axion B decay","Low-scale gluonic axion: B decay boosted 25% by spectator"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000265,"raw_usage":{"total_tokens":1723,"prompt_tokens":1174,"completion_tokens":549,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":790,"completion_tokens_details":{"reasoning_tokens":465}},"tokens_in":790,"tokens_out":549,"duration_ms":5064,"temperature":1.0,"reasoning_tokens":465,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:10:37.919701+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the soft-overlap matching of $O_{\\partial ag}$ using the momentum-dependent charm-loop form factor $E_0(q^2)$ evaluated at $q^2\\sim m_b^2$ instead of the constant $|E_0|=5$; if the resulting soft amplitude changes by more than the quoted uncertainties, the central $25$–$30\\%$ spectator-to-soft ratio and the rate enhancement $1.5$–$1.7$ shift accordingly. Alternatively, measure the $m_a$ dependence of the $B\\to K a$ rate: the low-scale gluonic scenario predicts a total SCET amplitude that decreases with $m_a$, whereas the conventional form-factor prediction increases.","supporting_citations":[{"cited_title":"Beneke and D","cited_arxiv_id":null,"evidence_quote":"It supplies the equation-of-motion elimination of the redundant operator $O_{Dg}$ that produces the new operator $O_{\\partial ag}$."},{"cited_title":"(27) are absorbed into the kaon projector, 6 while the soft spectator quark and the heavy-quark field combine to form theB-meson projector","cited_arxiv_id":null,"evidence_quote":"It defines the heavy QCD axion scenario in $b\\to s$ transitions and supplies the two-loop $b\\to s a$ Wilson coefficient for the above-electroweak realization."},{"cited_title":"ALPs Across Scales","cited_arxiv_id":null,"evidence_quote":"It supplies the conventional light-cone-sum-rule form-factor treatment and the branching-fraction estimates against which the SCET results are compared."}],"review_version":1}