Pith. sign in

REVIEW 3 major objections 4 minor 3 cited by

First lattice QCD determination of the Standard Model neutrino background for J/ψ → γ + invisible yields 1.00(9)(7)×10^-10.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 08:04 UTC pith:A7OJBJO5

load-bearing objection First lattice QCD number for the J/psi neutrino floor; useful and mostly sound, but the precision claim needs shoring up and the abstract/main-text mismatch is sloppy. the 3 major comments →

arxiv 2601.18209 v2 pith:A7OJBJO5 submitted 2026-01-26 hep-lat hep-exhep-ph

Lattice determination of the neutrino background for J/psi rightarrow γ + textrm{invisible}

classification hep-lat hep-exhep-ph PACS 12.38.Gc13.20.Gd
keywords lattice QCDJ/psi radiative invisible decayneutrino backgroundbranching fractionform factorcharmoniumtwisted-mass fermionsdark matter searches
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper attempts to establish, from first principles, the Standard Model decay rate for J/ψ → γ ν ν̄, which is the irreducible neutrino background that contaminates searches for J/ψ → γ + invisible. Using lattice QCD with three lattice spacings and a method that captures the full photon-energy distribution, it obtains a branching fraction of 1.00(9)(7)×10^-10. A sympathetic reader would care because next-generation charm-factory experiments aim to reach exactly this sensitivity level; having a parameter-free, ab initio number to subtract is what turns any remaining events into evidence for invisible new particles. The calculation also provides a template for computing similar neutrino backgrounds in other quarkonium channels.

Core claim

The paper claims that the branching fraction for the Standard Model decay J/ψ → γ ν ν̄ is Br = 1.00(9)(7)×10^-10, where the first uncertainty is statistical (including lattice-spacing and continuum-extrapolation errors) and the second is an estimate of systematic effects. This is presented as the first lattice QCD calculation of this quantity, obtained by computing the hadronic matrix element ⟨0|T{J_em^μ(x) J_Z^ν(0)}|J/ψ⟩ nonperturbatively, parameterizing it with a single form factor, and integrating over the full phase space. The paper contrasts this ab initio result with an earlier phenomenological estimate of 0.7×10^-10 and argues that the lattice value is model-independent and precise en

What carries the argument

The central object is the hadronic tensor H^{μνα}(q,p) = ε^{μναβ} q_β F_{γνν̄}(q^2), which reduces all strong-interaction physics to a single form factor F. The authors use the scalar function method: they multiply H by ε^{μναβ} p_β, average over photon directions, and integrate against a spherical Bessel function, obtaining a scalar I(E_γ, Δt) from which F is extracted directly at every photon energy. This avoids the model-dependent interpolation of discrete lattice momenta used previously. The form factor feeds a Monte Carlo phase-space integral for the decay width; the dimensionless ratio Γ/f_J/ψ is then continuum-extrapolated in a^2 and rescaled with the lattice decay constant and the wo

Load-bearing premise

The load-bearing assumption is that the charm-quark vector-axial hadronic matrix element in the J/ψ → γ transition is insensitive to the light-quark sea mass at the ~7% level, so using N_f=2 ensembles with pion masses of 300–365 MeV — without any chiral extrapolation or corresponding systematic error — is safe; if the sea-quark effect is larger than that, the central value moves outside the quoted error budget.

What would settle it

