⚛️ quantum physics

Efficient classical simulation of large-scale unitary cluster Jastrow circuits

Este artículo presenta un algoritmo clásico de tiempo polinomial capaz de simular eficientemente circuitos Jastrow de clúster unitario de una sola capa a gran escala, permitiendo que una computadora portátil reproduzca y supere un reciente experimento cuántico de 77 cúbits en menos de un minuto.

Autores originales: Hrishikesh Belagali, Thomas Van Camp, R. Pradeep, Sourin Das, Namit Anand, Ryan LaRose

Publicado 2026-07-24

📖 6 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Hrishikesh Belagali, Thomas Van Camp, R. Pradeep, Sourin Das, Namit Anand, Ryan LaRose

Artículo original bajo licencia CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Imagina que estás intentando resolver un nudo masivo y enredado de cuerda que representa una molécula. En el mundo de la química, averiguar cómo estos nudos se asientan en su forma más estable (su "estado fundamental") es como intentar predecir el clima en un huracán: es increíblemente difícil porque cada trozo de cuerda tira de todos los demás. Durante décadas, los científicos han esperado que las computadoras cuánticas —máquinas que utilizan las reglas extrañas de las partículas diminutas para hacer matemáticas— pudieran desatar estos nudos más rápido que una computadora normal. La idea es que, mientras una computadora regular tiene que revisar cada posible giro uno por uno, una computadora cuántica puede revisar muchos giros a la vez. Recientemente, investigadores han estado construyendo estas máquinas cuánticas para simular moléculas complejas, como los cúmulos de hierro-azufre que se encuentran en la naturaleza, con la esperanza de encontrar nuevos medicamentos o materiales. La gran pregunta ha sido: ¿Están estas máquinas cuánticas haciendo realmente algo que una computadora regular no pueda hacer, o simplemente estamos usando un mazo para romper una nuez que una computadora regular podría haber manejado con un poco más de ingenio?

Este artículo trata sobre un equipo de investigadores que decidió probar ese mazo. Analizaron un tipo específico de experimento cuántico llamado circuito "Unitary Cluster Jastrow" (UCJ), que se ha utilizado en algunos de los experimentos de química cuántica más grandes e impresionantes hasta la fecha. Estos experimentos involucraron computadoras cuánticas con hasta 77 qubits (la versión cuántica de los bits) y miles de puertas, ejecutándose en supercomputadoras con miles de nodos solo para procesar los resultados. Los investigadores se preguntaron: "¿Podemos hacer este mismo cálculo en una laptop regular?". La respuesta, sorprendentemente, es sí. Desarrollaron un truco matemático nuevo y súper rápido que permite a una computadora estándar calcular la energía de estos circuitos cuánticos específicos en tiempo polinomial, lo que significa que el tiempo que tarda crece razonablemente con el tamaño del problema, en lugar de explotar hacia el infinito.

El equipo descubrió que podían reproducir los resultados del experimento más grande jamás realizado sobre un cúmulo de hierro-azufre (que utilizó 77 qubits y 10,570 puertas) en menos de un minuto en una laptop. De hecho, debido a que su método era tan rápido, pudieron ajustar los parámetros del circuito para encontrar un estado de energía incluso más bajo que el alcanzado por el experimento cuántico, el cual había requerido 6,400 nodos de supercomputadora y horas de procesamiento de tiempo. Sin embargo, hay un inconveniente. El artículo establece explícitamente que este "triunfo" solo funciona para circuitos de una sola capa. Si se añaden más capas al circuito (haciéndolo más profundo y complejo), las matemáticas se rompen y el problema probablemente se vuelva demasiado difícil para las computadoras regulares nuevamente. Los autores concluyen que los circuitos UCJ de una sola capa no son suficientes para demostrar la ventaja cuántica (el punto donde las computadoras cuánticas superan a las clásicas), pero necesitaremos construir circuitos mucho más profundos y de múltiples capas para ver realmente qué pueden hacer las computadoras cuánticas que las clásicas no puedan.

Para entender cómo hicieron esto, piensa en el circuito cuántico como una máquina compleja que transforma un estado inicial (como una hoja de papel plana) en una forma final arrugada. Usualmente, para ver cómo es la forma final, tienes que ejecutar la máquina y luego intentar realizar ingeniería inversa de la arruga, lo cual es una pesadilla para las computadoras regulares porque el número de posibilidades es astronómico. El truco de los autores fue trabajar hacia atrás. En lugar de intentar predecir la forma arrugada final, tomaron las "reglas" de la molécula (el Hamiltoniano) y las hicieron correr hacia atrás a través de la máquina.

