Spectral densities from Euclidean correlators via integral transforms: theoretical framework

Resumen Técnico: Densidades espectrales a partir de correladores euclidianos mediante transformadas integrales: marco teórico

Planteamiento del problema
Las densidades espectrales son fundamentales en la teoría cuántica de campos (QFT), vinculando las propiedades dinámicas con observables experimentales como las secciones eficaces hadrónicas inclusivas y las tasas de decaimiento semileptónico. En regímenes no perturbativos, particularmente en la Cromodinámica Cuántica (QCD), estas densidades están codificadas en la dependencia temporal euclidiana de las funciones de correlación a través de la descomposición espectral de Källén-Lehmann, que las relaciona mediante una transformada de Laplace. La extracción de la densidad espectral a partir de datos de QCD en el retículo (lattice QCD) requiere resolver la Transformada de Laplace Inversa (ILT), un problema notoriamente mal planteado.

Desafíos actuales:

1. Efectos de discretización: Los correladores en el retículo se definen en pasos de tiempo discretos, lo que introduce artefactos de orden $O(a)$ o $O(a^2)$.
2. Extensión temporal finita: Las simulaciones prácticas en el retículo están limitadas por un tiempo euclidiano finito $t_{max}$, lo que conduce a una falta de conocimiento del correlador en tiempos grandes. Esta truncación impide la reconstrucción exacta de la densidad espectral.
3. Ruido y resolución: A medida que la distancia-tiempo aumenta, la relación señal-ruido se deteriora y el espectro de estados se vuelve más denso en simulaciones de "campo maestro" (master-field), complicando la extracción de funciones espectrales continuas a partir de datos discretos.

Metodología
Los autores proponen un marco teórico basado en transformadas integrales para invertir analíticamente la relación de la transformada de Laplace. La estrategia central consiste en:

1. Formulación de la transformada integral: En lugar de una inversión numérica directa, los autores mapean la función de correlación $C(t)$ y la densidad espectral $\rho(\omega)$ en un espacio conjugado utilizando transformadas integrales específicas (Mellin, Kontorovich–Lebedev y Mehler-Fock).
2. Relación algebraica: Al elegir funciones de base que son eigenfunciones de operadores de (cuasi-)Carleman, la ecuación integral que relaciona $C(t)$ con $\rho(\omega)$ se convierte en una relación algebraica simple en el dominio de la transformada: $C_s = \lambda_s \hat{\rho}_s$.
3. Reconstrucción inversa: La densidad espectral se recupera mediante una integral de la transformada inversa: $\rho(\omega) = \int (C_s / \lambda_s) \hat{u}_s(\omega) ds$. Los autores abordan explícitamente la "cancelación delicada" entre el numerador y el denominador en esta integral, que típicamente decae exponencialmente.
4. Gestión de la discretización y la truncación:
   • Muestreo discreto: El marco se adapta para datos de retículo mediante la discretización de las transformadas integrales. Los autores demuestran que si el correlador del retículo está mejorado de orden $O(a)$, la densidad espectral extraída también está mejorada de orden $O(a)$, apareciendo los efectos de discretización dominantes solo a orden $O(a^2)$.
   • Transformadas incompletas: Para correladores conocidos solo hasta un $t_{max}$, los autores introducen transformadas integrales "incompletas". Establecen rigurosamente límites para el error sistemático introducido por la cola desconocida del correlador ($t > t_{max}$).
5. Densidades espectrales suavizadas (smeared): Para mitigar la inestabilidad de la reconstrucción de la densidad exacta, los autores se centran en densidades espectrales suavizadas (integrales de la densidad contra un núcleo $\kappa(\omega)$, por ejemplo, Breit-Wigner). Demuestran que si la transformada integral de la función de suavizado decae lo suficientemente rápido con la variable conjugada $s$, la contribución de la cola desconocida puede limitarse y controlarse.

Contribuciones clave y resultados

• Fórmulas analíticas: El artículo deriva fórmulas analíticas explícitas para extraer densidades espectrales reguladas y/o suavizadas tanto de correladores euclidianos continuos como discretos.
• Límites de error rigurosos: Una contribución principal es la derivación de límites superiores rigurosos para la incertidumbre sistemática derivada de la extensión temporal finita del retículo. Se demuestra que el error está suprimido exponencialmente con $t_{max}$, siempre que la transformada de la función de suavizado decaiga lo suficientemente rápido.
• Análisis de discretización: Los autores prueban que la solución propuesta en el retículo converge a su contraparte en el continuo con efectos de $O(a^2)$, siempre que el correlador del retículo de entrada esté mejorado de orden $O(a)$.
• Validación numérica: El marco se prueba en correladores de una sola exponencial y de múltiples exponenciales. Los resultados demuestran que:
  • La densidad espectral suavizada reconstruida converge al resultado exacto conforme $t_{max}$ aumenta.
  • Los límites de error derivados teóricamente coinciden con las desviaciones observadas en las pruebas numéricas.
  • El preacondicionamiento (sustracción de estados exponenciales líderes conocidos) mejora significamente la reconstrucción, aumentando efectivamente la brecha de masa relevante para el correlador residual.
• Generalizabilidad: El formalismo se presenta en un contexto general, aplicable a las transformadas de Mellin, Kontorovich–Lebedev y Mehler-Fock, y se extiende a densidades espectrales reguladas (donde $\rho(\omega)$ puede divergir en $\omega$ bajos).

Significado y pretensiones
El artículo pretende proporcionar un marco robusto y teóricamente fundamentado para la extracción directa de información espectral a partir de datos de QCD en el retículo. Su importancia radica en:

1. Control sistemático: Ofrece un método para cuantificar y limitar rigurosamente los errores sistemáticos asociados con la extensión temporal finita de los retículos, una limitación común en las simulaciones actuales.
2. Precisión: Al aprovechar correladores mejorados de orden $O(a)$ y transformadas integrales, el método evita introducir errores de discretización adicionales de orden $O(a)$ durante el proceso de inversión.
3. Utilidad práctica: El enfoque en densidades espectrales suavizadas se alinea con los observables físicos (como las secciones eficaces) que son naturalmente inclusivos, lo que hace que el método sea directamente aplicable a la fenomenología.
4. Base para trabajos futuros: Los autores declaran que la implementación práctica de esta estrategia, incluyendo su aplicación a datos reales de Monte Carlo multinivel, se detallará en un artículo complementario. El presente trabajo establece la maquinaria teórica necesaria y el análisis de error.

Los autores mantienen la modestia respecto a las aplicaciones experimentales inmediatas, posicionando este trabajo como un marco teórico para permitir cálculos no perturbativos más precisos en el futuro, particularmente a medida que las técnicas de retículo (como la integración multinivel) continúen mejorando la calidad de los datos de los correladores euclidianos.