Realization of all-to-all fermion propagator for the first principle high accuracy strong interaction prediction

Planteamiento del Problema

La Cromodinámica Cuántica en el retículo (LQCD) es una herramienta primaria para calcular la estructura de los hadrones desde los primeros principios, como el acoplamiento axial vectorial del nucleón ($g_A$). Sin embargo, alcanzar una precisión sub-porcentual requiere abordar dos desafíos mayores:

1. Costo Computacional de los Propagadores de Fermiones All-to-All: El cálculo de funciones de correlación de $N$ puntos (donde $N \ge 4$), que son esenciales para las correcciones de QED, la ruptura de la simetría de isospín (ISB) y los bucles de quarks desconectados, requiere propagadores de fermiones all-to-all. El cálculo exacto escala como $O(V^2)$ con el volumen del retículo $V$, lo que resulta prohibitivamente costoso.
2. Contaminación por Estados Excitados (ESC): La extracción de elementos de matriz del estado fundamental se ve complicada por la contaminación de estados excitados (por ejemplo, estados de tipo Roper, estados multifotónicos como $\pi N$). Aunque el método de destilación (proyección sobre modos de bajo momento) suprime eficazmente la ESC para los espectros de hadrones, descarta los modos de alto momento que son esenciales para las inserciones de corrientes locales, limitando su utilidad para funciones de $N$ puntos y elementos de matriz.

Las técnicas existentes, como las fuentes secuenciales, los enfoques estocásticos con sustitución de modos bajos (LMS) y el esparcimiento de campos, ofrecen soluciones parciales, pero tienen dificultades para manejar simultáneamente funciones de $N$ puntos de alta precisión, diagramas desconectados y la supresión de la ESC sin incurrir en costos computacionales masivos.

Metodología: El Algoritmo de "Blending" (Mezcla)

Los autores proponen un nuevo método estocástico llamado algoritmo de "blending" para aproximar eficientemente y sin sesgo los propagadores de fermiones all-to-all.

• Construcción del Espacio de Blending: El método construye un "espacio de blending" para un corte de tiempo fijo combinando dos subespacios:
  1. Modos de Baja Frecuencia (LFM): El subespacio $L_1$ generado por los $N_e$ autovalores más bajos del Laplaciano espacial covariante con la gauge (similar a la destilación). Este captura la física de larga distancia y permite la construcción de estados de hadrones espacialmente extendidos.
  2. Modos de Alta Frecuencia (SHFM): El complemento ortogonal $L_2$ se aproxima estocásticamente utilizando $N_{st}$ vectores aleatorios ($\eta_j$) muestreados uniformemente de $L_2$.
• Estimación Sin Sesgo: El operador identidad en el espacio completo se aproxima mediante una suma ponderada de proyectores sobre estos vectores de base.
  • Los proyectores de modo bajo retienen su peso exacto (1).
  • Los proyectores estocásticos de modo alto se asignan un peso de $\dim(L_2)/N_{st}$ para asegurar una estimación sin sesgo de la identidad en el subespacio complementario.
  • Para operadores que involucran productos de campos (por ejemplo, en funciones de $N$ puntos), se aplica un tensor de peso $\Omega^{(2)}_{ij}$ para corregir el doble conteo y asegurar que el estimador permanezca sin sesgo, particularmente para los términos estocásticos diagonales.
• Compresión: El propagador de quarks all-to-all $S$ se comprime proyectándolo en este espacio de blending. Esto reduce los requisitos de almacenamiento de $O((4N_c N_L^3 N_T)^2)$ a $O((4N_T(N_e + N_{st}))^2)$, permitiendo el cálculo flexible de funciones de correlación con operadores interpolantes y de corriente arbitrarios.
• Control de Estados Excitados: El método facilita el uso del análisis de Problema de Autovalores Generalizados (GEVP) con una base enriquecida de interpoladores (incluyendo operadores "involucrados en la corriente" como estados de dispersión $N f_1$) sin costos de inversión adicionales, ya que los propagadores comprimidos pueden reutilizarse para múltiples combinaciones de operadores.

Contribuciones Claves y Resultados

Los autores demuestran la eficiencia y precisión del método de blending a través de diversas aplicaciones en dos conjuntos de CLQCD: C24P29 ($m_\pi \approx 292$ MeV, $a \approx 0.105$ fm) y F64P13 ($m_\pi \approx 134$ MeV, $a \approx 0.077$ fm).

