Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

Resumen Técnico: Sondeo del Estado Fundamental del Modelo de Heisenberg Antiferromagnético en la Red de Kagome mediante un Algoritmo de Eigensolver Cuántico Variacional Geométricamente Informado

Planteamiento del Problema
El artículo aborda el desafío de encontrar el estado fundamental del modelo de Heisenberg antiferromagnético en una red de kagome, un sistema conocido por su frustración geométrica y potenciales fases de líquido de espín cuántico. Determinar el estado fundamental de tales Hamiltonianos $k$-locales es un problema QMA-duro, lo que hace que la simulación clásica sea intratable para sistemas grandes debido al crecimiento exponencial del espacio de Hilbert. Aunque los Eigensolvers Cuánticos Variacionales (VQE) ofrecen un enfoque híbrido cuántico-clásico adecuado para dispositivos de Escala Intermedia con Ruido Cuántico (NISQ), enfrentan obstáculos significativos: la necesidad de circuitos profundos, la susceptibilidad al ruido (errores coherentes e incoherentes), los barrenos de planicie (barren plateaus) donde los gradientes se desvanecen, y la dificultad de diseñar ansätze que sean tanto expresivos como eficientes para el hardware. Específicamente, los autores pretenden caracterizar las propiedades del estado fundamental, como las correlaciones espín-espín y los factores de estructura estática, en hardware cuántico real sin depender de reducciones pesadas específicas del problema o de una mitigación de errores extensiva durante el bucle de optimización.

Metodología
Los autores proponen un enfoque VQE geométricamente informado adaptado para la red de kagome, centrándose en dos geometrías fundamentales: un triángulo individual y una estrella de kagome (seis triángulos que comparten esquinas).

1. Ansatz Geométricamente Informado:
   La innovación central es un ansatz de bajo nivel y eficiente para el hardware, diseñado a medida. A diferencia del ansatz de la Variación del Hamiltoniano (HVA), que puede ser demasiado profundo para los dispositivos NISQ, este ansatz se construye utilizando un patrón de rotaciones de un solo qubit parametrizadas ($R_y$) y puertas de entrelazamiento restringidas a la conectividad nativa del hardware (evitando las puertas SWAP).

   • Espacio de Parámetros Euclídeo: Una característica crítica de este diseño es que se ha demostrado analíticamente que el tensor métrico de Fubini-Study es diagonal y constante. Esto asegura que el espacio de parámetros sea naturalmente euclídeo, libre de singularidades y de sobreparametrización. Consecuentemente, el gradiente natural cuántico es equivalente al gradiente estándar, simplificando el paisaje de optimización.
   • Compatibilidad con el Hardware: La profundidad del circuito se minimiza (una capa de entrelazamiento) y el número de parámetros se reduce alternando la colocación de las puertas $R_y$.

2. Estrategia de Optimización (AQNGD/AGD):
   Los autores emplean un optimizador de Descenso de Gradiente Adaptativo (AGD). Debido a que el ansatz produce un métrico de Fubini-Study diagonal y constante, el Descenso de Gradiente Natural Cuántico Adaptativo (AQNGD) se reduce a AGD. La tasa de aprendizaje se adapta dinámicamente utilizando la regla de retroceso de Armijo. Esto permite pasos iniciales grandes para acelerar la convergencia, mientras reduce automáticamente el tamaño del paso a medida que el optimizador se acerca al mínimo, evitando el estancamiento en geometrías estrechas. Esto se compara con la Aproximación Estocástica de Perturbación Simultánea (SPSA).

3. Mitigación de Errores (Post-procesamiento):
   La mitigación de errores se aplica después de la optimización del VQE para refinar las estimaciones finales de energía y las propiedades estructurales, en lugar de durante el entrenamiento iterativo.

   • Mitigación de Error de Lectura (REM): Se utiliza un enfoque de matriz de respuesta particionada para revertir los errores de SPAM (Preparación de Estado y Medición), escalando de manera más eficiente que la inversión de matriz completa para mediciones dispersas.
   • Extrapolación de Ruido Cero (ZNE): Los autores aplican ZNE utilizando extrapolación cuadrática (y regresión polinómica bayesiana) para extrapolar los resultados al límite de ruido cero. Observan que, aunque ZNE no obedece estrictamente el principio variacional de Rayleigh-Ritz, se utiliza aquí para mejorar la precisión.

