Pith. sign in

REVIEW 3 major objections 4 minor 28 references

Di-$\pi^0$ Production and Generalized Distribution Amplitudes at Future Electron-Ion Colliders

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read This paper claims that exclusive di-$\pi^0$ production at the Electron-Ion Collider can measure pion generalized distribution amplitudes, with e-Au collisions producing about 10,000 times more events than the current e+e- measurement.

desk verdict The e-p part is plausible, but the e-Au '10^4' claim rests on a bare nuclear flux without a survival factor. read the letter →

arxiv 2501.01267 v1 pith:3J6FH7HS submitted 2025-01-02 hep-ph

Bing'ang Guo , Jing Han , Ya-Ping Xie , Xurong Chen This is my paper · ORCID
classification hep-ph
keywords generalizeddistributionamplitudesdi-pionproductionequivalentphotonapproximationElectron-IonColliderexclusivemesonpairheavy-ionfluxpi0structureultra-peripheralcollisions
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 hadron generalized distribution amplitudes (GDAs), which encode how a quark-antiquark pair turns into a pair of hadrons and reveal the three-dimensional structure of hadrons, can be measured at the future Electron-Ion Collider through exclusive di-$\pi^0$ production. In electron-hadron collisions the ion emits a quasi-real photon whose flux grows with the square of the nuclear charge, so the $\gamma^*\gamma\to\pi^0\pi^0$ cross section should be much larger than in electron-positron colliders. Using the equivalent photon approximation and leading-twist $\pi^0$ GDAs extracted from e+e- data, the authors predict that e-p collisions at the EIC would produce about as many di-$\pi^0$ events as the current e+e- experiment in the same kinematic region, and e-Au collisions about $10^4$ times more. If right, this turns EIC heavy-ion beams into a high-statistics laboratory for GDA physics and for related quantities such as hadron gravitational form factors.

What carries the argument

The machinery is the equivalent photon approximation (EPA) for the quasi-real photon emitted by the hadron or ion, combined with the leading-twist $\pi^0$ generalized distribution amplitude as the nonperturbative input for the $\gamma^*\gamma\to\pi^0\pi^0$ subprocess. The photon flux for a proton is $f_{\gamma/p}(x)\propto \alpha_{\rm em}/(2\pi x)\,[1+(1-x)^2](\ln\Omega-11/6+\dots)$, while for a nucleus of charge $Z$ it is $f_{\gamma/A}(x)\propto 2Z^2\alpha_{\rm em}/(\pi x)\,[\xi K_0(\xi)K_1(\xi)-\xi^2/2\,(K_1^2(\xi)-K_0^2(\xi))]$, with the $Z^2$ factor giving the large enhancement in e-Au collisions. This flux is convoluted with the differential $e\gamma$ cross section expressed through helicity amplitudes built from the $\pi^0$ GDAs, using GDAs extracted from a leading-twist analysis of the Belle measurement. The factorized formula $d\sigma(eX\to eX\pi^0\pi^0)=\int dx\, f_{\gamma/X}(x)\,d\sigma(e\gamma\to e\pi^0\pi^0)\,\delta(\dots)$ is what converts a known e+e- measurement into a prediction for e-p and e-Au rates.

What would settle it

At the EIC, measure the exclusive $e\text{Au}\to e\text{Au}\,\pi^0\pi^0$ cross section near $W\simeq 1$ GeV and $Q^2\simeq 16$ GeV$^2$; if the event rate is more than a factor of a few below the EPA prediction including the $Z^2$ enhancement, the central event-rate claim is falsified.

Watch

Extended reading notes

Core claim

The central claim is that exclusive two-pion production in electron-ion collisions, with the ion radiating a quasi-real photon and the electron a virtual photon, gives access to the same $\pi^0$ GDA that was measured in e+e- collisions, but with event rates boosted by the ion's large electric charge. Concretely, the paper predicts that within the same kinematic region e-p collisions at the EIC produce a number of di-$\pi^0$ events comparable to those generated by Belle II, while e-Au collisions produce roughly $10^4$ times more. The calculation is made by folding the equivalent-photon flux of the proton or nucleus with the $e\gamma\to e\pi^0\pi^0$ cross section computed from leading-twist $\pi^0$ GDAs, and includes the $f_2(1270)$ resonance peak in the invariant-mass dependence. The authors therefore propose electron-heavy-ion collisions at the EIC as a precision route to GDAs and the nonperturbative structure of QCD.

