Digitized counterdiabatic quantum critical dynamics

1. Planteamiento del Problema

El trabajo aborda un desafío fundamental en la computación cuántica y la simulación de materia condensada: la formación de defectos topológicos (como paredes de dominio o vórtices) cuando un sistema cuántico es impulsado a través de una transición de fase a una velocidad finita.

• Mecanismo de Kibble-Zurek (QKZM): En regímenes lentos (adiabáticos), la densidad de defectos sigue una ley de potencia universal determinada por la clase de universalidad de la transición. Sin embargo, en quench rápidos (cambios rápidos de parámetros), el sistema no tiene tiempo de relajarse, y la densidad de defectos se satura en un valor constante (un "plateau") que depende de la profundidad del quench, no de la velocidad.
• Limitaciones de la Computación Cuántica Adiabática: Para evitar defectos, se requiere un tiempo de evolución muy largo (lento), lo cual excede los tiempos de coherencia actuales de los procesadores cuánticos.
• El Desafío: ¿Es posible suprimir la formación de defectos en regímenes de quench rápido utilizando protocolos de control avanzados en hardware cuántico ruidoso (NISQ)?

2. Metodología

Los autores implementaron y probaron experimentalmente un protocolo de conducción contraadiabática (Counterdiabatic driving, CD) digitalizado en procesadores cuánticos superconductores de IBM.

• Protocolo CD: Se añade un término auxiliar al Hamiltoniano ($H_{CD} = \dot{\lambda} A_{\lambda}$) que cancela exactamente las transiciones a estados excitados. Dado que el potencial de gauge adiabático ($A_{\lambda}$) es no local y difícil de implementar, se utiliza una aproximación de conmutadores anidados (Nested-Commutator, NC) de primer orden para obtener un término local y tratable.
• Digitalización: El tiempo de evolución se discretiza utilizando la descomposición de Suzuki-Trotter de primer orden. El circuito cuántico se construye con puertas de un y dos qubits nativas de los procesadores.
• Sistemas Simulados:
  • Se estudiaron modelos de Ising con campo transversal (TFIM) en una y dos dimensiones.
  • Geometrías: Cadena 1D ($N=100$), red hexagonal pesada (heavy-hexagonal) de 156 qubits, escalera de tres patas y red cuadrada 2D.
• Hardware: Experimentos realizados en los procesadores cuánticos de IBM Fez y Marrakesh (hasta 156 qubits) a través de la nube.
• Análisis de Datos: Se midió la densidad de defectos (kinks en la magnetización) y se analizaron sus primeros tres cumulantes (media $\kappa_1$, varianza $\kappa_2$ y asimetría $\kappa_3$) para caracterizar la distribución completa de defectos.
• Validación: Los resultados experimentales se compararon con soluciones analíticas exactas (en 1D) y simulaciones numéricas avanzadas utilizando Estados de Producto Matricial (MPS) con el principio variacional dependiente del tiempo (TDVP) para las geometrías 2D.

3. Contribuciones Clave

1. Primera Evidencia Experimental: Demostración experimental de que los protocolos CD reducen la formación de defectos topológicos en quench rápidos a través de transiciones de fase cuánticas en sistemas de muchos cuerpos.
2. Supresión de Defectos en Regímenes Rápidos: Se logró una reducción de la densidad de defectos de hasta un 48% en el régimen de quench rápido, un resultado difícil de alcanzar con el recocido cuántico digitalizado estándar bajo los niveles actuales de ruido.
3. Caracterización de la Ruptura de Escala: Se confirmó experimentalmente y teóricamente la ruptura de la ley de escala de Kibble-Zurek en quench rápidos, donde la densidad de defectos se satura. El protocolo CD logra reducir el valor de este plateau.
4. Análisis de Cumulantes: Más allá de la media, el estudio caracterizó la varianza y la asimetría de la distribución de defectos, revelando comportamientos complejos dependientes de la geometría (por ejemplo, cambios en el signo de la asimetría entre geometrías 1D y 2D).
5. Solución Analítica en 1D: Derivación de una solución analítica exacta para la distribución de defectos en el límite de quench rápido para el modelo 1D, validando los resultados experimentales.

4. Resultados Principales

• Reducción de Defectos: En todas las geometrías probadas (1D, heavy-hexagonal, escalera, cuadrada), la aplicación del protocolo CD redujo la densidad media de defectos ($\kappa_1$) en comparación con el recocido digitalizado sin CD.
• Comportamiento Temporal:
  • En tiempos de evolución cortos ($T \lesssim 0.5/J$), la densidad de defectos sin CD se mantiene en un plateau alto (cerca del estado inicial desordenado).
  • Con CD, este plateau se reduce significativamente incluso para tiempos muy cortos ($T \to 0$), donde el operador de evolución unitaria se vuelve independiente del tiempo total.
• Impacto del Ruido: Los resultados experimentales coinciden estrechamente con las predicciones teóricas y numéricas en tiempos cortos. Sin embargo, a tiempos más largos, el ruido del hardware (decoherencia y errores de puertas) causa desviaciones, aunque la ventaja del protocolo CD se mantiene en el régimen de interés.
• Geometrías 2D: Para geometrías donde las simulaciones clásicas son difíciles, el uso de métodos MPS permitió validar que la reducción de defectos es un fenómeno robusto, aunque la implementación en hardware 2D (red cuadrada) sufrió más errores debido a la necesidad de operaciones de intercambio (SWAP) para acoplar qubits distantes.
• Asimetría ($\kappa_3$): Se observó que la asimetría de la distribución de defectos cambia de signo dependiendo de la geometría (positiva en 1D y heavy-hexagonal, negativa en la red cuadrada), indicando que el quench rápido favorece eventos raros con densidades de defectos por debajo de la media en ciertas configuraciones 2D.

5. Significado e Impacto

Este trabajo es significativo por varias razones:

• Validación de Protocolos de "Shortcuts to Adiabaticity": Confirma que las técnicas de atajo a la adiabaticidad, específicamente la conducción contraadiabática, son viables y efectivas en hardware cuántico real, no solo en simulaciones ideales.
• Mejora para la Optimización Cuántica: Dado que la formación de defectos reduce la fidelidad del estado final, este método ofrece una vía práctica para mejorar la precisión en algoritmos de optimización cuántica (como QAOA) y la preparación de estados, especialmente en problemas donde la adiabaticidad estricta es inalcanzable debido a limitaciones de tiempo.
• Simulación de Materiales: Facilita la simulación de dinámica crítica en materiales y sistemas cuánticos complejos en regímenes no explorados previamente, permitiendo estudiar la física de no equilibrio en condiciones más realistas.
• Escalabilidad: Demuestra la utilidad de los procesadores cuánticos de escala media (NISQ) para realizar simulaciones de sistemas de muchos cuerpos en 2D, superando las limitaciones de las simulaciones clásicas en ciertos regímenes de tiempo y precisión.

En resumen, el artículo establece que los protocolos contraadiabáticos digitalizados son una herramienta poderosa para controlar la dinámica crítica cuántica, mitigando la formación de defectos topológicos en regímenes rápidos y ruidosos, lo cual es crucial para el desarrollo de algoritmos cuánticos robustos y la simulación de materiales.