Resource Estimation for VQE on Small Molecules: Impact of Fermion Mappings and Hamiltonian Reductions

Título: Estimación de Recursos para VQE en Moléculas Pequeñas: Impacto de las Mapeos de Fermiones y Reducciones del Hamiltoniano

1. El Problema

La determinación precisa de las energías del estado fundamental de las moléculas es un desafío fundamental en la química cuántica, con implicaciones críticas para el descubrimiento de fármacos y el diseño de materiales. Aunque el algoritmo Variational Quantum Eigensolver (VQE) con el ansatz UCCSD (Unitary Coupled Cluster Singles and Doubles) es un paradigma híbrido prometedor para abordar este problema, su implementación práctica en hardware actual (era NISQ - Noisy Intermediate-Scale Quantum) y futura (FASQ - Fault-Tolerant Application-Scalable Quantum) está severamente limitada por:

• El ruido y la decoherencia en los dispositivos actuales.
• La escalabilidad restringida del hardware (número limitado de qubits y profundidad de circuitos).
• La falta de una comprensión detallada de cómo los recursos computacionales (número de qubits, puertas cuánticas, cadenas de Pauli) escalan con la complejidad molecular y las estrategias de codificación.

Existe una necesidad urgente de estimaciones de recursos a priori para determinar qué simulaciones son realistas y cómo optimizarlas mediante la elección de mapeos fermión-qubit y técnicas de reducción del Hamiltoniano.

2. Metodología

Los autores realizaron un análisis sistemático y agnóstico al hardware de los requisitos de recursos para VQE con UCCSD. El flujo de trabajo se dividió en cuatro etapas secuenciales:

1. Modelado del Hamiltoniano: Se definieron las geometrías moleculares y se calculó el Hamiltoniano electrónico en segunda cuantización ($\hat{H}_{sq}$) utilizando el paquete químico cuántico PySCF. Se aplicó la aproximación de núcleo congelado (Frozen-Core, FC) para eliminar orbitales internos no correlacionados.
2. Mapeo Fermión-Qubit (FTQM): El Hamiltoniano fermiónico se transformó en un Hamiltoniano de qubits ($\hat{H}_{qubit}$) utilizando tres esquemas de mapeo estándar:
   • Jordan-Wigner (JW): Mapeo directo pero con cadenas de Pauli largas.
   • Bravyi-Kitaev (BK): Mapeo que distribuye la información de paridad y ocupación, reduciendo la localidad del operador.
   • Paridad (Pa): Codifica la información de paridad localmente en cada qubit.
3. Reducción del Hamiltoniano: Se aplicaron estrategias de reducción para minimizar los recursos:
   • Aproximación de Núcleo Congelado (FC): Reduce el espacio de orbitales activos.
   • Poda de Qubits por Simetría $Z_2$ (Tapering): Elimina qubits redundantes explotando simetrías del grupo abeliano $Z_2$ (como la paridad del número de partículas) presentes en el Hamiltoniano.
4. Preparación de Circuitos y Compilación: Se construyó el ansatz UCCSD y se compiló para un simulador (Qiskit-Aer) con optimización de nivel 3. Se analizaron métricas físicas de nivel de circuito: ancho de circuito (número de qubits), profundidad, conteo total de puertas (unión de puertas de 1 y 2 qubits) y número de cadenas de Pauli únicas.

El estudio cubrió un conjunto representativo de 13 moléculas (desde $H_2$ hasta $C_2H_4$) utilizando la base STO-3G.

3. Contribuciones Clave

• Análisis Comparativo Exhaustivo: Se proporciona una evaluación sistemática de tres mapeos (JW, BK, Pa) bajo cuatro configuraciones de reducción (sin reducción, solo poda, solo núcleo congelado, y ambas combinadas) para un conjunto diverso de moléculas.
• Desacoplamiento de Efectos: Se aíslan los efectos dependientes del mapeo de las limitaciones específicas del hardware, proporcionando una línea base física clara.
• Guías Prácticas de Selección: Se establecen criterios claros para seleccionar el mapeo y la estrategia de reducción óptima según la simetría molecular y la disponibilidad de orbitales de núcleo.
• Datos de Referencia: Se publican métricas detalladas (conteo de puertas, cadenas de Pauli, qubits) que sirven como base para futuros análisis de recursos lógicos en sistemas tolerantes a fallos.