Load-bearing premise

The rate calculation assumes every quasi-real photon emitted by the proton or gold nucleus produces a pion pair with no extra suppression from absorption, survival probability, or higher-twist corrections; if those suppress the rate, the quoted event counts, especially for gold, will be too high.

Editorial extensions

If this is right

  • If the prediction holds, e-p collisions at the EIC with $15\,\text{fb}^{-1}\cdot\text{y}^{-1}$ will deliver a di-$\pi^0$ sample comparable to the current e+e- experiment, giving an independent and higher-energy handle on the $\pi^0$ GDA.
  • Electron-gold collisions at the EIC would produce roughly $10^4$ times more di-$\pi^0$ events than the current e+e- measurement in the same kinematic region, turning heavy-ion beams into a precision source for GDA extraction.
  • The $Z^2$ scaling of the nuclear photon flux means heavier isocurrent nuclei give larger rates, so the e-Au channel should be the preferred setting for high-precision GDA studies.
  • The predicted invariant-mass spectrum shows a clear $f_2(1270)$ resonance peak, which can be used as a cross-check of the GDA input and of the equivalent-photon treatment.
  • These cross sections provide event-rate estimates needed to design e-p and e-Au ultra-peripheral collision measurements at the EIC and EicC.

Reading between the lines

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

  • The quoted event counts assume no suppression from nuclear survival probability, absorptive corrections, or higher-twist effects; if those are sizable, the e-Au rate could be substantially lower than $10^4$ times the e+e- rate.
  • The same EPA-plus-GDA construction should apply to other meson pairs such as charged pions, kaons, or charm mesons, so the method, if validated on $\pi^0\pi^0$, would open a program of GDA spectroscopy at the EIC.
  • Because the cross section is computed from e+e--fitted GDAs, a high-statistics EIC measurement would effectively test the scale dependence and higher-twist sensitivity of the leading-twist GDA extraction, something a single e+e- data set cannot do.
  • A direct e-p versus e-Au comparison at fixed invariant mass and $Q^2$ would isolate the $Z^2$ photon-flux enhancement and provide a clean check of the equivalent-photon treatment independent of GDA modelling.
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 / 4 minor

Summary. The paper proposes that generalized distribution amplitudes (GDAs) of the pion can be studied at future Electron-Ion Colliders through the subprocess gamma* gamma -> pi0 pi0, where the quasi-real photon is emitted by the proton or nucleus and the virtual photon by the electron. Using the equivalent photon approximation (EPA) for the hadron/ion flux (Eqs. (4)-(5)) and a Belle-extracted leading-twist pi0 GDA in the e-gamma cross section (Eq. (7)), the authors compute differential cross sections for e-p and e-Au collisions at EIC and EicC and estimate event rates. They claim that e-p collisions at the EIC will produce a number of di-pi0 events comparable to Belle II, while e-Au collisions will produce about 10^4 times more, making the EIC a high-statistics facility for GDA studies.

Significance. If the central quantitative claims hold, the paper would establish a new and potentially powerful way to access pion GDAs, extending the GDA program from e+e- colliders to electron-hadron and electron-nucleus collisions. The Z^2 enhancement of the nuclear photon flux is kinematically robust and the EPA framework is a standard tool, so the basic idea is credible. The paper also makes concrete, falsifiable predictions for EIC and EicC kinematics, which is a strength. However, the advertised event-rate gain for e-Au depends on a bare coherent photon flux with no survival/absorption factor for nuclear breakup, and the event-rate comparisons are not backed by tabulated integrated yields or uncertainty estimates. These issues are load-bearing for the paper's headline claims.

