Optimal fermion-qubit mappings via quadratic assignment

Título: Mapeos óptimos fermión-qubit mediante asignación cuadrática

1. El Problema

La simulación de sistemas fermiónicos es una de las aplicaciones más prometedoras de la computación cuántica, abarcando desde la química cuántica hasta la física de altas energías. Sin embargo, existe un desafío fundamental: traducir las interacciones fermiónicas a operadores y estados de qubits (mapeo fermión-qubit).

• Dilema de Recursos: La literatura actual presenta una dicotomía entre dos paradigmas:
  • Mapeos sin áncilas (Ancilla-free): Minimizan la cantidad de qubits necesarios (igual al número de modos fermiónicos), pero a menudo resultan en operadores de Pauli no locales con pesos altos (ej. la transformación de Jordan-Wigner tiene un peso que escala linealmente con el tamaño del sistema).
  • Codificaciones locales: Minimizan la complejidad de las puertas (peso de Pauli constante), pero requieren un número de qubits áncila que escala linealmente con el tamaño del sistema, lo cual es prohibitivo en hardware actual donde los qubits lógicos son escasos.
• Objetivo: Encontrar un equilibrio que reduzca la complejidad de las puertas (peso de Pauli) sin incurrir en un costo excesivo de qubits, optimizando tanto el número total como el máximo de qubits involucrados en cada término del Hamiltoniano.

2. Metodología

Los autores proponen dos enfoques computacionales para construir mapeos generales bajo una restricción limitada de qubits:

A. Optimización del Orden Fermiónico como Problema de Asignación Cuadrática (QAP)

• Concepto: El orden en que se etiquetan los modos fermiónicos es arbitrario desde el punto de vista físico, pero afecta drásticamente el peso de Pauli en el Hamiltoniano de qubits resultante.
• Enfoque: Se formula la búsqueda del orden óptimo de los modos fermiónicos como una instancia del Problema de Asignación Cuadrática (QAP).
• Aplicación: Se aplica a una clase prominente de mapeos sin áncilas conocidos como "codificaciones lineales" (Jordan-Wigner, Bravyi-Kitaev, Ternary Tree y Parity Basis).
• Función de Costo: Se define una matriz de costos $D$ que calcula el peso de Pauli (total o máximo) de los términos de salto (hopping) e interacción basándose en la estructura de la matriz de transformación $U$ y el orden $\sigma$. El objetivo es minimizar $\sum D_{ij}$ o $\max D_{ij}$.

B. Adición Incremental de Qubits Áncila a Jordan-Wigner

• Concepto: Extender la estrategia de [29] para agregar áncilas de manera incremental a la transformación de Jordan-Wigner, en lugar de usar codificaciones locales completas que requieren $O(n)$ áncilas.
• Mecanismo:
  1. Se comienza con un Hamiltoniano de Jordan-Wigner optimizado por orden.
  2. Se introduce un qubit áncila adicional.
  3. Para cada término de salto, se elige entre dos opciones para modificar el operador de Pauli:
     • Opción 1: Multiplicar por una cadena de $Z$ ($Z_I$) en un subconjunto de qubits y un operador en el áncila.
     • Opción 2: Multiplicar solo por un operador en el áncila.
  4. La elección se optimiza para cancelar cadenas largas de operadores $Z$ (que causan no-localidad) reduciendo el peso de Pauli, manteniendo las relaciones de conmutación requeridas (Teorema 1).
• Optimización: Se resuelve un problema de optimización para encontrar los subconjuntos de qubits $I$ que minimizan el aumento del peso de Pauli al agregar cada áncila. Esto se repite iterativamente hasta agotar el presupuesto de áncilas o dejar de haber mejoras.

3. Contribuciones Clave

1. Generalización de la Optimización de Orden: Se extiende la optimización de la etiqueta fermiónica más allá de Jordan-Wigner para incluir Bravyi-Kitaev, Ternary Tree y Parity Basis, formulándola como un QAP aplicable a cualquier patrón de salto e interacción (no solo redes cuadradas).
2. Protocolo de Áncilas Constantes: Se desarrolla un método general para agregar cualquier número de qubits áncila a la transformación de Jordan-Wigner, creando mapeos de "ánclas constantes" que superan a los mapeos sin áncilas.
3. Reducción Significativa de Peso: Demostración de que con un número pequeño de áncilas (hasta 10), es posible reducir drásticamente el peso de Pauli, superando incluso a las mejores transformaciones sin áncilas conocidas.
4. Reproducibilidad: Se proporciona código fuente y datos para la asignación de qubits.

4. Resultados

• Comparación de Codificaciones (Sin Áncilas):
  • En sistemas pequeños, Jordan-Wigner con orden optimizado es competitivo.
  • Para sistemas más grandes (hasta 225 modos), la transformación Ternary Tree (y Bravyi-Kitaev) tiende a superar a Jordan-Wigner en términos de peso total y máximo de Pauli.
  • Se recuperó el ordenamiento óptimo conocido (Mitchison–Durbin) para redes cuadradas, validando el método.
• Impacto de las Áncilas (Jordan-Wigner + Áncilas):
  • En una red de $8 \times 8$ (64 modos), la adición de 10 qubits áncila redujo el peso total de Pauli de los términos de salto en un 39.3% comparado con Jordan-Wigner sin áncilas.
  • En diversos grafos de 64 modos (redes hexagonales, triangulares, aleatorias), el método con 10 áncilas redujo el peso de Pauli hasta en un 67%.
  • Comparación Directa: Los mapeos de Jordan-Wigner con 10 áncilas lograron un peso total de Pauli menor o comparable al de la transformación Ternary Tree (que es sin áncilas), demostrando que un pequeño costo en qubits ofrece una ventaja masiva en complejidad de puertas.

5. Significancia

Este trabajo es crucial para la computación cuántica de mediano plazo (NISQ) y el futuro inmediato:

• Equilibrio de Recursos: Demuestra que no es necesario elegir estrictamente entre "pocos qubits" o "baja complejidad de puertas". Se puede sacrificar un número muy pequeño de qubits (constante, no escalable con $n$) para obtener una reducción drástica en la complejidad de las puertas.
• Viabilidad de Simulaciones: Al reducir el peso de Pauli (y por ende la profundidad de circuito y el número de mediciones necesarias), se hacen factibles simulaciones de sistemas fermiónicos más grandes en hardware actual.
• Nueva Dirección de Investigación: Abre la puerta a optimizaciones computacionales asistidas para encontrar mapeos específicos para problemas concretos, superando las limitaciones de las transformaciones analíticas fijas.

En resumen, el paper establece que la optimización computacional del orden de los modos y la adición estratégica de un número limitado de qubits áncila pueden superar el estado del arte en la eficiencia de los mapeos fermión-qubit, haciendo que la simulación cuántica sea más accesible para dispositivos con recursos limitados.