Spectral functions on a quantum computer through system-environment interaction

Enunciado del Problema
Las funciones espectrales, medidas mediante espectroscopía de fotoemisión resuelta en ángulo (ARPES), son críticas para elucidar la estructura de bandas de los materiales. Teóricamente, estas funciones corresponden a funciones dinámicas de un punto en un estado de equilibrio. Aunque son difíciles de calcular clásicamente, su medición en computadoras cuánticas enfrenta obstáculos significativos. Los enfoques estándar, que implican calcular correlaciones dinámicas $\langle c^\dagger_j(t)c_0(0)\rangle$ y realizar transformadas de Fourier, sufren de una gran sobrecarga de muestreo. Específicamente, debido a que los operadores $c_j + c^\dagger_j$ para diferentes sitios $j$ no conmutan, solo se puede medir una correlación por disparo. Para reconstruir la función espectral para todos los $N$ momentos, esto requiere una sobrecarga de muestreo de $O(N)$. Esto es particularmente problemático para las computadoras cuánticas de trampas de iones, que poseen alta fidelidad de puertas y conectividad pero tiempos de disparo relativamente lentos, convirtiendo la cuenta de disparos de $O(N)$ en un cuello de botella para aplicaciones a escala de utilidad.

Metodología
Los autores proponen un método para medir funciones espectrales $A(k, \omega)$ directamente modelando la interacción física entre el sistema y el entorno tal como ocurre en los experimentos de ARPES. El enfoque involucra:

1. Acoplamiento Sistema-Entorno: El sistema de interés (Hamiltoniano $H_{sys}$) se acopla a un entorno de $N$ sitios (Hamiltoniano $H_{env}$) mediante un término de interacción $H_{int}$. El Hamiltoniano total es $H = H_{sys} + \epsilon H_{int} + \omega H_{env}$.
2. Preparación del Estado y Evolución: El sistema se inicializa en un autoestado $|E\rangle$ de $H_{sys}$, mientras que el entorno se inicializa en un estado de vacío $|0\rangle$ (para $A_+$) o en un estado completamente lleno $|all\rangle$ (para $A_-$). El sistema total evoluciona bajo $H$ durante un tiempo fijo $t$.
3. Medición: En lugar de medir correlaciones dinámicas, el protocolo mide el número de ocupación de momento $n(k)$ en el entorno. Esto requiere realizar una transformada de Fourier fermiónica (FFT) en los qubits del entorno al final de la evolución, seguida de la medición de matrices de Pauli locales $Z$.
4. Resultado Teórico: En el límite de acoplamiento débil ($\epsilon \to 0$) y tiempo largo ($t \to \infty$), el valor esperado $\langle n(k) \rangle$ es proporcional a la función espectral $A(k, \omega)$ (convolucionada con un núcleo que se aproxima a una función delta a medida que $t \to \infty$).
5. Implementación de FFT Fermiónica: Para manejar tamaños de sistema arbitrarios (específicamente $N=27$), los autores desarrollan una descomposición de puertas eficiente para una FFT fermiónica genérica de radix-$n$. Utilizan una estructura recursiva que involucra operaciones de "entrelazado" para hacer adyacentes modos fermiónicos no adyacentes para la transformada, optimizada mediante técnicas de decimación de grafos para arquitecturas con acoplamiento todo-con-todo.

Contribuciones Clave

• Muestreo Eficiente: La contribución principal es una reducción en la sobrecarga de muestreo. Al medir la ocupación de momento del entorno, el método obtiene datos espectrales para todos los $N$ momentos $k$ en un solo disparo, reduciendo el número de disparos requeridos por un factor de $O(N)$ en comparación con los enfoques estándar de correlación dinámica.
• Mejora del Tiempo de Ejecución: Esta reducción de disparos se traduce en una mejora de $O(N)$ en el tiempo de ejecución general, lo cual es crucial para plataformas como las computadoras cuánticas de trampas de iones donde el tiempo de disparo es el factor limitante.
• Robustez ante el Error de Trotter: Los autores argumentan que su enfoque es más robusto ante errores de evolución temporal (Trotterización) que los métodos estándar. Mientras que los métodos estándar pueden producir valores negativos no físicos para las funciones espectrales debido a errores, el método propuesto mide un número de ocupación, que permanece como una cantidad positiva incluso con una evolución imperfecta.
• Evaluación de Hardware: Los autores implementaron una FFT fermiónica de radix-3 en 27 qubits y evaluaron diversas estrategias de compilación (CZ locales, escaleras CX, decimación de grafos) en los sistemas de trampas de iones H2-1 y H2-2 de Quantinuum. Identificaron que la compilación con escalado $\log_2(N)$ proporciona la menor energía e infidelidad.
• Demostración del Algoritmo Completo: El algoritmo completo se demostró en un sistema de 54 qubits (27 sitios del sistema + 27 sitios del entorno) en el hardware H2-2 de Quantinuum, calculando con éxito la función espectral de un Hamiltoniano de fermiones libres unidimensional.

Resultados

• Rendimiento del Hardware: En el hardware H2-2, el método reprodujo con éxito la forma de la banda espectral para un sistema de 27 sitios. Aunque el ruido del hardware causó una atenuación general de la señal (un "piso de ruido" cerca de 0.5), las características físicas, como las posiciones y anchos de banda, permanecieron sin sesgo.
• Análisis de Errores: Las comparaciones con simulaciones exactas mostraron que el método propuesto alcanza la misma precisión con significativamente menos pasos de Trotter que el enfoque estándar de correlación dinámica, particularmente para acoplamientos pequeños $\epsilon$.
• Análisis de Costos: Para el experimento específico de 27 sitios, el método incurrió en una sobrecarga de puertas de dos qubits (debido al entorno y la FFT), pero ahorró un factor de $\sim 27$ en el número de disparos. Los autores estiman que sin este método, el costo total del tiempo de ejecución habría sido aproximadamente 11 veces mayor para la misma precisión.
• Efectos de Acoplamiento Finito: El análisis teórico y las simulaciones de fermiones libres indican que un acoplamiento finito $\epsilon$ introduce un ensanchamiento de bandas y "bandas fantasma". Sin embargo, para $\epsilon t \lesssim \pi$, estos efectos son despreciables y el método aproxima fielmente la función espectral.

Significado
El artículo afirma que este enfoque mejora fundamentalmente la escalabilidad del cálculo de funciones espectrales en computadoras cuánticas. Al intercambiar una sobrecarga de qubits y puertas por una reducción masiva en los requisitos de muestreo, el método aborda un cuello de botella crítico para las computadoras cuánticas de trampas de iones. Los autores postulan que a medida que los tamaños de sistema aumenten hacia la escala de utilidad (cientos o miles de sitios), la reducción de $O(N)$ en disparos se convertirá en un factor decisivo en la viabilidad práctica de simular propiedades de materiales en hardware cuántico. El trabajo también proporciona una implementación concreta y eficiente de transformadas de Fourier fermiónicas para tamaños de sistema que no son potencias de dos, un componente necesario para simular modelos de materiales realistas.