A concrete check: compute the same decay width on a lattice ensemble with a physical pion mass (or with N_f=2+1 dynamical quarks) and see whether the central value moves by more than the quoted total uncertainty of ±0.16×10^-10. Alternatively, include the OZI-suppressed disconnected diagrams and compare; if either change exceeds the error budget, the central claim is falsified. A direct experimental measurement at a future charm-tau factory reaching 10^-10 sensitivity that excludes 1.0×10^-10 would also settle it.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Searches for invisible new particles in J/ψ → γ + invisible can now subtract a known Standard Model background, so any remaining excess is attributable to new physics.
  • Next-generation charm-tau facilities with roughly 10^5 times more J/ψ events, targeting 10^-10 sensitivity, will land on this neutrino floor; the lattice value tells them exactly where that floor is.
  • The same method transfers directly to Υ and φ radiative invisible decays, providing ab initio background estimates for those channels.
  • The result, 1.00×10^-10, is about 1.4 times larger than the earlier nonrelativistic color-singlet estimate of 0.7×10^-10, demonstrating that phenomenological models can be off by tens of percent.
  • The full-q^2 scalar function method removes a source of model dependence, so the quoted precision is directly testable by finer lattice calculations rather than depending on fitting assumptions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The quoted 7% systematic uncertainty does not include any estimate for the unphysical sea-quark mass, since the ensembles use pion masses of 300–365 MeV with no chiral extrapolation; if the charm axial matrix element varies by more than ~7% when extrapolated to the physical pion mass, the central value would shift outside the stated error budget.
  • Using the ratio Γ/f_J/ψ before continuum extrapolation likely cancels some common renormalization and discretization effects; a direct computation without this ratio would clarify how much of the precision comes from that cancellation.
  • A natural test of the claim is a repeat calculation with N_f=2+1 or physical-mass sea quarks, or one that includes the disconnected diagrams currently dropped via OZI suppression; either could move the central value beyond the quoted uncertainty.
  • If a next-generation experiment directly measures J/ψ → γ ν ν̄, the comparison with this lattice prediction would double as a low-energy test of the Standard Model's neutrino couplings in charmonium decays.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper presents the first lattice QCD calculation of the Standard Model decay J/ψ → γνν̄, which constitutes an irreducible neutrino background to searches for J/ψ → γ + invisible. Using three N_f=2 twisted-mass ensembles with lattice spacings a ≃ 0.067, 0.085, 0.098 fm and pion masses 300–365 MeV, the authors compute the hadronic function via a scalar-function method, extract Γ(J/ψ→γνν̄)/f_{J/ψ} with a two-state excited-state fit, extrapolate linearly in a², and rescale by the continuum J/ψ decay constant. The final result is Br = 1.00(9)(7)×10⁻¹⁰ (Eq. 17), with the systematic error estimated from dropping the coarsest ensemble. The paper argues that this ab initio prediction is needed for upcoming STCF and other searches.

Significance. If the central result is correct, this is a valuable first non-perturbative benchmark for a background process that will soon be experimentally relevant. The scalar-function method, the explicit finite-volume check via spatial truncation, and the multi-state treatment of excited states are genuine strengths, and the comparison with the phenomenological NRQCD estimate is useful. The calculation is independent of the experimental bound and does not use the target branching fraction as input. However, the reliability of the quoted precision rests on systematic assumptions that are not quantitatively supported in the manuscript, so the result cannot currently be taken as a robust 10% benchmark.

major comments (3)
  1. [Numerical Setup / Table I; Eq. (17)] The calculation uses only three N_f=2 ensembles with pion masses 300, 315, 365 MeV and performs no chiral extrapolation or sea-quark-mass systematics. The continuum extrapolation in a² is necessarily correlated with m_π because all three ensembles lie on a single trajectory in (a², m_π) space. The claim that the heavy-charm matrix element is insensitive to the light sea quark at the 7% level is unsupported. Since a 7% shift would move the central value outside the quoted systematic, this is a load-bearing gap. The authors should either provide an estimate of the sea-quark dependence (e.g., from a partially quenched comparison, heavy-quark scaling, or a dedicated ensemble at a different m_π) or explicitly enlarge the systematic error.
  2. [Numerical Results, Fig. 4 and Table II] The continuum extrapolation rests on three points, and the a85 point is poorly described by the linear a² fit. From Table II, the values are 1.852(44), 1.503(34), 1.371(25) for a67, a85, a98. A straight line through a67 and a98 predicts ≈1.59 at the a85 lattice spacing, about 2.6σ above the measured 1.503(34). This contradicts the statement that 'the fitting curves ... describe the lattice data well.' The systematic error is estimated only from the difference between the three-point fit and the two-point fit (1.00 vs 1.07), which does not capture the poor internal consistency of the three-point fit. The a² extrapolation needs a more careful treatment or a larger uncertainty, especially because the central value is directly read off this fit.
  3. [Numerical Results, final paragraph; Br formula (Eq. 16)] The neglect of disconnected diagrams is justified only by OZI suppression, with references to other charmonium studies, but no estimate of the size of this effect in the present observable is provided. The weak current in Eq. (6) contains light-quark contributions, and the J/ψ state can couple to disconnected light-quark loops; OZI suppression is qualitative. The same applies to the connected sea-quark effects noted above. Since the paper claims a parameter-free, precision benchmark, these unquantified systematics should be bounded by a numerical check or explicitly folded into the second error.