major comments (3)
  1. [Eq. (5) and Sec. III] The e-Au event-rate claim rests on the bare coherent photon flux of Eq. (5), which is integrated over all impact parameters and counts photons even from configurations in which the Au nucleus breaks up or undergoes additional inelastic/electromagnetic interactions. For the exclusive final state e Au -> e Au pi0 pi0, each impact parameter must be weighted by a survival probability S(b) that the nucleus remains intact; no such factor appears in Eqs. (3)-(5) or in the event estimate. Since the advertised factor ~10^4 is the headline result, the omission is load-bearing; a survival factor of order 0.5-0.8 would change the central claim by factors of roughly 2 or more, and electromagnetic dissociation of Au in e-A collisions can be sizable. Please include a survival/absorption factor, or justify quantitatively why it is negligible in the selected kinematic region.
  2. [Sec. III, event estimate] The statement that e-p at EIC produces events 'comparable to' Belle II and e-Au 'about 10^4 times' more is not reproducible from the manuscript. No integrated yields are tabulated; the reader is not told the exact ranges of W, Q^2, se_gamma, and x over which the events are counted, the integrated luminosity per year used for each facility beyond the 15 fb^-1 y^-1 quoted for e-p, or the angular/kinematic acceptance implied by the uniform 70% efficiency. Please add a table with N(e-p), N(e-Au), and N(Belle II) for identical phase-space cuts, and state all luminosity and efficiency assumptions.
  3. [Eq. (7) and Sec. III] The subprocess cross section is evaluated with pi0 GDAs extracted from Belle data [11], but no uncertainty on the GDA parameters is propagated into the cross sections or event rates. Since the Belle extraction itself was limited by 'considerable statistical uncertainty' (Sec. III), the predictions inherit an unquantified uncertainty that is comparable in size to the proposed precision gain. Please show an error band from the GDA fit parameters and, ideally, a comparison with an alternative GDA model.
minor comments (4)
  1. [Section heading] The heading 'NUMERICAL RESUL TS' contains a typo and should read 'NUMERICAL RESULTS'.
  2. [Fig. 3 caption] The label 'f_gamma/Au,EicC x 10^9' is not explained in the caption; please state the scaling explicitly so the reader can interpret the curve.
  3. [Eq. (3)] The delta function in Eq. (3) is redundant; consider writing the convolution directly as dsigma/dQ^2 dW^2 = f_gamma/X(s_e_gamma/s_eX) dsigma_e_gamma/dQ^2 dW^2 after fixing x, or explicitly showing the integration over x.
  4. [Sec. III, EicC] The exclusion of e-Au events at EicC is stated without a numerical estimate; a one-line estimate of the expected number of events would make the decision transparent.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the EIC event-rate prediction is a forward application of Belle-extracted GDAs and standard EPA photon fluxes, not a fit to the claimed result.

full rationale

The derivation chain is: (i) parameterize the γ*γ → π0π0 amplitudes with π0 GDAs 'extracted from leading-twist analysis of the Belle measurements [11]' (Sec. II); (ii) fold the eγ cross section of Eq. (7) with the equivalent-photon fluxes of Eqs. (4)-(6); (iii) integrate to obtain cross sections and multiply by luminosity and detection efficiency to estimate events in Sec. III. The GDA input comes from external Belle data via Ref. [11] (Kumano, Song, Teryaev), not from the present authors, and the eγ cross-section formula is cited to Refs. [3,18]. The photon-flux formulas in Eqs. (4) and (5) are standard EPA results cited jointly to the external review [21] and the authors' own Ref. [22]; the nuclear radius formula in Eq. (5) is cited jointly to the authors' Ref. [23] and to Lappi-Mantysaari [24]. These self-citations are not load-bearing because the formulas are textbook-standard and independently corroborated by the external references. No parameter is fitted to the EIC event rate, and no predicted quantity reduces by construction to an input of the calculation. The reviewer's concern about missing survival or absorption factors for coherent e-Au events is a physics-validity risk that would affect the magnitude of the prediction, but it is not a circularity of the derivation. Overall, the paper's central feasibility claim is a straightforward forward application of existing external inputs, so the circularity burden is low.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central event-rate numbers rest entirely on the Belle-fitted GDA input and on the validity of EPA for nuclear beams; no new particles, forces, or entities are introduced by this paper beyond the assumed detection efficiency.

free parameters (2)
  • pi0 GDA parameters from Belle fit = from Kumano, Song, Teryaev, PRD 97, 014020 (2018)
    Input hadronic amplitude; all cross section magnitudes depend on it.
  • Detection efficiency = 70%
    Chosen uniformly across facilities to convert cross sections to event counts; does not affect cross-section ratios.
