Worldline-Susceptibility Scheduling for Quantum Annealing Beyond Local-Adiabatic Evolution

Resumen Técnico: Programación de Susceptibilidad de Línea de Mundo para Recocido Cuántico

Planteamiento del Problema
El rendimiento del recocido cuántico (QA) depende críticamente del programa de recocido (annealing schedule), el cual dicta cómo se distribuye el tiempo total a lo largo de la trayectoria de evolución. Si bien el programa local-adiabático de Roland–Cerf es teóricamente óptimo —asignando el tiempo de forma inversamente proporcional al cuadrado de la brecha espectral instantánea ($\Delta(s)^{-2}$)—, es prácticamente inviable para problemas de optimización de gran escala. La construcción de este programa requiere el conocimiento completo del espectro de baja energía instantáneo, lo que exige la diagonalización repetida de un Hamiltoniano exponencialmente grande. En consecuencia, existe la necesidad de un sustituto computacionalmente económico que pueda identificar las regiones críticas del recocido sin información espectral explícita.

Metodología
Los autores proponen un marco de programación sustituto basado en la susceptibilidad de la magnetización de la línea de mundo ($\chi_m$) medida durante el Recocido Cuántico Simulado (SQA). La metodología procede de la siguiente manera:

1. Simulación SQA: El proceso de recocido cuántico se mapea a un sistema equivalente de mecánica estadística clásica utilizando la descomposición de Suzuki–Trotter. Esto transforma el modelo de Ising de campo transversal en un sistema clásico de $(d+1)$ dimensiones con $M$ réplicas acopladas (rebanadas de Trotter).
2. Medición de Observables: Se muestrean configuraciones de equilibrio mediante técnicas de Monte Carlo clásico. A partir de estas, se calcula la magnetización de la línea de mundo $m$, y su susceptibilidad se deriva mediante el teorema de fluctuación–disipación: $\chi_m = NM(\langle m^2 \rangle - \langle |m| \rangle^2)$.
3. Construcción del Sustituto: La susceptibilidad medida $\chi_m(s)$ actúa como un sustituto de la inversa del cuadrado de la brecha espectral. Se define un peso de asignación de tiempo $w(s) = \chi_m(s) + \chi_0$ (donde $\chi_0$ es una constante de regularización). La función de asignación de tiempo acumulada $\tau(s)$ se construye integrando este peso, y el programa final $s(t)$ se obtiene invirtiendo $\tau(s)$.
4. Validación: El marco, implementado en el software de código abierto Qanneal, valida los programas sustitutos frente a resultados de diagonalización exacta (ED) para instancias de vidrio de espín de Sherrington–Kirkert (SK). El rendimiento se cuantifica mediante la probabilidad final del estado fundamental ($P_{GS}$) obtenida resolviendo la ecuación de Schrödinger dependiente del tiempo para programas lineales, de Roland–Cerf exactos y el programa sustituto propuesto.

Resultos Clave
El estudio demuestra que el programa de susceptibilidad de la línea de mundo supera consistentemente tanto al recocido lineal convencional como, de forma contraintuitiva, al programa exacto de Roland–Cerf para una fracción sustancial de las instancias de SK.

• Identificación de Regiones Críticas: El pico de la susceptibilidad de la línea de mundo identifica de manera fiable la región amplia donde la brecha espectral es mínima, a pesar de que no reproduce cuantitativamente el perfil detallado de la brecha.
• Desempeño Superior en Tiempo Finito: Para instancias representativas ($n=10$ y $n=12$), el programa sustituto produce probabilidades de estado fundamental más altas que el programa exacto de Roland–Cerf a través de varios tiempos totales de recocido.
• Descubrimiento de Modos de Fallo: Los autores identifican dos modos de fallo distintos de tiempo finito del programa local-adiabático exacto que explican su bajo rendimiento:
  1. Trampa de Brecha de Frontera (Boundary-Gap Trap): Cuando la brecha mínima ocurre al final del recocido ($s \approx 1$), el programa de Roland–Cerf asigna la mayor parte del tiempo a una región donde el impulsor de campo transversal prácticamente ha desaparecido. Esto hace que el tiempo adicional sea dinámicamente ineficaz, ya que las fluctuaciones cuánticas están suprimidas. El programa basado en la susceptibilidad evita esto al concentrar el tiempo en la región de cruce más amplia donde las fluctuaciones cuánticas permanecen activas.
  2. Inestabilidad Oscilatoria: Cuando la brecha mínima es interior, el programa exacto concentra el tiempo de forma tan aguda alrededor de la brecha que induce una interferencia de niveles múltiples coherente. Esto resulta en una dependencia no monótona y oscilatoria de la probabilidad de éxito respecto al tiempo total de recocido, algo que el programa basado en la susceptibilidad, al ser más suave, logra evitar.
• Robustez Estadística: Experimentos numéricos a gran escala a través de realizaciones de desorden y tamaños de sistema ($n=10$ a $n=20$) confirman que la superioridad del programa sustituto es una característica genérica y no un artefacto de instancias específicas. A medida que el tamaño del sistema aumenta, la prevalencia de instancias de brecha de frontera crece, degradando aún más el rendimiento del programa exacto de Roland–Cerf en relación con el sustituto.

Significancia y Reivindicaciones
El artículo argumenta que, bajo condiciones realistas de tiempo finito, la robustez en la programación es más crítica que la reproducción precisa del perfil de la brecha espectral instantánea. Los autores afirman que:

• Observables de equilibrio económicos (específicamente la susceptibilidad de la línea de mundo desde SQA) pueden servir como sustitutos efectivos para las estrategias basadas en la brecha espectral.
• La condición local-adiabática exacta, aunque teóricamente óptima en el límite asintótico, puede ser contraproducente en regímenes de tiempo finito debido a la "trampa de brecha de frontera" y la "inestabilidad oscilatoria".
• La metodología propuesta ofrece un camino escalable para la optimización de programas que evita la intratabilidad computacional de los cálculos espectrales exactos, haciéndola aplicable a tamaños de problema que superan el alcance de la diagonalización exacta.

El trabajo está implementado en la plataforma de código abierto Qanneal, proporcionando un entorno unificado para SQA, la medición de susceptibilidad y la construcción de programas, facilitando así la investigación reproducible en la dinámica del recocido cuántico.