Resultos Clave
El estudio se llevó a cabo en hardware cuántico real (IBM Oslo, IBM Kyoto, IBM Torino e IBM Algiers) y se validó mediante simulaciones de vector de estado.

• Convergencia y Optimización:

  • En un sistema de triángulo de 3 qubits, el VQE optimizado por AGD convergió suavemente al estado fundamental de energía ($E \approx -3$) con alta fidelidad (99.89% en simulación) utilizando solo 19 iteraciones. Los parámetros finales coincidieron estrechamente con los valores teóricos exactos.
  • En una estrella de kagome de 12 qubits, el AGD demostró una convergencia superior comparada con SPSA, requiriendo menos iteraciones para alcanzar la energía del estado fundamental ($E \approx -18$). La fidelidad final en la simulación fue de 96.37%.
  • Se demostró que la adaptación del tamaño del paso mediante la regla de Armijo permite navegar eficazmente el paisaje de costo, comenzando con pasos grandes y refinando a medida que la convergencia se aproxima.

• Resiliencia al Ruido de las Propiedades Estructurales:

  • Un hallazgo significativo es que las correlaciones espín-espín y el factor de estructura estática ($S(q)$) son resilientes al ruido. Incluso sin mitigación de errores, los mapas experimentales de $S(q)$ retuvieron la firma cualitativa del estado fundamental (por ejemplo, la supresión en el momento cero que indica la ausencia de orden ferromagnético).
  • El coeficiente de correlación de Pearson entre el $S(q)$ experimental ruidoso y el $S(q)$ teórico exacto se mantuvo extremadamente alto (>99.9%).
  • Los experimentos de simulación revelaron una dependencia "escalonada" del factor de estructura estática respecto a los niveles de ruido: la estructura de correlación es robusta hasta un umbral crítico de ruido, tras el cual se degrada rápidamente, asemejándose a una transición de fase. Esto contrasta con la energía, que decae monótonamente con el ruido.

• Eficacia de la Mitigación de Errores:

  • Aplicar REM y ZNE post-optimización mejoró significativamente las estimaciones de energía. Para la estrella de kagome, la energía bruta fue de $\approx -14.9$, mientras que el REM por sí solo la mejoró a $\approx -17.3$, y ZNE (cuadrática) la acercó al valor exacto de $-18$.
  • Los autores observaron que combinar REM y ZNE a veces provocaba sobrepasar la energía del estado fundamental (violando el límite variacional), un comportamiento consistente con la naturaleza no variacional de ZNE. ZNE por sí solo a menudo proporcionó la mejor precisión entre las configuraciones probadas.

Significancia y Reivindicaciones
El artículo sostiene que un diseño de ansatz informada geométricamente, que garantiza un espacio de parámetros naturalmente euclídeo, permite el entrenamiento eficiente de VQE en el hardware actual de la era NISQ sin la necesidad de reducciones complejas y específicas del problema.

• Robustez: El trabajo demuestra que propiedades físicas significativas, específicamente las correlaciones espín-espín y el factor de estructura estática, pueden extraerse de hardware cuántico ruidoso con suficiente precisión para caracterizar el estado fundamental, incluso sin mitigación de errores durante el bucle de optimización.
• Eficiencia de Optimización: El uso de AGD con tamaños de paso adaptativos demostró ser más eficiente que SPSA para este ansatz de eficiencia de hardware, logrando una convergencia más rápida en términos de iteraciones.
• Escalabilidad: Al evitar circuitos profundos y una mitigación de errores pesada durante el entrenamiento, el enfoque ofrece un camino escalable para sondear imanes cuánticos frustrados. Los autores sugieren que este marco puede extenderse a redes de kagome más grandes y otros sistemas cuánticos frustrados.
• Modestia: Los autores reconocen que, si bien su método caracteriza con éxito pequeñas celdas fundamentales (triángulo y estrella), el estado fundamental general de la red de kagome infinita sigue siendo una pregunta abierta en la literatura. No afirman haber resuelto el problema de la red infinita, sino que han demostrado una metodología viable para sondear estos sistemas en dispositivos actuales. También señalan que, aunque ZNE mejora las estimaciones de energía, no preserva el principio variacional, y que la combinación de REM y ZNE requiere una interpretación cuidadosa.