Pith. sign in

REVIEW 3 major objections 5 minor 19 references

The hadronic decay of vector charmonium

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read The paper extracts the hadronic decay width of ψ(3770) to a D-meson pair from lattice QCD using a narrow-width ratio method, obtaining 24.2(6.4) MeV, compatible with experiment.

desk verdict The ratio-method application is promising, but the printed decay-width formula contradicts the paper's own table, so the headline agreement with experiment is currently unverifiable. read the letter →

arxiv 2412.15915 v1 pith:HFSJN67X submitted 2024-12-20 hep-lat hep-ph

classification hep-lathep-ph
keywords latticeQCDhadronicdecaycharmoniumpsi(3770)narrow-widthapproximationfinitevolumeratiomethod3P0quarkmodel
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that a narrow-width ratio method, computationally cheaper than standard multi-volume scattering analyses, can extract the hadronic decay width of an excited quarkonium state directly from lattice QCD. The key step is to isolate the finite-volume mixing $\langle \bar D D | \psi(3770)\rangle$ from correlator ratios, then convert it to the decay width $\Gamma(\psi(3770)\to D\bar D)$ through Fermi's golden rule. On two ensembles at $m_\pi\sim 440\,\mathrm{MeV}$ the extracted width is compatible with experiment, with the larger volume giving $\Gamma=24.2(6.4)\,\mathrm{MeV}$. The paper further shows that the ${}^3P_0$ quark model reproduces the momentum dependence of the extracted mixing, suggesting the quark-rearrangement mechanism is qualitatively correct. If correct, the method gives a practical alternative for decays where the narrow-width approximation holds.

What carries the argument

The central object is the hadronic mixing element $x_{31}\equiv \langle \bar D D | \psi(3770)\rangle$ in a finite box. It is isolated from the long-time behavior of the ratio $R(t)=|\bar T_3(t)|/\sqrt{\bar P^\psi_{33}(t) P_{DD}(t)}$, whose spectral decomposition at large $t$ is $|x_{31}|/\big(\Delta\sinh(t\Delta)+A e^{-t\Delta}\big)$, where $\Delta=(m_\psi-E_{DD})/2$ is half the energy gap. The conversion to the decay width assumes non-interacting plane-wave $D$ mesons and uses Fermi's golden rule, $\Gamma = \frac{L^3}{24\pi p_i m_\psi}|x_{31}|^2$, with $p_i^2=m_\psi^2/4-m_D^2$. Partially twisted boundary conditions on the charm quarks set the $D$-meson momenta so the on-shell condition $m_\psi=2E_D$ can be reached, and the $\bar D D$ operator is projected to a $p$-wave.

What would settle it

Compute the width on a third lattice volume at the same pion mass, keeping the on-shell momentum fixed: if the width from Eq (6) changes with $L$ beyond statistical error, the assumed $L^3$ compensation is wrong. A complementary test is a full finite-volume scattering analysis on the same ensembles, which would give an independent width to compare against.

Watch

Extended reading notes

Core claim

The central claim is that the finite-volume hadronic mixing between the $\psi(3770)$ charmonium state and a $D\bar D$ pair in a $p$-wave, extracted from the ratio $R(t)$ of Eq (4), is the same quantity that enters Fermi's golden rule for the physical decay, so the decay width can be written as $\Gamma = \frac{L^3}{24\pi p_i m_\psi} |x_{31}|^2$ at the on-shell momentum. Inserting the lattice-determined mixing, the paper obtains $\Gamma(\psi(3770)\to D\bar D) = 24.2(6.4)\,\mathrm{MeV}$ on the larger ensemble, consistent with the experimental value. This is presented as evidence that the narrow-width ratio method is a viable alternative to standard finite-volume scattering analyses for this decay. The paper also extracts the mixing for off-shell momenta and compares the resulting width function with the ${}^3P_0$ quark-model prediction, finding qualitative agreement, particularly on the larger volume.

Load-bearing premise

The method assumes that the finite-volume mixing extracted with non-interacting plane-wave $D$ mesons, divided by the $L^3$ factor in the golden-rule formula, is exactly the same coupling that drives the decay in infinite volume; the paper explicitly notes it never establishes this connection.

Editorial extensions