minor comments (4)
  1. [Abstract vs. Eq. (17)] The abstract quotes Br = 1.04(7)(8)×10⁻¹⁰, while the main text and conclusion quote 1.00(9)(7)×10⁻¹⁰ (Eq. 17). These are inconsistent in both central value and errors; the manuscript must use one consistent set of numbers.
  2. [Fig. 4 caption] The caption refers to 'Γ_{η_c γγ}/f_{J/ψ}', but the plotted quantity is Γ(J/ψ→γνν̄)/f_{J/ψ}. Please correct the typo.
  3. [Introduction] The phrase 'futural experiments' should be 'future experiments'.
  4. [Supplementary, Eq. (S7)] The definition of Z_{0i} involves the point-source operator O_i = Z_A ar c γ_i c; it would be clearer to state explicitly that the renormalization constant Z_A is the same as that used for the weak current, and how the statistical correlation between f_{J/ψ} and R_f is handled in the ratio rescaled to the physical width.

Circularity Check

0 steps flagged

No significant circularity: the J/psi->gamma+invisible branching fraction is an independent lattice calculation, with only non-load-bearing self-citations.

full rationale

The central result Br = 1.00(9)(7)e-10 is produced by a direct lattice computation: the hadronic tensor is measured via three-point functions, the form factor is extracted through the scalar-function construction in Eqs. (12)-(13), and the width is obtained by phase-space integration (Eq. 14) after continuum extrapolation of R_f and f_J/psi (Eqs. 15-16). The target branching fraction never appears as an input; the only external inputs are standard electroweak couplings, the PDG total width, and renormalization constants Z_V (from previous work) and Z_A (from ETMC/RI-MOM), each with independent determinations. The ratio R_f = Gamma/f_J/psi is a standard variance-reduction device: R_f is measured on the lattice, and f_J/psi is measured and separately cross-checked against the PDG and HPQCD values, so rescaling R_f by f_J/psi is not fitting the answer. Self-citations [22-27] supply the scalar-function method and finite-volume checks, but the paper includes a self-contained derivation and a direct R-truncation check (Fig. 6), so these citations are not load-bearing. No uniqueness theorem is imported from the authors, no fitted quantity is renamed as a prediction, and no quantity is defined in terms of the target. The abstract/main-text numerical mismatch (1.04(7)(8) vs 1.00(9)(7)) and the unquantified N_f=2 sea-quark-mass dependence are correctness risks, not circularity.

Axiom & Free-Parameter Ledger

1 free parameters · 4 axioms · 0 invented entities

The calculation rests on standard SM electroweak inputs (G_F, α, sin²θ_W, m_Z), the single-form-factor parameterization of Ref. [21], OZI-suppression of disconnected diagrams, and the unquantified assumption that N_f=2 ensembles at m_π=300–365 MeV are adequate. The valence charm quark mass is tuned to the physical J/ψ mass. No new particles or forces are introduced.

free parameters (1)
  • valence charm quark mass = tuned per ensemble so m_J/ψ matches the physical value (details in Ref. [22])
    The charm quark mass is tuned to reproduce the physical J/ψ mass on each ensemble. This is a standard renormalization condition, but it is an input that affects the matrix element and is not a prediction.
