Investigadores demuestran que algoritmos cuánticos funcionan como protocolos de enfriamiento con recursos universales

Investigadores han formalizado la conexión entre dos métodos fundamentales de computación cuántica, demostrando que ambos actúan como protocolos de enfriamiento con propiedades universales, según un estudio publicado en Physical Review Research en 2026.

El trabajo establece que el Algoritmo de Optimización Aproximada Cuántica (QAOA) y el Recocido Cuántico (QA) generan distribuciones de probabilidad pseudo-Boltzmann cuya temperatura disminuye con los recursos invertidos, según los autores del estudio. En el caso de QAOA, la temperatura más fría es inversamente proporcional al número de capas del circuito cuántico.

Ambos métodos fueron propuestos para resolver problemas de optimización combinatoria que pueden mapearse a problemas de Optimización Binaria Cuadrática sin Restricciones (QUBO), demostrando que la tarea es al menos NP-difícil, según explica el artículo. El QA es un protocolo continuo basado en una deformación lenta de la dinámica del sistema, mientras que QAOA imita una versión trotterizada de QA mediante un circuito cuántico con capas que alternan evolución con hamiltonianos de mezcla y costo.

La investigación analizó 800 problemas aleatorios para sistemas de 8 a 16 qubits y 500 problemas para 18 y 20 qubits, con hasta 50 capas de compuertas, según detalla el estudio. Los resultados muestran que los ángulos óptimos de un circuito QAOA multicapa convergen hacia trayectorias universales de QA con interacción creciente monótona y hamiltoniano de mezcla progresivamente decreciente.

El análisis de las distribuciones de probabilidad revela una estructura bimodal en circuitos QAOA multicapa, según los investigadores. Esta estructura se caracteriza por dos temperaturas inversas distintas: una asociada al régimen de baja energía que se concentra cerca del estado fundamental, y otra correspondiente al régimen de alta energía que representa el ruido de fondo.

Los datos experimentales muestran que la temperatura inversa asociada al componente más frío de la función de onda crece linealmente con el número de capas, mientras que la temperatura del ruido se satura a un valor asintótico, según indica el estudio. En el límite de una sola capa, ambas temperaturas tienden al mismo valor, recuperando las distribuciones pseudo-Boltzmann reportadas previamente para QAOA de capa única.

El peso relativo de la cola de la distribución, que representa el ruido térmico, disminuye exponencialmente a medida que aumenta el número de capas, según los resultados. Para un número grande de capas, la contribución del componente más frío domina la distribución global, consistente con el límite de optimización perfecta con un número infinito de capas QAOA donde el sistema alcanza temperatura cero.

Los investigadores introdujeron una formulación que conecta la dinámica de QA con el protocolo QAOA mediante variables integradas que no dependen del tiempo sino de los pesos integrados de las contribuciones del hamiltoniano de mezcla e interacción. Esta formulación permite definir los valores finales de estas coordenadas como los recursos del protocolo de recocido, según explica el artículo.

Los ángulos de QAOA son las contrapartes discretas de la trayectoria hamiltoniana de QA, según establece el estudio. Esta identificación permite definir el recurso de un protocolo QAOA de manera similar como los ángulos integrados sobre las capas del circuito, conectando los recursos de QAOA con la calidad de la aproximación del estado fundamental o la temperatura efectiva de los estados producidos.

El análisis de las trayectorias hamiltonianas inferidas por el protocolo QAOA muestra que colapsan en curvas universales que aproximan trayectorias de recocido en el límite de grandes números de capas, según los resultados. Estas trayectorias convergen a círculos deformados de radio 1 tanto al aumentar el número de qubits como al incrementar el número de capas con qubits fijos.

Las trayectorias límite de recocido descritas por los ángulos QAOA corresponden a ecuaciones polares específicas con un parámetro pequeño que caracteriza la desviación de tamaño finito del paso desde el círculo de radio 1, según detalla el estudio. Los ángulos QAOA se distribuyen uniformemente a lo largo de esta trayectoria límite con desviaciones muy menores.

Los recursos integrados muestran un colapso excepcionalmente bueno de los valores promedio sobre diferentes tamaños de problema, con un aumento suave y monótono de los recursos con el número de capas, según indica la investigación. Este resultado proporciona evidencia de que los recursos demandados por los circuitos QAOA para lograr distribuciones con temperaturas progresivamente más bajas exhiben tendencias muy prometedoras.

El estudio se basó en una estadística extensa de problemas QUBO generados aleatoriamente definidos en grafos completamente conectados, donde todos los elementos no diagonales de las matrices de coeficientes son distintos de cero, según la metodología descrita. Los valores distintos de cero de cada muestra se extrajeron aleatoriamente de una distribución normal con media cero y varianza uno.

La propiedad de concentración de ángulos en trayectorias optimizadas de QAOA ha sido demostrada previamente para MaxCut en grafos regulares y el modelo Sherrington-Kirkpatrick, según reconocen los autores. Sin embargo, este trabajo es el primero en demostrar la universalidad de estas trayectorias concentradas, definida sobre conjuntos de vidrios de espín densos, y su conexión con trayectorias continuas de QA.

Los investigadores identificaron los errores de estos protocolos de recocido discretizados como distribuciones cuasitérmicas con propiedades ajustables, extendiendo la interpretación física de QAOA a ambos protocolos simultáneamente, según el artículo. Este descubrimiento permite la aplicación de protocolos QAOA y QA sin optimización previa y su uso como algoritmo de enfriamiento o simuladores de funciones de partición con temperatura ajustable.

La potencia de enfriamiento promedio de ambos métodos exhibe escalas algebraicas favorables con respecto a la temperatura objetivo y el tamaño del problema, según concluye el estudio. En QAOA, la temperatura más fría es inversamente proporcional al número de capas y a los ángulos integrados, o interacciones integradas en QA.

La investigación utilizó una divergencia de Kullback-Leibler modificada para calibrar los parámetros libres de la distribución bimodal de Boltzmann, ponderando la distribución con una función que mejora la sensibilidad del ajuste a la cola de la distribución, según la metodología. Este peso captura muy bien tanto las partes de alta como de baja energía del espectro.

El hamiltoniano de la red de espín cuántico puede mapearse a problemas de optimización clásica del tipo QUBO o de optimización binaria sin restricciones de orden superior, según explica el estudio. Muchos problemas de optimización del mundo real están contenidos en familias de problemas binarios NP-difíciles como QUBO.

Los resultados proporcionan evidencia de que tanto QAOA como QA son protocolos de enfriamiento y simuladores de funciones de partición cuya temperatura objetivo puede ajustarse reescalando la trayectoria universal, según los investigadores. Esta propiedad abre nuevas posibilidades para la aplicación práctica de estos métodos cuánticos en problemas de optimización combinatoria.

El estudio fue publicado en Physical Review Research, volumen 8, artículo 013211, en 2026, según la referencia proporcionada. La investigación representa un avance significativo en la comprensión de los recursos y el rendimiento de los protocolos de optimización cuántica, estableciendo una base teórica y empírica para su aplicación sin necesidad de optimización previa de parámetros.