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

Resumo Técnico: Amplificação Espectral para Estimativa de Energia do Estado Fundamental em Primeira Quantização

Enunciado do Problema
Simular a estrutura eletrônica de sistemas químicos e de materiais é um alvo primário para algoritmos quânticos. Embora a maior parte do trabalho se concentre em representações de segunda quantização, abordagens de primeira quantização oferecem uma vantagem exponencial de espaço, escalando a memória como $O(\eta \log N)$ (onde $\eta$ é o número de elétrons e $N$ é o número de funções de base), comparado a $O(N)$ para a segunda quantização. Isso torna a primeira quantização particularmente adequada para visar o limite de conjuntos de bases grandes e simulações de materiais que requerem convergência para o limite termodinâmico.

No entanto, apesar da escala de memória favorável, a complexidade de portas para a estimativa da energia do estado fundamental em primeira quantização permanece proibitivamente alta quando se utilizam conjuntos de bases de ondas planas (plane-wave) extensos. O principal gargalo é a constante de normalização do codificação de bloco ($\lambda$) necessária para a Estimativa de Fase Quântica (QPE). Trabalhos recentes utilizando conjuntos de bases compactos tentaram minimizar este fator, mas os custos de portas permanecem elevados. O desafio reside em reduzir a constante de normalização sem introduzir aproximações físicas adicionais ou sacrificar a vantagem de espaço da primeira quantização.

Metodologia
Os autores propõem o emprego do framework de Amplificação Espectral de Soma de Quadrados (SOSSA) para melhorar a complexidade assintótica de portas na estimativa da energia do estado fundamental. O núcleo da metodologia envolve a construção de uma representação analítica de Soma de Quadrados (SOS) do Hamiltoniano de estrutura eletrônica de primeira quantização em uma base de ondas planas.

1. Representação SOS: O Hamiltoniano $H$ é deslocado por uma constante $\beta$ para formar um operador não negativo $H + \beta I = \sum_\alpha O_\alpha^\dagger O_\alpha$. Os autores identificam um conjunto específico de operadores $O_\alpha$ gerados pelo operador de densidade de carga total.
   • O termo de energia cinética é reproduzido pelos operadores $T_{j,w}$.
   • A energia potencial (elétron-elétron e elétron-núcleo) é reproduzida pelos operadores $\tilde{V}_\nu$, que correspondem à transformada de Fourier da densidade de carga total (densidade eletrônica menos densidade nuclear).
   • Um termo de autointeração decorrente do quadrado de $\tilde{V}_\nu$ é convertido em um termo de identidade (um deslocamento constante) usando um gerador SOS adicional $\tilde{X}_{j,\nu}$.
2. Construção da Codificação de Bloco: Os autores compilam uma codificação de bloco para o operador $H_{SOSSA} = \sum_\alpha O_\alpha$. Isso envolve uma estrutura aninhada de oráculos de preparação e seleção "externos" e "internos".
   • Para gerenciar as diferentes estruturas de somatório dos termos SOS (onde os índices eletrônicos aparecem no exterior para a energia cinética, mas no interior para os termos potenciais), os autores introduzem um espaço de Hilbert de soma direta com um registro de seleção de ramo.
   • Eles utilizam amplificação de amplitude oblíqua e qubits de flag de overflow específicos para lidar com condições de contorno (deslocamentos de momento dentro vs. fora dos limites) de forma eficiente, evitando aritmética condicional complexa.
3. Amplificação Espectral: Ao utilizar a raiz quadrada da representação SOS, o protocolo constró-se um operador de caminhada com autofases $\phi_k = \cos^{-1}((2E_k + \beta)/\lambda_{SOS} - 1)$. O custo da QPE escala com um fator de normalização efetivo $\lambda_{eff} = \sqrt{E_{gap}(\lambda_{SOS} - E_{gap})}$, onde $E_{gap}$ é o gap entre o estado fundamental e o limite deslocado.

Principais Contribuições

• Construção Analítica SOS: O artigo fornece a primeira representação SOS analítica para um Hamiltoniano de estrutura eletrônica de primeira quantização em uma base de ondas planas que é implementada com custos comparáveis aos métodos anteriores de Combinação Linear de Unitárias (LCU).
• Melhoria Assintótica: Os autores derivam uma nova escala para a normalização da codificação de bloco. A normalização efetiva $\lambda_{eff}$ escala como $O(\eta \Delta^{-1.5} + \eta^{1.5} \Delta^{-1})$, onde $\Delta$ é o espaçamento da grade. Isso melhora a escala LCU do estado da arte anterior de $O(\eta \Delta^{-2} + \eta^2 \Delta^{-1})$.
• Compilação de Recursos: O artigo detalha a compilação da codificação de bloco SOSSA, incluindo otimizações específicas para oráculos de preparação de estado, seleção e reflexão, demonstrando que o overhead do framework SOS é mínimo comparado ao método LCU.

Resultos
Os autores fornecem estimativas detalhadas de recursos lógicos para vários sistemas de benchmark, incluindo sistemas moleculares (etileno carbonato, $\text{LiPF}_6$), supercélulas de lítio metálico e alvos de gás eletrônico uniforme (UEG).

• Magnitude do Speedup: O método SOSSA demonstra um speedup sobre o método LCU anterior variando de 2$\times$ a 44$\times$ nas contagens totais de portas Toffoli.
• Dependência do Sistema: O speedup é mais pronunciado para sistemas maiores e mais realistas (por exemplo, o alvo de deutério mostra uma melhoria de 44$\times$).
• Verificação de Escala: A análise numérica confirma a escala assintótica teórica. No limite contínuo (dominado pela energia cinética), a complexidade reduz de $O(\eta \Delta^{-2})$ para $O(\eta \Delta^{-1.5})$. No limite de sistema grande (dominado pela interação Coulomb), reduz de $O(\eta^2 \Delta^{-1})$ para $O(\eta^{1.5} \Delta^{-1})$.
• Custos Dominantes: A análise identifica que o custo dominante no algoritmo permanece o movimento de dados (troca de registros de elétrons para registros de workspace), consistente com outros algoritmos de primeira quantização. No entanto, a redução na constante de normalização efetiva traduz-se diretamente em menos chamadas ao operador de caminhada, resultando no speedup observado.

Significância
O artigo afirma que este trabalho resolve a dificuldade dos altos custos de portas em simulações de primeira quantização sem introduzir novas aproximações físicas. Ao alavancar o operador de densidade de carga total para gerar uma representação SOS de baixa frustração, os autores alcançam uma vantagem assintótica provável. Isso torna os algoritmos de primeira quantização uma escolha mais viável para simulações de materiais ab initio, particularmente para sistemas que exigem grandes conjuntos de bases e convergência para o limite termodinâmico. Os autores observam que, embora a implementação atual foque na aproximação de Born-Oppenheimer, o framework é generalizável para casos não-Born-Oppenheimer e trabalhos futuros o estenderão para incluir pseudopotenciais, um passo necessário para aplicações práticas.