Mass spectrum and decay widths of charmonium-like mesons: A diabatic approach with complex scaling

Resumen Técnico: Espectro de masas y anchuras de decaimiento de mesones tipo charmonio: un enfoque adiabático con escalamiento complejo

Planteamiento del problema
La observación experimental de numerosos estados tipo charmonio (estados XYZ) cerca de los umbrales de charm abierto desafía la descripción convencional de $q\bar{q}$. Estados como el $X(3872)$ exhiben masas y patrones de decaimiento inconsistentes con modelos de potencial simples. Si bien el potencial de tipo Cornell describe con éxito los estados de baja excitación, falla al proporcionar una descripción unificada para estados altamente excitados o cercanos al umbral, donde los efectos de canales acoplados y las contribuciones del continuo de mesón-mesón son significativos. Los marcos adiabáticos previos abordaron el espectro, pero a menudo descuidaron las contribuciones del continuo a las anchuras de decaimiento o trataron los estados ligados y resonantes por separado, requiriendo correcciones perturbativas o parámetros libres adicionales para estimar las anchuras.

Metodología
Este trabajo extiende un marco adiabático previo integrando el Método de Escalamiento Complejo (CSM, por sus siglas en inglés) para tratar los estados ligados y las resonancias sobre una misma base.

1. Formalismo Adiabático: El sistema se describe mediante una ecuación de Schrödinger multicanal con un potencial de matriz. El Hamiltoniano incluye un núcleo de quark-antiquark pesado ($c\bar{c}$) y canales de sabor abierto mesón-mesón. La matriz de potencial incorpora:

   • Términos Diagonales: Un potencial de Cornell dependiente del espín para el núcleo $c\bar{c}$ (que incluye términos de confinamiento, Coulomb, espín-espín, espín-órbita y tensor) y umbrales constantes para los canales mesón-mesón.
   • Términos Fuera de la Diagonal: Potenciales de mezcla ($V_{mix}$) que representan la ruptura de la cuerda, parametrizados como funciones gaussianas de la diferencia de energía entre el potencial $c\bar{c}$ y el umbral mesón-mesón.
   • Dependencia del Espín: Inclusión explícita de interacciones dependientes del espín derivadas del intercambio de un solo gluón, regularizadas mediante suavizado gaussiano para manejar las singularidades de corta distancia.

2. Método de Escalamiento Complejo (CSM): Para extraer las anchuras de decaimiento sin aproximaciones perturbativas, las coordenadas y los momentos se rotan hacia el plano complejo ($r \to re^{i\theta}$, $p \to pe^{-i\theta}$). Esto transforma el Hamiltoniano en una forma no hermítica donde los estados resonantes aparecen como autovalores discretos complejos ($E = M - i\Gamma/2$) separados del continuo rotado.

   • Solución Numérica: Las ecuaciones de canales acoplados se resuelven utilizando el Método de Expansión Gaussiana (GEM).
   • Observables: La parte real del autovalor proporciona la masa ($M$), y la parte imaginaria proporciona la anchura de decaimiento ($\Gamma$). El estudio también calcula radios cuadráticos medios complejos ($r_{RMS}$) y pesos de los componentes de los canales ($P_i$) utilizando el producto $c$ (producto biorortogonal) para analizar las estructuras internas.

3. Parámetros: El modelo utiliza parámetros fijados mediante el ajuste de estados de charmonio establecidos y resultados de QCD en la red (lattice QCD). La fuerza de la mezcla ($\Delta$) está restringida por la masa de la $\chi_{c1}(3872)$, mientras que el rango de mezcla ($\rho$) se fija basándose en estudios adiabáticos previos. No se introducen parámetros libres adicionales para el tratamiento del continuo.

Contribuciones Clave y Resultados
El estudio calcula el espectro de masas y las anchuras de decaimiento para estados tipo charmonio por debajo de 4.3 GeV en diversos sectores de $J^{PC}$:

