{"id":"efd4a3cc-f970-4d9a-90b2-cd302498d351","arxiv_id":"2608.11814","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The authors derive relative branching ratios for two- and three-gluon glueball decays with Fierz rearrangement and a QCD sum rule check, favoring f0(1710) and eta(2370) as glueball candidates and identifying K*(892) anti-K*(892), pi pi omega, and K Kbar phi as promising search channels.","lead":"A glueball is a hypothetical particle made only of gluons, and this paper calculates how glueballs would decay by treating each gluon as turning into a quark-antiquark pair. The predicted channels, such as kaon-star pairs and pi-pi-omega, give experimenters specific places to look for glueballs.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claims rest on the undetermined vector-only gluon-to-quark vertex in Eq. (16); one must show the branching-ratio rankings survive replacing that vertex with other allowed Lorentz structures.","rationale":"The concern is load-bearing because it targets the origin of the predicted branching-ratio patterns, not just their normalization. The authors are transparent about the effective nature of the framework and about omitted dynamical effects; I am not accusing them of overclaiming beyond their caveats. However, the abstract and conclusion present the f0(1710), eta(2370), and tensor-glueball statements as physical findings, and those findings are only as secure as Eq. (16). The proposed computation is a direct stress test: it asks whether the same algebraic machinery with a different, equally plausible Lorentz coupling preserves the rankings. If the rankings are stable, the reader's and my concern is resolved; if not, the paper should be revised to present the results as illustration of one ansatz rather than as guidance for experiment. The reader's conditional verdict already captures this, so no change to the verdict is needed.","tokens_in":20858,"tokens_out":8371,"duration_ms":88296,"concrete_test":"Recompute the relative branching ratios in Tables II and III with Eq. (16) replaced by a one-parameter family Gamma^mu = gamma^mu + kappa gamma^mu gamma^5 (and, in a second pass, by each independent Dirac basis element gamma^mu, gamma^mu gamma^5, sigma^{munu}, 1, gamma^5 separately), keeping all masses and Table I couplings fixed. Check whether the channel rankings that support the central claims—eta eta / pi pi approximately 0.42 for 0++, K*(892)Kbar*(892)/pi pi approximately 116 and VV-over-PP dominance for 2++, omega phi and omega omega dominance for 0−+, and pi pi omega as the leading 0++ three-gluon mode—are stable for kappa of order unity and across the basis. If any ranking flips, the conclusions are contingent on the unverified vector ansatz. As a secondary check, re-derive Eq.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's main physical identifications—f0(1710) with a sizable gluon component, eta(2370) as a 0−+ glueball, and K*(892)Kbar*(892) as the preferred tensor-glueball channel—all follow from ratios computed with the effective vertex A_i^mu -> lambda^i_ab qbar_a gamma^mu q_b (Eq. 16). That vertex is a constituent-model ansatz, not a QCD-derived current. Every Fierz amplitude in Secs. IV and V is a rearrangement of products of this vector coupling, so the relative weights of the scalar, pseudoscalar, vector, and tensor meson operators in Eqs. (20) and (A1) are fixed by the choice of gamma^mu. An axial-vector, tensor, or mixed Lorentz coupling would redistribute these weights and thus alter every entry in Tables II and III; the paper does not show that the quoted hierarchies (e.g., K*Kbar* dominance for 2++, omega-phi and omega-omega dominance for 0−+, the near-vanishing of the small three-body modes) are stable under such changes. The text explicitly acknowledges in Sec. III and Appendix A that the overall normalization and dynamical overlap/FSI corrections are undetermined, but that caveat does not address the Lorentz-structure dependence, which is a condition on the predictions rather than a common scale. The QCD sum rule check in Sec. VI covers only the scalar channel, is leading order, and relies on the OPE result Eq. (29) without derivation, so it cannot independently secure the pseudoscalar, tensor, or three-gluon claims. The same-flavor pairing rule stated in Sec. III is another unquantified bookkeeping assumption, and the absence of displayed amplitudes for many Table II/III channels (explicitly omitted for three-gluon states in Sec. V) prevents a reader from testing these dependencies from the paper alone.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends a previously proposed Fierz-rearrangement framework to compute relative branching ratios for two-gluon glueballs with J^PC = 0++, 0-+, 2++, and 2-+, and for three-gluon glueballs with 0++ and 1+-. The gluon fields are converted into quark-antiquark currents through the effective vertex A_i^mu -> lambda^i_ab qbar_a gamma^mu q_b, followed by color and Lorentz Fierz rearrangements; the resulting operators are matched to meson states using the couplings in Table I. The paper also performs a leading-order QCD sum rule analysis of the scalar glueball two-pion decay and compares the resulting Gamma(eta eta)/Gamma(pi pi) ratio with the Fierz prediction. Based on the decay patterns, the authors argue for a sizable gluon component in f0(1710), a 0-+ glueball interpretation of eta(2370), and K*(892) Kbar*(892) as a favorable search channel for the tensor glueball, and they identify pi pi omega and K Kbar phi as promising three-body modes.","tokens_in":21254,"tokens_out":3325,"duration_ms":36030,"significance":"If the effective-vertex assumption is accepted, the paper provides a systematic and transparent classification of the flavor and Lorentz structures of two- and three-gluon glueball decays, with explicit relative branching ratios that could guide experimental searches. The scalar-channel QCD sum rule check is a useful independent consistency test, and the algebraic Fierz decompositions in Appendix A are presented in enough detail to be checked. The significance of the main phenomenological claims is, however, conditional on the undetermined Lorentz structure of the gluon-to-quark conversion: the claimed hierarchies, such as K*Kbar* dominance for the tensor glueball and omega-phi/omega-omega dominance for the pseudoscalar glueball, are inherited from the vector-only vertex in Eq. (16), and the paper does not demonstrate their stability under alternative Lorentz couplings.","major_comments":[{"comment":"The central predictions all rest on the ad hoc effective vertex A_i^mu -> lambda^i_ab qbar_a gamma^mu q_b. Every Fierz amplitude in Secs. IV and V is obtained by rearranging products of this vector coupling, so the relative weights of the scalar, pseudoscalar, vector, and tensor meson operators in Eqs. (20) and (A1) are fixed by the choice of gamma^mu. An axial-vector, tensor, or mixed Lorentz structure would redistribute these weights and would change essentially every entry in Tables II and III. The manuscript does not show that the quoted hierarchies, including K*Kbar* dominance for 2++, omega-phi and omega-omega dominance for 0-+, and the suppression of the small three-body modes, survive such a replacement. Because the abstract and Sec. VII present these hierarchies as the main results, this is a load-bearing model assumption rather than an overall normalization that cancels in ratios.","section":"Sec. III, Eq. (16)"},{"comment":"The three-gluon section explicitly states that the explicit decay amplitudes and their squared forms are not listed, and only the final relative branching ratios in Table III are given. This means the 37 branching ratios for the 0++ three-gluon glueball and the 29 ratios for the 1+- state cannot be independently checked from the text, and no ancillary or machine-readable derivation is provided. Since the three-gluon predictions are part of the paper's claimed results, the omission blocks verification and should be remedied by including the amplitudes (or a reproducible derivation) or by clearly marking these results as provisional.","section":"Sec. V, three-gluon amplitudes"},{"comment":"The QCD sum rule analysis is limited to the scalar glueball, and the central OPE expression T_OPE = <g_s^2 GG>^2/(48 pi^2 k^2) is presented without derivation. The check therefore provides only one number, Gamma(eta eta)/Gamma(pi pi) ~ 0.44 versus 0.42, and it does not constrain the pseudoscalar, tensor, pseudotensor, or three-gluon predictions. As a result, the sum rule cannot serve as an independent validation of the main claims beyond the scalar ratio, and the paper should state this limitation more prominently rather than presenting the consistency as a general check of the Fierz framework.","section":"Sec. VI, Eq. (29)"},{"comment":"The paper adopts the masses of f0(1710), eta(2370), f2(2340), and X(2600) as inputs and then uses agreement between the computed branching patterns and the measured or assumed properties of these states to support their glueball interpretation. Using a candidate's mass as input is not circular by itself, but the argument would be much stronger if the framework also showed how a conventional q-qbar state with the same quantum numbers would differ in the same observables. Without such a benchmark, the conclusion that the patterns 'support a sizable gluon component' rests on a comparison against only the glueball candidate, which limits the discriminative power of the analysis.","section":"Sec. IV, mass inputs and interpretation"}],"minor_comments":[{"comment":"In the last displayed term of Eq. (20), the momentum factor is written as (q^nu_a + q^nu_b) with lowercase indices, whereas all other terms in the same equation use q_A and q_B; this appears to be a typographical error.","section":"Eq. (20)"},{"comment":"The Borel-transformed expression in Eq. (30) uses an undefined variable T; please specify the Borel mass parameter and the Borel window used in the numerical analysis.","section":"Sec. VI, Eq. (30)"},{"comment":"The formatting of Table I is very compressed (for example, the kaon-strange-scalar rows), which makes the quark-flavor assignments and quantum numbers harder to read than necessary.","section":"Table I"},{"comment":"The paper says the three-gluon squared amplitudes are 'not particularly useful to display in full,' but the analogous two-gluon squared amplitudes are provided in Appendix B; a brief statement of the phase-space integration method and the Dalitz variables used for Table III would improve reproducibility.","section":"Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the Fierz algebra appears to be carried out carefully, but the central claims are conditioned on an unvalidated effective vertex. I would encourage the editor to request a revision in which the authors either justify the vector-only coupling from QCD or demonstrate that the quoted branching-ratio hierarchies are stable under alternative Lorentz structures."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent extension of the authors' own Fierz-rearrangement framework to 2++ and 2-+ two-gluon glueballs and to a more complete treatment of 0++ and 1+- three-gluon states, plus a leading-order QCD sum rule check for the scalar. It is honest about what it fixes and what it doesn't. The main caveat is the one you flagged: Eq. (16) chooses a vector gluon-to-quark coupling, and every branching ratio in Tables II and III inherits that choice. The paper does not test whether the rankings survive an axial-vector or tensor coupling, so the headline predictions (K*Kbar* dominance, omega-omega dominance for 0-+, the small three-body rates) are conditioned on that ansatz.\n\nWhat is genuinely new and good: the tensor and pseudotensor tables and the three-gluon tables go beyond Ref. [1]; Appendix A gives the full scalar Fierz decomposition; Appendix B lists squared amplitudes for many two-body channels. The sum rule for the scalar gives Gamma(eta eta)/Gamma(pi pi) = 0.44, close to the Fierz value 0.42, and the paper correctly notes that rho rho and omega omega vanish at leading order. The vanishing of all scalar-scalar tensor decays (Eq. 22) is a sharp, testable consequence of the framework. The writing is candid: the authors repeatedly state that normalization and dynamical corrections are undetermined.\n\nSoft spots, in proportion. The biggest is the one above: Eq. (16) is an ansatz, not a QCD-derived current, and the scalar sum-rule check does not rescue the pseudoscalar, tensor, and three-gluon predictions, because that check covers only 0++ and only at leading order. Second, the three-gluon amplitudes are not displayed; a reader cannot check Table III without trusting the intermediate algebra. That is a transparency problem, not a correctness problem. Third, the OPE result in Eq. (29) is stated without derivation; a few lines showing how the double-gluon condensate term arises would help. Fourth, the mass inputs for f0(1710), eta(2370), f2(2340), and X(2600) are candidate identifications; using them as inputs means the support for those identifications is weaker than the abstract suggests. That said, this is not circular in a damaging way: the ratios are not fitted to the data they are supposed to explain.\n\nThe citation pattern looks fine: the self-citations are to the framework being extended, and the literature coverage is broad. I saw no sign of fabrication.\n\nWho this is for: experimentalists planning glueball searches (especially in VV channels) and phenomenologists working in the QCD sum rule / Fierz tradition. It deserves a serious referee. A good referee will ask the authors to test the Lorentz-structure dependence of Eq. (16), display at least representative three-gluon amplitudes, and tone down the interpretation claims for eta(2370) and the tensor. Those are addressable in revision. I would recommend conditional acceptance after such revisions.","headline":"A competent, honest extension of the authors' Fierz program to tensor and three-gluon glueballs, with a useful scalar sum-rule cross-check, but the headline predictions rest on an untested vector-only effective vertex.","tokens_in":21854,"tokens_out":3975,"would_cite":true,"duration_ms":41571,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.39.Mk","13.25.-k"],"model":"deepseek-v4-flash","headline":"By converting each constituent gluon into a quark–antiquark pair and applying exact color and Lorentz Fierz rearrangements, this paper derives glueball branching-ratio patterns that support a sizable gluon component in f0(1710), favor the…","keywords":["glueball","interpolating current","QCD sum rules","Fierz rearrangement","two-gluon glueball","three-gluon glueball","relative branching ratios","f0(1710)"],"falsifier":"Measure the $K^*(892)\\bar K^*(892)$ rate of the tensor glueball candidate relative to $\\pi\\pi$: the Fierz analysis predicts a ratio near 116, so finding that vector–vector modes are not enhanced, or observing a sizable scalar–scalar mode such as $f_0(980)f_0(980)$, would falsify the central decay pattern.","tokens_in":20636,"feed_emoji":"⚛️","tokens_out":11703,"duration_ms":112475,"temperature":0.7,"pith_summary":"The paper tries to show that the decay patterns of glueballs are fixed largely by the algebraic color and Lorentz structure of QCD, not by a detailed dynamical model. Its method converts each constituent gluon into a quark–antiquark pair through a vector current, then uses Fierz rearrangements to rewrite two- and three-gluon glueball currents as sums of ordinary meson operators. The predicted relative branching ratios single out f0(1710) as having a sizable scalar-glueball component, support identifying eta(2370) as the $0^{{-+}}$ glueball, and point to vector–vector channels, especially K*(892) Kbar*(892), as the best place to find the tensor glueball. A separate QCD sum rule calculation for the scalar glueball gives an eta-eta/pi-pi ratio close to the Fierz result, which the authors read as evidence that the algebraic pattern captures the dominant physics. If these identifications are right, future searches already know which specific meson pairs to look at.","feed_headline":"Decay patterns identify f0(1710) and eta(2370) as glueball states","feed_subtitle":"Relative branching ratios from Fierz and sum rules favor eta(2370) as a 0−+ glueball and point tensor searches to K*(892) Kbar*(892).","key_machinery":"The load-bearing device is the effective gluon-to-quark vertex of Eq. (16), $A_i^\\mu \\to \\lambda^i_{ab}\\bar q_a\\gamma^\\mu q_b$, together with the color identity $\\lambda^i_{ab}\\lambda^i_{cd}=2\\delta_{ad}\\delta_{cb}-\\frac{2}{3}\\delta_{ab}\\delta_{cd}$ and the standard Lorentz-space Fierz identities. These transformations convert a glueball current into a finite sum of color-singlet quark-bilinear currents, which are then matched to physical meson states through the current–meson couplings of Table I. Two-gluon currents need one Fierz rearrangement; three-gluon currents need two, producing three-meson final states. The overall normalization of the amplitudes is not fixed, so the framework predicts relative branching ratios rather than absolute widths.","core_discovery":"The central claim is that the relative strengths of glueball decay channels can be obtained by exact Fierz rearrangement of the glueball interpolating currents, once each gluon is converted into a quark–antiquark pair through the effective vertex $A_i^\\mu \\to \\lambda^i_{ab}\\bar q_a\\gamma^\\mu q_b$. For the scalar two-gluon glueball this gives $B(\\eta\\eta)/B(\\pi\\pi)=0.42$ and a pattern that differs from f0(1500), consistent with a sizable gluon component in f0(1710). For the pseudoscalar glueball the calculation produces large vector–vector modes (omega-phi, omega-omega, phi-phi) and a suppressed $K\\bar K^*(892)$ mode, matching the observed behavior of eta(2370) and favoring its $0^{-+}$ glueball interpretation. For the tensor glueball the vector–vector channels dominate, with $K^*(892)\\bar K^*(892)$ the most favorable, and scalar–scalar modes vanish identically at leading order. For three-gluon glueballs the paper identifies $\\pi\\pi\\omega$ and $K\\bar K\\phi$ as promising channels for the $0^{++}$ state and $\\pi\\pi\\omega$, $\\pi\\pi\\phi$, and $K\\bar K\\phi$ for the $1^{+-}$ state.","pith_inferences":["The same algebraic machinery could be applied to other gluonic candidates, such as hybrid mesons, with the vanishing scalar–scalar couplings of the tensor glueball serving as a sharp null test against conventional $q\\bar q$ tensor mesons.","Because every prediction flows from the vector-current vertex of Eq. (16), repeating the calculation with axial-vector or tensor couplings would map how much of the branching-ratio hierarchy is model-dependent.","The predicted $K^*(892)\\bar K^*(892)$ dominance is strong enough to be tested in high-statistics partial-wave analyses of radiative quarkonium decays once the tensor candidate's sample size grows.","If future measurements find sizable rho-rho or omega-omega decays of the scalar glueball, that would point to subleading or non-vector mechanisms that the leading-order Fierz and sum-rule treatments both omit."],"forward_implications":["If f0(1710) is a predominantly two-gluon scalar glueball, its eta-eta to pi-pi branching ratio should be about 0.42 and its decay pattern should differ clearly from that of f0(1500).","If eta(2370) is the $0^{-+}$ glueball, its dominant decay modes should be omega-phi, omega-omega, and phi-phi, with the $K\\bar K^*(892)$ mode suppressed, as recently observed.","A tensor glueball near 2.34 GeV should be searched for in $K^*(892)\\bar K^*(892)$ and other vector–vector channels, while scalar–scalar final states are predicted to be absent at leading order.","The $0^{++}$ three-gluon glueball should be searched for in $\\pi\\pi\\omega$ and $K\\bar K\\phi$, and the $1^{+-}$ state in $\\pi\\pi\\omega$, $\\pi\\pi\\phi$, and $K\\bar K\\phi$.","The QCD sum rule estimate for the scalar glueball gives a total width around 107 MeV with vanishing rho-rho and omega-omega contributions, consistent with the Fierz pattern and supporting the use of the Fierz method as a survey tool."],"supporting_citations":[{"why":"This is the earlier Fierz analysis of two- and three-gluon glueball decays that the present paper extends and develops.","marker":"[1]"},{"why":"It supplies the masses and experimental branching fractions used as inputs and as comparisons for f0(1710), eta(2370), and f2(2340).","marker":"[2]"},{"why":"It provides the interpolating currents adopted for the two-gluon tensor/pseudotensor states and the three-gluon states.","marker":"[24]"},{"why":"It supplies the X(2600) mass used as representative input for the pseudotensor $2^{-+}$ glueball.","marker":"[27]"},{"why":"It provides the experimental $0^{-+}$ assignment of the state near 2370 MeV that is interpreted here as the pseudoscalar glueball.","marker":"[28]"},{"why":"It reports the experimentally suppressed $K\\bar K^*(892)$ mode used as a consistency check for the eta(2370) glueball interpretation.","marker":"[30]"},{"why":"It is the source of the Lorentz-space Fierz identities used in the color and Dirac rearrangements.","marker":"[59]"},{"why":"It provides the gluon condensate value entering the QCD sum rule calculation of the scalar glueball decay.","marker":"[73–75]"}],"fun_headline_variants":["Fierz and sum rules agree: f0(1710) and eta(2370) are glueballs","Tensor glueball decays favor K*(892)Kbar*(892)","eta(2370) fits 0-+ glueball decay pattern","Three-gluon glueballs decay to pi pi omega","Fierz rearrangement predicts glueball branching"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each gluon creates a quark–antiquark pair through the vector current of Eq. (16); if the true gluon-to-meson conversion is not dominated by this coupling, every relative branching ratio in the paper changes.","fun_headline_variants_meta":{"raw":{"variants":["Fierz and sum rules agree: f0(1710) and eta(2370) are glueballs","Tensor glueball decays favor K*(892)Kbar*(892)","eta(2370) fits 0-+ glueball decay pattern","Three-gluon glueballs decay to pi pi omega","Fierz rearrangement predicts glueball branching"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001289,"raw_usage":{"total_tokens":5339,"prompt_tokens":1091,"completion_tokens":4248,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":707,"completion_tokens_details":{"reasoning_tokens":4155}},"tokens_in":707,"tokens_out":4248,"duration_ms":31337,"temperature":1.0,"reasoning_tokens":4155,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:26:28.083727+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the $K^*(892)\\bar K^*(892)$ rate of the tensor glueball candidate relative to $\\pi\\pi$: the Fierz analysis predicts a ratio near 116, so finding that vector–vector modes are not enhanced, or observing a sizable scalar–scalar mode such as $f_0(980)f_0(980)$, would falsify the central decay pattern.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"This is the earlier Fierz analysis of two- and three-gluon glueball decays that the present paper extends and develops."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the masses and experimental branching fractions used as inputs and as comparisons for f0(1710), eta(2370), and f2(2340)."},{"cited_title":"Yamanaka, H","cited_arxiv_id":null,"evidence_quote":"It provides the interpolating currents adopted for the two-gluon tensor/pseudotensor states and the three-gluon states."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the X(2600) mass used as representative input for the pseudotensor $2^{-+}$ glueball."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the experimental $0^{-+}$ assignment of the state near 2370 MeV that is interpreted here as the pseudoscalar glueball."},{"cited_title":"Ablikimet al.[BESIII], Observation of a State X(2600) in theπ +π−η′ System in the ProcessJ/ψ→ γπ +π−η′, Phys","cited_arxiv_id":null,"evidence_quote":"It reports the experimentally suppressed $K\\bar K^*(892)$ mode used as a consistency check for the eta(2370) glueball interpretation."},{"cited_title":"De Rujula, H","cited_arxiv_id":null,"evidence_quote":"It is the source of the Lorentz-space Fierz identities used in the color and Dirac rearrangements."}],"review_version":1}