{"id":"64b0e44e-46d5-49df-a569-e52713833f02","arxiv_id":"2506.17030","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A first-principles lattice QCD calculation gives Gamma(chi_c1 to J/psi gamma) = 0.3265(79) MeV and a2 = -0.0666(22), consistent with experiment and far more precise than prior lattice results.","lead":"Lattice QCD now predicts the chi_c1 to J/psi gamma decay width and the magnetic quadrupole amplitude, with both matching experiment. This is the first full QCD calculation of this radiative charmonium transition at physical quark masses with four lattice spacings.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central continuum values rest on the unverified neglect of disconnected Wick contractions; the paper computes only connected diagrams and provides no direct bound on the few-percent shifts that could alter Eq. (22).","rationale":"The manuscript is careful: four lattice spacings, OS regularization to avoid the twisted-mass 1++/1-- mixing (numerically demonstrated in Fig. 1), twisted boundary conditions to set q^2=0, and stability checks on the continuum extrapolation. I do not find an algebraic inconsistency in the extraction formulas or an obvious misestimate in the fits. The weakest point is the truncated contraction structure. Equations (14) and (15) define the quantities that are fitted, and both omit disconnected diagrams; the final continuum extrapolation is therefore an extrapolation of connected-only form factors. Since the paper's central claim is that E1 and M2/E1 are the physical QCD values at the quoted few-percent precision, the size of the omitted loops is not a peripheral issue. The two-point mass comparison gives some confidence, but it is not a direct bound on the three-point disconnected matrix elements, which are OZI-suppressed but could still be at the 1-3% level; the slight width excess in Eq. (23) is consistent with that possibility. This is exactly the limitation the authors flag in Section V, and it is not resolved in the paper. I would not reject; the connected calculation is a substantial improvement and the a2 agreement with experiment is reassuring. But the numerical claim should be conditional: either compute the disconnected diagrams on at least one ensemble, or phrase the quoted values as connected-only and enlarge the systematic error accordingly.","tokens_in":14805,"tokens_out":15059,"duration_ms":174923,"concrete_test":"On the finest physical ensemble (D96), compute the disconnected contributions with the same OS valence setup and kinematics |k| = (m_chi^2 - m_J/psi^2)/(2 m_chi), t_J ~ 1.7 fm: (i) annihilation diagrams in the J/psi and chi_c1 two-point functions of Eq. (15); (ii) three-point diagrams of Eq. (14) in which the electromagnetic current couples to a sea-quark loop (u,d,s,c), using stochastic all-to-all propagators and the same dilution scheme as the connected computation. Re-extract E1(t_chi;t_J) and M2(t_chi;t_J), form the plateaus, and apply the same fitting procedure. If the shifts in E1 and M2/E1 are below about 0.3 sigma of the current errors, the omission is benign and Eq. (22) stands; if they reach the 1-3% scale, the continuum results must be revised.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantities E1=0.879(11) and M2/E1=-0.0668(22) in Eq. (22) are extracted from three-point functions in Eq. (14) that include only the quark-connected Wick contractions; the disconnected diagrams in the three-point functions and in the two-point functions of Eq. (15) are set to zero. The physical QCD values are therefore recovered only if these omitted loops shift E1 and M2/E1 by less than the quoted 1% and 3% errors. The paper's evidence for smallness is indirect: the masses in Eq. (20) agree with experiment at the 0.3-0.4% level (Section III), which bounds disconnected two-point effects, and the generic Zweig/SU(3) suppression argument in Section V. This does not directly bound the disconnected three-point pieces, which have their own lattice-spacing dependence and enter the extracted amplitudes through the Z-factors in Eq. (18). The observed 1-2 sigma excess of Gamma in Eq. (23) over experiment is the size one would expect from a missing few-percent contribution. The authors acknowledge this limitation, but it is load-bearing for the claim that Eq. (22) are the physical QCD form factors.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a lattice QCD determination of the E1 and M2 form factors for the radiative transition chi_c1 -> J/psi gamma at q^2=0. The calculation uses four ETMC Nf=2+1+1 twisted-mass ensembles with lattice spacings from about 0.057 to 0.091 fm, with physical light, strange, and charm sea quark masses except for the coarsest ensemble (m_pi ~ 175 MeV). The charm valence mass is tuned through m_Ds, twisted boundary conditions impose the on-shell photon momentum, and the chi_c1 interpolating operator is built in the OS regularization to avoid the parity/C mixing that would otherwise contaminate the J^PC=1++ channel. Three-point functions are computed from quark-connected Wick contractions only; two values of t_J on one ensemble are used to assign a systematic error via Eq. (21). Linear-in-a^2 continuum extrapolations, with a stability check excluding the coarsest lattice and a BAIC average, give E1=0.879(11), M2/E1=-0.0668(22), and hence Gamma(chi_c1->J/psi gamma)=0.3265(79) MeV and a2=-0.0666(22). The results are compared with the existing Fermilab measurement, with the CLEO and BES-III angular analyses, and with earlier quenched and Nf=2 lattice determinations.","tokens_in":15010,"tokens_out":17220,"duration_ms":188026,"significance":"If correct, this is the first unquenched, physical-quark-mass lattice determination of both transition form factors for this decay, and the quoted M2/E1 accuracy improves on the only existing quenched lattice result by about a factor of 30. The agreement of a2 with experiment is a nontrivial check of the lattice approach. The paper has several concrete strengths: four lattice spacings with good control of the continuum extrapolation, an explicit OS-regularization solution to the spurious-mixing problem, the use of twisted boundary conditions to reach q^2=0, an explicit treatment of the t_J systematic through Eq. (21), a BAIC-based fit averaging, and transparent reporting of the connected-only approximation. The main caveat is that disconnected Wick contractions are not computed; the authors argue they are OZI/SU(3)-suppressed and support this with the good continuum-limit agreement of the charmonium masses in Eq. (20) with experiment, but no direct numerical bound is provided. This is an acknowledged limitation that should be kept in mind when quoting Eqs. (22)-(24), but I do not regard it as invalidating the central results.","major_comments":[],"minor_comments":[{"comment":"The phrase 'first full QCD computation' should be qualified, since only quark-connected diagrams are evaluated and the coarsest ensemble has m_pi ~ 175 MeV rather than the physical pion mass. I suggest rewording to something like 'first computation with Nf=2+1+1 dynamical quarks at physical light-quark masses for the quark-connected contribution' or an equivalent explicit definition of what 'full QCD' means here.","section":"Section V"},{"comment":"It is not completely clear how the systematic uncertainties Delta E1 and Delta M2 defined in Eq. (21) from the single B64 ensemble are propagated into the final errors. The text says they are included before the continuum extrapolation, but it does not state whether they are added in quadrature to the statistical errors in Table II or applied as a global shift to all ensembles. Please make this propagation explicit.","section":"Section II and Table II"},{"comment":"There is a typographical error in the J/psi two-point function: 'e^{-E_J/Psi (T-t)}' should be e^{-E_{J/psi}(T-t)}.","section":"Eq. (16)"},{"comment":"In Fig. 4, the data points at t_J ~ 2.4 fm are 'slightly shifted horizontally for easier comparison', but the shift is not specified in the caption. A brief statement of the shift would make the comparison more transparent.","section":"Section III and Fig. 4"}],"recommendation":"minor_revision","confidential_remarks":"The paper is technically solid and the central results are clearly presented. The only substantive concern is the disconnected-diagram caveat, which the authors themselves acknowledge; in my view this warrants a wording change rather than a new computation. I would not require a direct calculation of the disconnected contributions for acceptance, but the final wording should avoid overclaiming 'full QCD' without qualification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers the first Nf=2+1+1 physical-mass, four-spacing lattice determination of E1 and M2 for chi_c1 -> J/psi gamma, with a clean solution to the twisted-mass mixing problem via the OS valence regularization. The quoted continuum values E1=0.879(11), M2/E1=-0.0668(22) are credible, and a2 improves on the only prior quenched result by about a factor of thirty.\n\nThe computational setup is careful: the charm mass is tuned via m_Ds, twisted boundary conditions put the photon exactly on shell, two tJ values are compared and their spread folded in via Eq. (21), the continuum extrapolation is a linear a^2 fit with a stability check omitting the coarsest lattice, and the BAIC average is used. The chi^2/dof values are small, and the continuum masses agree with experiment at the 0.3-0.4% level. The paper is also honest in its comparison with experiment, flagging the ambiguity from the PDG chi_c1 width and the fact that only E835 has measured this channel directly.\n\nThe main soft spot is exactly what the authors acknowledge: the disconnected Wick contractions are set to zero in both the two- and three-point functions, and the evidence that they are negligible is indirect. The mass agreement bounds disconnected effects in the two-point functions, and the Zweig/SU(3) suppression argument is plausible, but there is no direct bound on the disconnected three-point pieces, which enter the extracted form factors through the Z-factors in Eq. (18). The 1-2 sigma excess of the predicted width over experiment is about the size one would expect from a few-percent missing contribution, so this caveat is not cosmetic; it is the main limitation on the claim that Eq. (22) are the physical QCD form factors. A secondary minor point: M2 shows 10-15% cutoff effects on the coarsest lattice, so the M2/E1 continuum value leans more heavily on the finer spacings, though the stability check helps.\n\nNo data or code is released, which is typical for the field, but it means the quoted errors cannot be re-derived from the text.\n\nThis paper deserves a serious referee. The analysis is thorough, the result is a real step forward in precision charmonium phenomenology, and the authors have clearly identified the main caveat. I would recommend sending it to peer review, with a referee who can scrutinize the OS regularization argument and the continuum extrapolation. The disconnected contribution should be quantified in a follow-up, but that does not block publication of this result.","headline":"First physical-mass multi-spacing lattice QCD determination of the chi_c1 -> J/psi gamma form factors; the central numbers are credible, and the main caveat is the unquantified size of the omitted disconnected contractions.","tokens_in":15659,"tokens_out":1953,"would_cite":true,"duration_ms":21293,"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":"Lattice QCD fixes the chi_c1 to J/psi gamma decay width at 2.4 percent.","keywords":["lattice QCD","charmonium","radiative decay","transition form factors","chi_c1 to J/psi gamma","twisted mass fermions","magnetic quadrupole amplitude","continuum extrapolation"],"falsifier":"Evaluate the disconnected three-point diagrams on the same four ensembles and compare the resulting $E_1$ and $M_2/E_1$ with the connected-only values; a shift in $E_1$ larger than about 0.01 or in $M_2/E_1$ larger than about 0.003 would mean the quoted continuum results are not the full QCD answer. A second decisive check is a new high-precision measurement of the $\\chi_{c1}\\to J/\\psi\\,\\gamma$ width, which would settle whether the one-to-two-$\\sigma$ excess over the current average is real.","tokens_in":14577,"feed_emoji":"⚛️","tokens_out":13382,"duration_ms":120481,"temperature":0.7,"pith_summary":"This paper computes, directly from lattice QCD, the two form factors that control the radiative decay $\\chi_{c1}\\to J/\\psi\\,\\gamma$: the electric dipole $E_1$ and the magnetic quadrupole $M_2$. Using gauge-field ensembles with $N_f=2+1+1$ dynamical quarks at four lattice spacings, with physical charm and strange quarks and physical light quarks except on the coarsest lattice, the continuum extrapolation gives $E_1=0.879(11)$ and $M_2/E_1=-0.0668(22)$. These imply a decay width $\\Gamma(\\chi_{c1}\\to J/\\psi\\,\\gamma)=0.3265(79)$~MeV and a magnetic quadrupole fraction $a_2=-0.0666(22)$. The width agrees with experiment at the one-to-two-$\\sigma$ level, and $a_2$ agrees with the measured angular-analysis value while improving the only previous lattice determination by roughly a factor of thirty. If correct, this is the first full QCD calculation of both form factors at physical quark masses, and it gives a high-precision check of QCD in charmonium transitions.","feed_headline":"Lattice QCD fixes the chi_c1 to J/psi gamma width at 2.4 percent","feed_subtitle":"Continuum-extrapolated form factors agree with experiment and pin the quadrupole amplitude to -0.0666(22).","key_machinery":"The argument rests on a decomposition of the transition matrix element $\\langle J/\\psi(k,\\varepsilon)|J^\\mu_{\\rm em}|\\chi_{c1}(p,\\eta)\\rangle$ into dimensionless form factors $E_1(q^2)$, $M_2(q^2)$, and a $C_1(q^2)$ that does not contribute to the physical amplitude. The photon is put on shell by giving the $J/\\psi$ a three-momentum $|\\mathbf{k}|\\simeq389.4$~MeV through twisted boundary conditions on one charm propagator. A second load-bearing element is the use of the OS regularization, a valence-quark discretization that preserves exact charge conjugation, which forbids the $\\chi_{c1}$ interpolating operator from mixing with the lighter $J/\\psi$ ($1^{--}$) state and with $1^{+-}$ states. The form factors are extracted from long plateaus of three-point correlation functions at two source-sink separations, and the continuum limit is taken with a linear fit in $a^2$, with the coarsest-lattice-excluded fit used to set the extrapolation systematic.","core_discovery":"The central claim is that the continuum limit of lattice QCD, with $N_f=2+1+1$ dynamical quarks at physical masses, determines the on-shell transition form factors $E_1(0)$ and $M_2(0)$ for $\\chi_{c1}\\to J/\\psi\\,\\gamma$ to be $E_1=0.879(11)$ and $M_2/E_1=-0.0668(22)$. From these the decay width is $\\Gamma(\\chi_{c1}\\to J/\\psi\\,\\gamma)=0.3265(79)$~MeV and the magnetic quadrupole fractional amplitude is $a_2=-0.0666(22)$. The paper presents this as the first full QCD computation of both form factors at physical quark masses: it matches the experimental $a_2$ and is compatible with the measured width at one to two $\\sigma$, while differing from the only previous unquenched lattice result, a discrepancy the authors trace to the earlier work's interpolating operator overlapping spuriously with the $J/\\psi$.","pith_inferences":["The close agreement of $a_2$ with experiment while the width runs one to two sigma high suggests that any missing contribution would have to affect $E_1$ more strongly than the $M_2/E_1$ ratio; a targeted computation of the charm-quark disconnected contribution to the electric form factor would test this directly.","The same ensembles and valence regularization could be used to predict the radiative transition $\\Upsilon\\to\\eta_b\\gamma$, whose width is not known to comparable precision, giving a few-percent first-principles target for bottomonium physics.","A new independent measurement of the branching fraction in a channel less entangled with the $\\chi_{c1}$ total width would discriminate between the two experimental width assignments currently used to interpret the one-to-two-sigma tension."],"forward_implications":["If the central values are right, $\\Gamma(\\chi_{c1}\\to J/\\psi\\,\\gamma)=0.3265(79)$~MeV becomes the first-principles reference for this charmonium transition, and the one-to-two-sigma excess over the experimental average marks a real tension that a new width measurement could sharpen or resolve.","The ratio $M_2/E_1=-0.0668(22)$ gives $a_2=-0.0666(22)$, a precision lattice result competitive with the best angular analyses; it confirms that the magnetic quadrupole component is small and negative.","The paper's way of isolating the $1^{++}$ state, by choosing the OS valence regularization so that charge conjugation stays exact, should carry over directly to other charmonium radiative transitions computed with twisted-mass ensembles.","Because the quoted uncertainties are 1.2% on $E_1$ and 3% on $M_2/E_1$, the next decisive step is a direct evaluation of the omitted disconnected diagrams; if they fall inside the quoted errors, these results are the complete QCD answer."],"supporting_citations":[{"why":"Supplies the experimental branching fraction, meson masses, width averages, and measured quadrupole amplitude used for the comparisons in Section IV.","marker":"[7]"},{"why":"Provides the decomposition of the transition matrix element into E1, M2, and C1 that the calculation adopts, along with the first lattice study of this decay.","marker":"[10]"},{"why":"Gives the quenched single-lattice result for the decay width this work improves on and compares against.","marker":"[11]"},{"why":"Reports the only previous unquenched lattice result; the paper argues its much smaller width stems from spurious J/psi mixing in the twisted-mass interpolator.","marker":"[12]"},{"why":"Introduces the mixed-action valence regularization whose exact charge-conjugation symmetry suppresses the unwanted state mixing.","marker":"[14]"},{"why":"Documents the Nf=2+1+1 ensembles, lattice spacings, and vector-current renormalization constants used for the calculation.","marker":"[15]"},{"why":"Establishes the companion computation for h_c and h_b radiative decays, supplying the smearing, plateau-analysis, and continuum-fit averaging procedures reused here.","marker":"[18]"}],"fun_headline_variants":["Lattice QCD fixes chi_c1 to J/psi gamma width at 2.4 percent","Lattice QCD pins chi_c1 radiative width to 0.3265(79) MeV","Continuum lattice QCD reproduces chi_c1 decay width","Physical-mass lattice QCD computes chi_c1 width and quadrupole","QCD lattice yields chi_c1 width with 2.4% precision"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation includes only the dominant connected Feynman diagram, and the sea-quark pair-creation ('disconnected') loops are set to zero rather than computed, so if those omitted loops shift $E_1$ or $M_2/E_1$ by more than about one to three percent, the quoted width and $a_2$ change by more than the stated errors.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD fixes chi_c1 to J/psi gamma width at 2.4 percent","Lattice QCD pins chi_c1 radiative width to 0.3265(79) MeV","Continuum lattice QCD reproduces chi_c1 decay width","Physical-mass lattice QCD computes chi_c1 width and quadrupole","QCD lattice yields chi_c1 width with 2.4% precision"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000751,"raw_usage":{"total_tokens":3393,"prompt_tokens":1047,"completion_tokens":2346,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":663,"completion_tokens_details":{"reasoning_tokens":2238}},"tokens_in":663,"tokens_out":2346,"duration_ms":20088,"temperature":1.0,"reasoning_tokens":2238,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:13:13.342561+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the disconnected three-point diagrams on the same four ensembles and compare the resulting $E_1$ and $M_2/E_1$ with the connected-only values; a shift in $E_1$ larger than about 0.01 or in $M_2/E_1$ larger than about 0.003 would mean the quoted continuum results are not the full QCD answer. A second decisive check is a new high-precision measurement of the $\\chi_{c1}\\to J/\\psi\\,\\gamma$ width, which would settle whether the one-to-two-$\\sigma$ excess over the current average is real.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the only previous unquenched lattice result; the paper argues its much smaller width stems from spurious J/psi mixing in the twisted-mass interpolator."}],"review_version":2}