4. Resultados Principales

• Reducción de Qubits:

  • La combinación de núcleo congelado (FC) y poda por simetría (T) puede reducir el número de qubits en hasta un 50% en comparación con la representación sin reducir.
  • El mapeo de Paridad (Pa) ofrece una ventaja inherente de 2 qubits en la codificación inicial para ciertos sistemas, pero esto reduce la cantidad de simetrías adicionales disponibles para la poda posterior. En moléculas como $NH_3$ y $CH_4$, la poda sola no redujo más qubits con el mapeo Pa.
  • JW y BK muestran reducciones idénticas en qubits cuando se aplica la poda, ya que la reducción depende de la simetría molecular y no del mapeo en sí.

• Reducción de Puertas Cuánticas:

  • La reducción de puertas es altamente variable dependiendo de la molécula. La combinación de FC y T logra reducciones de hasta 27.5 veces (ej. en $BeH_2$).
  • Para la mayoría de las moléculas, la reducción media es de aproximadamente 9.55x con JW, 8.53x con BK y 7.55x con Pa cuando se aplican ambas técnicas.
  • Hallazgo Crítico: En configuraciones de "solo poda" (sin núcleo congelado), el mapeo Paridad (Pa) puede resultar en un aumento marginal del conteo de puertas (factor < 1.0) para moléculas como $NH_3$ y $CH_4$. Esto se debe a que la codificación de paridad redistribuye la información de ocupación, aumentando la estructura de entrelazamiento de las cadenas de Pauli si no se compensa con la reducción de orbitales. Por lo tanto, se recomienda aplicar la aproximación de núcleo congelado como requisito mínimo al usar el mapeo Pa.

• Reducción de Cadenas de Pauli (Costo de Medición):

  • El conteo de cadenas de Pauli únicas (que determina el costo de medición) es independiente del mapeo (JW, BK o Pa producen el mismo número de cadenas únicas para una configuración dada).
  • La aproximación de núcleo congelado (FC) es el mecanismo dominante para reducir el costo de medición, logrando reducciones medias de 2.09x.
  • La poda por simetría ($T$) tiene un impacto mínimo en el número de cadenas de Pauli (reducción media de 1.03x), ya que principalmente simplifica los pesos de los operadores sin eliminar bases de medición únicas, excepto en sistemas de alta simetría como $H_2$.

5. Significado e Impacto

Este trabajo establece un marco de referencia esencial para la ejecución de simulaciones químicas en hardware cuántico actual y futuro:

1. Viabilidad en NISQ: Demuestra que mediante la combinación inteligente de mapeos y reducciones (especialmente núcleo congelado), es posible reducir drásticamente los requisitos de recursos, haciendo factibles simulaciones de moléculas más complejas en dispositivos con pocos qubits.
2. Guía de Diseño de Algoritmos: Proporciona una jerarquía práctica para la selección de estrategias:
   • Si solo es posible la poda: Preferir JW o BK.
   • Si está disponible el núcleo congelado: Los tres mapeos son competitivos.
   • Si se aplican ambas reducciones: JW y BK ofrecen los mayores ahorros agregados de puertas, con BK mostrando ventajas adicionales en el número total de puertas para moléculas más grandes.
3. Base para Futuras Investigaciones: Los resultados físicos de nivel de circuito sirven como punto de partida para estimaciones de recursos lógicos (contando qubits lógicos y puertas T) en sistemas tolerantes a fallos, y destacan la necesidad de integrar futuras técnicas de mapeo (como árboles ternarios) y formulaciones no Born-Oppenheimer en flujos de trabajo de estimación de recursos.

En resumen, el estudio confirma que la optimización de la estructura del Hamiltoniano y la elección del mapeo no son meros detalles técnicos, sino determinantes críticos para la viabilidad de la química cuántica en la era NISQ y más allá.