Handling the Cornell potential within the Lagrange-mesh method in momentum space

A continuación se presenta un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Manejo del potencial de Cornell dentro del método de malla de Lagrange en el espacio de momentos

Autores: Cyrille Chevalier y Joachim Viseur
Institución: Universidad de Mons, Bélgica.

---

1. El Problema

El método de malla de Lagrange (LMM, por sus siglas en inglés) es una herramienta numérica altamente eficiente y precisa para resolver ecuaciones de tipo Schrödinger en problemas de dos cuerpos, especialmente en física hadrónica. Históricamente, el LMM se ha desarrollado y aplicado con gran éxito en el espacio de configuraciones (posición).

Sin embargo, para muchos sistemas físicos, particularmente aquellos que involucran cinemática exótica, interacciones dependientes del momento o descripciones basadas en teorías de campo cuántico, es más natural y ventajoso trabajar en el espacio de momentos.

El desafío principal identificado en este trabajo es que la metodología previa para aplicar el LMM en el espacio de momentos (desarrollada por G. Lacroix et al., Ref. [11]) es numéricamente inadecuada para manejar potenciales de largo alcance, específicamente:

• El potencial de Coulomb ($1/r$).
• El potencial lineal ($r$).

Estos potenciales son fundamentales en la fenomenología de hadrones, especialmente a través de interacciones tipo Cornell ($V(r) = -\kappa/r + ar + C$), que modelan la interacción entre quarks y gluones. El método anterior falla porque la transformada de Fourier de estos potenciales conduce a funciones que divergen en los puntos de la malla (diagonal $i=j$), haciendo imposible el cálculo de los elementos de matriz.

2. Metodología Propuesta

Los autores proponen una metodología alternativa para calcular los elementos de matriz del potencial en el espacio de momentos, inspirada en cómo se maneja la cinemática exótica en el LMM de espacio de configuraciones.

Puntos clave de la nueva metodología:

1. Inversión de la estrategia: En lugar de intentar calcular directamente la transformada de Fourier del potencial (lo cual falla para Coulomb/lineal), el método calcula primero los elementos de matriz del operador $r^2$ en la base de funciones de Lagrange en el espacio de momentos.
2. Cálculo de $r^2$: Dado que $r^2$ actúa como un operador diferencial en el espacio de momentos (análogo a $p^2$ en el espacio de posición), se utilizan fórmulas ya establecidas para calcular la matriz $R^2_{ij} = \langle f_i | r^2 | f_j \rangle$ mediante cuadraturas de Gauss-Laguerre.
3. Diagonalización: La matriz $R^2$ se diagonaliza para obtener sus autovalores ($D^2$) y la matriz de transición ($S$).
4. Evaluación del Potencial: En la base de autoestados de $r^2$, el operador potencial $V(r^2)$ es diagonal. Se calcula simplemente aplicando la función $V$ a los autovalores de $r^2$: $V_{D^2} = V(D^2)$.
5. Transformación de vuelta: Finalmente, se transforma la matriz de potencial de vuelta a la base original de Lagrange utilizando la relación de similitud: $V_{R^2} = S V_{D^2} S^{-1}$.

Este enfoque evita las divergencias directas asociadas con la transformada de Fourier de potenciales singulares, permitiendo un tratamiento numérico estable.

3. Contribuciones Clave

• Extensión del rango de potenciales: Se ha logrado adaptar el LMM en el espacio de momentos para incluir interacciones de Coulomb y lineales, que anteriormente eran inaccesibles con este método en dicha representación.
• Tratamiento eficiente del potencial de Cornell: Se demuestra que es posible resolver ecuaciones de tipo Schrödinger con el potencial de Cornell (crucial para modelos de mesones) manteniendo una descripción estrictamente en el espacio de momentos.
• Cálculo de densidades de probabilidad: Se desarrollaron fórmulas explícitas para obtener las densidades de probabilidad tanto en el espacio de momentos como en el de configuraciones (mediante transformada de Fourier numérica) a partir de los coeficientes de la solución en el espacio de momentos.
• Validación rigurosa: Se compararon los resultados con soluciones analíticas exactas (sistema Coulombiano) y con resultados de la literatura (modelo de mesones), demostrando la precisión del método.

4. Resultados

• Validación con Potencial de Coulomb:
  • Se resolvió el sistema de dos partículas no relativistas con potencial de Coulomb.
  • Con solo 50 puntos de malla, el método reproduce al menos tres dígitos significativos de las energías exactas. Con 300 puntos, la precisión supera los cinco dígitos.
  • Las densidades de probabilidad en el espacio de momentos convergen rápidamente a la solución analítica a medida que aumenta el tamaño de la malla ($N$).
• Aplicación al Modelo de Mesones (Potencial de Cornell):
  • Se aplicó el método a un modelo de mesones semirrelativistas con el potencial de Cornell.
  • Los autovalores de energía obtenidos coinciden con los reportados en la literatura original (Ref. [3]) y con cálculos de LMM en espacio de configuraciones (Ref. [10]), alcanzando una precisión de hasta cuatro dígitos.
  • Aunque la convergencia es ligeramente más lenta que en el espacio de configuraciones (debido a la elección de un parámetro de escala $h$ fijo en lugar de optimizado), los resultados son consistentes y válidos.
• Comparación con métodos previos:
  • En el caso de potenciales de corto alcance (Gaussiano), el nuevo método ofrece una precisión comparable a la del método anterior (Ref. [11]), confirmando que la nueva aproximación no degrada el rendimiento en casos donde el método antiguo funcionaba.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

1. Versatilidad: Permite a los físicos utilizar la representación en el espacio de momentos para una gama mucho más amplia de problemas, incluyendo aquellos dominados por el confinamiento (potencial lineal) y la interacción electromagnética (Coulomb).
2. Física Hadrónica: Facilita el estudio de modelos de quarks y gluones donde las interacciones dependientes del momento o la cinemática semirrelativista son críticas, sin sacrificar la precisión numérica.
3. Complementariedad: Ofrece una alternativa robusta al LMM en espacio de configuraciones. Dependiendo de la naturaleza específica del problema (por ejemplo, si se prioriza la interpretación de momentos o la facilidad de manejo de ciertos operadores), los investigadores ahora tienen la libertad de elegir la representación más conveniente sin perder capacidad de cálculo.
4. Resolución de Divergencias: Proporciona una solución elegante a un problema numérico de larga data (la divergencia de los elementos diagonales en la representación de momentos para potenciales singulares), utilizando la diagonalización de un operador auxiliar ($r^2$) como puente matemático.

En resumen, el artículo establece un nuevo estándar para la aplicación del método de malla de Lagrange en el espacio de momentos, eliminando una barrera técnica importante y abriendo nuevas vías para la investigación en física de partículas y nuclear.