If this is right

  • The narrow-width ratio method offers a computationally cheaper path to hadronic decay widths of narrow excited states, avoiding the need for several lattice volumes and group irreps required by standard scattering analyses.
  • The extraction is not restricted to on-shell kinematics: the lattice provides the mixing as a function of $D$-meson momentum, enabling direct comparison with quark-model predictions over a range of momenta.
  • Because the method works at unphysical pion mass ($m_\pi\sim 440\,\mathrm{MeV}$), a continuum extrapolation and a move to physical quark masses would test whether the agreement with experiment persists away from this single point.
  • The ${}^3P_0$ quark model, with parameters fixed by the lattice spectroscopy, reproduces the momentum dependence of the mixing, supporting the quark-pair-creation picture as a useful effective description.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The method's validity hinges on the narrow-width approximation; for broader resonances or decays further from threshold, the assumed two-level system would be a poorer approximation, and the extracted width would need to be cross-checked against a full finite-volume analysis.
  • A natural next test is to apply the same ratio method to a decay with a known but wider resonance to see how the quality of the width degrades as the width-to-mass gap grows.
  • The explicit caveat that no infinite-volume connection is established suggests a possible follow-up derivation: a rigorous relation between this ratio-extracted mixing and the infinite-volume coupling, analogous to the quantization condition for two-particle energies, would place the method on firmer theoretical ground.
  • If the $L^3$ compensation is exact, the ratio method could also be used to extract decay couplings for states above multiple thresholds, where standard single-channel finite-volume analyses become complicated.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The proceedings paper applies the Michael–Pennanen ratio method to extract the finite-volume hadronic mixing x31 = <DD|psi(3770)> on two Nf=2 CLS ensembles at a single lattice spacing (a ~ 0.066 fm) and an unphysical pion mass (m_pi ~ 440 MeV). Using partially twisted boundary conditions to tune the D-meson pair to near on-shell kinematics, the authors fit the ratio R(t) of Eq. (4) to a parabola in the D-meson momentum, evaluate the on-shell point, and convert the mixing to a decay width with Eq. (6). The main result is Gamma(psi(3770) -> DD) = 24.2(6.4) MeV on the larger E5 ensemble, quoted as compatible with the experimental width, and the paper also compares the lattice mixing with the ^3P_0 quark model.

Significance. If the method is quantitatively sound, it offers a computationally cheaper alternative to Lüscher-type finite-volume analyses for narrow hadronic resonances, and the paper is explicit about the practical ingredients: GEVP spectroscopy, PTBC kinematics, and the ratio estimators of Eqs. (4) and (5). The strength of the paper is its concrete presentation of the lattice setup and of the Wick contractions, which makes the numerical procedure reproducible in principle. The explicit admission that no connection to the infinite-volume decay is established (Section 2) is an honest statement of the main limitation. However, the central numerical claim is currently undermined by an apparent inconsistency between the printed conversion formula and the tabulated decay width, and the absence of continuum, chiral, and systematic uncertainties prevents the quoted compatibility with experiment from being conclusive.

major comments (3)
  1. [Section 2, Eq. (6), and Table 4] The printed formula Gamma = L^3/(24 pi p_i m_psi) |x31|^2 is dimensionally inconsistent (|x31| has mass dimension one, so the right-hand side has dimension E^-3) and does not reproduce the quoted width. With the E5 numbers L = 10.67 GeV^-1, p = 0.468 GeV, m_psi = 3.946 GeV, and x31 = 0.0284 GeV, Eq. (6) gives about 7 MeV, while Table 4 quotes Gamma = 24.2(6.4) MeV; in lattice units the same formula gives Gamma/a ~ 0.19 against the tabulated Gamma/a = 0.0081(21). The tabulated value is instead reproduced by Gamma = L^3 p_i m_psi/(24 pi) |x31|^2. The authors must correct Eq. (6) or, if x31 carries an undocumented normalization, must state that normalization explicitly. Because this conversion is the only bridge between the lattice extraction and the physical width, the manuscript as written does not allow the central claim to be checked.
  2. [Section 2, after Eq. (6)] The paper states that 'at no point we establish an explicit connection to the infinite volume decay'. This is a load-bearing gap: the method assumes that the finite-volume mixing extracted with non-interacting plane-wave D mesons is the coupling entering Fermi's golden rule, and that the L^3 factor provides the correct compensation. Since the key result is the decay width, the authors should either provide a derivation of Eq. (6) in the context of their non-relativistic two-level system, or perform a cross-check (for instance, comparing with a Lüscher-type analysis at the same lattice parameters or at least demonstrating volume independence with more than the two ensembles used here). Without such a check, the compatibility with experiment could be coincidental.
  3. [Sections 3 and 4, Table 4] The quoted result is obtained at a single lattice spacing and a single unphysical pion mass, and only statistical errors are reported. The central claim of compatibility with experiment would require an estimate of systematic uncertainties from, at minimum, the GEVP fit ranges, the smearing choices, the twist-angle interpolation, finite-volume effects, and the dependence on m_pi and a. As it stands, the agreement of the E5 central value with experiment is suggestive but not yet a quantitative validation of the method.
