{"id":"0b31f5c9-7447-4f56-840e-20a7356892eb","arxiv_id":"2608.10570","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"The two gluons from chi_b0 decay carry correlated linear polarizations, so the angular correlation between two pion pairs can directly probe the previously unmeasured linearly polarized gluon dihadron fragmentation function.","lead":"Bottomonium chi_b0 decays into two gluons whose polarizations are linked, and this paper shows the two pion pairs they produce should carry a characteristic angular correlation. Measuring that correlation would give the first direct handle on how linearly polarized gluons fragment into hadrons.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerical predictions integrate over |P_i| ~ M_i where Eq. (1) is invalid; the percent-level asymmetry and Belle-reach claim are unestablished without a collinear cut.","rationale":"The reader correctly identifies the kinematic validity of Eq. (1) as the weakest point. The paper states the condition |P_i| >> M_i immediately before Eq. (1), but the phase-space integration in Fig. 4 uses only z_i >= 2M_i/m_chi, which allows |P_i| = 0 and includes large regions with |P_i| <= M_i. In those regions the pair is slow in the chi_b0 rest frame, so the collinear factorization into independent gluon DiFFs is not justified. Since all numerical statements, including the percent-level asymmetry and the Belle/Belle II reach estimates, come from these integrals, the central quantitative claim is not yet supported. This is a testable, fixable issue rather than a fundamental flaw: imposing a collinearity cut and recomputing A12 would show whether the effect survives. The spectator model is admittedly a benchmark, so model uncertainty alone would not change my verdict. The lack of a full derivation of Eq. (1) is a related weakness, but the kinematic inconsistency is the more concrete, load-bearing problem. I therefore keep the CONDITIONAL verdict, in agreement with the reader.","tokens_in":9813,"tokens_out":16611,"duration_ms":175827,"concrete_test":"Recompute A12(z1,M1), the z1- and M1-projected asymmetries, and the Belle/Belle II sensitivity estimates with an explicit collinearity cut |P_i| >= 2 M_i (equivalently z_i >= 4 M_i/m_chi) imposed on both pion pairs in Eq. (3), keeping the same spectator-model parameters and JAM inputs. If the integrated asymmetry changes by more than a factor of two or drops below about 1%, the 'within reach of Belle' claim is not supported by the factorized calculation as presented.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing concern is an internal inconsistency between the stated domain of Eq. (1) and the kinematic integrations used for the numerical results. Eq. (1) is justified only for |P_i| >> M_i, i.e., z_i m_chi/2 >> M_i for each pion pair. Yet the integrations in Fig. 4 run over z_i >= 2M_i/m_chi, where |P_i| = 0 at the boundary and is comparable to M_i over a substantial part of the allowed phase space (for M_i ~ 1 GeV, |P_i| < M_i when z_i < 0.406). In that region the pion pair is slow in the chi_b0 rest frame, so collinear fragmentation of a single energetic gluon is not justified and the factorized product of two independent DiFFs, together with the cos(2phi1-2phi2) correlation, is undefined. Because A12 and the Belle/Belle II statistical projections are obtained by integrating over these regions, the percent-level benchmark is not a prediction of the stated factorization. No cut or systematic check is provided to isolate the collinear region, so the headline quantitative claim is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new observable for accessing the linearly polarized gluon dihadron fragmentation function H_1^{\\sphericalangle,g}. Specifically, it argues that in the color-singlet leading-order decay \\chi_{b0}\\to gg, the two gluons carry correlated linear polarizations, and that collinear fragmentation of each gluon into a \\pi^+\\pi^- pair generates an Artru--Collins-type angular correlation proportional to \\cos(2\\phi_1-2\\phi_2). The corresponding azimuthally averaged rate constrains the unpolarized gluon DiFF D_1^g. A scalar-spectator model fit to the JAM D_1^g is used to estimate the size of the asymmetry, yielding percent-level values that the authors argue could be accessible with existing Belle data and more precisely with Belle II.","tokens_in":10079,"tokens_out":31228,"duration_ms":311748,"significance":"If the factorization in Eq. (1) is valid, this is a genuinely new and interesting proposal: it would be the first direct experimental access to the linearly polarized gluon dihadron fragmentation function, and it simultaneously provides a new constraint on the poorly known unpolarized gluon DiFF. The paper is honest about the model dependence of its numerical benchmark and about the weak constraints on the JAM gluon DiFF, and the observable is in principle falsifiable at B factories. The central idea deserves serious consideration. However, the numerical results and the experimental-reach claims rest on several load-bearing assumptions that are either not derived or not checked against the stated kinematics, so the present manuscript does not yet fully establish its quantitative conclusions.","major_comments":[{"comment":"The kinematic domain used for the numerical integration violates the stated validity condition of Eq. (1). The factorization is introduced with the requirement |P_i| \\gg M_i, but the Fig. 4 integrations run over z_i\\in[0.19,0.99] and M_i\\in[0.28,2.05] GeV subject only to the threshold constraint z_i\\ge 2M_i/m_{\\chi}. At the boundary |P_i|=0, and for a substantial part of the allowed region |P_i|<M_i (for example, for M_i=2 GeV, |P_i|<M_i for z_i\\lesssim 0.57). In this region the pion pair is not collinear with the fragmenting gluon, so the factored product of two independent DiFFs and the \\cos(2\\phi_1-2\\phi_2) correlation are not defined. The numerical values of A_{12} and the Belle/Belle II statistical projections in Fig. 4 are therefore not predictions of the stated factorization. The authors should impose an explicit collinearity cut, for example |P_i| \\ge \\kappa M_i with \\kappa\\gg 1, and demonstrate that the asymmetry remains after restricting to the region where Eq. (1) is valid.","section":"§4, Fig. 4 and Eq. (1)"},{"comment":"The conversion from the unintegrated DiFFs to the (z,M_h)-dependent functions is incomplete. In Eq. (12), the right-hand sides still depend on |R_T|, but the relation between |R_T| and the integration variables (\\xi, M_h) is never stated. Without this relation and the corresponding Jacobian, the definitions of D_1^g(z,M_h) and H_1^{\\sphericalangle,g}(z,M_h) are mathematically ill-defined, and the spectator-model results in Figs. 3 and 4 cannot be reproduced. The authors should give the explicit mapping, e.g. M_h^2=4(m_\\pi^2+R_T^2)/(1-\\xi^2) if that is the intended relation, and the resulting measure.","section":"§3, Eq. (12)"},{"comment":"The statistical sensitivity estimates are not reproducible because the event number N is never evaluated. The text defines N as the number of selected events after kinematic cuts, but no estimate is provided for the expected \\chi_{b0} yield at Belle with \\mathcal{L}=24.7\\,\\text{fb}^{-1}, nor are the relevant branching ratios (\\Upsilon(2S)\\to\\gamma\\chi_{b0}, \\chi_{b0}\\to 2\\pi^+2\\pi^-+X) or selection efficiencies given. Without these inputs, the plotted statistical bands and the claim that existing Belle data may already be sensitive to the asymmetry are unsupported. The authors should provide a concrete event-count estimate and state the cuts used.","section":"§4, Eq. (14) and Fig. 4"},{"comment":"The factorization formula itself is asserted rather than derived. The text states the result with a citation to collinear factorization, but does not show how the NRQCD decay amplitude for \\chi_{b0}\\to gg is combined with the gluon fragmentation functions, how the thrust axis replaces the partonic axis in the presence of hadron transverse momenta, or how the hard coefficients C^{q,g} in Eq. (2) are obtained. Since Eq. (1) is the foundation of the proposed observable, a more explicit derivation or a reference in which the same double-factorization is established should be provided.","section":"§2, Eq. (1)"},{"comment":"The numerical benchmark is model-dependent in a way that is not reflected in the quoted uncertainty bands. The 10-parameter spectator-model fit is anchored to the JAM D_1^g, which the paper itself describes as weakly constrained and obtained under indirect assumptions, and no theory uncertainty from the fit is propagated to A_{12}. The text correctly says the best-fit values are used only as a benchmark, but the abstract and conclusion state 'percent-level asymmetries, potentially within reach of existing Belle data' without the accompanying caveat. The authors should either present a range of benchmark values reflecting the fit uncertainties and model sensitivity, or explicitly reframe the reach claim as conditional on the spectator model.","section":"§4, Fig. 4 and Eq. (3)"}],"minor_comments":[{"comment":"Equation (4) contains a garbled factor '(k -)^2' that appears to be a typo; it should presumably be (k-P_h)^2 or an explicit k^2 factor. Please correct the notation.","section":"§3, Eq. (4)"},{"comment":"The statement that the scalar \\chi_{b0} decay 'does not induce a transverse-spin correlation between the produced quark and antiquark' is too strong. A J=0 decay can produce a spin-singlet q\\bar q correlation, and the color-octet ^3S_1^{[8]} channel is a spin-triplet. The relevant point is that the quark-channel spin correlations do not project onto the particular cos(2\\phi_1-2\\phi_2) moment, because quark interference DiFFs enter at first order in the azimuthal angle. Please clarify the argument.","section":"§2, text after Eq. (2)"},{"comment":"The fitted spectator-model parameters are quoted without uncertainties or a fit-quality indicator. Given that the fit is used as the numerical basis for the benchmark, a \\chi^2/dof and parameter covariance would be needed for reproducibility.","section":"§3, Eq. (13)"},{"comment":"The notation H_b^1 for the NRQCD matrix element is nonstandard and can be confused with a fragmentation function; consider using the conventional \\langle O(^3P_0^{[1]})\\rangle notation.","section":"Throughout"},{"comment":"The qualitative discussion of \\eta_b and \\chi_{b2} is useful, but the title and abstract promise a treatment of '\\chi_b decays' in general; the actual quantitative analysis is only for \\chi_{b0}. It would help to state explicitly that \\eta_b and \\chi_{b2} are left for future work.","section":"Introduction and Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The central observable is well motivated and potentially important, but the quantitative claims need the kinematic, definitional, and statistical fixes listed above. I am also somewhat concerned that the paper leans on an unpublished preprint [35] for a key input (the independent determination of \\rho_8); the authors should either make that work openly available or soften the dependence. The benchmark claim in the abstract should be tempered until the spectator-model uncertainty and the event-count estimate are provided."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What's new: this is the first proposal I've seen that gets at the linearly polarized gluon DiFF H_1^{angular,g} through two-gluon decay of chi_b0 and a cos(2phi1-2phi2) correlation between dihadron pairs. The combination of NRQCD and collinear factorization is natural, and the paper is honest that the quark CO channel dilutes but cannot fake the asymmetry. The spectator-model benchmark is clearly labeled as such, with the authors admitting the fitted parameters carry large uncertainties. That is the right attitude for a first estimate.\n\nThe soft spot that matters: Eq. (1) is stated to hold for |P_i| >> M_i, but the numerical integrations in Fig. 4 run down to z_i = 2M_i/m_chi, where |P_i| = 0. For M_i around 1 GeV, a substantial part of the integrated phase space has |P_i| < M_i. In that region, collinear fragmentation of a single gluon into the pair is not justified, and the factorized form with two independent DiFFs is not defined. The stress-test note is right: the percent-level asymmetry and the Belle/Belle II reach claim are not yet predictions of the stated factorization. This is fixable by imposing a cut like |P_i| > 2 M_i (or z_i > 4 M_i/m_chi) and showing the result is stable, but it has to be done before the numbers mean anything.\n\nMinor points: the derivation of Eq. (1) is only sketched, which is normal for a Letter, but a reader needs to see why the CS and CO terms add incoherently and how the thrust axis maps to the partonic gg axis. The choice mu = m_chi with no scale variation is acceptable for a benchmark, but it would be good to state that the prediction is LO only. The paper also does not propagate any theory uncertainty from the spectator model; again fine for a Letter if the kinematic cut issue is resolved.\n\nOverall, the central idea holds up. The factorization structure is standard, the spin correlation is derived from the P-wave spin state, and the proposed observable is genuinely novel. The numerics need work, not the concept.\n\nI would send this to peer review. A serious referee can sort out the kinematic cut and push for a more careful statement of the factorization domain. The paper will be useful to anyone working on gluon TMDs, DiFFs, or quarkonium decays; it is a solid pointer for the community, even if the current numbers are not final.","headline":"A genuinely new observable for the linearly polarized gluon dihadron fragmentation function, with an honest but model-dependent benchmark and a kinematic-region inconsistency in the numerics that should be fixed before publication.","tokens_in":10639,"tokens_out":2236,"would_cite":true,"duration_ms":25339,"reading_group":"maybe","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 argues that $\\chi_{b0}$ decays give a first direct probe of the linearly polarized gluon dihadron fragmentation function, through an Artru–Collins-type angular correlation between the two pion pairs.","keywords":["gluon linear polarization","dihadron fragmentation function","Artru-Collins asymmetry","chi_b0 decay","NRQCD factorization","bottomonium","spectator model","unpolarized gluon dihadron fragmentation function"],"falsifier":"A measurement of the $\\cos(2\\phi_1-2\\phi_2)$ asymmetry in $\\chi_{b0}\\to\\pi^+\\pi^-\\pi^+\\pi^-X$ over the same $z$ and $M$ bins as Fig. 4, restricted to high-thrust events where the collinear condition $|P_i|\\gg M_i$ holds; if the modulation is statistically consistent with zero where the spectator-model benchmark predicts percent-level values, while the azimuthally averaged rate matches the predicted $D_1^g$, the central claim is falsified. A second check is the sign and the $z_1$-dependence: the model predicts a specific growing trend with $z_1$, so a flat or opposite trend would rule out the proposed linear-polarization mechanism.","tokens_in":9601,"feed_emoji":"⚛️","tokens_out":10975,"duration_ms":101104,"temperature":0.7,"pith_summary":"This paper proposes a way to measure, for the first time, the dihadron fragmentation function of a linearly polarized gluon. The idea is to use the two-gluon decay of the $P$-wave bottomonium state $\\chi_{b0}$: at leading order the two gluons carry correlated linear polarizations, and when each fragments into a $\\pi^+\\pi^-$ pair, the two pairs inherit an angular correlation of the Artru–Collins type, specifically a $\\cos(2\\phi_1-2\\phi_2)$ modulation. If the prediction holds, measuring that correlation directly accesses the linearly polarized gluon dihadron fragmentation function $H_1^{\\angle,g}$, while the azimuthally averaged rate constrains the unpolarized gluon dihadron fragmentation function $D_1^g$. A spectator-model estimate puts the asymmetry at the percent level, within reach of existing lepton-collider data, with much better precision at a future high-luminosity run. This would open a new window on the spin-dependent hadronization of gluons, which has so far remained experimentally unconstrained.","feed_headline":"χ_b0 decays may reveal the gluon's hidden linear polarization","feed_subtitle":"An azimuthal angle between the two pion pairs carries the signal at the percent level—within reach of current data.","key_machinery":"The load-bearing object is the linearly polarized gluon dihadron fragmentation function $H_1^{\\angle,g}$, the gluon analogue of the quark interference dihadron fragmentation function that appears in Collins-type asymmetries. It enters the decay distribution through the product $H_1^{\\angle,g}(z_1,M_1)H_1^{\\angle,g}(z_2,M_2)$ multiplying $\\cos(2\\phi_1-2\\phi_2)$, the Artru–Collins-type angular correlation between the two dihadron planes. The derivation combines NRQCD factorization, which separates the $\\chi_{b0}$ annihilation into short-distance coefficients and long-distance matrix elements, with collinear factorization for the fragmentation of each gluon into a $\\pi^+\\pi^-$ pair; the hard two-gluon production from the color-singlet channel is what supplies the correlated linear polarizations. To make a numerical prediction, the gluon dihadron fragmentation functions are modeled with a scalar-spectator model for $g\\to\\pi^+\\pi^-X$, with parameters fitted against the global estimate of $D_1^g$ and constrained by the $\\rho$ and $\\omega$ resonance structure in the pion-pair invariant mass.","core_discovery":"At leading order in the NRQCD velocity expansion, the color-singlet component of the $P$-wave bottomonium state $\\chi_{b0}$ annihilates into two gluons; because the $\\chi_{b0}$ is a scalar, the two gluons emerge with correlated linear polarizations. In collinear factorization, the fragmentation of those gluons into two dihadron pairs gives a differential decay distribution with an azimuthally independent term proportional to $D_1^g(z_1,M_1)D_1^g(z_2,M_2)$ and a modulation proportional to $H_1^{\\angle,g}(z_1,M_1)H_1^{\\angle,g}(z_2,M_2)\\,\\cos(2\\phi_1-2\\phi_2)$. The paper shows that this modulation is the first direct experimental probe of the linearly polarized gluon dihadron fragmentation function $H_1^{\\angle,g}$, and that the same decay rate provides direct access to the poorly constrained unpolarized gluon dihadron fragmentation function $D_1^g$. The color-octet subprocess, which produces a light quark-antiquark pair instead of two gluons, dilutes the asymmetry but cannot generate it, because the scalar decay does not correlate the quark and antiquark transverse spins. Using a spectator model for $g\\to\\pi^+\\pi^-X$ tuned to the current global estimate of $D_1^g$, the asymmetry is predicted at the percent level, suggesting existing data could already be sensitive to this observable.","pith_inferences":["If the asymmetry is confirmed, the same observable could be used to test the universality of $H_1^{\\angle,g}$ by comparing with future measurements in other hard processes, since fragmentation functions are universal by the QCD factorization theorems.","The extraction of the LDME ratio $\\rho_8$ from different input data feeds into the denominator of the asymmetry; comparing the $z_1$-dependence of $A_{12}$ at large $z_1$, where quark dilution is strongest, could help resolve the current discrepancy between different determinations of $\\rho_8$.","A null result would not by itself rule out gluon linear polarization, since the spectator-model normalization is a benchmark rather than a rigorous prediction; the decisive test is whether the $\\cos(2\\phi_1-2\\phi_2)$ modulation appears with the predicted sign and kinematic shape, not just its integrated size.","The thrust-axis reconstruction is an experimental proxy for the partonic $gg$ axis; a dedicated comparison of the measured modulation at different thrust values would test the factorization assumption underlying the prediction."],"forward_implications":["A measurement of the $\\cos(2\\phi_1-2\\phi_2)$ asymmetry in $\\chi_{b0}\\to\\pi^+\\pi^-\\pi^+\\pi^-X$ would give the first direct experimental determination of the linearly polarized gluon dihadron fragmentation function.","The azimuthally averaged semi-inclusive decay rate directly constrains the unpolarized gluon dihadron fragmentation function $D_1^g$, which current global analyses constrain only indirectly or not at all.","At the benchmark estimate the asymmetry is at the percent level, so existing data collected at the $\\Upsilon(2S)$ resonance may already have statistical sensitivity to the signal.","A dedicated high-luminosity data set at the same energy would substantially improve the statistical precision, allowing the kinematic dependence of both gluon dihadron fragmentation functions to be mapped in $z$ and invariant mass.","The same framework extends to $\\eta_b$ decays, where the leading-order color-octet quark channel is absent, providing a cleaner probe of gluon fragmentation once enough data become available."],"supporting_citations":[{"why":"Supplies the collinear factorization of dihadron fragmentation and the definition of the interference dihadron fragmentation functions used in Eq. (1).","marker":"[10]"},{"why":"Provides NRQCD factorization, separating the short-distance annihilation coefficients from the long-distance matrix elements in the decay formula.","marker":"[40]"},{"why":"Provides the global extraction of the unpolarized gluon dihadron fragmentation function used as the benchmark for the spectator-model fit and numerical estimates.","marker":"[19]"},{"why":"Supplies the spectator-model parameterization for gluon fragmentation into a pion pair, including the resonant s- and p-wave structures used in Eq. (7).","marker":"[51]"},{"why":"Gives the lattice-NRQCD value of the LDME ratio $\\rho_8=0.044$ used as the reference value in the benchmark asymmetry calculation.","marker":"[49]"},{"why":"Provides the alternative extraction $\\rho_8=0.160$ used to test how sensitive the predicted asymmetry is to the poorly known LDME ratio.","marker":"[47]"},{"why":"Shows a previous independent determination of $\\rho_8$ from $\\chi_{b2}$ dihadron fragmentation observables, which the current paper builds on to motivate the observable's discriminating power.","marker":"[35]"}],"fun_headline_variants":["χ_b0 decays reveal gluon's hidden spin via pion pairs","Gluon polarization glimpsed in χ_b0 two-gluon decay","First direct probe of gluon linear polarization in χ_b0","χ_b0 decay yields azimuthal clue to gluon fragmentation","Percent-level asymmetry in χ_b0 decay exposes gluon spin"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction stands on the assumption that the collinear factorization formula, which requires each pion-pair momentum to be much larger than the pair's invariant mass, remains valid over the full kinematic ranges integrated in Fig. 4, including regions near $z_i = 2M_i/m_\\chi$ where the pair momentum can vanish.","fun_headline_variants_meta":{"raw":{"variants":["χ_b0 decays reveal gluon's hidden spin via pion pairs","Gluon polarization glimpsed in χ_b0 two-gluon decay","First direct probe of gluon linear polarization in χ_b0","χ_b0 decay yields azimuthal clue to gluon fragmentation","Percent-level asymmetry in χ_b0 decay exposes gluon spin"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000653,"raw_usage":{"total_tokens":3031,"prompt_tokens":1023,"completion_tokens":2008,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":639,"completion_tokens_details":{"reasoning_tokens":1918}},"tokens_in":639,"tokens_out":2008,"duration_ms":15096,"temperature":1.0,"reasoning_tokens":1918,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:37:06.504518+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the $\\cos(2\\phi_1-2\\phi_2)$ asymmetry in $\\chi_{b0}\\to\\pi^+\\pi^-\\pi^+\\pi^-X$ over the same $z$ and $M$ bins as Fig. 4, restricted to high-thrust events where the collinear condition $|P_i|\\gg M_i$ holds; if the modulation is statistically consistent with zero where the spectator-model benchmark predicts percent-level values, while the azimuthally averaged rate matches the predicted $D_1^g$, the central claim is falsified. A second check is the sign and the $z_1$-dependence: the model predicts a specific growing trend with $z_1$, so a flat or opposite trend would rule out the proposed linear-polarization mechanism.","supporting_citations":[{"cited_title":"Inclusive Charm Production in chi_b Decays","cited_arxiv_id":"0704.2599","evidence_quote":"Supplies the spectator-model parameterization for gluon fragmentation into a pion pair, including the resonant s- and p-wave structures used in Eq. (7)."},{"cited_title":"Inclusive chi_bJ(nP) Decays to D0 X","cited_arxiv_id":"0807.3757","evidence_quote":"Gives the lattice-NRQCD value of the LDME ratio $\\rho_8=0.044$ used as the reference value in the benchmark asymmetry calculation."},{"cited_title":"Unveiling Light-Quark Yukawa Flavor Structure via Dihadron Fragmentation at Lepton Colliders","cited_arxiv_id":"2512.16492","evidence_quote":"Shows a previous independent determination of $\\rho_8$ from $\\chi_{b2}$ dihadron fragmentation observables, which the current paper builds on to motivate the observable's discriminating power."}],"review_version":1}