Multiconfiguration Pair-Density Functional Theory Calculations of Low-lying States of Complex Chemical Systems with Quantum Computers

A continuación presento un resumen técnico detallado del artículo "Multiconfiguration Pair-Density Functional Theory Calculations of Low-lying States of Complex Chemical Systems with Quantum Computers" en español.

1. El Problema

La simulación precisa de la estructura electrónica molecular es fundamental para predecir reactividad y propiedades espectroscópicas, pero enfrenta un desafío central: la descripción simultánea de la correlación electrónica estática (asociada a la deslocalización de electrones y estados de múltiples referencias) y la correlación dinámica (interacciones de corto alcance entre electrones).

• Limitaciones Clásicas: Los métodos de un solo determinante (como Hartree-Fock o DFT estándar) fallan en sistemas fuertemente correlacionados. Los métodos multireferencia clásicos (como CASSCF) capturan la correlación estática pero carecen de precisión en la dinámica, mientras que métodos de alto nivel (como CCSD(T)) son computacionalmente prohibitivos para sistemas grandes.
• Limitaciones Cuánticas (NISQ): Los algoritmos cuánticos actuales, como el Variational Quantum Eigensolver (VQE) con ansatz de tipo UCCSD, requieren circuitos excesivamente profundos para capturar la correlación dinámica completa, lo que excede las capacidades de los dispositivos cuánticos ruidosos de escala intermedia (NISQ) debido a la decoherencia y el ruido.

2. Metodología: VQE-MC-PDFT Híbrido

Los autores proponen una estrategia híbrida que combina el VQE (cuántico) con la Teoría de Funcional de Densidad de Pareja Multiconfiguracional (MC-PDFT) (clásica) para desacoplar y tratar eficientemente ambos tipos de correlación.

• División del Trabajo:
  • Correlación Estática (Cuántica): Se trata dentro de un espacio activo químico utilizando un circuito cuántico compacto basado en un ansatz CASCI (Configuración Interacción de Espacio Activo Completo). El VQE prepara el estado multireferencia y mide las matrices de densidad reducida (RDMs) de 1 y 2 cuerpos.
  • Correlación Dinámica (Clásica): Se recupera mediante un funcional de densidad de pareja "on-top" (sobre la densidad) evaluado clásicamente a partir de las RDMs medidas. Se utiliza el funcional ftPBE (fully-translated PBE), que depende de la densidad electrónica $\rho(r)$ y la densidad de pareja on-top $\Pi(r)$.
• Optimización Autoconsistente: A diferencia de enfoques previos, este método introduce una optimización autoconsistente de los orbitales moleculares. Se calcula el gradiente de rotación orbital derivado del funcional MC-PDFT y se codifica como una unidad cuántica ($\hat{U}_{orb}$) que se aplica al circuito, permitiendo que tanto los coeficientes de configuración como la forma de los orbitales se optimicen variacionalmente.
• Mitigación de Errores: Se implementa un protocolo avanzado de mitigación de errores de lectura inspirado en el Método de Elementos Finitos (FEM). Este divide el sistema de qubits en clústeres para modelar correlaciones locales de ruido, superando las limitaciones de los modelos de tensor product estándar. Además, se emplean técnicas de Extrapolación a Ruido Cero (ZNE) y Ajuste de Clifford (CF).
• Estrategia de Partición de Circuitos: Para escalar a espacios activos grandes, se utiliza una estrategia de partición de circuitos que descompone cúmulos débilmente entrelazados, permitiendo reconstrucción clásica y simulación de sistemas que requieren más qubits de los disponibles físicamente.

3. Contribuciones Clave

1. Algoritmo VQE-MC-PDFT: Desarrollo de un flujo de trabajo híbrido que reduce la profundidad del circuito cuántico al evitar la necesidad de circuitos de correlación dinámica profundos (como UCCSD), confinando la parte cuántica a la correlación estática.
2. Optimización de Orbitales Variacional: Integración exitosa de la optimización de orbitales basada en gradientes de MC-PDFT dentro del bucle VQE, mejorando la precisión sin aumentar significativamente la complejidad del hardware.
3. Escalabilidad mediante Partición: Demostración de cálculos en espacios activos masivos (hasta 48 electrones en 42 orbitales) mediante la partición de circuitos y reconstrucción, superando las limitaciones físicas de los qubits disponibles.
4. Validación Experimental: Ejecución exitosa en el dispositivo superconductor de 13 qubits Tianji-S2, demostrando la viabilidad práctica en hardware real ruidoso.

4. Resultados Principales

El método fue validado en tres sistemas de referencia desafiantes:

• Dímero de Carbono ($C_2$):

  • Se calcularon curvas de energía potencial (PECs) para el estado fundamental y estados excitados.
  • Precisión: Se logró un error absoluto medio (MAE) de 0.006 Å en longitudes de enlace y 0.048 eV en energías de excitación vertical/adiabática, comparables a la precisión experimental y superando significativamente a CASSCF puro.
  • El método reprodujo correctamente la intersección cónica y la estructura de pozos de potencial, mostrando robustez frente al ruido del hardware.

• Dímero de Cromo ($Cr_2$):

  • Este sistema presenta un enlace sextuple y requiere un espacio activo enorme (48e, 42o) para una descripción correcta.
  • Superación de Limitaciones: Mientras que VQE-CASSCF (sin funcional dinámico) predijo curvas repulsivas incorrectas, VQE-MC-PDFT recuperó una curva de energía potencial unida (bound).
  • Escalabilidad: Mediante partición de circuitos, se simuló un espacio equivalente a 84 qubits. El resultado de energía de enlace (-0.465 eV) convergió sistemáticamente hacia los valores de referencia clásicos de alta precisión (MPS-LCC y SHCI), situándose entre ellos.

• Benceno ($C_6H_6$):

  • Se calcularon energías de excitación $\pi \to \pi^*$.
  • El método obtuvo un MAE de 0.048 eV frente a estimaciones teóricas de mejor calidad (QUEST), superando al TDDFT estándar (MAE 0.37 eV) y siendo competitivo con métodos Variational Monte Carlo avanzados (NES-VMC).

5. Significado e Impacto

Este trabajo establece un camino práctico para realizar predicciones químicas confiables en hardware cuántico de la era NISQ. Al separar estratégicamente la correlación estática (tratada cuánticamente) de la dinámica (tratada clásicamente mediante funcionales de densidad), el método:

• Reduce drásticamente la profundidad de los circuitos y el número de puertas de entrelazamiento necesarias.
• Mantiene el rigor físico necesario para sistemas fuertemente correlacionados.
• Demuestra que la combinación de algoritmos cuánticos con teoría de funcionales de densidad clásica es una ruta viable para superar las limitaciones actuales de los ordenadores cuánticos, permitiendo el estudio de moléculas complejas que antes eran inaccesibles.

En resumen, el estudio valida que la separación de tipos de correlación dentro de un flujo de trabajo variacional cuántico-clásico es una solución robusta y escalable para la química cuántica en dispositivos reales.