{"id":"a86fbf14-1e63-457f-807f-49c27523d124","arxiv_id":"2606.24455","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"NNLO perturbative calculation with PMC scale setting for η_c2 and η_b2 → γγ decay widths and branching ratios in NRQCD effective theory.","lead":"This paper computes NNLO QCD corrections to the two-photon decay of D-wave spin-singlet heavy quarkonia in NRQCD, applying the Principle of Maximum Conformality to remove renormalization scale dependence. A smart generalist might read it to understand how higher-order perturbative techniques and scale-setting methods refine predictions for rare decays that test strong-interaction bound states.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"LDME uncertainties omitted from quoted errors, despite being the dominant nonperturbative input for the width","rationale":"The reader’s weakest_assumption directly flags the LDME accuracy and NRQCD validity; the abstract’s explicit uncertainty statement confirms those errors are not included, making this the single most load-bearing gap for the numerical claim. Full-text details on LDME sources would be needed to close it.","tokens_in":2044,"tokens_out":382,"duration_ms":17635,"concrete_test":"Locate the numerical inputs section; extract the central LDME values and any quoted uncertainties used for η_c2 and η_b2; recompute both Γ^PMC values after varying each LDME by its stated error (added in quadrature to the existing sources); if the total relative uncertainty increases by >25% the headline error bands are incomplete.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim gives Γ_ηc2→γγ^PMC = 3.322^{+0.899}_{-0.828} eV (and analog for bottomonium) with uncertainties stated to arise only from Δm_Q and uncalculated higher orders. The width factorizes as Γ ~ |SDC_NNLO(Q_*)|^2 × LDME, where the D-wave LDME (e.g. ⟨O^η_Q2(1D0)⟩ or equivalent) is taken from external sources. No propagation of LDME uncertainty appears in the error budget. For D-wave states these matrix elements are typically extracted from potential models or fits and carry 15–40% relative errors; including them would enlarge the total uncertainty and weaken the quoted precision. The paper’s PMC procedure also invokes RG running of the LDMEs themselves, which presupposes a factorization scale choice whose consistency with the nonperturbative definition is not independently verified.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper computes NNLO QCD corrections to the two-photon decays of D-wave spin-singlet heavy quarkonia η_{Q2} within NRQCD factorization. It decomposes the width into short-distance coefficients (SDCs) and D-wave LDMEs, applies the Principle of Maximum Conformality (PMC) to eliminate renormalization- and factorization-scale dependence by recursively using RGEs for α_s and the LDMEs, determines optimal scales Q*=1.483 GeV (charm) and 4.246 GeV (bottom), and reports scale-invariant numerical results Γ_ηc2→γγ^PMC = 3.322^{+0.899}_{-0.828} eV and Γ_ηb2→γγ^PMC = 0.0188^{+0.0014}_{-0.0013} eV together with the corresponding branching ratios.","tokens_in":2265,"tokens_out":641,"duration_ms":16893,"significance":"If the central results hold after proper error propagation, the work supplies the first PMC-improved NNLO predictions for these rare decays, removing conventional scale ambiguities and potentially improving convergence of the perturbative series. Such predictions could be compared with future experimental searches at e+e- or hadron colliders and help test the applicability of NRQCD factorization at NNLO for D-wave states.","major_comments":[{"comment":"Abstract: the quoted uncertainties for Γ_ηc2→γγ^PMC = 3.322^{+0.899}_{-0.828} eV (and the analogous bottomonium result) are stated to arise only from Δm_Q and uncalculated higher orders. However, the width factorizes as Γ ~ |SDC_NNLO(Q*)|^2 × LDME; the D-wave LDMEs are external nonperturbative inputs whose typical 15–40% relative uncertainties are omitted from the error budget. This omission is load-bearing for the claimed precision.","section":"abstract"},{"comment":"Abstract (and numerical-results section): no explicit numerical values, sources, or references are supplied for the D-wave LDMEs (e.g. ⟨O^η_Q2(1D0)⟩ or equivalent) that enter the final widths. Without these inputs the numerical claims cannot be independently verified or reproduced.","section":"abstract"},{"comment":"Abstract: the claim that PMC “naturally improves the convergence of the perturbative series for SDCs” is asserted but not demonstrated by any explicit comparison of the conventional versus PMC series (e.g., term-by-term ratios or partial sums) in the provided text.","section":"abstract"},{"comment":"Factorization and PMC implementation: the recursive application of RGEs to the D-wave LDMEs presupposes a consistent factorization-scale choice whose matching to the nonperturbative definition of the LDMEs is not independently verified; any mismatch would introduce an additional systematic uncertainty not quantified in the error budget.","section":null}],"minor_comments":[],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the thorough review and constructive comments on our manuscript. We address each major comment point by point below, agreeing where revisions are needed to enhance clarity and completeness while defending our approach on substantive grounds.","responses":[{"response":"We agree that a complete error budget should incorporate the uncertainties from the input D-wave LDMEs. The abstract emphasized the perturbative improvements from PMC, but the full manuscript uses LDME values from the literature. We will revise the abstract and numerical-results section to propagate and quote the LDME uncertainties (typically 15-40%) alongside the existing ones, providing a more comprehensive total uncertainty.","revision_made":"yes","referee_comment":"[abstract] Abstract: the quoted uncertainties for Γ_ηc2→γγ^PMC = 3.322^{+0.899}_{-0.828} eV (and the analogous bottomonium result) are stated to arise only from Δm_Q and uncalculated higher orders. However, the width factorizes as Γ ~ |SDC_NNLO(Q*)|^2 × LDME; the D-wave LDMEs are external nonperturbative inputs whose typical 15–40% relative uncertainties are omitted from the error budget. This omission is load-bearing for the claimed precision."},{"response":"The full manuscript provides the LDME values, sources, and references in the numerical-results section. To improve accessibility, we will add the explicit LDME numerical inputs and their references directly into the abstract (or a footnote) and ensure the numerical-results section cross-references them clearly for reproducibility.","revision_made":"yes","referee_comment":"[abstract] Abstract (and numerical-results section): no explicit numerical values, sources, or references are supplied for the D-wave LDMEs (e.g. ⟨O^η_Q2(1D0)⟩ or equivalent) that enter the final widths. Without these inputs the numerical claims cannot be independently verified or reproduced."},{"response":"The manuscript text discusses the removal of renormalon contributions via PMC, but we acknowledge an explicit side-by-side comparison would strengthen the claim. We will add a dedicated paragraph or table in the revised manuscript showing term-by-term ratios and partial sums for both conventional and PMC series to demonstrate the improved convergence.","revision_made":"yes","referee_comment":"[abstract] Abstract: the claim that PMC “naturally improves the convergence of the perturbative series for SDCs” is asserted but not demonstrated by any explicit comparison of the conventional versus PMC series (e.g., term-by-term ratios or partial sums) in the provided text."},{"response":"The PMC scale choice is determined self-consistently within the NRQCD factorization framework, with RGEs applied to both α_s and LDMEs to ensure scale invariance. This follows standard practice in PMC applications to effective theories. We will expand the discussion in the revised manuscript to clarify the matching procedure and note that any residual mismatch is absorbed into the LDME definition, but we do not introduce an additional ad-hoc systematic beyond the quoted higher-order errors.","revision_made":"partial","referee_comment":"[—] Factorization and PMC implementation: the recursive application of RGEs to the D-wave LDMEs presupposes a consistent factorization-scale choice whose matching to the nonperturbative definition of the LDMEs is not independently verified; any mismatch would introduce an additional systematic uncertainty not quantified in the error budget."}],"tokens_in":1845,"tokens_out":739,"duration_ms":18113,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is a clean NNLO calculation of the two-photon widths for the D-wave states eta_c2 and eta_b2 inside NRQCD, using PMC to remove renormalization and factorization scale dependence. The authors run both alpha_s and the LDMEs with the RGE to obtain fixed PMC scales (1.48 GeV and 4.25 GeV) and report widths of 3.322 eV and 0.0188 eV with branching ratios around 7.5e-6 and 6.5e-7.\n\nThe work does what it sets out to do: it extends the PMC procedure to this previously uncomputed channel at NNLO, shows the series convergence improves once the renormalon pieces are absorbed, and produces explicit numerical predictions. The factorization setup and the recursive application of the RGE look standard and internally consistent.\n\nThe soft spot is the error budget. The quoted uncertainties come only from quark-mass variations and an estimate of higher orders; the D-wave LDMEs, which multiply the short-distance coefficients, are taken from outside sources and carry their own 15-40% uncertainties that are not propagated. For D-wave states those matrix elements are the leading nonperturbative input, so the total theoretical error is larger than presented. The paper also assumes NRQCD factorization remains reliable at this order without large velocity corrections, which is plausible but not re-examined here.\n\nThis is a narrow but competent piece for people who work on quarkonium decays and scale-setting methods. A reader who already follows PMC applications or needs updated numbers for these channels will get value from it. It is worth sending to a serious referee; the calculation is non-trivial and the method is applied without obvious internal contradictions, even though the nonperturbative error treatment needs checking.","headline":"This paper gives the first NNLO PMC results for eta_Q2 to gamma gamma decays, with cleaner scale dependence, but the error budget skips the dominant LDME uncertainties.","tokens_in":2752,"tokens_out":444,"would_cite":false,"duration_ms":15710,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"PMC removes scale ambiguities from NNLO predictions for D-wave eta_Q2 to two-photon decays.","keywords":["NNLO QCD corrections","Principle of Maximum Conformality","heavy quarkonia","D-wave states","two-photon decay","NRQCD factorization","eta_c2","eta_b2"],"falsifier":"A measurement of Br(eta_c2 -> gamma gamma) lying outside the interval 5.6 x 10^{-6} to 9.5 x 10^{-6} would falsify the PMC prediction.","tokens_in":2936,"feed_emoji":"","tokens_out":870,"duration_ms":15411,"temperature":0.7,"pith_summary":"This paper computes the two-photon decays of D-wave spin-singlet heavy quarkonia eta_Q2 at NNLO in QCD using NRQCD factorization. The short-distance coefficients originally carry large renormalization and factorization scale uncertainties. The Principle of Maximum Conformality is applied by recursively solving the renormalization group equations for alpha_s and the D-wave long-distance matrix elements, producing an effective scale Q* that absorbs all nonconformal beta terms. The resulting series becomes scale-invariant and more convergent, giving explicit numerical widths with uncertainties dominated by quark-mass variations. A reader cares because the method supplies cleaner theoretical values for these rare processes that experiments can target.","feed_headline":"PMC fixes scale dependence in NNLO eta_Q2 gamma gamma widths","feed_subtitle":"Charmonium width 3.322 eV and bottomonium width 0.0188 eV after determining optimal scales Q* = 1.483 GeV and 4.246 GeV.","key_machinery":"The Principle of Maximum Conformality, which sets the renormalization scale by absorbing all beta-function dependent terms into the running coupling to obtain a conformal, scale-invariant perturbative series for the short-distance coefficients.","core_discovery":"The total decay width factors into perturbatively calculable short-distance coefficients and nonperturbative D-wave long-distance matrix elements. The original NNLO series for the coefficients suffers from sizable scale dependence. Recursively applying the renormalization group equations inside the PMC framework determines an effective strong coupling alpha_s(Q*) consistent with the expansion coefficients, removing divergent renormalon contributions and yielding a scale-invariant perturbative series. The PMC scales are fixed at Q*=1.483 GeV for eta_c2 and Q*=4.246 GeV for eta_b2, producing Gamma(eta_c2->gamma gamma)^PMC = 3.322^{+0.899}_{-0.828} eV and Gamma(eta_b2->gamma gamma)^PMC = 0.0188","pith_inferences":["These widths can be compared directly with future measurements at e+e- or hadron colliders to test the combined NRQCD+PMC framework.","The same PMC procedure can be applied to other rare decays of D-wave or higher-L states once the corresponding LDMEs are available.","Agreement with data would support extending the method to N3LO or to processes involving different final states."],"forward_implications":["The decay widths receive their dominant remaining uncertainties from variations in the heavy-quark masses.","The branching ratios are (7.463^{+2.020}_{-1.860})x10^{-6} for eta_c2 and (6.460^{+0.481}_{-0.447})x10^{-7} for eta_b2.","The perturbative series for the short-distance coefficients converges better after removal of renormalon contributions.","NRQCD factorization separates short-distance coefficients from long-distance matrix elements up to NNLO for this channel."],"fun_headline_variants":["PMC resolves scale ambiguity in NNLO eta_Q2 to gamma gamma","NNLO D-wave quarkonia gamma decay with PMC scale setting","Effective scales set by PMC for eta_c2 and eta_b2 gamma decays","PMC yields scale-invariant NNLO series for heavy quarkonia decays","PMC determines Q* for NNLO eta_c2 and eta_b2 gamma gamma widths"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The nonperturbative D-wave long-distance matrix elements are known to sufficient accuracy from independent sources and NRQCD factorization remains valid without large velocity-suppressed corrections at NNLO.","fun_headline_variants_meta":{"raw":{"variants":["PMC resolves scale ambiguity in NNLO eta_Q2 to gamma gamma","NNLO D-wave quarkonia gamma decay with PMC scale setting","Effective scales set by PMC for eta_c2 and eta_b2 gamma decays","PMC yields scale-invariant NNLO series for heavy quarkonia decays","PMC determines Q* for NNLO eta_c2 and eta_b2 gamma gamma widths"]},"model":"grok-4.3","cost_usd":0.006065,"raw_usage":{"total_tokens":3036,"prompt_tokens":1005,"num_sources_used":0,"completion_tokens":93,"cost_in_usd_ticks":60649500,"prompt_tokens_details":{"text_tokens":1005,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1938,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":1005,"tokens_out":93,"duration_ms":12732,"temperature":1.0,"reasoning_tokens":1938,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-25T23:39:38.551119+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A measurement of Br(eta_c2 -> gamma gamma) lying outside the interval 5.6 x 10^{-6} to 9.5 x 10^{-6} would falsify the PMC prediction.","supporting_citations":[],"review_version":1}