Imagina que tienes una receta para un pastel, pero en lugar de hornearlo y luego probarlo, comienzas con el pastel terminado y trabajas hacia atrás a través de la receta para ver exactamente qué ingredientes se usaron. Los autores demostraron que, para estos circuitos específicos de una sola capa, puedes ejecutar la receta hacia atrás sin que el número de ingredientes explote hacia el infinito. Utilizaron una herramienta matemática llamada fórmula de Löwdin, que es como una calculadora especial que puede determinar rápidamente el valor de una forma compleja sin tener que medir cada punto individualmente. Al combinar este método de ejecución hacia atrás con una forma ingeniosa de manejar las "fases" (el tiempo y los ángulos de los movimientos cuánticos), mantuvieron la eficiencia del cálculo.

Los resultados fueron impactantes. Cuando aplicaron su método al experimento del cúmulo de hierro-azufre, obtuvieron un valor de energía de -326.796 Hartrees, que es más bajo (mejor) que los -326.645 Hartrees alcanzados por el experimento cuántico utilizando el método de Diagonalización Cuántica Basada en Muestreo (SQD). El experimento cuántico tuvo que usar una supercomputadora masiva (Fugaku) con 6,400 nodos para obtener su resultado, mientras que los autores obtuvieron un mejor resultado en una laptop en menos de un minuto. También probaron esto en cadenas de hidrógeno con hasta 160 qubits, mostrando que su método escala bien, aunque el tiempo que toma crece a medida que el sistema se hace más grande (siguiendo aproximadamente una potencia de 4.448, lo cual es mucho mejor que la explosión exponencial que ocurriría con métodos más antiguos).

Sin embargo, el artículo es muy claro sobre lo que esto NO significa. No significa que las computadoras cuánticas sean inútiles. Los autores señalan que su método solo funciona para la "simulación débil", lo que significa calcular la energía, pero no puede realizar una "simulación fuerte" del proceso de muestreo de cadenas de bits aleatorias (el dato bruto de salida) que produce la computadora cuántica. El experimento original utilizó esas cadenas de bits aleatorias para limpiar errores y mejorar el resultado. Debido a que el método de los autores no genera esas cadenas aleatorias, no puede realizar ese tipo específico de corrección de errores. Pero, al ser tan rápidos, simplemente pudieron optimizar la configuración del circuito mejor que los experimentadores, superando el resultado final de todos modos.

La conclusión más importante es una línea divisoria. Los autores argumentan que los circuitos de una sola capa no son el "santo grial" de la ventaja cuántica. Si una computadora regular puede simularlos en un minuto, no están haciendo nada especial. Para realmente superar a las computadoras clásicas, necesitamos construir circuitos con dos o más capas (L ≥ 2). Las matemáticas se vuelven demasiado complicadas para que su truque funcione una vez que se añaden esas capas extra, y es ahí donde probablemente reside el verdadero poder cuántico. Así que, aunque este artículo muestra que podemos simular los experimentos actuales de "estado del arte" en química cuántica en una laptop, también nos dice que la verdadera carrera acaba de comenzar: necesitamos construir circuitos más profundos y complejos para encontrar los problemas que las computadoras clásicas verdaderamente no pueden resolver.

Resumen Técnico: Simulación Clásica Eficiente de Circuitos de Clúster Jastrow Unitarios de Gran Escala

Planteamiento del Problema
Experimentos recientes en química cuántica han utilizado el ansatz de Clúster Jastrow Unitario (UCJ), particularmente su variante local (LUCJ), para simular estados fundamentales de moléculas fuertemente correlacionadas en hardware cuántico. Experimentos notables, como los reportados en Sci. Adv. 11, 25 (2025), han ejecutado circuitos con hasta 77 cúbits y más de 10,000 puertas, requiriendo recursos masivos de supercomputación clásica (por ejemplo, 6400 nodos en Fugaku) para el post-procesamiento mediante la Diagonalización Cuántica Basada en Muestreo (SQD). Si bien se ha establecido la dificultad clásica de muestrear de los circuitos UCJ (vinculada a las computaciones IQP), la tarea específica de calcular el valor de la esperanza de la energía para circuitos UCJ de una sola capa no ha sido caracterizada completamente. Las estrategias de simulación clásica existentes para estos circuitos enfrentan un escalamiento exponencial o dependen de truncamientos que comprometen la precisión química. El problema central abordado es si la energía de los circuitos UCJ de una sola capa puede ser calculada exactamente y de manera eficiente en hardware clásico, independientemente de las restricciones de localidad impuestas al hardware cuántico.

Metodología
Los autores presentan un algoritmo clásico de tiempo polinomial para el cálculo exacto de la energía de cualquier circuito UCJ de una sola capa ( $L=1$ ). El algoritmo opera en el formalismo de segundo cuantizado y evita el crecimiento exponencial de términos típicamente asociado con la retropropagación en otros métodos de simulación. El enfoque procede en tres pasos distintos:

Retropropagación a través de Rotaciones Orbitales: El algoritmo primero propaga el Hamiltoniano fermiónico hacia atrás a través de los operadores de rotación orbital final ( $e^{K_\mu}$ ). Utilizando el Lema 1, las integrales de un electrón ( $h$ ) y de dos electrones ( $g$ ) se actualizan mediante transformaciones unitarias. Este paso se realiza exactamente en tiempo $O(N^3)$ y $O(N^5)$ respectivamente, donde $N$ es el número de orbitales, sin aumentar el número de términos en el Hamiltoniano.
Propagación a través de Operadores Jastrow: El Hamiltoniano es entonces propagado a través de los operadores Jastrow ( $e^{iJ_\mu}$ ). El Lema 2 demuestra que esta operación también puede realizarse exactamente en $O(N^5)$ de tiempo. Crucialmente, este paso no aumenta el conteo de términos; en su lugar, cada término en el Hamiltoniano acumula una fase escalar y una fase vectorial dependiente de los números de operador. Esto se mantiene válido tanto para los circuitos UCJ generales como para los circuitos LUCJ dispersos.
Evaluación de la Energía vía la Fórmula de Löwdin: La energía final se calcula como el valor de la esperanza del Hamiltoniano modificado contra el determinante de Slater de referencia. Debido a que la propagación de Jastrow introduce fases no ortogonales, los elementos de la matriz se evalúan entre determinantes de Slater no ortogonales. Los autores utilizan la fórmula de Löwdin (Lema 3) para calcular estos elementos de la matriz de manera eficiente. El cuello de botella computacional implica la formación y factorización de una matriz $S$ (relacionada con el solapamiento de los determinantes), lo que conduce a una complejidad algorítmica global de $O(N^7)$ . Los autores señalan que, en la práctica, la memorización de elementos de la matriz reduce el escalamiento observado a aproximadamente $O(N^{4.45})$ .

Contribuciones Clave

Algoritmo de Tiempo Polinomial: El artículo establece que calcular la energía de un circuito UCJ de una sola capa es clásicamente tratable en tiempo $O(N^7)$ , independiente de las restricciones de conectividad del hardware cuántico.
Simulación Exacta: A diferencia de métodos previos que dependen de truncamientos o aproximaciones heurísticas, este algoritmo proporciona un cálculo exacto del valor de la esperanza de la energía.
Distinción de la Simulación Fuerte: Los autores aclaran que su método constituye una "simulación débil" (calcular observables) en lugar de una "simulación fuerte" (muestrear cadenas de bits o bitstrings). Por consiguiente, el algoritmo no puede realizar el post-procesamiento basado en el muestreo de cadenas de bits como la SQD, que se utiliza en experimentos para mitigar el ruido del hardware.

Resultados
Los autores validan su algoritmo a través de dos casos de estudio primarios:

Clúster de Hierro-Azufre (72 Cúbits): El algoritmo fue aplicado al experimento del clúster de hierro-azufre de 77 cúbits (72 cúbits activos) de Sci. Adv. 11, 2025. Mientras que el experimento original requirió 6400 nodos en Fugaku durante 1.44 horas para lograr una energía de $-326.645$ Ha vía SQD, el algoritmo clásico de los autores simuló el mismo circuito en una computadora portátil en menos de un minuto. Además, optimizando los parámetros del UCJ de una sola capa utilizando el evaluador de energía clásico rápido (L-BFGS-B), lograron una energía de estado fundamental de $-326.796$ Ha, superando el resultado experimental obtenido mediante SQD.
Cadenas de Hidrógeno (Hasta 160 Cúbits): El algoritmo se escaló para simular cadenas de Hidrógeno de hasta 160 cúbits, representando los límites actuales de la computación cuántica de estado sólido de última generación. Los resultados mostraron que el ansatz UCJ de una sola capa, inicializado con parámetros de Coupled Cluster Singles and Doubles (CCSD), sigue de cerca las energías de Hartree-Fock (diferencia promedio de ~11 mHa). La optimización de los parámetros para sistemas de hasta 40 cúbits produjo energías comparables a CCSD.

Significancia y Reivindicaciones
El artículo concluye que los circuitos UCJ de una sola capa son insuficientes para demostrar la ventaja cuántica en química. La capacidad de simular estos circuitos de manera exacta y eficiente en hardware clásico implica que la "ventaja cuántica" reclamada en experimentos recientes de gran escala depende del post-procesamiento específico (SQD) o de la necesidad potencial de circuitos más profundos.

Los autores afirman explícitamente que lograr una computación más allá de lo clásico probablemente requiere $L \ge 2$ capas del ansatz (L)UCJ. Argumentan que su algoritmo falla para $L \ge 2$ porque las rotaciones orbitales entre las capas de Jastrow rompen la estructura diagonal requerida para la fórmula de Löwdin. Dado que la alternancia de capas de rotaciones de Givens y puertas CZ es conocida por ser universal para la química cuántica, los autores sugieren que, si bien $L \ge 2$ es necesario para la ventaja cuántica, puede no ser suficiente sin abordar la significativa sobrecarga de puertas (gate overhead) en el hardware actual. El trabajo sirve como un punto de referencia, demostrando que los experimentos actuales de una sola capa aún no exceden las capacidades de los algoritmos clásicos optimizados.

Más como este