Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

Resumen Técnico: Simulación Cuántica de Distribuciones de Partones Generalizadas en el Modelo de Schwinger

Planteamiento del Problema
Las Distribuciones de Partones Generalizadas (GPDs, por sus siglas en inglés) son fundamentales para comprender la estructura interna de los hadrones, codificando correlaciones entre los parámetros de impacto transversales y el momento longitudinal. Sin embargo, el cálculo de las GPDs desde los primeros principios en Cromodinámica Cuántica (QCD) se ve obstaculizado por el problema de la signatura en la QCD de red euclídea, lo que impide el cómputo directo de las funciones de correlación en el tiempo real en la luz (light-cone). Mientras que los enfoques clásicos como la Teoría de Momento Efectivo de Gran Momento (LMET) ofrecen soluciones alternativas mediante quasi-PDFs, estos requieren un impulso de momento extremo y sufren de correcciones de potencia. La computación cuántica ofrece un cambio de paradigma al permitir la simulación directa de la dinámica hamiltoniana y la evolución en el tiempo real, superando potencialmente el problema de la signatura. A pesar de este potencial, la simulación de las GPDs sigue siendo un desafío formidable, particularmente en lo que respecta a la preservación de las simetrías discretas (Conjugación de Carga $C$, Paridad $P$, Reversión Temporal $T$) y la preparación eficiente de estados hadrónicos con momento no nulo.

Metodología
Los autores proponen un algoritmo cuántico integral para simular las GPDs en el modelo de Schwinger (electrodinámica cuántica) en (1+1) dimensiones, que sirve como un sustituto (proxy) de la QCD completa. La metodología se distingue por varias elecciones técnicas clave:

1. Fermiones de Wilson: A diferencia de las simulaciones cuánticas actuales que utilizan predominantemente fermiones escalonados (staggered fermions), este trabajo emplea fermiones de Wilson para la discretización en red. Esta elección es crítica para la preservación estricta de la simetría exacta de conjugación de carga ($C$) en la red. Los autores argumentan que la simetría $C$ es vital para estudiar mesones de carga neutra, ya que elimina contaminaciones no físicas de la parte real en las funciones de correlación y asegura que las GPDs sigan siendo funciones impares de la fracción de momento $x$.
2. Formulación Hamiltoniana: El modelo de Schwinger se formula en el calibre axial. El Hamiltoniano continuo se discretiza utilizando fermiones de Wilson, aumentados por un término de Wilson para evitar el problema del doble de fermiones (fermion doubling) mientras se preservan las simetrías $C$, $P$ y $T$. El campo de gauge temporal se elimina mediante la ley de Gauss, lo que resulta en un término de interacción no local en el Hamiltoniano.
3. Preparación de Estados: El algoritmo utiliza un Algoritmo Cuántico Variacional de Resolución de Número de Estado Cuántico (VQE). Para preparar estados hadrónicos con momentos no nulos específicos ($|h(k)\rangle$) y el vacío ($|\Omega\rangle$), los autores construyen un conjunto de estados de referencia ortogonales. Un circuito de inyección de momento, $D(k)$, mapea un estado de referencia a un estado con momento definido. El ansatz del VQE, $U(\theta)$, está diseñado para preservar la simetría, optimizando una función de pérdida multi-estado para extraer simultáneamente el vacío y los estados excitados hadrónicos dentro del sector sin carga ($Q=0$).
4. Medición de la Función de Correlación: El núcleo del algoritmo evalúa la función de correlación de la luz invariante de gauge, que incluye líneas de Wilson. Los autores implementan una versión controlada del operador de traslación de momento y enlaces de gauge controlados ($W_0$) utilizando un qubit auxiliar. Al medir los valores de expectativa fuera de la diagonal ($\sigma_x$ y $\sigma_y$) del qubit auxiliar después de la evolución en el tiempo real, se extraen las funciones de correlación espacial $\tilde{H}_{q/h}(N, \xi)$.
5. Análisis de Complejidad: Los requisitos de recursos se analizan utilizando la fórmula de producto de Suzuki-Trotter de segundo orden para la evolución temporal. La complejidad total escala polinomialmente con el tamaño del sistema $N$ y la precisión inversa $\epsilon$, específicamente como $O(N^p + N^{5/2}/\epsilon^{1/2})$.

Contribuciones Clave

• Marco Unificado: El artículo establece un marco algorítmico cuántico unificado capaz de computar tanto las Funciones de Distribución de Partones (PDFs) como las GPDs, en lugar de tratarlas como evaluaciones disjuntas.
• Preservación de Simetría: Demuestra que la implementación de fermiones de Wilson en un marco de computación cuántica no conlleva una mayor complejidad computacional en comparación con los fermiones escalonados, pero ofrece propiedades de simetría superiores. Específicamente, preserva rigurosamente la simetría $C$ exacta, la cual se muestra esencial para obtener resultados físicamente significativos para mesones neutros.
• Construcción Algorítmica: El trabajo detalla la construcción de circuitos cuánticos específicos para la inyección de momento ($D(k)$), la implementación del enlace de gauge y el ansatz VQE que preserva la simetría requerido para el modelo de Schwinger.

Resultidos
Los autores evaluaron su enfoque utilizando Diagonalización Exacta (ED) en una red espacial de tamaño $N=18$ (y hasta $N=22$ para verificaciones específicas):

• Espectro de Masa: El espectro de masa calculado del estado ligado $q\bar{q}$ más ligero reprodujo el límite continuo analítico conocido ($m_h/g \approx 0.564$) mediante extrapolación lineal, validando la formulación hamiltoniana.
• Funciones de Correlación: Se extrajeron las funciones de correlación espacial de la luz $\tilde{H}_{q/h}(N, \xi)$. Los resultados confirmaron que la parte real de la función de correlación desaparece asintóticamente a medida que el volumen de la red aumenta, lo cual es consistente con la simetría $C$.
• Comparación con Fermiones Escalonados: Una comparación directa mostró que la razón entre las partes real e imaginaria de la función de correlación está significativamente más suprimida para los fermiones de Wilson en comparación con los fermiones escalonados. Esto confirma que la formulación de Wilson elimina eficazmente las contaminaciones no físicas derivadas de los efectos de volumen finito y las violaciones de simetría.
• GPDs y PDFs: Las transformadas de Fourier de las funciones de correlación produjeron las GPDs en el espacio de momentos. Las PDFs resultantes (límite $\xi=0$) presentan un pico en $x=0.5$, reflejando correctamente la naturaleza de dos quarks de valencia del estado ligado en 1+1 dimensiones. Se encontró que las GPDs son funciones impares de $x$, satisfaciendo el requisito teórico de que el factor de forma de carga elástica para un estado de carga neutra debe ser nulo.

Significancia y Reivindicaciones
El artículo afirma que este trabajo constituye un precursor vital para las simulaciones de QCD en (3+1) dimensiones. Al demostrar que los fermiones de Wilson pueden implementarse eficientemente en hardware cuántico sin aumentar la complejidad, los autores argumentan que este enfoque ofrece un camino robusto hacia la simulación de fenómenos no perturbativos mientras se adhiere estrictamente a las simetrías fundamentales. La escala polinomial del algoritmo con respecto al tamaño del sistema y la precisión resalta la ventaja teórica de las computadoras cuánticas para simular directamente las correlaciones de la luz. Los autores posicionan este trabajo como un paso fundacional hacia la generalización de estos algoritmos a teorías de gauge no abelianas (por ejemplo, Yang-Mills SU(2) o SU(3)) y, finalmente, a la QCD completa en futuros hardware cuánticos tolerantes a fallos.