⚛️ phenomenology

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

This paper presents a relativistic study using the Bethe-Salpeter equation to predict the exclusive production cross sections of P-wave charmonia ( $\chi_{c1}+\eta_c$ and $h_c+h_c$ ) via two virtual photons in electron-positron annihilation at $\sqrt{s}=10.6$ GeV, yielding benchmark results of order $10^{-3}$ fb that exhibit a smooth power-law energy dependence and a hierarchy determined by meson wave function overlaps.

Original authors: Avinash Okram, Shashank Bhatnagar

Published 2026-07-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Avinash Okram, Shashank Bhatnagar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, high-speed dance floor where particles are constantly colliding and creating new partners. This paper is about a very specific, very rare dance move that happens when an electron and a positron (its antimatter twin) crash into each other at a particle collider called a "B-factory."

Here is the story of what the researchers found, explained without the heavy math.

The Setup: A Forbidden Dance

Usually, when these particles crash, they create new heavy particles (called "charmonia") by swapping a single messenger particle (a photon), kind of like two dancers passing a ball to create a new couple.

However, the researchers are looking at a special case where the new couples are made of P-wave particles. Think of these as dancers spinning in a very specific, complex way. Because of the rules of physics (specifically "charge parity"), a single ball-pass (one photon) is strictly forbidden for creating these specific couples.

Instead, the collision must happen through a double-pass mechanism. Two virtual photons (messengers that exist only for a split second) have to work together to rearrange the quarks inside the particles. It's like two dancers trying to swap partners simultaneously without ever letting go of each other's hands first. This is a "quark-rearrangement" dance, and it's incredibly difficult to pull off.

The Method: The "Relativistic Blueprint"

To predict how often this rare dance happens, the authors used a sophisticated mathematical tool called the Bethe-Salpeter equation.

The Analogy: Imagine trying to predict the shape of a spinning top. A simple model might just look at the top as a solid ball. But this team used a "relativistic blueprint" that accounts for the top spinning at near-light speeds, wobbling, and stretching. They didn't just guess the rules; they built their model using data from how these particles vibrate and weigh (spectroscopy), ensuring their blueprint was accurate before trying to predict the dance.

The Results: Tiny Numbers, Big Patterns

The team calculated the "cross-section," which is just a fancy way of saying "how likely is this dance to happen?"

The Numbers are Tiny: The chance of this happening is incredibly small—about 0.001 femtobarns. To put that in perspective, if you watched this collision happen a trillion times, you might only see this specific dance move a few times. It is so rare that current detectors might not even see it yet.
The Hierarchy (The Ranking): Even though the numbers are tiny, the team found a clear pattern in which dances are more likely than others.
- The "Ground State" Wins: Dances involving the simplest, most stable versions of the particles happen more often than dances involving excited, wobbly versions.
- The Surprise: For one specific pair (two identical "hc" particles), the team found that the "excited" version actually danced more often than the simple version. Why? Because the "wobbly" excited particles have a specific internal structure (a node, or a zero-point in their wave) that actually helps them grab onto the double-photon mechanism better in this specific scenario. It's like a dancer with a specific limp who, by accident, is perfectly suited for a specific complex step.

The Sensitivity Test: How Steady is the Prediction?

The researchers asked, "What if our measurements of the particle's weight or the strength of the force holding them together are slightly off?"

The "Spring" Matters Most: They found that their predictions were very sensitive to one specific parameter: the "spring constant" (a measure of how tightly the quarks are bound together). If you change how tight that spring is, the predicted number of dances changes dramatically.
The Other Factors: Changing the weight of the quark or the scale of the force had much less effect.
The Takeaway: This tells us that if we ever observe this dance in a real experiment, it will be a very precise test of how tightly these particles are bound together.

Comparison with Other Theories

The team compared their "relativistic blueprint" results with another popular theory called NRQCD (which treats the particles more like slow-moving balls).

Agreement: Both theories agree on the general ranking (Ground state is usually king).
Disagreement: The actual numbers differ. The "blueprint" method predicts slightly different rates because it accounts for the fact that these particles are moving so fast that they stretch and warp. The authors suggest their method might be more accurate for these high-speed, complex spins.

The Bottom Line

This paper provides a theoretical map for a dance that hasn't been seen yet. It predicts that:

These specific particle pairs are extremely rare.
They follow a predictable pattern based on their internal "shape" and how tightly they are bound.
Future, more powerful colliders (like the upgraded Belle II) might finally be sensitive enough to catch a glimpse of this rare event.

If they do catch it, it will confirm that our understanding of how heavy particles behave at near-light speeds is correct, specifically regarding how their internal "springiness" dictates their behavior.

Technical Summary: Exclusive Production of P-Wave Charmonia through Two Virtual Photons in $e^- + e^+$ Annihilation

Problem Statement
The paper addresses the theoretical prediction of exclusive double-charmonium production in electron-positron annihilation at a center-of-mass energy of $\sqrt{s} = 10.6$ GeV. Specifically, it investigates the processes $e^-e^+ \to \gamma^*\gamma^* \to \chi_{c1} + \eta_c$ and $e^-e^+ \to \gamma^*\gamma^* \to h_c + h_c$ . These channels are of particular interest because they are forbidden in single-photon annihilation due to charge conjugation ( $C$ -parity) and spin-parity constraints. Consequently, they proceed purely through two-virtual-photon exchange mechanisms involving quark-rearrangement diagrams, rather than the photon-fragmentation mechanisms that dominate other double-charmonium channels (such as $J/\Psi + J/\Psi$ ). The absence of fragmentation enhancements makes these processes a "clean" probe of relativistic quarkonium dynamics and the interplay between QED and QCD mechanisms at intermediate energies.