• Sector $1^{++}$:

  • $\chi_{c1}(3872)$: Identificado como un polo cerca del umbral $D\bar{D}^*$ con una anchura de decaimiento ínfima. El análisis revela un peso dominante del componente $D\bar{D}^*$ (88.4%) y un radio RMS complejo grande, lo que respalda su interpretación como una molécula hadrónica débilmente ligada.
  • $X(3940)$: Se encuentra un polo en $M \approx 3941$ MeV con $\Gamma \approx 40$ MeV. Los autores lo proponen como un candidato para el $X(3940)$ observado experimentalmente con $J^{PC}=1^{++}$. El estado muestra una mezcla significativa entre configuraciones $c\bar{c}$ y $D\bar{D}^*$, y la asignación $1^{++}$ es consistente con la supresión observada del modo de decaimiento $D\bar{D}$.

• Sector $0^{++}$:

  • $\chi_{c0}(3860)$: Se identifica un polo en $M \approx 3865$ MeV con una anchura estrecha ($\Gamma \approx 4$ MeV). Aunque la masa se alinea bien con la $\chi_{c0}(3860)$ experimental, la anchura calculada es significativamente más estrecha que el valor central experimental ($\sim 200$ MeV). El estado es predominantemente $c\bar{c}$ (77%). El artículo señala que esta tensión sugiere que la naturaleza experimental de la $\chi_{c0}(3860)$ puede ser más compleja.
  • $X(3915)$: No se encuentra ningún polo correspondiente al $X(3915)$, lo cual es consistente con la visión de que no puede describirse como un estado $c\bar{c}$ convencional acoplado principalmente a canales de charm abierto.

• Sector $2^{++}$:

  • $\chi_{c2}(3930)$: La masa calculada ($\approx 3900$ MeV) y la anchura ($\approx 36$ MeV) concuerdan bien con la $\chi_{c2}(3930)$ experimental. El análisis de componentes indica una mezcla sustancial, con un peso de $c\bar{c}$ de $\sim 48\%$, resaltando la importancia de las nubes mesón-mesón.
  • Predicción: Se predice un segundo polo en $M \approx 4009$ MeV ($\Gamma \approx 17$ MeV), asignado tentativamente al estado $\chi_{c2}(1F)$.

• Sector $1^{--}$ (Vectorial):

  • El modelo reproduce la multiplicidad de los estados vectoriales observados ($\psi(2S)$, $\psi(3770)$, $\psi(4040)$, $\psi(4160)$, $\psi(4230)$).
  • $\psi(4040)$: Identificado con un componente molecular $D^*\bar{D}^*$ significativo (peso de $c\bar{c} \sim 34\%$).
  • $\psi(4230)$: Identificado con un componente molecular $D_s\bar{D}_s^*$ dominante, a pesar de la naturaleza no hermítica de los pesos.

• Otros Sectores: Los estados ligados para $\eta_c(1S, 2S)$ y $h_c(1P)$ se reproducen con alta precisión. Se realizan predicciones para resonancias en los sectores $0^{-+}$ ($\sim 4085$ MeV) y $1^{+-}$ ($\sim 3937$ MeV).

Significancia y Reivindicaciones
El artículo afirma que la integración del Método de Escalamiento Complejo en el marco adiabático proporciona una descripción unificada y libre de parámetros tanto de los estados ligados como de los resonantes tipo charmonio. Los logros clave son:

1. Extracción Directa de Anchuras: Las anchuras de decaimiento se obtienen directamente de los autovalores complejos sin correcciones perturbativas ni suposiciones ad-hoc.
2. Información Estructural: El cálculo de los pesos de los componentes complejos y los radios RMS elucida la estructura interna de los candidatos exóticos, confirmando la naturaleza molecular de la $\chi_{c1}(3872)$ y la naturaleza mixta de la $\chi_{c2}(3930)$.
3. Identificación de Candidatos: El trabajo proporciona un fuerte apoyo teórico para la asignación $1^{++}$ del $X(3940)$ y ofrece una identificación consistente en masa para la $\chi_{c0}(3860)$ como la $\chi_{c0}(2P)$, al tiempo que destaca las discrepancias en la anchura que requieren mayor investigación.
4. Limitaciones y Direcciones Futuras: Los autores señalan explícitamente que el modelo actual está restringido a canales de charm abierto. La ausencia de estados como el $X(3915)$, $\chi_{c1}(4140)$ y $\chi_{c1}(4274)$ en el espectro sugiere que los mecanismos que involucran canales de charm oculto (por ejemplo, $J/\psi\phi$, $\omega J/\psi$) son necesarios para describir completamente estos candidatos exóticos. El artículo concluye que la incorporación de estos canales de charm oculto es una extensión necesaria para una imagen física completa.