Pith. sign in

REVIEW 2 major objections 4 minor 42 references

Relativistic quark-model calculation predicts ~10^{-3} fb cross sections for two C-even P-wave double-charmonium channels that only two virtual photons can produce.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-11 21:29 UTC pith:VT2TLOKI

load-bearing objection Solid, parameter-free LO predictions for two unobserved two-photon channels; the heavy-quark freeze of photon virtualities is the main uncontrolled step but does not erase the hierarchy or the usefulness of the numbers. the 2 major comments →

arxiv 2607.04126 v1 pith:VT2TLOKI submitted 2026-07-05 hep-ph

Exclusive production of P-wave charmonia through two virtual photons in electron-positron annihilation

Avinash Okram , Shashank Bhatnagar This is my paper
classification hep-ph
keywords Bethe-Salpeter equationdouble charmonium productionP-wave charmoniatwo-photon annihilationcross sectionsquark rearrangement
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.

The paper calculates exclusive production of two C-even P-wave charmonium pairs, χ_c1 + η_c and h_c + h_c, in electron-positron collisions at B-factory energy. Because charge conjugation and spin-parity forbid single-photon or photon-fragmentation paths, the reactions proceed only through two-virtual-photon quark-rearrangement diagrams. Using a fully relativistic Bethe-Salpeter equation whose parameters were already fixed by charmonium mass spectra, the authors obtain leading-order cross sections of order 10^{-3} fb and a definite hierarchy among ground and radially excited channels that is set by the radial nodes of the bound-state wave functions. The energy dependence is a smooth power-law fall-off. The absolute rates differ from non-relativistic estimates, yet the qualitative ordering agrees, furnishing clean benchmarks for future high-luminosity searches.

Core claim

With parameters fixed solely from spectroscopy, the 4×4 Bethe-Salpeter equation under the covariant instantaneous ansatz yields leading-order cross sections of order 10^{-3} fb for e^{-}e^{+} → γ*γ* → χ_c1(nP)+η_c(nS) and h_c(nP)+h_c(nP) at √s = 10.6 GeV. The hierarchy is σ(χ_c1(1P)η_c(1S)) > σ(χ_c1(2P)η_c(1S)) > σ(χ_c1(1P)η_c(2S)) > σ(χ_c1(2P)η_c(2S)), while for the identical-P-wave channel σ(h_c(2P)h_c(2P)) > σ(h_c(1P)h_c(1P)), both driven by the momentum-space structure of the radial wave functions that enter the overlap integrals.

What carries the argument

The 4×4 Bethe-Salpeter wave functions of the P-wave and S-wave charmonia, reduced under the covariant instantaneous ansatz to three-dimensional Salpeter amplitudes whose radial shapes and normalizers are taken unchanged from earlier spectroscopic fits; these amplitudes supply the relativistic vertices that enter the two-photon quark-rearrangement diagrams.

Load-bearing premise

The calculation assumes that the internal relative momenta of the heavy quarks can be neglected, so each virtual photon simply carries half the collision energy; if that approximation fails, the amplitude factorization and all numerical rates change.

What would settle it

A high-luminosity e^{+}e^{-} measurement (or a sufficiently tight upper limit) of the χ_c1(1P)η_c(1S) or h_c(1P)h_c(1P) cross section near 10.6 GeV that lies well outside the predicted 10^{-3} fb range and its quoted parameter uncertainty.

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

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

2 major / 4 minor

