Spectral Gaps via Imaginary Time

Resumo Técnico: Lacunas Espectrais via Tempo Imaginário

Definição do Problema
A lacuna espectral, definida como $\Delta E = E_1 - E_0$ (a diferença de energia entre o estado fundamental e o primeiro estado excitado), é um parâmetro crítico em diversos contextos físicos, variando desde comprimentos de onda de fluorescência e transferência de energia excitônica em materiais até tempos de execução adiabática e fidelidades de transferência de estado em informação quântica. Determinar essa lacuna é computacionalmente desafiador; no pior caso, decidir se um Hamiltoniano local é lacunado é indecidível, e estimar o espectro de baixas energias é geralmente NP-difícil. Métodos numéricos existentes, como o Monte Carlo de Fermiôni de Tempo Quântico (que sofre do problema do sinal) e métodos de Redes de Tensores (limitados pela dimensão de ligação), enfrentam compromissos significativos entre tamanho do sistema, tamanho da base e sobrecarga de medição. Embora a Evolução no Tempo Imaginário (ITE) seja uma ferramenta clássica robusta para encontrar estados fundamentais, extrair lacunas espectrais específicas sem resolver a densidade total de estados ou realizar a estimação de fase quântica permanece um desafio específico.

Metodologia
Os autores propõem um método para estimar a lacuna espectral $\Delta E$ diretamente a partir da razão entre dois valores esperados de comutadores aninhados avaliados ao longo de uma trajetória de tempo imaginário.

1. Estrutura Teórica:
   Para um Hamiltoniano $H$ de dimensão finita com decomposição espectral $H = \sum E_n \Pi_n$, e um estado inicial $|\phi_0\rangle$ com suporte nos subespaços do estado fundamental ($\Pi_0$) e do primeiro estado excitado ($\Pi_1$), o sistema evolui no tempo imaginário $\tau$ como $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$.
   O método utiliza um observável "coordenador" $O$ que acopla os subespaços do estado fundamental e do primeiro estado excitado. Os autores definem o $M$-ésimo comutador aninhado $[H, O]_M$.
   O resultado central (Teorema 1) demonstra que conforme $\tau \to \infty$, a razão dos valores esperados converge para o quadrado da lacuna de energia:
   $$\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 é agnóstico a degenerações em níveis excitados superiores porque se baseia em formulações de projetores, em vez de autovetores individuais.

2. Estimadores Logarítmicos:
   O artigo também deriva estimadores logarítmicos para extrair $E_0 + E_1$ (a partir da inclinação de $\ln|\langle [H, O]_M \rangle_\tau|$ vs. $\tau$) e a segunda lacuna $E_2 - E_1$ (a partir do resíduo do ajuste da razão), fornecendo alternativas quando a divisão numérica de números pequenos é instável.

3. Implementação Quântica (ITQDE):
   Para implementar isso em hardware quântico, onde a evolução unitária é nativa mas o tempo imaginário não é, os autores utilizam a Emulação Dinâmica Quântica de Tempo Imaginário (ITQDE). Eles empregam uma transformação de Hubbard-Stratonovich para expressar o núcleo euclidiano $e^{-\tau H^2}$ como uma média Gaussiana de unitárias de tempo real. Este núcleo é discretizado usando quadratura de Gauss-Hermite, aproximando a evolução como uma soma ponderada de unitárias.
   Os valores esperados são então reduzidos a somas ponderadas de partes reais de sobreposições da forma $\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$, onde $P_\alpha$ são strings de Pauli. Essas sobreposições podem ser medidas usando testes de Hadamard padrão. Os autores observam que, devido à simetria do núcleo, apenas as partes reais dessas sobreposições são necessárias.

Resultados Principais
Os autores testaram o método usando simulações clássicas (via biblioteca QuSpin) em dois modelos padrão:

• Modelo de Ising de Campo Transverso (TFIM) 1D: Para tamanhos de sistema $L=4$, o erro relativo na lacuna estimada diminuiu exponencialmente com $\tau$ (atingindo $\sim 10^{-8}$) antes de aumentar novamente em $\tau$ muito grande devido à supressão do primeiro estado excitado abaixo do ruído numérico.
• Modelo de Fermi-Hubbard 1D: Um decaimento exponencial de erro semelhante foi observado para uma cadeia com meio preenchimento ($L=4$), com erros relativos atingindo $\sim 10^{-4}$. O método lidou com sucesso com o ambiente fermiônico interagente com estruturas de spin e carga não triviais.
• Simulação Quântica: Uma simulação de prova de conceito foi conduzida no TFIM usando o backend Aer do Qiskit com $10^6$ shots. Embora o método fosse compatível com o framework ITQDE, o erro mínimo alcançado foi de $\sim 10^{-1}$. Os autores atribuem isso ao "compromisso de largura de banda": o filtro Gaussiano necessário para suprimir estados superiores também leva as sobreposições relevantes para o nível do ruído de disparo (shot-noise floor) se o tempo imaginário $\tau$ e a ordem da quadratura $K$ não forem cuidadosamente ajustados.

Significância e Alegações
O artigo afirma oferecer uma "vantagem estrutural modesta, mas genuína" para a estimação espectral.

• Simplicidade e Escopo: O método oferece uma maneira de estimar lacunas de baixas energias sem resolver autovetores individuais, realizar estimação de fase quântica ou reconstruir a densidade total de estados. Ele converte questões espectrais em questões sobre relaxação ao longo de um fluxo de um único parâmetro.
• Robustez: O estimador é robusto contra degenerações em níveis de energia superiores e aplica-se tanto a álgebra linear esparsa clássica quanto a Hamiltonianos de soma de termos quânticos.
• Modéstia: Os autores afirmam explicitamente que esta abordagem não contorna a complexidade de pior caso do problema do Hamiltoniano local. Não é um substituto para a reconstrução espectral completa ou estimação de fase quando uma densidade de estados detalhada é necessária.
• Perspectiva Futura: O trabalho serve como uma prova de conceito. Os autores observam que tornar o método competitivo em hardware quântico requer importar a análise completa de largura de banda de trabalhos anteriores para otimizar parâmetros ($\tau, K$), projetar melhores observáveis coordenadores e empregar estimadores de sobreposição mais eficientes em termos de medição. O artigo posiciona essa perspectiva focada em lacunas como uma ferramenta que se integra a programas mais amplos que utilizam Hamiltonianos efetivos para resolver equações diferenciais, onde taxas de relaxação e lacunas são frequentemente as quantidades de interesse primário.