minor comments (5)
  1. [Section 4, Figure 5] The ^3P_0 quark-model comparison fits the coupling beta to the same lattice data (as stated in the text), so the 'remarkable agreement' is a fitted reproduction rather than an independent prediction; the wording should make this explicit.
  2. [Equation (8)] The notation is inconsistent: the text defines p = sqrt(3) theta/L, while Eq. (8) writes sqrt(3) a theta_0/L; the factors of the lattice spacing should be made uniform.
  3. [Figure 1] The text describing the figure refers to 'red lines' while the caption says 'solid line'; the color and line-style descriptions should be aligned.
  4. [Abstract and Conclusions] The phrase 'compatible with experiment' should be qualified by noting that the comparison is made at an unphysical pion mass, one lattice spacing, and with statistical errors only, so the compatibility does not yet establish a precision prediction.
  5. [Table 2] The level-2 mass, m_cc = 3859(52) MeV, is close to the DD threshold and has a large uncertainty; a comment on how the level-3 identification with psi(3770) is justified would help the reader assess the GEVP analysis.

Circularity Check

1 steps flagged · score 3.0 of 10

Fitted 3P0 β is called a prediction; the central lattice extraction is independent, though Eq. (6) as printed does not yield the tabulated width.

  1. fitted input called prediction [Section 4 (Analysis), 3P0 quark-model paragraph (page 6); echoed in Conclusions (page 7)]
    "The remaining free parameter β, which gives the quark-pair creation strength, is fixed fitting the model to the lattice data, see [8] for more details. We observe that the quark model describes our data on a qualitative level, and the agreement is especially good on ensemble E5."

    The quark-model curve shown in Fig. 5 is not an independent prediction: its only free parameter β is fitted to the same lattice mixing data against which the model is then said to agree. Calling the resulting qualitative match a “prediction by the 3P0 quark model” (Conclusions) is therefore a fitted reproduction renamed as a prediction. The central lattice extraction itself remains self-contained because x31 is obtained from correlator data and the final width is benchmarked against experiment, so this is a secondary, partial circularity rather than a collapse of the main derivation.

full rationale

The central derivation chain—GEVP spectroscopy, the ratio R(t) of Eq. (4), extraction of the hadronic mixing x31, and the final width compared with experiment—is self-contained in the sense that x31 is obtained from lattice correlator data and not from the experimental width it is meant to predict. The ratio method is attributed to the independent earlier work [9] and [15], with the authors' companion paper [8] used mainly for details of this project; that self-citation is descriptive rather than load-bearing. The paper's explicit limitation, “at no point we establish an explicit connection to the infinite volume decay” (Section 2), is an acknowledged assumption about the finite-volume-to-infinite-volume conversion, not a tautology. The one genuine circular step is the 3P0 comparison: β is fixed by fitting the model to the same lattice data that the comparison then “confirms,” so the claimed qualitative agreement is a fitted reproduction. This affects only the interpretive part of the paper, not the lattice extraction itself. Separately, there is a correctness defect that is not circularity: Eq. (6) as printed, Γ = L^3/(24π p_i m_ψ) |x31|^2, has inverse-energy dimension and, with the Table 4 values on E5, gives roughly 0.19/a rather than the quoted 0.0081/a; the tabulated value instead matches Γ = L^3 p_i m_ψ/(24π) |x31|^2. The paper's admission that no infinite-volume connection is established compounds this concern, but neither fact makes the central result equivalent to its own input. Overall, the main extraction is not circular; the 3P0 “prediction” is partially circular because it is fit to the data it claims to reproduce.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The central lattice extraction uses no invented entities and has one fitted parameter pair (c1,c2) defining the momentum dependence of the mixing; the quark model comparison adds the fitted parameter beta. The key axioms are the narrow-width approximation and the unproven infinite-volume conversion, both acknowledged in the text.