Summary. The manuscript computes leading-order cross sections for the C-even exclusive channels e^{-}e^{+} → γ*γ* → χ_c1(nP)+η_c(nS) and e^{-}e^{+} → γ*γ* → h_c(nP)+h_c(nP) at √s = 10.6 GeV in the 4×4 Bethe–Salpeter equation under the covariant instantaneous ansatz. Photon-fragmentation topologies are forbidden by C-parity and spin-parity, so only quark-rearrangement diagrams contribute. Using a kernel and parameters previously fixed from charmonium spectroscopy, the authors reduce the amplitudes under a heavy-quark approximation (k^{2} ≃ k'^{2} ≃ s/4), evaluate the Dirac traces and 3D overlap integrals G1,G2 (χ_c1η_c) and F1,F2 (h_c h_c), and obtain LO rates of order 10^{-3} fb together with a definite radial hierarchy and a smooth power-law fall-off with energy. Parameter variations are performed and uncertainties quoted; qualitative ordering is compared with existing NRQCD results.

Significance. If the absolute rates and the reported hierarchies survive a controlled estimate of the O(v^{2}) corrections that were neglected in the photon kinematics, the work supplies genuine, parameter-free predictions for two experimentally unobserved two-photon channels that cleanly isolate relativistic quark-rearrangement dynamics. The unified spectroscopy-to-production pipeline and the explicit contrast between unweighted (G1G2) and momentum-weighted (F1F2) overlaps are useful benchmarks for Belle II and future e^{+}e^{-} facilities, and they complement NRQCD by exposing the role of full Dirac structure in P-wave production.

major comments (2)
  1. Sections 3 and 4 (after Eqs. (5)–(6) and again before Eq. (22)): the entire numerical content of Tables 3 and 5 rests on the heavy-quark approximation that freezes the photon virtualities at k^{2} ≃ k'^{2} ≃ s/4, allowing the constant 16/s^{2} to be pulled outside every loop integral. Because the P-wave BS vertices already carry explicit factors of ˆq (or ˆq·ϵ), the same momentum region that was declared negligible for the photon kinematics is precisely the region sampled by the production overlaps. No estimate of the size of the resulting O(v^{2}) shift, nor any recomputation with a momentum-dependent photon propagator, is provided. Absolute rates of order 10^{-3} fb and the precise numerical hierarchy are therefore controlled by an uncontrolled approximation that must be quantified before the central claims can be regarded as robust.
  2. Section 4.1 and the scaling arguments following Table 5: the claimed enhancement σ(h_c(2P)+h_c(2P)) > σ(h_c(1P)+h_c(1P)) is attributed to the larger value of |N_A F1| for the 2P state. That quantity is itself proportional to the same high-ˆq region whose neglect underpins the factorization. Without a controlled O(v^{2}) assessment, it remains unclear whether the reported 2.6-fold enhancement is a genuine dynamical prediction or an artifact of the approximation that freezes the photon momenta.
minor comments (4)
  1. Abstract and Introduction: the phrase “order 10^{-3} fb” is used for both channels, yet Table 5 shows the h_c(1P)+h_c(1P) rate is an order of magnitude smaller; a more precise statement would avoid overstating the common magnitude.
  2. Eq. (20) and Eq. (28): the lengthy |M|^{2} expressions contain several nested square roots and trigonometric factors that are never reduced or cross-checked against a simpler kinematic limit; a short appendix verifying the high-s or θ=0 limits would improve transparency.
  3. Figures 2 and 3: the energy-dependence plots lack error bands corresponding to the parameter variations already performed in Tables 1 and 4; adding them would make the claimed power-law fall-off more informative.
  4. References: several recent NRQCD and light-cone calculations of related two-photon double-charmonium channels (post-2020) are not cited; a brief comparison would place the present LO BSE results in clearer context.

Circularity Check

0 steps flagged

No significant circularity: parameters and 3D wave functions fixed once from mass spectroscopy yield genuine LO predictions for unobserved two-photon channels.

full rationale

The free parameters (ω₀ = 0.22 GeV, m_c = 1.490 GeV, Λ = 0.250 GeV, C₀, A₀) and the radial Salpeter wave functions ϕ_A(ˆq), ϕ_P(ˆq) are taken unchanged from earlier spectroscopic fits to experimental charmonium masses (refs. [21,22] of the same group). Those fits do not involve the production amplitudes of the present paper. The LO cross sections for χ_c1(nP)+η_c(nS) and h_c(nP)+h_c(nP) are then obtained by inserting the already-normalized 3D BS vertices into the quark-rearrangement diagrams, evaluating the overlap integrals G1G2 or F1F2, and integrating the resulting |M̄_fi|^{2}. Because no production data exist and no parameters are re-tuned, the numerical results (Tables 3 and 5) and the reported hierarchy are genuine predictions of the previously constrained wave functions. The only residual self-citation is the reuse of the authors’ own spectroscopic solutions and the CIA reduction; this is ordinary sequential model-building, not a circular reduction of the claimed cross sections to their inputs. External comparison with NRQCD further anchors the qualitative ordering. The heavy-quark approximation k^{2} ≃ k'^{2} ≃ s/4 is an uncontrolled modeling assumption, not a circularity.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The central numerical claims rest on a small set of parameters previously fitted to spectroscopy, on the standard BSE under the covariant instantaneous ansatz, and on a handful of kinematic approximations that simplify the two-photon amplitudes. No new dynamical entities are introduced.

free parameters (4)
  • ω_{0} (flavour-independent spring constant) = 0.22 GeV
    Controls the strength of the confining piece of the BSE kernel; fixed to 0.22 GeV by earlier mass-spectrum fits and varied here only for sensitivity.
  • m_c (constituent charm-quark mass) = 1.490 GeV
    Appears in the quark propagators and in the inverse-range parameters β; fixed to 1.490 GeV by spectroscopy.
  • Λ (QCD scale in α_s) = 0.250 GeV
    Enters the running coupling that multiplies both the one-gluon-exchange and confining kernels; fixed to 0.250 GeV.
  • C_{0}, A_{0} (dimensionless constants in confining potential) = C_{0}=0.69, A_{0}=0.01
    Additional shape parameters of the confining kernel, taken from the same spectroscopic fits (C_{0}=0.69, A_{0}=0.01).
axioms (4)
  • domain assumption Covariant Instantaneous Ansatz: the interaction kernel depends only on the transverse relative momentum q̂, allowing exact reduction of the 4D BSE to 3D Salpeter equations while preserving Lorentz covariance of the wave functions.
    Stated in Section 2 and used throughout the amplitude reductions of Sections 3 and 4.
  • domain assumption Heavy-quark approximation: internal relative momenta of the quarkonia may be neglected relative to the hard scale, so that the two photon virtualities satisfy k^{2} eq k'² eq s/4.
    Explicitly invoked in the opening paragraphs of Sections 3 and 4 to factor the photon propagators out of the loop integrals.
  • domain assumption Constituent (constant) quark mass is an adequate approximation for the heavy-quark propagators inside the BSE.
    Defended in Section 2 by reference to earlier literature; the authors note that a future upgrade to dressed propagators is planned.
  • domain assumption Color-octet and higher-order QCD corrections are power-suppressed and may be omitted at leading order in α_em^{4}.
    Stated in the introduction and again for the h_c+h_c channel in Section 4.

pith-pipeline@v1.1.0-grok45 · 24004 in / 3357 out tokens · 23451 ms · 2026-07-11T21:29:41.779798+00:00 · methodology

0 comments
read the original abstract

We present a relativistic study of the exclusive double-charmonium processes $e^-e^+ \to \gamma^*\gamma^* \to \chi_{c1}+\eta_c$ and $e^-e^+ \to \gamma^*\gamma^* \to h_c+h_c$ at $\sqrt{s}=10.6$ GeV in the framework of the $4\times4$ Bethe-Salpeter equation. Since these channels are forbidden in single-photon annihilation, they proceed purely through two-photon quark-rearrangement mechanisms and thus provide a clean probe of relativistic quarkonium dynamics. Using parameters fixed from mass spectroscopy in our earlier work, we obtain leading-order predictions for the total cross sections and their energy dependence. The predicted cross sections are of order $10^{-3}$ fb and exhibit a hierarchy that originates from the radial structure of the meson wave functions, which determines the overlap integrals contributing to the production amplitude. While the absolute magnitudes of the cross sections differ from NRQCD estimates, the qualitative ordering of the channels is consistent with NRQCD expectations. We find a smooth power-law falloff of the cross section for both $\chi_{c1}+\eta_c$ and $h_c+h_c$. These results provide benchmark predictions for future searches at high-luminosity $e^-e^+$ colliders.

Figures

Figures reproduced from arXiv: 2607.04126 by Avinash Okram, Shashank Bhatnagar.

Figure 1
Figure 1. Figure 1: Feynman diagrams involved in e −e + → γ ∗γ ∗ → H + H′ at leading order. ∼ O(v 2 ), implying ˆq 2/m2 << 1. In this limit, the photon momenta are effectively insensitive to the internal quark motion, and the photon virtualities satisfy k 2 ≈ k ′2 ≈ s/4. The product of the photon propagators then satisfies, 1/(k 2k ′2 ) ≈ 16/s2 . This is further justified since the photon virtualities are dominated by the ext… view at source ↗
Figure 2
Figure 2. Figure 2: Plot of cross section for e −e + → γ ∗γ ∗ → χc1(1P) + ηc(1S) (in fb) versus √ s (in GeV) Further from the cross section results it can be checked that the hierarchy is physically reasonable and consistent with expectations, and it can be justified quite cleanly within both BSE-type models and general QCD intuition. For the channels where NRQCD results are available (1P+1S and 1P+2S), both approaches agree … view at source ↗
Figure 3
Figure 3. Figure 3: Plot of cross section for e −e + → γ ∗γ ∗ → hc(1P) + hc(1P) (in fb) versus √ s (in GeV) The plot of the cross section for e −e + → hc(1P) + hc(1P) versus √ s in Fig.3 displays a peak near the threshold followed by the asymptotic fall at higher energies. 13 [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

42 extracted references · 12 linked inside Pith

  1. [1]

    Abe et al.(Belle), Phys

    K. Abe et al.(Belle), Phys. Rev. Lett. 89, 142001 (2002);[hep-ex/0205104]

    Pith/arXiv arXiv 2002

  2. [2]

    Abe et al.(Belle), Phys

    K. Abe et al.(Belle), Phys. Rev. D 70, 071102 (2004); [hep-ex/0407009]

    Pith/arXiv arXiv 2004

  3. [3]

    Aubert et al.(BABAR),Phys

    B. Aubert et al.(BABAR),Phys. Rev. D 72, 031101(2005); [hep-ex/0506062]

    Pith/arXiv arXiv 2005

  4. [4]

    G. T. Bodwin, J. Lee, and E. Braaten, Phys. Rev. D 67, 054023(2003), [hep-ph/0212352].[Erratum: Phys.Rev.D 72, 099904 (2005)]

    Pith/arXiv arXiv 2003

  5. [5]

    Bodwin, E

    G.T. Bodwin, E. Braaten, G.P. Lepage, Phys. Rev. D 51 (1995) 1125, Phys. Rev. D 55 (1997) 5853

    1995

  6. [6]

    G. T. Bodwin, J. Lee, and E. Braaten, Phys. Rev. Lett.90. 162001(2003), [hep-ph/0212181]

    Pith/arXiv arXiv 2003

  7. [7]

    G. T. Bodwin, E. Braaten, J. Lee, and C. Yu, Phys. Rev. D 74, 074014(2006), [hep-ph/0608200]

    Pith/arXiv arXiv 2006

  8. [8]

    Gong and J.-X

    B. Gong and J.-X. Wang, Phys. Rev. Lett. 100, 181803(2008), [arXiv:0801.0648]

    Pith/arXiv arXiv 2008

  9. [9]

    Gong, R-C Niu, H-M Yu, J-X Wang, JHEP 02, 055 (2024); arxiv: 2311.04751[hep-ph]

    X-D Huang, B. Gong, R-C Niu, H-M Yu, J-X Wang, JHEP 02, 055 (2024); arxiv: 2311.04751[hep-ph]

    Pith/arXiv arXiv 2024

  10. [10]

    Y.Fan, J.Lee, C.Yu, Phys. Rev. D 87, 094032 (2013); arxiv:1211.4111[hep-ph]

    Pith/arXiv arXiv 2013

  11. [11]

    Lee, Phys

    E.Braaten, J. Lee, Phys. Rev. D67, 054007,(2003)

    2003

  12. [12]

    A.V.Luchinsky,Phys. Atom. Nucl. 67, 1338 (2004); hep-ph/0301190

    Pith/arXiv arXiv 2004

  13. [13]

    Davier, M

    M. Davier, M. E. Peskin, and A. Snyder, hep-ph/0606155

    Pith/arXiv arXiv

  14. [14]

    S.Bhatnagar, H.Negash, Nucl. Phys. A1053, 122969 (2025)

    2025

  15. [15]

    V.V.Braguta, Phys. Rev. D78, 054025 (2008); arxiv:0712.1475[hep-ph]

    Pith/arXiv arXiv 2008

  16. [16]

    S.Bhatnagar, V.Guleria, Nucl. Phys. A1041, 122783(2024)

    2024

  17. [17]

    64, 82(2023)

    M.Narang, S.Bhatnagar, Few-Body Systs. 64, 82(2023)

    2023

  18. [18]

    S.Bhatnagar, Phys. Rev. D112, 054011 (2025). 16

    2025

  19. [19]

    S.Bhatnagar, L.Alemu, Phys. Rev. D97, 034021 (2018)

    2018

  20. [20]

    Physics 99, 169 (2025)

    V.Guleria, E.Gebrehana, S.Bhatnagar, Pramana J. Physics 99, 169 (2025)

    2025

  21. [21]

    E.Gebrehana, S.Bhatnagar, H.Negash, Phys. Rev. D100, 054034 (2019)

    2019

  22. [22]

    V.Guleria, S.Bhatnagar, Intl. J. Theor. Phys. 60, 3143 (2021)

    2021

  23. [23]

    V.Guleria, E.Gebrehana, S.Bhatnagar, Phys. Rev. D104, 094045 (2021)

    2021

  24. [24]

    S.Bhatnagar, E.Gebrehana, Phys. Rev. D102, 094024 (2020)

    2020

  25. [25]

    A.N.Mitra, S.Bhatnagar, Intl. J. Mod. Phys. A07, 121 (1992)

    1992

  26. [26]

    T. Wang, Y. Jiang, W.L. Ju, H. Yuan, G.L. Wang, J. High Energy Phys. 03, 209 (2016)

    2016

  27. [27]

    J.K. He, Eur. Phys. J. C 79, 393 (2019)

    2019

  28. [28]

    J.K. He, C.J. Fan, Phys. Rev. D 103, 114006 (2021)

    2021

  29. [29]

    Z.H.Wang, G.L.Wang, Phys. Rev. D106, 054037 (2022)

    2022

  30. [30]

    H.Negash, S.Bhatnagar, Intl. J. Mod. Phys. E25, 1650059 (2016)

    2016

  31. [31]

    Bhagwat, A

    M.S. Bhagwat, A. Krassnigg, P. Maris, C.D. Roberts, arXiv :nucl -th /0612027

  32. [32]

    Bender, C.D

    A. Bender, C.D. Roberts, L.V. Smekal, Phys. Lett. B 380 (1996) 7

    1996

  33. [33]

    Maris, A

    P. Maris, A. Raya, C.D. Roberts, S.M. Schmidt, Eur. Phys. J. A 18 (2003) 231

    2003

  34. [34]

    Maris, C.D

    P. Maris, C.D. Roberts, Int. J. Mod. Phys. E 12 (2003) 297

    2003

  35. [35]

    El-Bennich, M.A

    B. El-Bennich, M.A. Ivanov, C.D. Roberts, Phys. Rev. C 83 (2011) 025205

    2011

  36. [36]

    Chang, J.X

    C.H. Chang, J.X. Wang, X.G. Wu, Sci. China, Phys. Mech. Astron. 53 (2010) 2031, arXiv :1005 .4723 [hep -ph]

    2010

  37. [37]

    Wang, W.M

    Z.G. Wang, W.M. Yang, S.L. Wan, Phys. Lett. B 615 (2005) 75, arXiv :hep -ph /0411142

    2005

  38. [38]

    J. Yi, S-Y. Li, Y-R. Liu, Z-X. Meng, Z-G. Si, T. Yao, Phys. Rev. C 102 (2020) 015201

    2020

  39. [39]

    C.H.L.Smith, Ann. Phys. (N.Y.) 53, 521 (1969)

    1969

  40. [40]

    Alkofer and L

    R. Alkofer and L. V. Smekel, Phys. Rep. 353, 281 (2002)

    2002

  41. [41]

    Q-L.Liao, Y. Yu, Eur. Phys. J. C84, 1000 (2024)

    2024

  42. [42]

    Q.L. Liao, J. Jiang, Y.H. Zhao, Eur. Phys. J. C83, 22 (2022). 17

    2022