Reshaping quantum annealing landscapes with diagonal catalysts

Resumen Técnico: Remodelación de los Paisajes de Recocido Cuántico mediante Catalizadores Diagonales

Declaración del Problema
El recocido cuántico (QA, por sus siglas en inglés) es un algoritmo de optimización diseñado para encontrar el estado fundamental de un Hamiltoniano de costo ($H_P$) que representa un problema de Optimización Binaria No Restringida Cuadrática (QUBO). Sin embargo, el rendimiento del QA se ve frecuicamente limitado porque la población queda atrapada en mínimos locales que están lejos de la solución global en términos de distancia de Hamming. Este problema se exacerba en instancias difíciles donde la brecha de energía entre el estado fundamental y los estados excitados dinámicamente accesibles se vuelve exponencialmente pequeña. En consecuencia, los recocidos de tiempo finito suelen operar de forma diabática, lo que conduce a la redistribución de la población a través de cruces evitados entre estados localizados en distintos pozos del paisaje de energía. Un desafío crítico es que las configuraciones de baja energía no corresponden necesariamente a configuraciones cercanas a la solución en términos de distancia de Hamming; las configuraciones distantes pueden "socavar" a las cercanas, proporcionando candidatos de baja energía para pozos competidores.

Metodología y Marco Teórico
Los autores introducen un marco matemático que vincula la energía y la distancia de Hamming para diseñar catalizadores diagonales de dos cuerpos. Estos catalizadores se construyen únicamente a partir de la estructura de acoplamiento del problema (signos, magnitudes y topología del grafo) sin requerir conocimiento previo de la solución.

1. Análisis de Capas de Hamming: Los autores definen "capas de Hamming" ($S_{d, z^\star}$) como conjuntos de configuraciones a una distancia de Hamming $d$ específica de un estado fundamental de referencia $z^\star$. Derivan fórmulas exactas para la media y la varianza de la energía sobre estas capas para cualquier operador de Ising de dos cuerpos libre de campo.

   • Teorema 1 establece que la media de la energía de una capa sigue una curva parabólica determinada por factores combinatorios, mientras que la varianza (ancho de la capa) depende de la matriz de acoplamiento específica.
   • En el problema base, las medias de las capas están correctamente ordenadas, pero la amplia varianza causa un solapamiento de las capas, permitiendo que configuraciones distantes tengan energías más bajas que las cercanas.

2. Diseño del Catalizador ($C$): El objetivo es construir un catalizador diagonal $C$ que aumente la separación entre las medias de las capas en relación con sus fluctuaciones dentro de la capa, efectivamente "canalizando" el paisaje.

   • El catalizador ideal sería proporcional a $-z^\star_i z^\star_j$, pero esto requiere el conocimiento del estado fundamental desconocido.
   • Construcción de Camino Abierto (Open-Path): Los autores proponen estimar el patrón del estado fundamental utilizando caminos auto-evitantes en el grafo de interacción. Al propagar patrones de espín a lo largo de estos caminos basados en los signos de acoplamiento ($g(p)_{v_{a+1}} = -\text{sgn}(J_{v_a v_{a+1}}) g(p)_{v_a}$), generan patrones localmente satisfactorios.
   • El catalizador $C^{(m)}$ se forma sumando las contribuciones de todos los caminos auto-evitantes de $m$ aristas. Cada camino contribuye con una estimación de rango uno ponderada de la estructura de acoplamiento ideal. El operador final es proyectado ortogonal a los acoplamientos originales y normalizado.
   • Esta construcción crea un grafo de interacción más denso que el problema original, penalizando las configuraciones que frustran los enlaces fuertes y reforzando aquellas que se alinean con los patrones propagados.

Resultados Clave
Los autores evaluaron su método en 200 instancias aleatorias 3-regulares con $n=20$ espines y acoplamientos extraídos de $[-1, 1]$. Las simulaciones se realizaron con tiempos de barrido $T = 2, 6, 10$.

• Remodelación del Paisaje de Energía: El catalizador logra remodelar el paisaje estático, reduciendo el solapamiento entre las capas de Hamming y profundizando el embudo de energía hacia el estado fundamental.
• Redistribución de la Probabilidad:
  • Espacio de Energía: El catalizador desplaza la masa de probabilidad desde las bandas de energía intermedia hacia la banda del estado fundamental ($E/E_{GS} \ge 0.95$). Con $T=6$, la masa mediana de la banda superior aumentó de 0.067 (sin catalizador) a 0.324, lo que representa una ganancia mediana pareada de +351%.
  • Espacio de Hamming: La distribución se desplaza hacia la variedad del estado fundamental. Con $T=6$, la masa de probabilidad en la distancia de Hamming $\delta \le 2$ aumentó de 0.10 a 0.35 (+239%).
• Robustez y Comparación:
  • Las ganancias persisten a través de diferentes conectividades de grafos (grafos 3-, 4- y 5-regulares) y en el modelo Sherrington–Kirkler de conectividad completa, aunque la magnitud de la mejora disminuye a medida que aumenta la conectividad.
  • El catalizador de camino abierto supera significativamente a una línea base de "problema reforzado" ($C = H_P$), demostrando que el beneficio proviene de la redistribución estructural específica de los acoplamientos y no simplemente de aumentar la escala de energía general.
  • Existe un compromiso entre el orden del camino ($m$) y la densidad del grafo; órdenes moderados (por ejemplo, $m=4$ para grafos 3-regulares) ofrecen el mejor equilibrio entre resolver la estructura local y mantener la fidelidad al orden del verdadero estado fundamental.

Significancia y Reivindicaciones
El artículo afirma proporcionar un método práctico para mejorar el recocido cuántico en hardware de próxima generación utilizando controles longitudinales nativos (operadores diagonales) que son de dos cuerpos y no requieren los complejos gadgets necesarios para los catalizadores no estocásticos.

• Diseño Agnóstico a la Solución: Los catalizadores se construyen sin conocimiento de la solución, basándose únicamente en la matriz de acoplamiento del problema.
• Criterios Estáticos para el Diseño: El teorema de momentos de capa derivado ofrece criterios estáticos (profundidad del embudo y solapamiento de capas) para filtrar candidatos a catalizadores antes de ejecutar simulaciones de recocido dinámico.
• Compatibilidad de Hardware: Al ser diagonales y de dos cuerpos, los catalizadores preservan la estructura del campo transversal y son compatibles con las restricciones de conectividad y rango de acoplamiento de los recocedores actuales.
• Direcciones Futuras: Los autores sugieren que el marco puede extenderse a entornos digitales (por ejemplo, como un separador de fase en QAOA), combinarse con protocolos diabáticos optimizados o aplicarse al recocido inverso. También señalan que caracterizar el escalamiento con el tamaño del sistema y generalizar a funciones de costo que incluyan campos sigue siendo una tarea abierta.