axioms (4)
  • domain assumption Single-form-factor parameterization H^{μνα} = ε^{μναβ} q_β F(q²) (Eq. 7)
    Taken from Ref. [21]; justified by gauge invariance and C-parity, but not derived in this paper. The lattice extraction projects onto this Lorentz structure and would miss any additional form factors.
  • domain assumption Disconnected diagrams are negligible (OZI suppression)
    Section 'Numerical Results': 'Neglected disconnected diagrams are known to yield only small contributions in charmonium system [36–39] due to OZI suppression.' No numerical bound is estimated in this work.
  • domain assumption Wick rotation convergence; η_c(q=0) contamination vanishes (Eqs. S4–S6)
    Supplementary Eqs. (S4)-(S6) argue that the η_c(0) intermediate state decouples from the projected hadronic function. This is load-bearing for extracting the Minkowski form factor from the Euclidean correlator.
  • domain assumption N_f=2 ensembles with m_π=300–365 MeV are sufficient; no chiral extrapolation needed
    Table I lists pion masses well above physical; the paper does not discuss pion-mass dependence or assign an uncertainty for this approximation.

pith-pipeline@v1.3.0-alltime-deepseek · 12020 in / 29658 out tokens · 310494 ms · 2026-08-03T08:04:50.519454+00:00 · methodology

0 comments
read the original abstract

Searching for dark matter is a primary goal of modern astronomy and particle physics. Invisible decays of heavy quarkonia are particularly promising for probing light dark matter, attracting broad interest due to their unique sensitivity. Experiments searching for radiative invisible decays of the $J/\psi$ have steadily improved upper limits, and upcoming facilities will push sensitivity further--making the precise determination and subtraction of the neutrino background indispensable. Here, we present the first lattice QCD calculation of the Standard Model decay $J/\psi \to \gamma\nu\bar{\nu}$, an irreducible background to $J/\psi \rightarrow \gamma + \textrm{invisible}$. Our result for the branching fraction is $\operatorname{Br}(J/\psi \rightarrow \gamma\nu\bar{\nu})=1.04(7)(8)\times 10^{-10}$, where the first uncertainty is statistical and the second is our systematic estimate. This work advances lattice-based determinations of neutrino backgrounds to quarkonium invisible decays, delivering an ab initio benchmark for $J/\psi \rightarrow \gamma + \textrm{invisible}$. Our approach generalizes to other quarkonium channels (e.g., $\Upsilon/\phi \rightarrow \gamma+\textrm{invisible}$) and provides critical theoretical support for dark matter searches at colliders.

Figures

Figures reproduced from arXiv: 2601.18209 by Chuan Liu, Haobo Yan, Ke-Long Zhang, Ning Li, Xue-Ze Zhang, Yu Meng.