assumptions (4)
  • domain assumption Leading-twist factorization of gamma* gamma to pi pi in terms of GDAs is valid at Q^2 = 9-25 GeV^2 and W = 0.5-2 GeV, including the f2(1270) resonance region.
    The Belle-fitted GDA is leading-twist; at low W near the resonance, higher-twist contributions may be significant. The paper cites higher-twist effects [5,6,18] but neglects them.
  • domain assumption The equivalent photon approximation accurately describes the photon flux from protons and nuclei in the considered kinematic region.
    Eqs. (4) and (5) from Refs. [21,22] are used without modification or validation at the EIC kinematics.
  • domain assumption The gluon GDA contribution is negligible compared to the quark GDA contribution.
    Stated in Sec. II as suppressed by higher orders of alpha_s, but the suppression is not quantified numerically.
  • ad hoc to paper Exclusive e-A di-pion events can be modeled without a survival probability or veto for additional hadronic interactions and nuclear breakup.
    The calculation multiplies the coherent nuclear photon flux from Eq. (5) by the full gamma* gamma cross section with no absorptive correction; this likely overestimates the exclusive e-Au rate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Di-$\pi^0$ Production and Generalized Distribution Amplitudes at Future Electron-Ion Colliders." pith.science (2026). https://pith.science/paper/3J6FH7HS

@misc{pith2026250101267,
  author       = {Pith},
  title        = {Pith review of: Di-$\pi^0$ Production and Generalized Distribution Amplitudes at Future Electron-Ion Colliders},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3J6FH7HS}},
  note         = {Machine review of arXiv:2501.01267}
}
abstract

Generalized distribution amplitudes (GDAs) offer valuable insights into the three-dimensional structure of hadrons, delineating the amplitudes associated with the transition from a quark-antiquark pair to a hadron pair. Currently, hadron GDAs can be probed in electron-positron collisions, with experimental feasibility demonstrated at facilities such as Belle and BESIII. In this study, we put forth the proposition that hadron GDAs can also be investigated in electron-hadron collisions at forthcoming Electron-Ion Colliders (EICs), specifically through the subprocess $\gamma^*\gamma \to h_1 h_2$. In this framework, a quasi-real photon, emitted by the ion, exhibits a photon flux proportional to the square of the ion's electric charge. Consequently, we anticipate that the cross sections in EICs will be substantially larger than those in electron-positron collisions. We present numerical calculations pertaining to di-$\pi^0$ production employing the equivalent photon approximation (EPA). Our findings suggest that, within the same kinematic region, electron-proton ($e$-$p$) collisions at the EIC could yield an event rate comparable to that of Belle II, while electron-gold ($e$-Au) collisions are expected to generate an even greater number of events. This enhanced event rate facilitates a high-precision examination of di-$\pi^0$ GDAs at the EIC.

Figures

Figures reproduced from arXiv: 2501.01267 by the authors.

