Lattice Quantum Chromodynamics for Quantum Simulations

Resumen Técnico: Cromodinámica Cuántica de Red para Simulaciones Cuánticas

Planteamiento del Problema
La Cromodinámica Cuántica (QCD) es la teoría fundamental de la interacción fuerte, pero su régimen no perturbativo —esencial para comprender la estructura de los hadrones y la materia de alta densidad— sigue siendo computacionalmente desafiante. Los enfoques estándar de QCD en red dependen de integrales de trayectoria euclidianas y muestreo por importancia, los cuales fallan en presencia de un "problema de signo", como ocurre en la dinámica en tiempo real o cuando se incluye un ángulo $\theta$. Si bien las formulaciones hamiltonianas ofrecen una vía para la simulación en tiempo real, el crecimiento exponencial del espacio de Hilbert con el volumen de la red hace que la simulación clásica sea intratable para sistemas grandes. En consecuencia, existe una necesidad crítica de marcos que mapeen las teorías de gauge en la red sobre computadoras cuánticas, abordando específicamente la codificación de los grados de libertad de gauge y fermiónicos, la construcción de estados físicos (invariantes de gauge) y el desarrollo de algoritmos eficientes de evolución temporal.

Metodología
Los autores desarrollan un marco integral para simulaciones cuánticas de la teoría de gauge $SU(N_c)$ con quarks, centrándose en $N_c=3$ (QCD) en dos y tres dimensiones espaciales. La metodología se basa en los siguientes pilares:

1. Formulación Hamiltoniana y Discretización:
   El trabajo utiliza la formulación del Hamiltoniano de Kogut-Susskind. Los autores consideran tanto fermiones staggered como de Wilson para abordar el problema del doble de fermiones (fermion doubling). La teoría se discretiza en una red espacial cúbica con tiempo continuo. El Hamiltoniano se descompone en términos eléctricos ($H_E$), magnéticos ($H_B$), de theta ($H_\theta$), cinéticos ($H_{kin}$) y de masa ($H_{mass}$). El término $\theta$ se incluye como un acoplamiento a la densidad de carga topológica, introduciendo coeficientes complejos en el Hamiltoniano.

2. Base de Representación y Espacio de Hilbert Físico:
   Una contribución central es la formulación del espacio de Hilbert físico en una base de representación. En lugar de mapear elementos de grupo directamente a qubits, los grados de libertad de gauge en los enlaces y los grados de libertad fermiónicos en los sitios se codifican utilizando representaciones irreducibles (irreps) de $SU(N_c)$.

   • Espacio de Hilbert de los Enlaces: Los enlaces se describen mediante estados $|(R, r_R, r_L)\rangle$, donde $R$ es una irrep y $r_{R,L}$ son vectores de base para la mitad derecha e izquierda de los enlaces. Esta base diagonaliza los operadores de Casimir cuadráticos.
   • Espacio de Hilbert de los Sitios: Los fermiones se tratan mediante una base de número de ocupación, que luego se transforma en una base de representación utilizando coeficientes de Clebsch-Gordan (CGC). Esto permite que los estados de partícula y antipartícula sean clasificados por sus propiedades de transformación bajo $SU(N_c)$.
   • Invariancia de Gauge: Los estados físicos se construyen como singletes (representaciones triviales) formados por el producto tensorial de todas las irreps que convergen en un sitio de la red; los autores definen explícitamente la convención de "orden de sitio" necesaria para construir de forma inequívoca estos estados singlete, asegurando que se cumpla la restricción de la ley de Gauss.

3. Truncamiento y Codificación:
   Dado que el espacio de Hilbert de $SU(N_c)$ es de dimensión infinita, los autores implementan esquemas de truncamiento. Estos incluyen cortes en el peso más alto de las irreps ($n_1$) o en el autovalor del Casimir cuadrático ($C_2$) para los enlaces de gauge, y límites en los números de ocupación de color para los fermiones. Los estados físicos se codifican en qubits utilizando un esquema de codificación local, donde los qubits representan irreps específicas e índices de multiplicidad de singlete. Se señala que este enfoque es más eficiente en términos de qubits que las codificaciones directas para espacios de Hilbert grandes.

4. Algoritmos de Simulación Cuántica:
   Los autores proponen un algoritmo de Trotterización para la evolución temporal. El operador de evolución temporal $e^{-iHt}$ se aproxima como un producto de unitarias correspondientes a los términos descompuestos del Hamiltoniano.

   • Construcción de Puertas: La evolución de los operadores de plaqueta y de enlace se mapea a rotaciones de Givens (unitarias de dos niveles) y desplazamientos de fase. Los autores detallan la construcción de circuitos cuánticos para estas operaciones, incluyendo el manejo de las cadenas de Jordan-Wigner (cambios de fase) requeridos para las relaciones de anticonmutación fermiónica.
   • Elementos de Matriz: El artículo proporciona fórmulas explícitas para los elementos de matriz de todos los operadores del Hamiltoniano, expresados como productos de "factores de sitio" que involucran CGCs. Estas fórmulas son esenciales para calcular los ángulos de rotación y las fases requeridas para las puertas cuánticas.

