Digital-Analog Quantum Computing with Qudits

1. El Problema

La computación cuántica actual opera principalmente en el régimen de Ruido Intermedio a Escala (NISQ). Existen dos paradigmas principales con ventajas y desventajas complementarias:

• Computación Cuántica Analógica (AQC): Es robusta frente al ruido pero carece de versatilidad para la computación universal.
• Computación Cuántica Digital (DQC): Es universal y flexible, pero extremadamente susceptible al ruido, requiriendo corrección de errores que impone una sobrecarga masiva de recursos físicos.

El paradigma Digital-Analog (DAQC) busca combinar la robustez de la AQC (usando la Hamiltoniana natural del sistema como recurso de entrelazamiento) con la flexibilidad de la DQC (usando puertas locales). Sin embargo, la mayoría de los protocolos DAQC existentes están diseñados para qubits (sistemas de 2 niveles).

El desafío abordado en este trabajo es generalizar el marco DAQC a sistemas de $d$-niveles (qudits). Los qudits ofrecen ventajas significativas, como una mayor capacidad computacional por partícula, una reducción en el número de operaciones de entrelazamiento necesarias y una plataforma natural para simular fenómenos de espín de alta dimensión y teorías de gauge. El problema específico es cómo simular Hamiltonianas objetivo arbitrarias de dos cuerpos en qudits utilizando una Hamiltoniana fuente disponible, manteniendo la eficiencia en la complejidad del circuito y la robustez ante el ruido.

2. Metodología

Los autores proponen una extensión del protocolo DAQC para qudits basada en los siguientes pilares:

• Base de Operadores Weyl-Heisenberg: En lugar de la base de Pauli utilizada para qubits, el protocolo utiliza la base completa de operadores Weyl-Heisenberg ($W_{nm} = Z^n X^m$) para definir las operaciones locales. Estos operadores generalizan las matrices de Pauli mediante el formalismo de reloj y desplazamiento.
• Conjugación de Bloques Analógicos: El objetivo es sintetizar un propagador $U_P = e^{-iT H_P}$ (donde $H_P$ es la Hamiltoniana objetivo) a partir de una Hamiltoniana fuente $H_S$. Esto se logra mediante una secuencia de conjugaciones:
  $$U_P = \prod_q G_q^\dagger e^{-i t_q H_S} G_q$$
  Donde $G_q$ son puertas unitarias locales de un solo qudit seleccionadas de la base Weyl-Heisenberg.
• Formulación del Sistema Lineal: El problema se reduce a resolver un sistema de ecuaciones lineales $M \vec{t} = T \vec{h}_{p/h_s}$, donde:
  • $\vec{t}$ son las duraciones de los bloques analógicos.
  • La matriz $M$ codifica los cambios de fase (signos o factores de fase complejos) introducidos por las conjugaciones en los acoplamientos de la Hamiltoniana.
  • El número de filas de $M$ corresponde a los acoplamientos (hasta $O(d^4 n^2)$) y las columnas a las configuraciones de conjugación.
• Prueba de Existencia de Solución: Mediante el Teorema de Carathéodory y un análisis algebraico de las submatrices de $M$ (específicamente el bloque $M[2]$ que representa conjugaciones simultáneas en dos qudits), los autores demuestran que el determinante de la submatriz es no nulo y que existe una solución con tiempos positivos ($t_q > 0$).
• Esquema de Implementación (bDAQC): Se utiliza el enfoque "banged-DAQC" (bDAQC), donde la Hamiltoniana permanece activa continuamente mientras se aplican las puertas digitales sobre la dinámica continua, lo cual es más factible experimentalmente que encender/apagar la interacción.

3. Contribuciones Clave

1. Generalización a Qudits: Se extiende el marco DAQC desde sistemas de 2 niveles ($d=2$) a sistemas de dimensión arbitraria $d$, utilizando la base Weyl-Heisenberg como conjunto natural de operaciones.
2. Protocolo Constructivo: Se proporciona un método explícito para diseñar circuitos DAQC que simulan Hamiltonianas de dos cuerpos arbitrarias en qudits.
3. Eficiencia en Recursos: El método requiere resolver un sistema lineal, pero la implementación óptima necesita solo $O(d^4 n^2)$ bloques analógicos, lo cual es polinomialmente más eficiente que protocolos anteriores que intentaban aislar cada acoplamiento individualmente.
4. Aplicabilidad a Espines SU(2): Se demuestra que para Hamiltonianas expresadas en la base de operadores de espín $\{S_x, S_y, S_z\}$ (dinámica SU(2)), la matriz $M$ puede simplificarse, reduciendo la complejidad a $O(n^2)$ y el número de columnas a $4n$, similar al caso de qubits, siempre que las operaciones $e^{-i\pi S_\nu}$ estén disponibles.
5. Generación de Nuevos Operadores: El formalismo permite transformar elementos de una base en otros que comparten componentes Weyl idénticos, permitiendo generar operadores efectivos (por ejemplo, transformar operadores cuadráticos de espín en lineales).

4. Resultados

• Simulación de Modelos de Espín: Se ilustra la validez del protocolo simulando un modelo Ising bilineal-bicuadrático (BLBQ) para espines $s=1$ (qutrits, $d=3$).
• Comparación de Fidelidad: Se comparó el rendimiento del protocolo bDAQC frente a circuitos puramente digitales (DQC) para $n=6$ qutrits y diferentes ángulos $\theta$.
  • El protocolo bDAQC mostró mayores fidelidades en el rango $\theta < \pi/2$.
  • La fidelidad disminuye en ciertos rangos debido a la necesidad de más puertas de un solo qudit (que introducen errores), pero se recupera cerca de $\theta = \pi$ cuando algunos bloques son demasiado cortos y se descartan, reduciendo el número de puertas.
  • Se consideraron errores de puerta (fidelidad del 95% para puertas de dos qutrits y ~99.4% para puertas de un qutrit) y decoherencia.
• Viabilidad Experimental: El enfoque es compatible con plataformas de vanguardia como trampas de iones de $d$-niveles y circuitos superconductores, donde ya se han demostrado procesadores basados en qudits.

5. Significado e Impacto

Este trabajo es fundamental por varias razones:

• Ampliación del Régimen NISQ: Ofrece una vía para ejecutar algoritmos cuánticos complejos en dispositivos ruidosos actuales utilizando la robustez de la simulación analógica, pero con la flexibilidad de la programación digital.
• Optimización de Recursos: Al utilizar qudits, se reduce la necesidad de mapeos de múltiples qubits para representar estados de alta dimensión (crucial en simulaciones bosónicas y teorías de gauge), disminuyendo la profundidad del circuito y la exposición al ruido.
• Nuevas Capacidades de Simulación: Habilita la simulación eficiente de sistemas con términos de cuadrupolo magnético y otras interacciones de alta dimensión que son difíciles de modelar con qubits.
• Marco Experimental Sólido: Proporciona un marco teórico completo para implementar circuitos DAQC en plataformas reales que ya soportan qudits, sentando las bases para futuros algoritmos escalables y la simulación de teorías de gauge con menos recursos cuánticos.

En resumen, el artículo establece un puente teórico y práctico entre la computación analógica robusta y la digital flexible, adaptado específicamente para explotar las ventajas de los sistemas de alta dimensión (qudits) en la era post-NISQ.