1. Validación de la Ausencia de Sesgo:

   • El método fue probado utilizando el elemento de matriz de la corriente vectorial conservada del pion ($\pi^+$). Los resultados mostraron que el uso de solo modos bajos ($N_{st}=0$) producía un resultado sesgado ($\approx 0.64$), mientras que la inclusión de modos estocásticos ($N_{st} \ge 50$) restauró el resultado a la expectativa teórica de 1.000(29).
   • Se realizaron comprobaciones de consistencia para el factor de forma eléctrico del pion ($f_{\pi\pi}(Q^2)$) y el radio de carga utilizando tanto funciones de correlación de 3 puntos como de 4 puntos. Los resultados son consistentes dentro de las incertidumbres, validando el cálculo de la función de 4 puntos.

2. Acoplamiento Axial del Nucleón ($g_A$):

   • Iso-vector ($g_A^{u-d}$): Utilizando un ajuste conjunto de 3 estados a través de varios conteos de modos bajos ($\bar{N}_e$) y separaciones de fuente-sumidero, los autores obtuvieron $g_A^{u-d} = 1.2337(84)$ en la masa física del pion. Este resultado elimina la ambigüedad asociada con la selección de un $\bar{N}_e$ específico y muestra una excelente concordancia con otros estudios de alta precisión.
   • Descomposición de Singlete y de Sabor: El método permitió el cálculo de bucles de quarks desconectados para la carga axial de singlete ($g_A^{u+d+s}$) y las contribuciones de sabor individuales.
     • Utilizando un análisis GEVP con una base que incluye interpoladores $N$, $N f_1$ y $N f'_1$, suprimieron eficazmente la contaminación de estados excitados.
     • Los resultados finales en el punto físico ($a=0.077$ fm) son:
       • $g_A^u = 0.8408(86)$
       • $g_A^d = -0.3929(86)$
       • $g_A^s = -0.0381(57)$
       • $g_A^{u+d+s} = 0.410(20)$

3. Eficiencia Computacional:

   • El método de blending reduce significativamente el número de inversiones del operador de Dirac requeridas en comparación con los métodos tradicionales de fuente secuencial o estocástica.
   • Para el isovector $g_A$, el costo se redujo aproximadamente en un orden de magnitud en comparación con estudios recientes de alta precisión (por ejemplo, ETM, RQCD, PNDME) que apuntan a una precisión similar.
   • Para el $g_A$ de singlete (que involucra bucles desconectados), el método alcanzó una precisión del 4% con 0.70 millones de inversiones, comparado con los 13.7–34.1 millones de inversiones en trabajos previos comparables. Esta eficiencia proviene de la capacidad de reutilizar el propagador comprimido para múltiples interpoladores e inserciones de corriente sin nuevas inversiones.

Significado y Reivindicaciones

El artículo afirma que el método de blending proporciona un marco robusto para cálculos de alta precisión de la estructura de los hadrones en LQCD, abordando específicamente los dobles desafíos del costo computacional y la contaminación de estados excitados.

• Control Sistemático: Al permitir la construcción de funciones de correlación con múltiples operadores interpolantes que comparten los mismos números cuánticos, el método facilita la extracción robusta de elementos de matriz del estado fundamental mediante el análisis variacional (GEVP), suprimiendo eficazmente la ESC.
• Escalabilidad: En la masa física del pion, el método hace que las predicciones de precisión sub-porcentual para los elementos de matriz del nucleón sean alcanzables con recursos computacionales significativamente reducidos.
• Versatilidad: El enfoque es aplicable a diversas funciones de $N$ puntos, incluyendo diagramas desconectados, correcciones de QED y efectos de ISB, los cuales son cruciales para probar el Modelo Estándar.
• Limitaciones: Los autores señalan que el costo computacional sigue siendo prohibitivo para los fermiones Overlap debido a sus mayores costos de inversión, lo que sugiere la necesidad de más optimizaciones para una aplicabilidad más amplia.

En resumen, el método de blending representa un paso significativo para hacer factibles los cálculos de alta precisión de observables hadrónicos complejos desde los primeros principios, permitiendo futuros estudios de dispersión de muchos cuerpos y física de partones con reducidas incertidumbres estadísticas y sistemáticas.