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

Resumen Técnico: Producción Exclusiva de Charmonia de Onda P mediante Dos Fotones Virtuales en la Aniquilación $e^- + e^+$

Planteamiento del Problema
El artículo aborda la predicción teórica de la producción de doble charmonio exclusivo en la aniquilación electrón-positrón a una energía de centro de masa de $\sqrt{s} = 10.6$ GeV. Específicamente, investiga los procesos $e^-e^+ \to \gamma^*\gamma^* \to \chi_{c1} + \eta_c$ y $e^-e^+ \to \gamma^*\gamma^* \to h_c + h_c$. Estos canales son de particular interés porque están prohibidos en la aniquilación de un solo fotón debido a las restricciones de paridad de carga conjugación ($C$) y de espín-paridad. En consecuencia, proceden puramente a través de mecanismos de intercambio de dos fotones virtuales que involucran diagramas de reordenamiento de quarks, en lugar de los mecanismos de fragmentación de fotones que dominan otros canales de doble charmonio (como $J/\Psi + J/\Psi$). La ausencia de aumentos por fragmentación convierte a estos procesos en una sonda "limpia" de la dinámica de quarks relativistas y de la interacción entre los mecanismos de QED y QCD a energías intermedias.

Metodología
Los autores emplean un marco de la Ecuación de Bethe-Salpeter (BSE) de $4 \times 4$ totalmente relativista bajo la Ansatz Instantánea Covariante (CIA). Este enfoque ofrece un vínculo consistente entre la espectroscopía de baja energía y las amplitudes de transición de alta energía sin introducir parámetros ad-hoc.

Características metodológicas clave incluyen:

• Formalismo Relativista: Los estados ligados de los mesones se describen mediante funciones de onda de la BSE 4D, reducidas a formas covariantes 3D utilizando la CIA. El núcleo de interacción incluye un término de intercambio de un solo gluón y un potencial de confinamiento con una estructura de tipo oscilador armónico.
• Fijación de Parámetros: Todos los parámetros del modelo (masa del quark charm $m$, escala QCD $\Lambda$ y la constante de muelle independiente del sabor $\omega_0$) se fijan exclusivamente a partir de la espectroscopía de masas de los estados de charmonio fundamentales y excitados ($0^{++}, 1^{--}, 0^{-+}, 1^{+-}, 1^{++}$) en trabajos previos. No se introducen nuevos parámetros para los cálculos de producción.
• Cálculo de Amplitudes: Las amplitudes de orden líder ($O(\alpha^4)$) se calculan para los dos diagramas de reordenamiento de quarks (Directo y de Intercambio). Los autores utilizan aproximaciones de quarks pesados donde los efectos del momento relativo interno están suprimidos ($\hat{q}^2/m^2 \ll 1$), permitiendo que las virtualidades de los fotones se aproximen como $k^2 \approx k'^2 \approx s/4$.
• Funciones de Onda: El cálculo utiliza funciones de onda de Salpeter 3D radiales para los estados de onda $P$ ($\chi_{c1}, h_c$) y onda $S$ ($\eta_c$), incluyendo tanto estados fundamentales ($1P, 1S$) como excitados ($2P, 2S$). La normalización de estas funciones de onda está estrictamente determinada por las condiciones de conservación de la corriente.

Contribuciones Clave y Resultados
El artículo proporciona predicciones de orden líder para las secciones eficaces totales y su dependencia energética para los canales especificados:

1. Magnitudes de la Sección Eficaz: Las secciones eficaces predichas son extremadamente pequeñas, del orden de $10^{-3}$ fb.

   • Para $e^-e^+ \to \chi_{c1}(1P) + \eta_c(1S)$, la sección eficaz es $1.315 \times 10^{-3}$ fb.
   • Para $e^-e^+ \to h_c(1P) + h_c(1P)$, la sección eficaz es $0.113 \times 10^{-3}$ fb.
   • El estudio también proporciona la primera predicción teórica para el canal del estado excitado $e^-e^+ \to h_c(2P) + h_c(2P)$, arrojando $0.298 \times 10^{-3}$ fb.

