Spectral amplification for ground-state energy estimation of electronic structure in first quantization

Resumen Técnico: Amplificación Espectral para la Estimación de la Energía del Estado Fundamental en la Primera Cuantización

Planteamiento del Problema
La simulación de la estructura electrónica de sistemas químicos y de materiales es un objetivo primordial para los algoritmos cuánticos. Mientras que la mayor parte del trabajo se centra en representaciones de segunda cuantización, los enfoques de primera cuantización ofrecen una ventaja exponencial de espacio al escalar la memoria como $O(\eta \log N)$ (donde $\eta$ es el número de electrones y $N$ es el número de funciones de base), comparado con $O(N)$ para la segunda cuantización. Esto hace que la primera cuantización sea particularmente adecuada para abordar el límite de conjuntos de bases extensos y simulaciones de materiales que requieren la convergencia al límite termodinámico.

Sin embargo, a pesar de tener un escalado de memoria favorable, la complejidad de las compuertas para la estimación de la energía del estado fundamental en la primera cuantización sigue siendo prohibitivamente alta cuando se utilizan conjuntos de bases de ondas planas de gran tamaño. El cuello de botella principal es la constante de normalización del codificado de bloques ($\lambda$) requerido para la Estimación de la Fase Cuántica (QPE). Trabajos recientes que utilizan conjuntos de bases compactos han intentado minimizar este factor, pero los costos de las compuertas siguen siendo elevados. El desafío radica en reducir la constante de normalización sin introducir aproximaciones físicas adicionales ni sacrificar la ventaja de espacio de la primera cuantización.

Metodología
Los autores proponen emplear el marco de Amplificación Espectral de Suma de Cuadrados (SOSSA, por sus siglas en inglés) para mejorar la complejidad de compuertas asintótica de la estimación de la energía del estado fundamental. El núcleo de la metodología consiste en construir una representación analítica de Suma de Cuadrados (SOS) del Hamiltoniano de la estructura electrónica en primera cuantización en una base de ondas planas.

1. Representación SOS: El Hamiltoniano $H$ se desplaza mediante una constante $\beta$ para formar un operador no negativo $H + \beta I = \sum_\alpha O_\alpha^\dagger O_\alpha$. Los autores identifican un conjunto específico de operadores $O_\alpha$ generados por el operador de densidad de carga total.
   • El término de energía cinética se reproduce mediante los operadores $T_{j,w}$.
   • La energía potencial (electrón-electrón y electrón-núcleo) se reproduce mediante los operadores $\tilde{V}_\nu$, que corresponden a la transformada de Fourier de la densidad de carga total (densidad electrónica menos densidad nuclear).
   • Un término de auto-interacción que surge de elevar al cuadrado $\tilde{V}_\nu$ se convierte en un término de identidad (un desplazamiento constante) utilizando un generador SOS adicional $\tilde{X}_{j,\nu}$.
2. Construcción del Codificado de Bloques: Los autores compilan un codificado de bloques para el operador $H_{SOSSA} = \sum_\alpha O_\alpha$. Esto implica una estructura anidada de oráculos de preparación y selección "externos" e "internos".
   • Para gestionar las diferentes estructuras de sumatoria de los términos SOS (donde los índices de los electrones aparecen en el exterior para la energía cinética pero en el interior para los términos de potencial), los autores introducen un espacio de Hilbert de suma directa con un registro de selección de rama.
   • Utilizan amplificación de amplitud de forma oblicua y qubits de bandera de desbordamiento (overflow-flag) específicos para manejar eficientemente las condiciones de contorno (desplazamientos de momento dentro y fuera de los límites), evitando aritmética condicional compleja.
3. Amplificación Espectral: Al utilizar la raíz cuadrada de la representación SOS, el protocolo construye un operador de caminata con autofases $\phi_k = \cos^{-1}((2E_k + \beta)/\lambda_{SOS} - 1)$. El costo de QPE escala con un factor de normalización efectivo $\lambda_{eff} = \sqrt{E_{gap}(\lambda_{SOS} - E_{gap})}$, donde $E_{gap}$ es la brecha entre el estado fundamental y el límite desplazado.

Contribuciones Clave

• Construcción Analítica de SOS: El artículo proporciona la primera representación SOS analítica para un Hamiltoniano de estructura electrónica en primera cuantización en una base de ondas planas que se implementa con costos comparables a los métodos previos de Combinación Lineal de Unidades (LCU).
• Mejora Asintótica: Los autores derivan un nuevo escalado para la normalización del codificado de bloques. La normalización efectiva $\lambda_{eff}$ escala como $O(\eta \Delta^{-1.5} + \eta^{1.5} \Delta^{-1})$, donde $\Delta$ es el espaciado de la rejilla. Esto mejora el escalado LCU del estado del arte previo de $O(\eta \Delta^{-2} + \eta^2 \Delta^{-1})$.
• Compilación de Recursos: El artículo detalla la compilación del codificado de bloques SOSSA, incluyendo optimizaciones específicas para la preparación de estados, oráculos de selección y reflexiones, demostrando que la sobrecarga del marco SOS es mínima en comparación con el método LCU.

Resultados
Los autores proporcionan estimaciones detalladas de recursos lógicos para varios sistemas de referencia, incluyendo sistemas moleculares (carbonato de etileno, $\text{LiPF}_6$), superceldas de metal de litio y objetivos de gas de electrones uniformes (UEG).

• Magnitud de la Aceleración: El método SOSSA demuestra una aceleración respecto al método LCU previo que oscila entre 2$\times$ y 44$\times$ en el conteo total de compuertas Toffoli.
• Dependencia del Sistema: La aceleración es más pronunciada para sistemas más grandes y realistas (por ejemplo, el objetivo de deuterio muestra una mejora de 44$\times$).
• Verificación de Escalado: El análisis numérico confirma el escalado asintótico teórico. En el límite continuo (dominado por la energía cinética), la complejidad se reduce de $O(\eta \Delta^{-2})$ a $O(\eta \Delta^{-1.5})$. En el límite de sistema grande (dominado por la interacción de Coulomb), se reduce de $O(\eta^2 \Delta^{-1})$ a $O(\eta^{1.5} \Delta^{-1})$.
• Costos Dominantes: El análisis identifica que el costo dominante en el algoritmo sigue siendo el movimiento de datos (intercambio de registros de electrones a registros de trabajo), consistente con otros algoritmos de primera cuantización. Sin embargo, la reducción en la constante de normalización efectiva se traduce directamente en menos llamadas al operador de caminata, produciendo la aceleración observada.

Significancia
El artículo afirma que este trabajo resuelve la dificultad de los altos costos de compuertas en las simulaciones de primera cuantización sin introducir nuevas aproximaciones físicas. Al aprovechar el operador de densidad de carga total para generar una representación SOS de baja frustración, los autores logran una ventaja asintótica demostrable. Esto hace que los algoritmos de primera cuantización sean una opción más viable para la simulación de materiales ab initio, particularmente para sistemas que requieren bases extensas y convergencia al límite termodinámico. Los autores señalan que, si bien la implementación actual se centra en la aproximación de Born-Oppenheimer, el marco es generalizable a casos no de Born-Oppenheimer y el trabajo futuro lo extenderá para incluir pseudopotenciales, un paso necesario para aplicaciones prácticas.