free parameters (3)
  • c1 (parabola coefficient) = 0.01869(32) (D5), 0.01276(56) (E5), in lattice units
    Coefficient of the parabola fit to |x31|(p) in Eq (9); the on-shell value used in the decay width comes from this fit, so the final result depends on it.
  • c2 (parabola curvature) = c2/a = 1.169(44) (D5), 1.183(52) (E5)
    Parabola curvature in Eq (9); together with c1 it sets the on-shell hadronic mixing.
  • beta (3P0 quark-pair creation strength) = Not quoted in this proceedings; fixed by fitting the 3P0 model to the lattice data in [8]
    The quark-model prediction shown in Fig. 5 uses beta fitted to the same lattice data, so the comparison is not a parameter-free test.
assumptions (5)
  • domain assumption Narrow-width approximation holds for psi(3770).
    The method relies on this approximation; the paper states it is satisfied for this decay (Section 2).
  • domain assumption The final-state D D pair is non-interacting and described by plane waves with PBCs.
    Used in Eq (6) to convert the mixing to a decay width; interactions near threshold are neglected.
  • domain assumption Eq (4) and (5) give the correct spectral decomposition of the correlator ratios for alpha=3 and large t.
    Requires the GEVP to isolate psi(3770) and a single dominant state in the fit window; deviations are absorbed in the constant A.
  • domain assumption The finite-volume mixing x31 can be converted to the physical infinite-volume width via Fermi's golden rule, Eq (6), with an L^3 compensation.
    Explicitly flagged in the text as an unproven connection ('at no point we establish an explicit connection to the infinite volume decay').
  • domain assumption The third GEVP level on the D5 ensemble is the psi(3770) state.
    Level 3 mass is 3946(13) MeV at these unphysical masses; identification is based on the charmonium spectrum pattern.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The hadronic decay of vector charmonium." pith.science (2026). https://pith.science/paper/HFSJN67X

@misc{pith2026241215915,
  author       = {Pith},
  title        = {Pith review of: The hadronic decay of vector charmonium},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HFSJN67X}},
  note         = {Machine review of arXiv:2412.15915}
}
abstract

The extraction of decay parameters using lattice techniques is a computationally expensive task, requiring several volumes and group irreps to relate the spectrum on a lattice simulation to the infinite volume scattering. In this project we employ an alternative method based on a narrow-width approximation to extract the hadronic mixing $<\bar{D}D|\psi(3770)>$, which is needed to compute the decay $\Gamma(\psi(3770)\to\bar{D}D)$ between the second excited state of vector charmonium and a pair of $D$-mesons in a $p$-wave. We carry out our lattice simulations on two CLS ensembles at $m_\pi \sim 440~\text{MeV}$ and $a\sim 0.066~\text{fm}$ and obtain results compatible with experiment. Furthermore, we interpret our results analytically using the ${}^3P_0$ quark model.

Figures

Figures reproduced from arXiv: 2412.15915 by the authors.

