Low-Depth Initial-State Preparation for Ground-State Energy Estimation of Two-Dimensional Strongly Correlated Systems

Resumen Técnico: Preparación de Estados Iniciales de Baja Profundidad para la Estimación de la Energía del Estado Fundamental de Sistemas Fuertemente Correlacionados Bidimensionales

Planteamiento del Problema
Los algoritmos cuánticos para la estimación de la energía del estado fundamental, como la Estimación de la Fase Cuántica (QPE), requieren un estado inicial con una fidelidad no despreciable respecto al verdadero estado fundamental. La probabilidad de extraer exitosamente la energía del estado fundamental escala con esta fidelidad ($F_0$); si $F_0$ decae exponencialmente con el tamaño del sistema, el algoritmo se vuelve intratable debido al número de repeticiones requerido ($1/F_0$). Aunque los Solucionadores Cuánticos Variacionales (VQE) ofrecen un enfoque híbrido, enfrentan desafíos en los costos de evaluación de gradientes y la simulabilidad clásica de circuitos poco profundos. Además, compilar aproximaciones clásicas precisas (por ejemplo, Estados de Producto de Matrices) en circuitos cuánticos suele resultar en una profundidad de circuito excesiva, lo cual es prohibitivo para las computadoras cuánticas tolerantes a fallos (FTQC) tempranas, donde los recursos como los qubits lógicos y los estados mágicos son limitados. El desafío central es diseñar circuitos cuánticos poco profundos y eficientes en recursos que preparen estados iniciales de alta fidelidad para sistemas fuertemente correlacionados, específicamente el modelo de Hubbard bidimensional (2D), sin depender de bucles de optimización profundos o de la intratabilidad clásica.

Metodología
Los autores proponen una estrategia para construir circuitos cuánticos poco profundos aprovechando la conexión de acoplamiento fuerte entre el modelo de Heisenberg 2D y el modelo de Hubbard en medio llenado mediante la transformación de Schrieffer–Wolff (SW). La metodología procede en tres etapas principales:

1. Construcción del Circuito de Heisenberg:

   • Ansatz: El circuito apunta al subespacio de singlete $S=0$ del modelo de Heisenberg antiferromagnético. El estado inicial es un producto de dímeros de singlete ($| \Phi \rangle$) definidos sobre un emparejamiento perfecto de la red.
   • Dinámica: El estado evoluciona utilizando capas de puertas SWAP exponencializadas ($e$SWAP). Estas puertas respetan la simetría SU(2), asegurando que el estado permanezca en el subespacio de singlete.
   • Estructura: Los enlaces de la red se particionan en cuatro subconjuntos disjuntos ($B_{x0}, B_{x1}, B_{y0}, B_{y1}$). El circuito aplica puertas $e$SWAP en una secuencia específica ($B_{y1} \to B_{y0} \to B_{x1} \to B_{x0}$) a través de múltiples bloques ($N_b$). Esta estructura imita la física de enlace de valencia resonante (RVB) del estado fundamental de Heisenberg.
   • Mapeo Fermiónico: Los estados de espín se mapean al subespacio de ocupación simple mediante un embebido unitario local ($V$) que introduce qubits ancilares.

2. Extensión de Hubbard mediante Puertas SW:

   • Para capturar las fluctuaciones de carga (dobletes y huecos) presentes en el modelo de Hubbard a un $U/t$ finito, el circuito de Heisenberg se "viste" con puertas de fluctuación de carga.
   • Estas puertas se derivan directamente del generador de la primera orden de SW ($\Omega$), aproximado mediante una descomposición de Suzuki–Trotter en unitarias de dos sitios ($W_{\langle i,j \rangle}$).
   • El ansatz variacional final combina la evolución de Heisenberg ($U_B$) y el bloque de fluctuación de carga ($W_B$), parametrizado por ángulos de rotación $\theta$ y $\phi$.