Figure 1
Figure 1. FIG. 1. Feynman diagram for di- [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. shows the eγ → eπ0π 0 process in the c.m. frame of the di-photon, corre￾3 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Photon fluxes under different facilities’ kinematics as functions of [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Differential cross sections for di- [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Differential cross sections for di- [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 11 canonical work pages

  1. [11]

    Kumano, Q.-T

    S. Kumano, Q.-T. Song, and O. V. Teryaev, Phys. Rev. D 97, 014020 (2018), arXiv:1711.08088 [hep-ph]

    arXiv 2018

  2. [1]

    M¨ uller, D

    D. M¨ uller, D. Robaschik, B. Geyer, F. M. Dittes, and J. Hoˇ rejˇ si, Fortsch. Phys.42, 101 (1994), arXiv:hep-ph/9812448

    arXiv 1994

  3. [2]

    Diehl, T

    M. Diehl, T. Gousset, B. Pire, and O. Teryaev, Phys. Rev. Lett. 81, 1782 (1998), arXiv:hep-ph/9805380. 9

    arXiv 1998

  4. [3]

    Diehl, T

    M. Diehl, T. Gousset, and B. Pire, Phys. Rev. D 62, 073014 (2000), arXiv:hep- ph/0003233

    work page arXiv 2000

  5. [4]

    Pion pair production in $e^+e^-$ annihilation

    Z. Lu and I. Schmidt, Phys. Rev. D 73, 094021 (2006), [Erratum: Phys.Rev.D 75, 099902 (2007)], arXiv:hep-ph/0603151

    work page Pith review arXiv 2006

  6. [5]

    Kinematical higher-twist corrections in $\gamma^* \to M_1 M_2 \gamma$: Neutral meson production

    B. Pire and Q.-T. Song, Phys. Rev. D 107, 114014 (2023), arXiv:2304.06389 [hep-ph]

    work page Pith review arXiv 2023

  7. [6]

    Kinematical higher-twist corrections in $\gamma^* \to M_1 M_2 \gamma$: Charged meson production

    B. Pire and Q.-T. Song, Phys. Rev. D 109, 074016 (2024), arXiv:2311.06005 [hep-ph]

    work page Pith review arXiv 2024

  8. [7]

    M. V. Polyakov and P. Schweitzer, Int. J. Mod. Phys. A 33, 1830025 (2018), arXiv:1805.06596 [hep-ph]

    arXiv 2018

Show all 28 references
  1. [8]

    V. D. Burkert, L. Elouadrhiri, F. X. Girod, C. Lorc´ e, P. Schweitzer, and P. E. Shanahan, Rev. Mod. Phys. 95, 041002 (2023), arXiv:2303.08347 [hep-ph]

    2023 arXiv

  2. [9]

    M. V. Polyakov, Nucl. Phys. B 555, 231 (1999), arXiv:hep-ph/9809483

    1999 arXiv

  3. [10]

    Goeke, M

    K. Goeke, M. V. Polyakov, and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 (2001), arXiv:hep-ph/0106012

    2001 arXiv

  4. [12]

    Chen and H.-n

    C.-H. Chen and H.-n. Li, Phys. Lett. B 561, 258 (2003), arXiv:hep-ph/0209043

    2003 arXiv

  5. [13]

    Wang, H.-n

    W.-F. Wang, H.-n. Li, W. Wang, and C.-D. L¨ u, Phys. Rev. D 91, 094024 (2015), arXiv:1502.05483 [hep-ph]

    2015 arXiv

  6. [14]

    Li, A.-J

    Y. Li, A.-J. Ma, W.-F. Wang, and Z.-J. Xiao, Phys. Rev. D 95, 056008 (2017), arXiv:1612.05934 [hep-ph]

    2017 arXiv

  7. [15]

    Kawamura and S

    H. Kawamura and S. Kumano, Phys. Rev. D 89, 054007 (2014), arXiv:1312.1596 [hep-ph]

    2014 arXiv

  8. [16]

    Kawamura, S

    H. Kawamura, S. Kumano, and T. Sekihara, Phys. Rev. D 88, 034010 (2013), arXiv:1307.0362 [hep-ph]

    2013 arXiv

  9. [17]

    Chang, S

    W.-C. Chang, S. Kumano, and T. Sekihara, Phys. Rev. D 93, 034006 (2016), arXiv:1512.06647 [hep-ph]

    2016 arXiv

  10. [18]

    Lorc´ e, B

    C. Lorc´ e, B. Pire, and Q.-T. Song, Phys. Rev. D106, 094030 (2022), arXiv:2209.11140 [hep-ph]

    2022 arXiv

  11. [19]

    Accardi et al., Eur

    A. Accardi et al., Eur. Phys. J. A 52, 268 (2016), arXiv:1212.1701 [nucl-ex]

    2016 arXiv

  12. [20]

    D. P. Anderle et al., Front. Phys. (Beijing) 16, 64701 (2021), arXiv:2102.09222 [nucl- ex]

    2021 arXiv

  13. [21]

    C. A. Bertulani, S. R. Klein, and J. Nystrand, Ann. Rev. Nucl. Part. Sci. 55, 271 (2005), arXiv:nucl-ex/0502005

    2005 arXiv

  14. [22]

    Xie and X

    Y.-P. Xie and X. Chen, Nucl. Phys. A 970, 316 (2018), arXiv:1805.06210 [hep-ph]

    2018 arXiv

  15. [23]

    X. Wang, Y. Xie, Y. Huang, and X. Chen, Phys. Rev. D 109, 016007 (2024), arXiv:2310.00980 [hep-ph]

    2024 arXiv

  16. [24]

    Lappi and H

    T. Lappi and H. Mantysaari, Phys. Rev. C 83, 065202 (2011), arXiv:1011.1988 [hep- ph]. 10

    2011 arXiv

  17. [25]

    Berger and W

    C. Berger and W. Wagner, Phys. Rept. 146, 1 (1987)

    1987

  18. [26]

    Abe et al

    T. Abe et al. (Belle-II), (2010), arXiv:1011.0352 [physics.ins-det]

    2010 arXiv

  19. [27]

    Abdul Khalek et al., Nucl

    R. Abdul Khalek et al., Nucl. Phys. A 1026, 122447 (2022), arXiv:2103.05419 [physics.ins-det]

    2022 arXiv

  20. [28]

    Aihara et al., (2024), arXiv:2406.19421 [hep-ex]

    H. Aihara et al., (2024), arXiv:2406.19421 [hep-ex]. 11

    2024

Pith tools

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