Figure 1
Figure 1. The Wick contractions considered in this project. The charm quark appears as a solid line and the light-quarks as a dashed line. We do not consider diagrams with charm-quark annihilation. O(𝑎)-improved Wilson quarks [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Ratio R (𝑡) of equation (4) as a function of time, together with fit to the right-hand side (RHS) of equation (4). Ensemble D5 at twist angle 𝜃 = 1.5 rad (LEFT), and ensemble E5 at 𝜃 = 3 rad (RIGHT). These twist angles lie close to the on-shell condition on their respective datasets, see table 4. the correlator into two pieces, one that can be computed using PTBC, and another containing the momentum exchange that is… view at source ↗
Figure 3
Figure 3. Dispersion relation for the D-meson on ensembles D5 (LEFT) and E5 (RIGHT) as a function of the D-meson momentum modulus, 𝑝 ≡ |𝑝|. The lattice data is fully compatible with a non-relativistic particle in the continuum. Ensemble 𝜒 2 /dof 𝑎𝑐1 𝑐2 /𝑎 D5 34.39/39 0.01869(32) 1.169(44) E5 42.51/35 0.01276(56) 1.183(52) [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: The hadronic mixing between 𝜓(3770) and DD on ensembles D5 (LEFT) and E5 (RIGHT) as a function of the 3-momentum modulus of each 𝐷-meson. The darker (brighter) points are obtained fitting R (𝑡) (𝑥𝑇 (𝑡)) to equation (4) (equation (5)). We parametrize the darker points w…
Figure 5
Figure 5. Figure 5: The decay width Γ(𝜓(3770 → DD) as a function of the 3-momentum modulus of each 𝐷-meson. The slashed line shows the 3 𝑃0 quark-model prediction, and the full line the prediction given by our lattice determination. The decay may only occur on-shell, marked by the vertica…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 10 canonical work pages

  1. [1]

    Particle Data Groupcollaboration, Review of particle physics, Phys. Rev. D110 (2024) 030001

  2. [2]

    Lüscher M.,Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 1. Stable Particle States,Commun. Math. Phys.104 (1986) 177

  3. [3]

    Lüscher M.,Volume Dependence of the Energy Spectrum in Massive Quantum Field Theories. 2. Scattering States,Commun. Math. Phys.105(1986) 153

  4. [4]

    Lüscher Martin,Two particle states on a torus and their relation to the scattering matrix, Nucl. Phys. B354(1991) 531

  5. [5]

    Lüscher Martin,Signatures of unstable particles in finite volume, Nucl. Phys. B364 (1991) 237

  6. [6]

    C.B. Lang, L. Leskovec, D. Mohler and S. Prelovsek,Vector and scalar charmonium resonances with lattice QCD,JHEP 09 (2015) 089 [1503.05363]. 8 The hadronic decay of vector charmonium Teseo San José

  7. [7]

    Piemonte, S

    S. Piemonte, S. Collins, D. Mohler, M. Padmanath and S. Prelovsek,Charmonium resonances with𝐽𝑃𝐶 = 1−− and 3−− from ¯𝐷𝐷 scattering on the lattice,Phys. Rev. D100 (2019) 074505 [1905.03506]

  8. [8]

    Blossier, J

    B. Blossier, J. Heitger, J. Neuendorf and T. San José,Hadronic decay of vector charmonium from the lattice, 2411.10123

Show all 19 references
  1. [9]

    Pennanen and C

    P. Pennanen and C. Michael,String breaking in zero temperature lattice QCD, hep-lat/0001015

  2. [10]

    Chung, J

    S.U. Chung, J. Brose, R. Hackmann, E. Klempt, S. Spanier and C. Strassburger,Partial wave analysis in k-matrix formalism,Annalen der Physik507 (1995) 404 [https://onlinelibrary.wiley.com/doi/pdf/10.1002/andp.19955070504]

  3. [11]

    Hanhart, S

    C. Hanhart, S. Kürten, M. Reboud and D. van Dyk,Analysis of the𝜓(3770) resonance in line with unitarity and analyticity constraints,Eur. Phys. J. C84(2024) 483 [2312.00619]

  4. [12]

    Micu,Decay rates of meson resonances in a quark model,Nuclear Physics B10(1969) 521

    L. Micu,Decay rates of meson resonances in a quark model,Nuclear Physics B10(1969) 521

  5. [13]

    Le Yaouanc, L

    A. Le Yaouanc, L. Oliver, O. Pene and J.C. Raynal,Naive quark pair creation model of strong interaction vertices, Phys. Rev. D8 (1973) 2223

  6. [14]

    Weisskopf and E

    V. Weisskopf and E. Wigner,Berechnung der natürlichen Linienbreite auf Grund der Diracschen Lichttheorie, Zeitschrift für Physik63 (1930) 54

  7. [15]

    McNeile and C

    C. McNeile and C. Michael,Hadronic decay of a vector meson from the lattice, Phys. Lett. B 556 (2003) 177 [hep-lat/0212020]

  8. [16]

    Tayduganov,Electroweak radiative B-decays as a test of the Standard Model and beyond, theses, Université Paris Sud - Paris XI, Oct., 2011

    A. Tayduganov,Electroweak radiative B-decays as a test of the Standard Model and beyond, theses, Université Paris Sud - Paris XI, Oct., 2011

  9. [17]

    Fritzsch, F

    P. Fritzsch, F. Knechtli, B. Leder, M. Marinkovic, S. Schaefer, R. Sommer et al.,The strange quark mass and Lambda parameter of two flavor QCD, Nucl. Phys. B865 (2012) 397 [1205.5380]

  10. [18]

    Heitger, G.M

    J. Heitger, G.M. von Hippel, S. Schaefer and F. Virotta,Charm quark mass and D-meson decay constants from two-flavour lattice QCD,PoS LATTICE2013(2014) 475 [1312.7693]

  11. [19]

    Sachrajda and G

    C.T. Sachrajda and G. Villadoro,Twisted boundary conditions in lattice simulations, Phys. Lett. B609 (2005) 73 [hep-lat/0411033]. 9

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.