Fault tolerant computation of the static structure factor and finite size effects

Resumo Técnico: Computação Tolerante a Falhas do Fator de Estrutura Estático e Efeitos de Tamanho Finito

Enunciado do Problema
Algoritmos quânticos tolerantes a falhas oferecem um caminho para estimar energias de estado fundamental de materiais periódicos além do alcance dos métodos clássicos de estrutura eletrônica. No entanto, um desafio significativo permanece: simulações quânticas tipicamente avaliam supercélulas finitas ou malhas discretas da zona de Brillouin (BZ), enquanto as propriedades dos materiais são definidas no limite termodinâmico. Essa discrepância introduz erros de tamanho finito com dois componentes primários:

1. Erros de um corpo (shell): Resultantes da amostragem discreta de estados de partícula única, particularmente severos em metais próximos à superfície de Fermi.
2. Erros de dois corpos: Resultantes da distorção das correlações de Coulomb de longo comprimento de onda em uma célula periódica finita.

Embora os erros de um corpo possam ser mitigados via média de torção (twist averaging), o erro residual de dois corpos é controlado por flutuações de densidade de longo comprimento de onda. Estratégias de mitigação tradicionais, como o "subamostragem" (down-sampling) — que consiste em repetir cálculos de energia em supercélulas progressivamente maiores ou malhas $k$ mais densas — são proibitivamente caras no regime tolerante a falhas, pois a estimativa da energia do estado fundamental já é a primitiva quântica dominante. O artigo aborda a necessidade de uma estratégia de pós-processamento mais direcionada e quântica-eficiente para corrigir esses erros de tamanho finito sem repetir a rotina completa de estimativa de energia.

Metodologia
Os autores propõem uma estratégia híbrida que combina a média de torção com a computação quântica tolerante a falhas do Fator de Estrutura Estático (SSF), $S(q)$, em baixos momentos.

1. Média de Torção (Twist Averaging): Utilizada para suprimir efeitos de shell de um corpo. Isso envolve a média das energias sobre um conjunto de condições de contorno torcidas ($\theta$) em uma célula de simulação fixa, evitando a necessidade de aumentar o tamanho do sistema.
2. Reconstrução de Infravermelho via SSF: O erro de dois corpos residual dominante é determinado pelo comportamento de baixo momento (infravermelho) do SSF. Os autores formulam o operador SSF $\hat{S}_q = \hat{\rho}^\dagger_q \hat{\rho}_q$ em uma base de orbitais de Bloch, preservando a estrutura de momento cristalino.
3. Codificação de Bloco e Estimativa:
   • O operador de densidade $\hat{\rho}_q$ é decomposto em uma Combinação Linear de Unitárias (LCU) envolvendo fatores de forma e rotações de Givens.
   • Uma codificação de bloco de $\hat{S}_q$ é construída via o produto refletido da codificação de bloco do operador de densidade: $U_{\hat{S}} = [U_{\hat{\rho}}]^\dagger \cdot \text{REF} \cdot U_{\hat{\rho}}$.
   • O valor esperado do estado fundamental $\langle \psi_0 | \hat{S}_q | \psi_0 \rangle$ é estimado usando um teste de Hadamard amplificado, combinando amplificação de amplitude com estimativa de fase para alcançar o escalonamento limitado por Heisenberg.
4. Seleção Adaptativa de Janela de Infravermelho: Para reconstruir a correção de tamanho finito, o algoritmo deve identificar o corte de momento $q_t$ onde $S(q)$ entra em seu regime assintótico de infravermelho ($S(q) \sim \alpha q^\beta$). O artigo introduz dois procedimentos de busca adaptativa:
   • Busca Adaptativa Global: Uma busca binária sobre todas as cascas (shells) de baixo-$q$ amostradas para encontrar o maior corte estável.
   • Busca de Janela Local: Uma abordagem de janela deslizante que reutiliza dados de cascas previamente computados para reduzir o número de novas avaliações quânticas, trocando armazenamento clássico por eficiência de consulta quântica.

Principais Contribuições

• Formulação do Operador: Os autores definem o operador SSF em uma base de orbitais de Bloch e constroem sua codificação de bloco. Esta representação separa as regras de seleção de momento da física intracelular, comprimindo informações específicas do material em fatores de forma pré-computados.
• Algoritmo Tolerante a Falhas: Eles fornecem uma arquitetura de circuito completa para estimar $S(q)$ usando um teste de Hadamard amplificado, incluindo estimativas exatas de recursos para portas Toffoli e qubits lógicos.
• Protocolos Adaptativos: A introdução de procedimentos de busca binária global e local permite a determinação eficiente da janela de ajuste de infravermelho sem conhecimento prévio do regime assintótico.
• Análise de Recursos: O artigo demonstra que o custo desta correção é subdominante em relação à rotina principal de estimativa de energia do estado fundamental.

Resultados e Escalonamento
A análise de recursos revela os seguintes comportamentos de escalonamento, onde $N_b$ é o número de orbitais de Bloch e $N_k$ é o número de pontos $k$:

• Qubits Lógicos: A correção requer $O(N_b N_k)$ qubits lógicos, o que é subdominante comparado aos $O(N_{pw}^{1/2} N_b N_k)$ necessários para a simulação completa do Hamiltoniano (onde $N_{pw}$ é o tamanho da base de ondas planas).
• Contagem de Toffoli: O escalonamento de Toffoli dominante para a correção é $O((N_b N_k)^3)$. Embora compartilhe a mesma dependência cúbica do tamanho da base como a estimativa de energia, ela evita o grande preefator de ondas planas ($N_{pw}^{3/2}$) inerente à simulação completa do Hamiltoniano.
• Requisitos de Precisão: A rotina de correção é significativamente mais tolerante ao ruído do que a própria estimativa de energia. O ruído permitido em $S(q)$ cresce com o tamanho do sistema ($p \leq p N_k^{1/3}$), o que significa que $S(q)$ não precisa ser computado com a mesma precisão absoluta que a energia total.
• Comparação com Subamostragem: Diferente da subamostragem, que escala como $(2(N_a N_k)^{1/3} + 1)(N_b N_k)^3 N_{pw}^{3/2}$ (exigindo repetidas simulações completas em células maiores), o método proposto adiciona um termo de correção que é assintoticamente mais barato e evita o crescimento exponencial no tamanho do sistema.

Significância
O artigo afirma que a mitigação de tamanho finito em simulações quânticas periódicas não precisa ser alcançada forçando o cálculo principal de muitos corpos para malhas $k$ excessivamente densas ou supercélulas maiores. Em vez disso, ao tratar o erro de tamanho finito como um problema de correlação de densidade de infravermelho, pode-se extrair a correção necessária de um pequeno conjunto de medições de SSF de baixo momento.

Esta abordagem é particularmente significativa para materiais metálicos e de condições extremas, onde é necessária uma resolução densa de $k$-espaço para controlar os efeitos de shell e o blindamento de longo alcance. Ao combinar a média de torção (para efeitos de um corpo) com uma avaliação quântica direcionada do SSF (para efeitos de dois corpos), o método fornece uma rota para correções do limite termodinâmico que é tanto fisicamente bem fundamentada quanto assintoticamente compatível com as restrições de recursos tolerantes a falhas. O trabalho estabelece que o SSF é um observável alvo útil no cenário tolerante a falhas porque captura a física dominante de erro de tamanho finito residual, permanecendo computacionalmente subdominante à tarefa primária de estimativa de energia.