Methodology
The authors employ a fully relativistic $4 \times 4$ Bethe-Salpeter Equation (BSE) framework under the Covariant Instantaneous Ansatz (CIA). This approach offers a consistent linkage between low-energy spectroscopy and high-energy transition amplitudes without introducing ad-hoc parameters.

Key methodological features include:

Relativistic Formalism: The meson bound states are described by 4D BSE wave functions, reduced to 3D covariant forms using the CIA. The interaction kernel includes a one-gluon-exchange term and a confining potential with a harmonic oscillator-like structure.
Parameter Fixing: All model parameters (charm quark mass $m$ , QCD scale $\Lambda$ , and the flavor-independent spring constant $\omega_0$ ) are fixed exclusively from the mass spectroscopy of ground and excited charmonium states ( $0^{++}, 1^{--}, 0^{-+}, 1^{+-}, 1^{++}$ ) in previous works. No new parameters are introduced for the production calculations.
Amplitude Calculation: The leading-order amplitudes ( $O(\alpha^4)$ ) are calculated for the two quark-rearrangement diagrams (Direct and Exchange). The authors utilize heavy-quark approximations where internal relative momentum effects are suppressed ( $\hat{q}^2/m^2 \ll 1$ ), allowing photon virtualities to be approximated as $k^2 \approx k'^2 \approx s/4$ .
Wave Functions: The calculation utilizes 3D radial Salpeter wave functions for the $P$ -wave ( $\chi_{c1}, h_c$ ) and $S$ -wave ( $\eta_c$ ) states, including both ground ( $1P, 1S$ ) and excited ( $2P, 2S$ ) states. The normalization of these wave functions is strictly determined by current conservation conditions.

Key Contributions and Results
The paper provides leading-order predictions for the total cross sections and their energy dependence for the specified channels:

Cross Section Magnitudes: The predicted cross sections are extremely small, of the order of $10^{-3}$ fb.
- For $e^-e^+ \to \chi_{c1}(1P) + \eta_c(1S)$ , the cross section is $1.315 \times 10^{-3}$ fb.
- For $e^-e^+ \to h_c(1P) + h_c(1P)$ , the cross section is $0.113 \times 10^{-3}$ fb.
- The study also provides the first theoretical prediction for the excited state channel $e^-e^+ \to h_c(2P) + h_c(2P)$ , yielding $0.298 \times 10^{-3}$ fb.
Hierarchies and Wave Function Structure:
- $\chi_{c1} + \eta_c$ : The ground-state channel dominates ( $\sigma_{1P+1S} > \sigma_{2P+2S}$ ). This hierarchy arises because the production amplitude involves unweighted integrals of the radial wave functions ( $G_1, G_2$ ). The nodal structure of the excited states ( $2P, 2S$ ) leads to cancellations that reduce the overlap integrals compared to the ground states.
- $h_c + h_c$ : In contrast, the excited state channel is enhanced relative to the ground state ( $\sigma_{2P+2P} > \sigma_{1P+1P}$ ). This is attributed to the momentum-weighted nature of the $h_c$ production amplitude, which involves integrals of the form $\int \hat{q}^2 \phi(\hat{q}) d^3\hat{q}$ . The momentum weighting ( $\hat{q}^2$ ) preferentially samples the high-momentum region of the wave function, where the $2P$ state has significant support beyond its radial node, overcoming the nodal suppression seen in the unweighted case.
Parameter Sensitivity:
- The results are relatively stable with respect to variations in the charm quark mass ( $m$ ) and the QCD scale ( $\Lambda$ ), showing changes of only a few percent.
- The cross sections exhibit strong sensitivity to the confinement parameter $\omega_0$ . Specifically, $\sigma(h_c + h_c) \sim \omega_0^7$ , reflecting the high power of internal momentum in the double $P$ -wave amplitude. In contrast, $\sigma(\chi_{c1} + \eta_c) \sim \omega_0^3$ .
Energy Dependence: Both channels exhibit a smooth, monotonic power-law falloff ( $\sigma \sim 1/s^n$ ) as energy increases, consistent with dimensional counting rules for exclusive processes. This contrasts with some other models that predict non-monotonic behavior at lower energies.
Comparison with NRQCD:
- The qualitative hierarchy of channels (ground states dominating over excited states for $\chi_{c1}\eta_c$ ) is consistent with Non-Relativistic QCD (NRQCD) expectations.
- However, the absolute magnitudes differ. For $h_c + h_c$ , the BSE-CIA prediction ( $0.113 \times 10^{-3}$ fb) is lower than the NRQCD estimate ( $0.19 \times 10^{-3}$ fb). The authors attribute this suppression to the relativistic vertex structure in the BSE framework, which reduces the overlap integrals compared to the non-relativistic approximations used in NRQCD.

Significance and Claims
The paper positions these results as "benchmark predictions" for future searches at high-luminosity $e^-e^+$ colliders (such as Belle II). The authors claim that:

These channels provide a unique, clean test of relativistic quark-rearrangement dynamics, free from the complications of photon fragmentation.
The $h_c + h_c$ channel serves as a sensitive probe of spin-singlet $P$ -wave dynamics and the momentum-space structure of Bethe-Salpeter wave functions.
The study demonstrates that a unified BSE framework, with parameters fixed solely from spectroscopy, can yield self-consistent predictions for highly suppressed exclusive production channels.
While the absolute cross sections are below current experimental sensitivity, the internal hierarchies and energy dependencies derived from the relativistic bound-state structure offer distinct signatures that can be tested as experimental luminosity increases. The paper does not claim to resolve existing experimental discrepancies but rather provides a theoretical baseline for future observations.