Contribuciones Clave

• Marco para $SU(3)$ con Ángulo $\theta$: El artículo presenta la primera formulación y simulación de QCD en red con un ángulo $\theta$ en tres dimensiones espaciales utilizando un enfoque hamiltoniano en una computadora cuántica.
• Formalismo de Base de Representación: Proporciona una construcción rigurosa del espacio de Hilbert físico utilizando irreps y singletes de CGC, ofreciendo una alternativa escalable a las bases de elementos de grupo.
• Detalles Algorítmicos: El trabajo detalla los circuitos cuánticos específicos para la evolución Trotterizada, incluyendo la descomposición de términos complejos del Hamiltoniano (como el término $\theta$ y los términos cinéticos de Wilson) en puertas elementales.
• Cálculos de Elementos de Matriz: Los autores derivan y tabulan los elementos de matriz para los operadores cinéticos de plaqueta, theta y fermión, proporcionando los datos necesarios para implementar los circuitos de simulación.

Resultados
Los autores realizaron simulaciones sin ruido en redes pequeñas (hasta una red periódica de $2\times2$ y un cubo de frontera abierta de $2\times2\times2$) utilizando hasta 32 qubits en simuladores clásicos (Qiskit Aer y CUDA-Q). Las simulaciones demostraron:

• Efectos de Truncamiento del Espacio de Hilbert: Análisis de cómo diferentes truncamientos de gauge y fermiones afectan observables como el campo eléctrico promedio ($\langle \Pi^2 \rangle$) y el número de ocupación de color ($\langle N_q \rangle$). Los resultados indicaron que el acoplamiento fuerte suprime los estados de alta población, haciendo que los truncamientos modestos sean efectivos para los fermiones staggered, mientras que los fermiones de Wilson mostraron una mayor sensibilidad al truncamiento.
• Dependencia de Parámetros: Se estudió la evolución de los observables en función de la masa fermiónica y el parámetro de Wilson $r$. Se encontró que las masas más pesadas suprimen las oscilaciones hacia estados ricos en fermiones.
• Dinámica del Ángulo Theta: Las simulaciones variando el ángulo $\theta$ revelaron un comportamiento aperiódico en el régimen de acoplamiento fuerte, con dinámicas distintas para $\theta$ positivo y negativo.
• Estados Hadrónicos y Dinámica de Cuerdas: Los autores demostraron la preparación de estados de "pion cargado" y la evolución de estados de cuerda (tanto de bucle como de cuerda fermiónica). Observaron la dinámica de ruptura de la cuerda, donde los tubos de flujo decaen en pares fermión-antifermión, y rastrearon la probabilidad de ruptura de la cuerda frente a configuraciones no rotas.
• Potencial Químico Bariónico: Utilizando diagonalización exacta en un solo cubo, los autores computaron el costo energético de añadir un barión, sirviendo como prueba de concepto para el estudio de la materia de alta densidad.

Significado y Reivindicaciones
El artículo se posiciona como una contribución fundacional para el desarrollo de simulaciones cuánticas prácticas de la física fundamental. Los autores afirman que su trabajo proporciona los "ingredientes" necesarios para diseñar simulaciones de QCD en red en la formulación hamiltoniana.

• Prueba de Concepto: La principal significación radica en demostrar que la compleja maquinaria de la teoría de gauge $SU(3)$ en la red, incluyendo el ángulo $\theta$ y la materia fermiónica, puede codificarse y evolucionarse en hardware cuántico (simulado aquí).
• Vía de Escalabilidad: La codificación local y la base de representación se presentan como una vía viable hacia simulaciones más grandes, que potencialmente requieren menos qubits que las codificaciones globales.
• Perspectiva Futera: Los autores declaran modestamente que los resultados actuales son exploratorios y aún no realizan declaraciones físicas serias sobre la QCD en tiempo real debido a los tamaños de red pequeños, el régimen de acoplamiento fuerte y los truncamientos agresivos utilizados. Identifican la necesidad de mejores Hamiltonianos (contra-términos), mejores algoritmos de preparación de estados y hardware tolerante a fallos para alcanzar el límite continuo y simular regímenes físicamente relevantes. El trabajo se enmarca como un paso hacia la simulación de materia densa (por ejemplo, estrellas de neutrones) y la dinámica de axión-QCD en el futuro.