{"id":"542ab194-038d-47fd-bf68-a1f755689598","arxiv_id":"2501.04188","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Analytic two-loop form factors for chi_{Q,J} production and decay are presented, including a new NRQCD singularity that couples P-wave states to the 3S1[8] configuration.","lead":"This paper computes the two-loop QCD corrections (form factors) needed for production and decay of P-wave heavy quarkonium states (chi_{Q,J}) at next-to-next-to-leading order. It finds a new infrared singularity in gluon-fusion channels that requires mixing with the 3S1[8] state, and provides high-precision hard functions for phenomenology.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"P-wave IBP reduction to the same 76 S-wave master integrals is asserted (footnote 9) but not demonstrated; if incomplete, every new two-loop form factor is invalid.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the P-wave two-loop integrals are asserted, not demonstrated, to be spanned by the 76 master integrals of the S-wave companion paper. This is the foundation of all analytic results in the paper; an incomplete IBP reduction would invalidate the gg, γg, and colour-octet channels that are the main new claims. The paper provides meaningful external cross-checks for the γγ channel and the one-loop 3S1[8] form factor, but those checks do not cover the reduction premise for the new channels. The concern is therefore real, concrete, and testable, and it justifies keeping the reader's conditional verdict rather than upgrading to acceptance. No internal inconsistency was found in the presented pole-cancellation structure or in the finite remainders; the identified gap is one of missing verification rather than demonstrated error.","tokens_in":71750,"tokens_out":12896,"duration_ms":130652,"concrete_test":"Independently reduce, with a different IBP solver (e.g., Kira or FIRE), the full set of two-loop P-wave integrals generated after the Eq. (2.6) derivative for all channels, including all numerator insertions and raised propagator powers, and verify that the master-integral sector count matches the 76 MIs of ref. [41] at every order in epsilon; if any additional master appears, the analytic form factors and hard functions in Section 6 are incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The entire calculation rests on the structural premise, stated in Section 2 and footnote 9, that after the relative-momentum derivative in Eq. (2.6), all two-loop P-wave integrals reduce via IBP to the same 76 master integrals as the 1S0 case of ref. [32]. The derivative is taken before q is set to zero and produces numerator insertions and raised propagator powers; in principle these can generate new master integrals at a given order in epsilon. The paper gives only a confirmation ('we confirm') and a speculative aside about D waves; no reduction tables, no explicit master-count comparison, and no proof of closure of the reduction identities are provided. Because every bare amplitude, renormalized finite remainder, and hard function in Section 6 is expressed as a linear combination of those 76 MIs, a single missing master would make all four channels' two-loop results incomplete, not just one helicity configuration. The cross-checks against ref. [38] (γγ channel) and refs. [82,84] (one-loop 3S1[8]) are genuine external support, but they do not exercise the new gg and γg P-wave channels where the reduction claim is most needed. This is a load-bearing but checkable gap rather than a demonstrated error.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript computes analytic two-loop form factors and helicity amplitudes for the P-wave NRQCD channels gamma-gamma -> 3P_J[1], g-g -> 3P_J[1], gamma-g -> 3P_J[8], and g-g -> 3P_J[8]. The calculation uses dimensional regularization, projects the heavy-quark pair onto a P-wave state via a derivative in the relative momentum q, performs IBP reduction, and expresses all bare amplitudes in terms of the 76 master integrals of the companion papers [32,41]. After on-shell and MS-bar renormalization, the authors subtract standard QCD IR singularities and an NRQCD subtraction factor Z_NRQCD. A central new claim is that the gg channels contain, in addition to the Coulomb singularity, a new NRQCD pole that is cancelled only after including the one-loop form factor gg -> 3S_1[8] through a helicity-space matrix in Z_NRQCD. The paper provides finite remainders, hard functions with scale dependence, 20-digit numerical bases, and analytic expressions in Appendix F and in electronic form.","tokens_in":71970,"tokens_out":5190,"duration_ms":59716,"significance":"If the results are correct, this is a substantial step for quarkonium physics: it supplies the first analytic two-loop results for the P-wave channels gg -> 3P_J[1], gamma-g -> 3P_J[8] and gg -> 3P_J[8], improves the previously numerical gamma-gamma channel, and identifies a previously unnoticed NRQCD pole structure whose cancellation requires the gg -> 3S_1[8] form factor. The paper ships no fitted parameters, agrees with the independent numerical calculation of ref. [38] for gamma-gamma, and its one-loop gg -> 3S_1[8] result agrees with refs. [82,84] up to the reported order. The new P-wave channels lack independent external checks, and the central technical assumption is the closure of the IBP reduction on the same 76 master integrals as the S-wave case; that assumption is currently asserted rather than demonstrated.","major_comments":[{"comment":"The statement that all two-loop P-wave integrals reduce via IBP to the same 76 master integrals as the 1S0 case is load-bearing but is only asserted. Equation (2.6) differentiates in q before setting q=0, which produces raised propagator powers and numerator insertions; after partial-fraction decomposition and IBP those objects can, in principle, generate new master integrals at a given order in epsilon. Since every bare amplitude in Section 3 and every finite remainder in Section 6 is expressed as a linear combination of the 76 master integrals, a single missing master would invalidate all four channels, not just one helicity configuration. Please supply the reduction evidence: for example, per-topology integral counts before and after IBP, the master-integral count at each epsilon order, and a check that the reduction identities close. Uploading reduction tables to the ancillary directory would be sufficient.","section":"§2, footnote 9; §7"},{"comment":"The cancellation of the new NRQCD singularity is the main conceptual claim of the paper, but the demonstration is implicit. Please show explicitly that, after applying the matrix Z_NRQCD defined in Eqs. (5.17)-(5.20), the pole part of Z_IR^{-1} F_{p,A_i} in the gg channels is exactly reproduced by the Coulomb and cross anomalous dimensions of Eqs. (5.7), (5.12), (5.14) and (5.15), and that no 1/epsilon terms remain in F^fin at order alpha_s^{q+2} for each colour structure. This check is the minimal evidence that the new pole has been correctly assigned to the |3S_1[8] g> Fock-state contribution rather than subtracted by an accidental choice of Z_NRQCD.","section":"§5, Eqs. (5.6)-(5.20)"},{"comment":"The agreement with ref. [38] is the only independent two-loop check of a P-wave channel, but it is reported only as 'overall agreement' without quantitative detail. Since the new gg and gamma-g channels have no external verification, this cross-check carries substantial weight. Please provide a side-by-side numerical table of the three helicity amplitudes with the results of ref. [38], or state exactly how many digits agree for each amplitude and, where relevant, for each colour contribution. This is necessary for readers to verify the claimed cross-check and to assess the precision statement made in Section 7.","section":"§6.1 and §6.2"}],"minor_comments":[{"comment":"The term Gfin,(1) D(1)_muLambda in Eq. (6.3) is introduced after the finite remainder Ffin,(1)_{3S1[8],A4} has been defined in Appendix D; please state explicitly at which scale Ffin,(1)_{3S1[8],A4} is to be evaluated and how the one-loop scale dependence is transferred to the D(1)_muLambda term.","section":"§6, Eq. (6.3)"},{"comment":"The 20-digit numerical bases in Eqs. (6.18)-(6.68) would be easier to use if the inline identifications such as a(2)_1 = a(2)_{γγ,[1],A1;FF} were collected into a table that also gives the corresponding variables in the electronic file, since the reader may otherwise need to search through several paragraphs to assemble the mapping.","section":"§6.1"},{"comment":"The speculative sentence about D-wave states ('This may apply presumably also to C-even quarkonia with higher L states...') is not needed for the present results and should be removed or explicitly labelled as a conjecture.","section":"§2, footnote 9"},{"comment":"The statement that the light-by-light contribution proportional to n_h in ref. [38] 'has uncertainties in the last two digits' would be more useful if it identified which amplitude and coefficient this refers to and gave the corresponding numbers from both calculations.","section":"§6.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically strong and within the scope of JHEP. The main risk is the undocumented closure of the IBP reduction to the 76 master integrals of the companion paper; this is checkable and, if supplied, would remove the main structural uncertainty. I would also ask for a quantitative comparison with ref. [38] and an explicit demonstration of the pole cancellation in the gg channels before accepting."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper delivers the first analytic two-loop form factors for the P-wave channels and, more importantly, identifies a new NRQCD pole that requires the gg -> 3S1[8] form factor. That singularity is a real finding, not a technical footnote, and the accompanying d-dimensional helicity decomposition handles it carefully. The agreement with the earlier numerical gamma-gamma result and with the known one-loop 3S1[8] amplitude gives me reasonable confidence that the calculation is correct.\n\nThe strong points: analytic results for four new form factors, a clear explanation of the extra pole via the higher Fock state, no fitted parameters, and a clean separation of colour and flavour structures. The ancillary files with the analytical expressions are a plus.\n\nThe soft spot is real but checkable. Footnote 9 asserts that after the q-derivative in eq. (2.6), all two-loop integrals reduce to the same 76 master integrals as the 1S0 case, and this is a load-bearing premise. The cross-check against ref. [38] does exercise the reduction for the gamma-gamma channel, so the claim is not baseless, but the new gg and gamma-g channels have no independent check of that closure. A referee should ask for either a reduction table, a count of masters before and after, or an independent code run. This is a gap in documentation, not evidence of an error.\n\nMinor point: the conclusions advertise 1000-digit numerics while the tables show 20 digits; presumably the electronic files have more, but the text should be precise.\n\nWho gets value from this: quarkonium phenomenologists wanting NNLO hard functions, and amplitude people interested in NRQCD factorization subtleties. The paper deserves a serious referee; the IBP closure question is exactly what the referee should push on.","headline":"A genuinely new two-loop calculation for P-wave quarkonium with a real structural insight; the main soft spot is an asserted but undemonstrated IBP closure that partial cross-checks mitigate.","tokens_in":72535,"tokens_out":2292,"would_cite":true,"duration_ms":27122,"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":"Analytic two-loop form factors for P-wave quarkonium are computed for the first time in three of the four channels, and the fourth is upgraded from numerical to fully analytic.","keywords":["two-loop form factors","P-wave quarkonium","NRQCD factorisation","helicity amplitudes","master integrals","Coulomb singularity","3S1[8] form factor","chi_QJ production and decay"],"falsifier":"Independently IBP-reduce the P-wave two-loop integral families and check for any additional master integral at any order in $\\epsilon$, or recompute one helicity amplitude, for example the two-loop $A_3$ coefficient of $\\gamma\\gamma \\leftrightarrow {}^3P_0^{[1]}$ with value $-9.1475107737851465282$, by an independent numerical method and compare to the 20-digit value.","tokens_in":71536,"feed_emoji":"⚛️","tokens_out":7927,"duration_ms":72688,"temperature":0.7,"pith_summary":"This paper presents analytic two-loop form factors for the production and decay of spin-triplet P-wave quarkonium states, the $\\chi_{Q,J}$ family. The author computes the helicity amplitudes for $\\gamma\\gamma \\leftrightarrow {}^3P_J^{[1]}$, $gg \\leftrightarrow {}^3P_J^{[1]}$, $\\gamma g \\leftrightarrow {}^3P_J^{[8]}$ and $gg \\leftrightarrow {}^3P_J^{[8]}$, where $[1]$ and $[8]$ denote colour-singlet and colour-octet; the last three channels were previously unknown and the first was known only numerically. The two-loop results are expressed through the same 76 master integrals used in the S-wave companion calculation and are evaluated to 20-digit accuracy, with hard functions that can feed NNLO predictions for charmonium and bottomonium. The paper also uncovers a new NRQCD infrared pole in the $gg$ channel whose cancellation requires including the $gg \\leftrightarrow {}^3S_1^{[8]}$ form factor, restoring factorisation through the higher Fock state $|{}^3S_1^{[8]} g\\rangle$.","feed_headline":"Two-loop P-wave quarkonium form factors computed analytically","feed_subtitle":"Covers gamma-gamma, gg, and gamma-g channels; a new pole is cured by the 3S1[8] form factor.","key_machinery":"The load-bearing object is the set of 76 two-loop master integrals shared with the pseudo-scalar S-wave calculation, combined with a helicity projection onto the three Lorentz structures $T_1,T_2,T_3$ that define the P-wave amplitude. The projection is done in $d$ dimensions, which matters because the $c_1$ coefficient vanishes at tree level but carries a two-loop pole, so a four-dimensional reduction would miss finite pieces in the $A_3$ helicity amplitude. The second key element is the extended NRQCD subtraction factor $Z_{\\rm NRQCD} = Z_{\\rm Coul.} + Z_{{}^3S_1^{[8]}}$, a matrix in helicity space whose off-diagonal entries link the $A_2$ and $A_3$ helicity amplitudes of $gg \\leftrightarrow {}^3P_J$ to the one-loop $gg \\leftrightarrow {}^3S_1^{[8]}$ amplitude $A_4$. The new cross anomalous dimensions $\\gamma_{{}^3P_J^{[1]} {}^3S_1^{[8]}} = -C_A/3$ and $\\gamma_{{}^3P_J^{[8]} {}^3S_1^{[8]}} = -C_A/6$ encode the scale dependence $\\mu_\\Lambda$ of this mixing.","core_discovery":"On its own terms, the paper claims to complete the two-loop short-distance input for $\\chi_{Q,J}$ phenomenology: analytic results for all three helicity form factors of $\\gamma\\gamma \\leftrightarrow {}^3P_J^{[1]}$, validated against the earlier numerical computation, and first-time results for $gg \\leftrightarrow {}^3P_J^{[1]}$, $\\gamma g \\leftrightarrow {}^3P_J^{[8]}$ and $gg \\leftrightarrow {}^3P_J^{[8]}$. It finds that the bare two-loop amplitudes reduce, via integration-by-parts identities, to the 76 master integrals already evaluated analytically in the companion paper, whose special functions include multiple polylogarithms and elliptic generalisations. After UV renormalisation and subtraction of standard QCD IR singularities, the remaining NRQCD pole is the expected Coulomb singularity in the $\\gamma\\gamma$ and $\\gamma g$ channels; in the $gg$ channels, however, an additional pole appears only in the helicity configurations with $J_z = 0$ and $\\lambda_1 = \\lambda_2 = \\pm 1$. The paper traces this pole to the transition ${}^3P_J^{[1,8]} \\to {}^3S_1^{[8]} + g$ with an ultra-soft gluon, and shows that including the one-loop $gg \\leftrightarrow {}^3S_1^{[8]}$ form factor, computed here with its $O(\\epsilon)$ term, cancels it and restores NRQCD factorisation.","pith_inferences":["The same pole-cancellation mechanism should recur in other C-even quarkonium channels with $L \\ge 1$ and gluon initial states, since it only uses the derivative expansion in the relative momentum and the intermediate ${}^3S_1^{[8]}$ state; D-wave channels present a natural place to look.","The $d$-dimensional helicity projection is not optional: the paper's own comparison with the earlier numerical gamma-gamma computation shows that the $d$ versus four-dimensional treatment only matters for $A_3$, so future implementations of these hard functions should keep the $\\epsilon$-dependent projectors.","The newly computed $O(\\epsilon)$ term of the $gg \\leftrightarrow {}^3S_1^{[8]}$ form factor, needed for the pole cancellation here, is likely needed in the same way for NNLO predictions of $J/\\psi$ and $\\Upsilon$ polarisation."],"forward_implications":["The $gg \\leftrightarrow {}^3P_J^{[1]}$ hard functions can now be combined with existing LDMEs and parton distributions to build the virtual NNLO contribution to inclusive $\\chi_{Q,J}$ hadroproduction.","The $\\gamma\\gamma \\leftrightarrow {}^3P_J^{[1]}$ analytic results upgrade the di-photon decay of $\\chi_{c0,2}$ from numerical to fully analytic NNLO, with 20-digit numerics.","The new one-loop $gg \\leftrightarrow {}^3S_1^{[8]}$ form factor supplies the other NNLO ingredient needed for $\\chi_{Q,J}$ production and decay in the S-wave channel that shares the same velocity power counting.","The colour-octet Coulomb anomalous dimensions, obtained here for the first time, provide the infrared boundary condition for octet P-wave production, relevant for $J/\\psi$ and $\\Upsilon$ phenomenology at higher orders.","The restored factorisation via the $|{}^3S_1^{[8]} g\\rangle$ Fock state gives a concrete mechanism for the previously unseen pole, so the finite remainders are scheme-consistent and scale-dependent in the expected way."],"supporting_citations":[{"why":"Supplies the 76 two-loop master integrals, computed analytically and numerically, that all P-wave amplitudes are reduced to.","marker":"[41]"},{"why":"Companion S-wave calculation providing the IBP setup, renormalisation constants, Z_IR construction, and shorthand notation reused throughout.","marker":"[32]"},{"why":"Previous numerical two-loop result for gamma gamma to ^3P_J[1] used as the cross-check for the new analytic computation.","marker":"[38]"},{"why":"Establishes the NRQCD factorisation framework that the form factors are built to serve.","marker":"[33]"},{"why":"Shows the Landau-Yang theorem can fail with antisymmetric colour structure, which is what allows gg to 3S1[8] not to vanish.","marker":"[76]"},{"why":"Earlier one-loop calculation of gg to 3S1[8] up to the finite term, extended here by the O(epsilon) term and full colour decomposition.","marker":"[84]"},{"why":"Independent earlier one-loop result for gg to 3S1[8], used as an additional cross-check in Appendix D.","marker":"[82]"},{"why":"Previous Coulomb anomalous dimensions for colour-singlet 3P_J states, reproduced by this calculation for the singlet channels.","marker":"[68]"},{"why":"First observation of binding-energy logarithms in the leptonic decay of P-wave quarkonia, attributed to the E1 transition to 3S1 plus a photon and identified as the origin of the new pole structure.","marker":"[72]"}],"fun_headline_variants":["Two-loop P-wave quarkonium form factors go analytic","First analytic two-loop results for chi_QJ production and decay","P-wave quarkonium two-loop form factors: new pole tamed","Analytic two-loop form factors complete for P-wave quarkonia","Two-loop P-wave quarkonium: all channels now analytic"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the claim that the P-wave two-loop integrals are spanned, after integration-by-parts reduction, by the same 76 master integrals already computed in the companion paper, a statement the paper makes as a confirmed observation rather than showing the full reduction.","fun_headline_variants_meta":{"raw":{"variants":["Two-loop P-wave quarkonium form factors go analytic","First analytic two-loop results for chi_QJ production and decay","P-wave quarkonium two-loop form factors: new pole tamed","Analytic two-loop form factors complete for P-wave quarkonia","Two-loop P-wave quarkonium: all channels now analytic"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000816,"raw_usage":{"total_tokens":3644,"prompt_tokens":1081,"completion_tokens":2563,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":697,"completion_tokens_details":{"reasoning_tokens":2479}},"tokens_in":697,"tokens_out":2563,"duration_ms":18975,"temperature":1.0,"reasoning_tokens":2479,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:39:08.203221+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Independently IBP-reduce the P-wave two-loop integral families and check for any additional master integral at any order in $\\epsilon$, or recompute one helicity amplitude, for example the two-loop $A_3$ coefficient of $\\gamma\\gamma \\leftrightarrow {}^3P_0^{[1]}$ with value $-9.1475107737851465282$, by an independent numerical method and compare to the 20-digit value.","supporting_citations":[{"cited_title":"A note on the fate of the Landau-Yang theorem in non-Abelian gauge theories","cited_arxiv_id":"1509.07853","evidence_quote":"Shows the Landau-Yang theorem can fail with antisymmetric colour structure, which is what allows gg to 3S1[8] not to vanish."},{"cited_title":"Breakdown of QCD Factorization for P-Wave Quarkonium Production at Low Transverse Momentum","cited_arxiv_id":"1405.3373","evidence_quote":"Earlier one-loop calculation of gg to 3S1[8] up to the finite term, extended here by the O(epsilon) term and full colour decomposition."},{"cited_title":"No Landau-Yang in QCD","cited_arxiv_id":"1508.07115","evidence_quote":"Independent earlier one-loop result for gg to 3S1[8], used as an additional cross-check in Appendix D."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First observation of binding-energy logarithms in the leptonic decay of P-wave quarkonia, attributed to the E1 transition to 3S1 plus a photon and identified as the origin of the new pole structure."}],"review_version":1}