3. Estimación de Fidelidad y Benchmarking:

   • Transferencia de Parámetros: Los parámetros optimizados en una red pequeña de $4 \times 4$ se transfieren a redes más grandes ($6 \times 6$ a $10 \times 10$) sin optimización adicional.
   • Simulación: Los estados del circuito se simulan utilizando Estados de Pares Entrelazados Proyectados (PEPS) con actualizaciones simples guiadas por propagación de creencias (BP).
   • Límites de Fidelidad: Dado que el estado fundamental exacto es desconocido para sistemas grandes, los autores estiman un límite inferior de la fidelidad con el estado fundamental ($F_0$). Esto se logra mediante:
     • El cálculo de la fidelidad ($F_{ref}$) entre el estado del circuito y un estado de referencia de Monte Carlo Variacional (VMC).
     • El cálculo del límite de la fidelidad del estado VMC respecto al estado fundamental exacto ($F_{ref,0}$) usando la energía de referencia, la energía del estado fundamental exacto (de la Expansión de Series Estocásticas) y la brecha de excitación de singlete ($\Delta_0$).
     • La aplicación de la desigualdad triangular para los ángulos de Fubini–Study para derivar un límite inferior riguroso para $F_0$.

Contribuciones Clave y Resultados

• Transferencia Efectiva de Parámetros: Para el modelo de Heisenberg, los parámetros optimizados en una red de $4 \times 4$ producen límites inferiores de fidelidad robustos para redes de hasta $10 \times 10$ sin re-optimización. Específicamente, para $N_b=3$, el límite inferior estimado de la fidelidad del estado fundamental permanece positivo ($\approx 0.028$) incluso para la red de $10 \times 10$, lo que sugiere que el estado inicial es adecuado para la amplificación de amplitud.
• Mejora de la Fidelidad de Hubbard: Los cálculos exactos en una red de $4 \times 4$ demuestran que la incorporación de las puertas de fluctuación de carga derivadas de SW mejora significamente la fidelidad con el estado fundamental de Hubbard. Para $N_b=4$, la fidelidad aumenta de $\approx 0.44$ (solo Heisenberg) a $\approx 0.89$ (con ángulos de carga optimizados).
• Eficiencia de Recursos: Los autores estiman el costo de preparación del estado para una red de $10 \times 10$. Asumiendo la síntesis de Clifford+$T$ con un error por rotación de $10^{-5}$, el costo total es del orden de $10^5$ puertas $T$. Esto sitúa el costo de preparación dentro del régimen de "megaquops" de las FTQC tempranas y significativamente por debajo de los costos asociados con los algoritmos posteriores de estimación de la energía del estado fundamental.
• Escalabilidad: La profundidad del circuito escala linealmente con el número de enlaces de la red ($|B|$), y la estrategia de compartir parámetros permite la construcción de circuitos para sistemas grandes basados en la optimización de sistemas pequeños.

Significancia y Reivindicaciones
El artículo afirma ofrecer una ruta viable para la preparación de estados iniciales de bajo costo para la estimación de la energía del estado fundamental en sistemas fuertemente correlacionados. Al utilizar la conexión física entre los modelos de Heisenberg y de Hubbard, los autores demuestran que los circuitos poco profundos pueden alcanzar una fidelidad no despreciable con el estado fundamental en redes que superan el alcance de la diagonalización clásica exacta. Los resultados sugieren que estos circuitos son candidatos atractivos para el régimen de FTQC temprana, donde minimizar la profundidad del circuito y el conteo de puertas $T$ es crítico. El trabajo respalda la hipótesis de que las simulaciones en redes pequeñas pueden guiar eficazmente el diseño de circuitos cuánticos para sistemas más grandes y clásicamente intratables, siempre que los parámetros sean transferibles y la profundidad del circuito se mantenga baja. Los autores señalan explícitamente que, aunque sus límites son rigurosos, las fidelidades reales pueden ser mayores, y que el trabajo futuro podría involucrar ansatzes variacionales de Hamiltoniano más profundos o técnicas de simulación más precisas (por ejemplo, PEPS de actualización completa) para refinar aún más estas estimaciones.