Spectral Gaps via Imaginary Time

Resumen Técnico: Brechas Espectrales mediante Tiempo Imaginario

Planteamiento del Problema
La brecha espectral, definida como $\Delta E = E_1 - E_0$ (la diferencia de energía entre el estado fundamental y el primer estado excitado), es un parámetro crítico en diversos contextos físicos, que van desde las longitudes de onda de fluorescencia y la transferencia de energía excitónica en materiales, hasta los tiempos de ejecución adiabática y las fidelidades de transferencia de estado en la información cuántica. Determinar esta brecha es computacionalmente desafiante; en el peor de los casos, decidir si un Hamiltoniano local tiene una brecha es indecidible, y estimar los espectros de baja energía es generalmente NP-duro. Los métodos numéricos existentes, como el Monte Carlo de Cuántica de Fermiones (que sufre el problema del signo) y los métodos de Redes de Tensores (limitados por la dimensión de enlace), enfrentan compromisos significativos entre el tamaño del sistema, el tamaño de la base y la sobrecarga de medición. Si bien la Evolución en Tiempo Imaginario (ITE) es una herramienta clásica robusta para encontrar estados fundamentales, extraer brechas espectrales específicas sin resolver la densidad completa de estados o realizar la estimación de fase cuántica sigue siendo un desafío específico.

Metodología
Los autores proponen un método para estimar la brecha espectral $\Delta E$ directamente a partir de la razón de dos valores de expectativa de conmutadores anidados evaluados a lo largo de una trayectoria de tiempo imaginario.

1. Marco Teórico:
   Para un Hamiltoniano de dimensión finita $H$ con descomposición espectral $H = \sum E_n \Pi_n$, y un estado inicial $|\phi_0\rangle$ con soporte en los subespacios del estado fundamental ($\Pi_0$) y del primer estado excitado ($\Pi_1$), el sistema evoluciona en el tiempo imaginario $\tau$ como $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$.
   El método utiliza un observable "coordinador" $O$ que acopla los subespacios del estado fundamental y del primer estado excitado. Los autores definen el $M$-ésimo conmutador anidado $[H, O]_M$.
   El resultado central (Teorema 1) demuestra que cuando $\tau \to \infty$, la razón de los valores de expectativa converge al cuadrado de la brecha de energía:
   $$\frac{\langle [H, O]_{M+2} \rangle_\tau}{\langle [H, O]_M \rangle_\tau} = (E_0 - E_1)^2 + \mathcal{O}(e^{-\tau(E_2 - E_1)})$$
   Este estimador es agnóstico a las degeneraciones en los niveles excitados superiores porque se basa en formulaciones de proyectores en lugar de autoestados individuales.

2. Estimadores Logarítmicos:
   El artículo también deriva estimadores logarítmicos para extraer $E_0 + E_1$ (a partir de la pendiente de $\ln|\langle [H, O]_M \rangle_\tau|$ frente a $\tau$) y la segunda brecha $E_2 - E_1$ (a partir del residuo del ajuste de la razón), proporcionando alternativas cuando la división numérica de números pequeños es inestable.

3. Implementación Cuántica (ITQDE):
   Para implementar esto en hardware cuántico, donde la evolución unitaria es nativa pero el tiempo imaginario no lo es, los autores utilizan la Emulación Dinámica de Tiempo Cuántico de Tiempo Imaginario (ITQDE). Utilizan una transformación de Hubbard-Stratonovich para expresar el núcleo euclidiano $e^{-\tau H^2}$ como un promedio gaussiano de unitarias en tiempo real. Este núcleo se discretiza mediante cuadratura de Gauss-Hermite, aproximando la evolución como una suma ponderada de unitarias.
   Los valores de expectativa se reducen entonces a sumas ponderadas de las partes reales de los solapamientos de la forma $\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$, donde $P_\alpha$ son cadenas de Pauli. Estos solapamientos pueden medirse utilizando pruebas de Hadamard estándar. Los autores señalan que, debido a la simetría del núcleo, solo se requieren las partes reales de estos solapamientos.

Resultados Clave
Los autores evaluaron el método utilizando simulaciones clásicas (vía la biblioteca QuSpin) en dos modelos estándar:

• Modelo de Ising de Campo Transverso (TFIM) 1D: Para tamaños de sistema $L=4$, el error relativo en la brecha estimada disminuyó exponencialmente con $\tau$ (alcanzando $\sim 10^{-8}$) antes de aumentar nuevamente a un $\tau$ muy grande debido a la supresión del primer estado excitado por debajo del ruido numérico.
• Modelo de Fermi-Hubbard 1D: Se observó un decaimiento exponencial similar en el error para una cadena de llenado medio ($L=4$), con errores relativos alcanzando $\sim 10^{-4}$. El método manejó con éxito el entorno de fermiones interactuantes con estructuras de espín y carga no triviales.
• Simulación Cuántica: Se realizó una simulación de prueba de principio utilizando el backend Aer de Qiskit para el TFIM con $10^6$ disparos (shots). Aunque el método era compatible con el marco ITQDE, el error mínimo alcanzado fue $\sim 10^{-1}$. Los autores atribuyen esto al "compromiso de ancho de banda" (bandwidth tradeoff): el filtro gaussiano requerido para suprimir estados superiores también impulsa los solapamientos relevantes hacia el umbral del ruido de disparo (shot-noise floor) si el tiempo imaginario $\tau$ y el orden de la cuadratura $K$ no se ajustan cuidadosamente.

Significancia y Reivindicaciones
El artículo afirma ofrecer una "ventaja estructural modesta pero genuina" para la estimación espectral.

• Simplicidad y Alcance: El método proporciona una forma de estimar brechas bajas sin resolver autoestados individuales, realizar la estimación de fase cuántica o reconstruir la densidad completa de estados. Convierte las preguntas espectrales en preguntas sobre la relajación a lo largo de un flujo de un solo parámetro.
• Robustez: El estimador es robusto contra las degeneraciones en los niveles de energía superiores y se aplica tanto a álgebra lineal dispersa clásica como a Hamiltonianos de suma de términos en sistemas cuánticos.
• Modestia: Los autores declaran explícitamente que este enfoque no evita la complejidad del peor de los casos del problema del Hamiltoniano local. No es un reemplazo para la reconstrucción espectral completa o la estimación de fase cuando se requiere la densidad detallada de estados.
• Perspectiva Futura: El trabajo sirve como una prueba de principio. Los autores señalan que hacer que el método sea competitivo en hardware cuántico requiere importar el análisis completo de ancho de banda de trabajos previos para optimizar los parámetros ($\tau, K$), diseñar mejores observables coordinadores y emplear estimadores de solapamiento más eficientes en mediciones. El artículo posiciona esta perspectiva centrada en la brecha como una herramienta que se integra con programas más amplios que utilizan Hamiltonianos efectivos para resolver ecuaciones diferenciales, donde las tasas de relajación y las brechas son a menudo las cantidades de interés primario.