A Quantum Algorithm with Polylogarithmic Depth per Trotter Step for the Extended Hubbard Model

Resumen Técnico: Q2FMM para el Modelo de Hubbard Extendido

Planteamiento del Problema
La simulación del modelo de Fermi-Hubbard extendido en una red bidimensional es un desafío central en la computación cuántica, particularmente para la física de la materia condensada y las estructuras electrónicas fuertemente correlacionadas. Mientras que los términos locales (salto e interacciones en el sitio) pueden simularse eficientemente con circuitos de profundidad constante o baja, la inclusión de interacciones de Coulomb de largo alcance presenta un cuello de botella significativo. Estas interacciones generan $O(N^2)$ términos (donde $N$ es el número de sitios de la red), lo que conduce a circuitos cuánticos profundos. Las estrategias estándar, como las redes de SWAP fermiónico, resultan en una profundidad de circuito por paso de Trotter que escala como $O(N)$. Las estrategias de truncamiento suelen ser ineficaces para los potenciales de Coulomb debido a su lento decaimiento $1/r$, lo que requiere radios de corte grandes para mantener la precisión.

Metodología: El Algoritmo Q2FMM
Los autores proponen Q2FMM, un algoritmo cuántico inspirado en el Método de Multipolos Rápidos (FMM) clásico para simular la evolución temporal del modelo de Hubbard extendido dentro de la fórmula de producto de Trotter. La metodología central consiste en reemplazar las interacciones directas sitio-sitio por interacciones entre "cajas" jerárquicas y de grano grueso a través de múltiples escalas de longitud.

1. Estructura Jerárquica: El algoritmo organiza la red 2D en una estructura de quadtree. Distingue entre interacciones de "campo cercano" (gestionadas directamente en el nivel más fino) e interacciones de "campo lejano" (aproximadas mediante expansiones multipolares).
2. Expansiones Multipolares:
   • Q2FMM de Orden 0: Utiliza una aproximación simplificada donde las cajas distantes interactúan basándose en sus números de ocupación totales ($N_A N_B$).
   • Q2FMM de Orden Superior: Refina la precisión reteniendo los coeficientes multipolares ($M_{\ell m}$) en lugar de solo los números de ocupación totales. Esto permite un compromiso sistemático entre precisión y costo computacional mediante un orden de truncamiento $p$.
3. Diseño del Circuito Cuántico:
   • Fase de Cómputo: El algoritmo procede desde el nivel más fino (sitios individuales) hacia el nivel más grueso. En cada nivel, agrega información de las cajas hijas a las cajas padres utilizando sumadores cuánticos. Para las versiones de orden superior, esto implica traducir y sumar coeficientes multipolares.
   • Aplicación de Fase: La energía de interacción para pares de cajas se calcula utilizando multiplicadores cuánticos para calcular productos de números de ocupación o términos multipolares. La fase resultante se aplica mediante puertas de fase de un solo qubit.
   • Fase de Descomputación (Uncompute): Para evitar la acumulación de estados basura y el entrelazamiento no deseado, el algoritmo emplea una fase de "descomputación" reversible. Este proceso revierte los pasos de agregación (descomponiendo la información del padre de nuevo en la información de los hijos) después de que la fase ha sido impresa, asegurando que el sistema regrese a un estado limpio mientras preserva la fase dinámica.
4. Optimizaciones: El artículo discute optimizaciones opcionales utilizando operaciones COPY (CNOTs transversales) para paralelizar los cálculos de fase y puertas de fan-out sin límites para reducir la profundidad de los circuitos aritméticos a $O(\log^* n)$, aunque el análisis de escalamiento principal se basa en circuitos reversibles estándar.

Resultados Clave y Complejidad
El circuito cuántico resultante logra las siguientes propiedades de escalamiento por paso de Trotter y para el tiempo total de simulación $t$ con un error objetivo $\epsilon$:

• Profundidad del Circuito: La profundidad por paso de Trotter escala polilogarítmicamente con el tamaño del sistema $N$. Específicamente, escala como $O(\log N \cdot \log Q)$, donde $Q$ es el número de electrones.
• Complejidad de Puertas: La complejidad de puertas por paso de Trotter escala cuasi-linealmente con el tamaño del sistema, denotado como $\tilde{O}(N)$.
• Complejidad Total: Al contabilizar el número de pasos de Trotter requeridos para simular el tiempo $t$ con un error $\epsilon$, la profundidad total del circuito escala como $\tilde{O}(N^{1/2} t^{3/2} / \epsilon^{1/2})$, y la complejidad total de puertas escala como $\tilde{O}(N^{3/2} t^{3/2} / \epsilon^{1/2})$.
• Qubits Ancilla: El algoritmo requiere un número de qubits ancilla que escala cuasi-linealmente con el tamaño del sistema, $\tilde{O}(N)$.

Significancia y Reivindicaciones
El artículo afirma que Q2FMM ofrece una ventaja significativa sobre los métodos existentes para simular interacciones de largo alcance en sistemas grandes.

• Comparación con el Cuadro de Interacción: Los autores comparan Q2FMM con el método del cuadro de interacción (Low y Wiebe, Ref. [82]), el cual tiene una complejidad de puertas de $\tilde{O}(N^2 t)$. Q2FMM resulta ventajoso en el régimen de sistemas grandes (N grande), mientras que el método del cuadro de interacción es más favorable para simulaciones de largo tiempo y alta precisión (t grande, $\epsilon$ pequeño). El punto de cruce se estima en $N \sim \tilde{O}(t/\epsilon)$.
• Idoneidad del Hardware: Los autores sugieren que Q2FMM es particularmente adecuado para computadoras cuánticas de átomos neutros en dos dimensiones, que soportan nativamente el transporte de átomos (atom shuttling) y puertas de entrelazamiento de largo alcance. También es compatible con arquitecturas de código de superficie donde las operaciones de largo alcance pueden realizarse en profundidad constante.
• Generalidad: Aunque se demuestra en el modelo de Hubbard extendido, el marco se presenta como aplicable a otros núcleos (por ejemplo, Yukawa, Helmholtz) y potencialmente a Hamiltonianos moleculares ab initio discretizados en mallas de espacio real.

Distinción de Trabajos Relacionados
Los autores distinguen Q2FMM de los enfoques recientes basados en FMM para Hamiltonianos moleculares (por ejemplo, Berry et al., Ref. [80]). En esos trabajos, el estado cuántico está en primera cuantización, lo que requiere una clasificación compleja de los electrones en cajas espaciales basadas en datos de posición cuántica. En contraste, Q2FMM opera en una red donde los qubits están pre-dispuestos por posición espacial (segunda cuantización), eliminando la necesidad de movimiento de datos cuánticos basados en la posición. En consecuencia, Q2FMM se centra en el uso eficiente y jerárquico de los resultados intermedios y en la gestión de los estados basura mediante la descomputación.