2. Jerarquías y Estructura de la Función de Onda:

   • $\chi_{c1} + \eta_c$: El canal del estado fundamental domina ($\sigma_{1P+1S} > \sigma_{2P+2S}$). Esta jerarquía surge porque la amplitud de producción involucra integrales no ponderadas de las funciones de onda radiales ($G_1, G_2$). La estructura de nodos de los estados excitados ($2P, 2S$) conduce a cancelaciones que reducen las integrales de traslape en comparación con los estados fundamentales.
   • $h_c + h_c$: En contraste, el canal del estado excitado es aumentado respecto al estado fundamental ($\sigma_{2P+2P} > \sigma_{1P+1P}$). Esto se atribuye a la naturaleza ponderada por el momento de la amplitud de producción de $h_c$, que involucra integrales de la forma $\int \hat{q}^2 \phi(\hat{q}) d^3\hat{q}$. La ponderación por el momento ($\hat{q}^2$) muestrea preferencialmente la región de alto momento de la función de onda, donde el estado $2P$ tiene un soporte significativo más allá de su nodo radial, superando la supresión por nodos vista en el caso no ponderado.

3. Sensibilidad de los Parámetros:

   • Los resultados son relativamente estables con respecto a las variaciones en la masa del quark charm ($m$) y la escala QCD ($\Lambda$), mostrando cambios de solo unos pocos por ciento.
   • Las secciones eficaces exhiben una fuerte sensibilidad al parámetro de confinamiento $\omega_0$. Específicamente, $\sigma(h_c + h_c) \sim \omega_0^7$, lo que refleja el alto exponente del momento interno en la amplitud doble de onda $P$. En contraste, $\sigma(\chi_{c1} + \eta_c) \sim \omega_0^3$.

4. Dependencia Energética: Ambos canales exhiben una caída de ley de potencia suave y monotónica ($\sigma \sim 1/s^n$) a medida que aumenta la energía, consistente con las reglas de conteo dimensional para procesos exclusivos. Esto contrasta con algunos otros modelos que predicen un comportamiento no monotónico a energías más bajas.

5. Comparación con NRQCD:

   • La jerarquía cualitativa de los canales (estados fundamentales dominando sobre los excitados para $\chi_{c1}\eta_c$) es consistente con las expectativas de la QCD no relativista (NRQCD).
   • Sin embargo, las magnitudes absolutas difieren. Para $h_c + h_c$, la predicción BSE-CIA ($0.113 \times 10^{-3}$ fb) es inferior a la estimación de NRQCD ($0.19 \times 10^{-3}$ fb). Los autores atribuyen esta supresión a la estructura del vértice relativista en el marco de la BSE, que reduce las integrales de traslape en comparación con las aproximaciones no relativistas utilizadas en NRQCD.

Significancia y Reivindicaciones
El artículo posiciona estos resultados como "predicciones de referencia" para futuras búsquedas en colisionadores $e^-e^+$ de alta luminosidad (como Belle II). Los autores afirman que:

• Estos canales proporcionan una prueba única y limpia de la dinámica de reordenamiento de quarks relativistas, libre de las complicaciones de la fragmentación de fotones.
• El canal $h_c + h_c$ sirve como una sonda sensible de la dinámica de onda $P$ de espín singlete y de la estructura en el espacio de momentos de las funciones de onda de Bethe-Salpeter.
• El estudio demuestra que un marco de BSE unificado, con parámetros fijados únicamente a partir de la espectroscopía, puede generar predicciones consistentes para canales de producción exclusiva altamente suprimidos.
• Aunque las secciones eficaces absolutas están por debajo de la sensibilidad experimental actual, las jerarquías internas y las dependencias energéticas derivadas de la estructura de estado ligado relativista ofrecen firmas distintas que pueden ser probadas a medida que la luminosidad experimental aumente. El artículo no pretende resolver las discrepancias experimentales existentes, sino que proporciona una línea base teórica para observaciones futuras.