Figure 1
Figure 1. Figure 1: FIG. 1. An example of the feynman diagram for the decay [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The lattice results of Γ [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The lattice results of Γ [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Lattice values of Γ [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Lattice results of [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. For ensemble a67, Γ [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Probing the dark axion portal via $J/\psi$ decays at BESIII and STCF

    hep-ph 2026-06 unverdicted novelty 4.0

    Monte Carlo studies show BESIII data already reaches previously unexplored regions of the dark axion portal via J/ψ → a γ' decays, with STCF projected to improve sensitivity by about an order of magnitude.

  2. Searching for dark photons in $J/\psi$ decays

    hep-ph 2026-04 unverdicted novelty 4.0

    NRQCD-based study projects BESIII and STCF sensitivity to light dark photons (m_U < 3 GeV) via visible and invisible J/ψ decay channels.

  3. Searching for dark photons in $J/\psi$ decays

    hep-ph 2026-04 unverdicted novelty 3.0

    NRQCD calculations predict low event yields (0-172) and poor significances (10^{-6} to 10^{-1}) for dark photon mediated J/ψ decays at BESIII, with four-body invisible channels offering relatively better prospects bel...

Reference graph

Works this paper leans on

40 extracted references · 26 linked inside Pith · cited by 2 Pith papers

  1. [1]

    Bertone, D

    G. Bertone, D. Hooper, and J. Silk, Phys. Rept.405, 279 (2005), arXiv:hep-ph/0404175

  2. [2]

    Arkani-Hamed, D

    N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, Phys. Rev. D79, 015014 (2009), arXiv:0810.0713 [hep-ph]

  3. [3]

    Cadamuro and J

    D. Cadamuro and J. Redondo, JCAP02, 032 (2012), arXiv:1110.2895 [hep-ph]

  4. [4]

    P. W. Graham, I. G. Irastorza, S. K. Lamoreaux, A. Lind- ner, and K. A. van Bibber, Ann. Rev. Nucl. Part. Sci. 65, 485 (2015), arXiv:1602.00039 [hep-ex]

  5. [5]

    Bauer, M

    M. Bauer, M. Neubert, S. Renner, M. Schnubel, and A. Thamm, JHEP09, 056 (2022), arXiv:2110.10698 [hep- ph]

  6. [6]

    Navaset al.(Particle Data Group), Phys

    S. Navaset al.(Particle Data Group), Phys. Rev. D110, 030001 (2024)

  7. [7]

    Edwardset al., Phys

    C. Edwardset al., Phys. Rev. Lett.48, 903 (1982)

  8. [8]

    V. P. Druzhininet al., Z. Phys. C37, 1 (1987)

  9. [9]

    Insleret al.(CLEO), Phys

    J. Insleret al.(CLEO), Phys. Rev. D81, 091101 (2010), arXiv:1003.0417 [hep-ex]

  10. [10]

    del Amo Sanchezet al.(BaBar), Phys

    P. del Amo Sanchezet al.(BaBar), Phys. Rev. Lett.107, 021804 (2011), arXiv:1007.4646 [hep-ex]

  11. [11]

    I. S. Seonget al.(Belle), Phys. Rev. Lett.122, 011801 (2019), arXiv:1809.05222 [hep-ex]

  12. [12]

    Ablikimet al.(BESIII), Phys

    M. Ablikimet al.(BESIII), Phys. Rev. D101, 112005 (2020), arXiv:2003.05594 [hep-ex]

  13. [13]

    Ablikimet al.(BESIII), Phys

    M. Ablikimet al.(BESIII), Phys. Lett. B838, 137698 (2023), arXiv:2211.12699 [hep-ex]

  14. [14]

    Ablikimet al.(BESIII), Phys

    M. Ablikimet al.(BESIII), Phys. Rev. D105, 012008 (2022), arXiv:2109.12625 [hep-ex]

  15. [15]

    Achasovet al., (2023), arXiv:2303.15790 [hep-ex]

    M. Achasovet al., (2023), arXiv:2303.15790 [hep-ex]

  16. [16]

    Aiet al., Nucl

    X.-C. Aiet al., Nucl. Sci. Tech.36, 242 (2025), arXiv:2509.11522 [physics.acc-ph]

  17. [17]

    Altmannshoferet al.(Belle-II), PTEP2019, 123C01 (2019), [Erratum: PTEP 2020, 029201 (2020)], arXiv:1808.10567 [hep-ex]

    W. Altmannshoferet al.(Belle-II), PTEP2019, 123C01 (2019), [Erratum: PTEP 2020, 029201 (2020)], arXiv:1808.10567 [hep-ex]

  18. [18]

    Aaijet al.(LHCb), (2018), arXiv:1808.08865 [hep-ex]

    R. Aaijet al.(LHCb), (2018), arXiv:1808.08865 [hep-ex]

  19. [19]

    S. P. Adhyaet al.(CEPC Study Group), (2025), arXiv:2510.05260 [hep-ex]

  20. [20]

    Lyu (STCF Working Group), PoS BEAUTY2020, 060 (2021)

    X.-R. Lyu (STCF Working Group), PoS BEAUTY2020, 060 (2021)

  21. [21]

    Gao, Phys

    D.-N. Gao, Phys. Rev. D90, 077501 (2014), arXiv:1408.4552 [hep-ph]

  22. [22]

    Y. Meng, X. Feng, C. Liu, T. Wang, and Z. Zou, Sci. Bull.68, 1880 (2023), arXiv:2109.09381 [hep-lat]

  23. [23]

    Z. Zou, Y. Meng, and C. Liu, Chin. Phys. C46, 053102 (2022), arXiv:2111.00768 [hep-lat]

  24. [24]

    Meng, J.-L

    Y. Meng, J.-L. Dang, C. Liu, Z. Liu, T. Shen, H. Yan, and K.-L. Zhang, Phys. Rev. D109, 074511 (2024), arXiv:2401.13475 [hep-lat]

  25. [25]

    Meng, J.-L

    Y. Meng, J.-L. Dang, C. Liu, X.-Y. Tuo, H. Yan, Y.- B. Yang, and K.-L. Zhang, Phys. Rev. D110, 074510 (2024), arXiv:2407.13568 [hep-lat]

  26. [26]

    Y. Meng, C. Liu, T. Wang, and H. Yan, Phys. Rev. D 111, 014508 (2025), arXiv:2411.04415 [hep-lat]

  27. [27]

    G. Fan, Y. Meng, C. Liu, Z. Liu, T. Shen, T.-X. Wang, K.-L. Zhang, and L. Zhang, (2025), arXiv:2510.14478 [hep-lat]

  28. [28]

    Blossieret al.(ETM), JHEP07, 043 (2009), arXiv:0904.0954 [hep-lat]

    B. Blossieret al.(ETM), JHEP07, 043 (2009), arXiv:0904.0954 [hep-lat]

  29. [29]

    Becirevic and F

    D. Becirevic and F. Sanfilippo, JHEP01, 028 (2013), arXiv:1206.1445 [hep-lat]

  30. [30]

    Albaneseet al.(APE), Phys

    M. Albaneseet al.(APE), Phys. Lett. B192, 163 (1987)

  31. [31]

    G¨ usken, Nucl

    S. G¨ usken, Nucl. Phys. B Proc. Suppl.17, 361 (1990)

  32. [32]

    Constantinouet al.(ETM), JHEP08, 068 (2010), arXiv:1004.1115 [hep-lat]

    M. Constantinouet al.(ETM), JHEP08, 068 (2010), arXiv:1004.1115 [hep-lat]

  33. [33]

    Y. Meng, C. Liu, X.-Y. Tuo, H. Yan, and Z. Zhang, Eur. Phys. J. C85, 458 (2025), arXiv:2411.11533 [hep-lat]

  34. [34]

    Alexandrou, R

    C. Alexandrou, R. Baron, J. Carbonell, V. Drach, P. Gui- chon, K. Jansen, T. Korzec, and O. Pene (ETM), Phys. Rev. D80, 114503 (2009), arXiv:0910.2419 [hep-lat]

  35. [35]

    Baronet al.(ETM), JHEP08, 097 (2010), arXiv:0911.5061 [hep-lat]

    R. Baronet al.(ETM), JHEP08, 097 (2010), arXiv:0911.5061 [hep-lat]

  36. [36]

    McNeile and C

    C. McNeile and C. Michael (UKQCD), Phys. Rev. D70, 034506 (2004), arXiv:hep-lat/0402012

  37. [37]

    de Forcrand, M

    P. de Forcrand, M. Garcia Perez, H. Matsufuru, A. Naka- mura, I. Pushkina, I.-O. Stamatescu, T. Takaishi, and T. Umeda (QCD-TARO), JHEP08, 004 (2004), arXiv:hep-lat/0404016

  38. [38]

    Levkova and C

    L. Levkova and C. DeTar, Phys. Rev. D83, 074504 (2011), arXiv:1012.1837 [hep-lat]

  39. [39]

    Hatton, C

    D. Hatton, C. T. H. Davies, B. Galloway, J. Koponen, G. P. Lepage, and A. T. Lytle (HPQCD), Phys. Rev. D 102, 054511 (2020), arXiv:2005.01845 [hep-lat]

  40. [40]

    Ellis, Comput

    J. Ellis, Comput. Phys. Commun.210, 103 (2017), arXiv:1601.05437 [hep-ph]. Supplementary Information – S1 SUPPLEMENT AR Y MA TERIAL In this supplementary material, we expand on a se- lection of technical details of the calculations on the J/ψ→γν¯νbranching fraction. Relationship of hadronic function in Minkowski and Euclidean space In this section, we der...