Quantum algorithm for solving differential equations using SLAC derivatives

Resumo Técnico: Algoritmo Quântico para Resolver Equações Diferenciais Usando Derivadas SLAC

Enunciado do Problema
Simular teorias de campo contínuas e sistemas quânticos de muitos corpos em um reticulado requer a discretização de operadores diferenciais. Discretizações locais padrão, como diferenças finitas simétricas, frequentemente falham em preservar estruturas contínuas essenciais. Especificamente, elas introduzem artefatos de reticulado não físicos, como duplicação de férmions (modos zero espúrios) e relações de dispersão incorretas, o que impacta cálculos de quantidades físicas como o calor específico. O teorema de Nielsen–Ninomiya estabelece que nenhum operador de Dirac em reticulado local e invariante por translação pode satisfazer simultaneamente a localidade espacial, o limite contínuo correto, a ausência de duplicação de férmions e a invariância quiral.

Para abordar isso, foi proposta a Ação de Reticulado Degenerada com Derivada Contínua (SLAC). A derivada SLAC é definida no espaço de momento para corresponder exatamente ao momento contínuo $p$ dentro da primeira zona de Brillouin, evitando assim a duplicação de férmions e preservando a hermiticidade. No entanto, essa fidelidade exata no espaço de momento tem o custo de uma não-localidade extrema no espaço de posição, resultando em matrizes densas $N \times N$. Em computadores clássicos, lidar com essas matrizes densas é computacionalmente caro. Em computadores quânticos, o desafio reside em codificar em blocos (block-encoding) eficientemente esses operadores densos e não unitários para utilizar Algoritmos de Sistemas Lineares Quânticos (QLSA) como HHL ou QSVT. Além disso, a resolução dos sistemas lineares resultantes é impedida pelo grande número de condição dos operadores discretizados, que escala linearmente com o tamanho do sistema $N$.

Metodologia
Os autores propõem uma estrutura para implementar eficientemente operadores de derivada SLAC em um computador quântico usando a Combinação Linear de Unitários (LCU) e técnicas de codificação em blocos. A metodologia consiste em quatro componentes principais:

1. Preparação de Estado via Testes de Desigualdade de Caixas Aninhadas:
   O principal gargalo na codificação em blocos de operadores densos é preparar o estado auxiliar com as amplitudes densas requeridas. Os autores adaptam uma técnica de preparação de estado de "teste de desigualdade de caixas aninhadas". Em vez de usar métodos de Grover-Rudolph, que podem sofrer com probabilidades de sucesso exponencialmente decrescentes, eles particionam a base computacional em intervalos diádicos (caixas). Ao usar testes de desigualdade contra um registro de referência, eles preparam estados com amplitudes proporcionais a $1/j$ (para a derivada de primeira ordem) e $1/j^2$ (para o Laplaciano) com uma probabilidade de sucesso constante e profundidade de circuito polilogarítmica.

2. Codificação em Blocos LCU de Operadores SLAC:
   O Laplaciano SLAC e a derivada de primeira ordem são matrizes circulares na base de posição. Os autores constroem codificações em blocos expressando esses operadores como combinações lineares de matrizes de deslocamento cíclico (permutação).

   • Laplaciano SLAC: Construído usando coeficientes derivados do limite de reticulado infinito, truncados para um reticulado finito. A codificação em blocos utiliza as amplitudes preparadas e um subtrator modular para a operação SELECT.
   • Derivada SLAC de Primeira Ordem: Similarmente construída, mas requer o manuseio de fases complexas e sinais alternados. Os autores mostram que o componente $j=0$ está ausente, simplificando ligeiramente a preparação do estado.
   • Manuseio de Invariância de Translação Quebrada: A estrutura é estendida para lidar com operadores com simetria de translação quebrada (por exemplo, na presença de potenciais externos ou cortes dependentes de linha) introduzindo oráculos de predicado (KEEP) que mascaram elementos matriciais específicos, tudo operando inteiramente dentro da base de posição para evitar transformações de base repetidas.

3. Transformada de Wavelet de Shannon Quântica (QSWT):
   Para abordar o problema do número de condição, os autores implementam a Transformada de Wavelet de Shannon Quântica. Diferentemente de wavelets locais, a wavelet de Shannon alcança uma separação de momento perfeita, dividindo o espaço de Hilbert em setores Infravermelho (IR) e Ultravioleta (UV) sem sobreposição. A QSWT é implementada eficientemente usando Transformadas de Fourier Quânticas (QFT), somadores modulares e unitários esparsos. Isso permite que o operador SLAC seja transformado em uma representação multiescala e em blocos diagonais, onde o operador é decomposto em blocos IR e UV independentes em escalas sucessivas.

4. Pré-condicionamento Diagonal:
   Na base de wavelet multiescala, o número de condição do operador SLAC é reduzido aplicando um pré-condicionador diagonal. O pré-condicionador escala cada bloco por $1/2^r$, onde $r$ é o índice de escala. Isso reduz o número de condição do sistema pré-condicionado a uma constante pequena ($\kappa \approx 4$) após projetar o modo do espaço nulo.

Contribuições e Resultados Chave

• Codificação em Blocos Ótima do Laplaciano SLAC: O artigo apresenta uma codificação em blocos $\epsilon$-aproximada para o Laplaciano SLAC de $n$ qubits com um fator de subnormalização $\alpha \approx (\pi^2 + 24)/3 \approx 11,29$ (uma constante independente de $N$). A construção usa $O(n)$ qubits auxiliares e tem complexidade de portas de $O(n^2)$. Esta é uma codificação em blocos ótima.
• Codificação em Blocos Eficiente da Derivada de Primeira Ordem: Uma codificação em blocos $\epsilon$-aproximada para a derivada SLAC de primeira ordem é construída com um fator de subnormalização $\alpha = 2(n-1) = O(n)$. Também requer $O(n)$ auxiliares e $O(n^2)$ portas. Isso é classificado como uma codificação em blocos "boa".
• Representação Multiescala: Os autores demonstram que aplicar a QSWT recursivamente ao operador SLAC codificado em blocos produz uma representação multiescala. Isso requer uma consulta à codificação em blocos e $O(n)$ consultas à rotina QSWT, mantendo uma complexidade total de portas de $O(n^2)$ por consulta.
• Redução do Número de Condição: Ao combinar a representação multiescala com o pré-condicionador de wavelet diagonal e projetar o espaço nulo, o número de condição do sistema é reduzido a uma constante $\kappa = O(1)$.
• Resolução de EDPs: O artigo estabelece que equações diferenciais parciais (EDPs) $d$-dimensionais discretizadas via o formalismo SLAC podem ser resolvidas usando QLSA. O estado solução pode ser preparado com $O(\alpha(k) \log(1/\epsilon))$ consultas à codificação em blocos, onde $\alpha(k)$ é o fator de subnormalização ($O(1)$ para Laplaciano, $O(n)$ para primeira ordem). A complexidade total de portas por consulta é $O(dn^2)$.

Significado e Alegações
O artigo afirma fornecer uma estrutura prática para simular sistemas físicos onde preservar propriedades de dispersão contínuas é crítico, como em modelos de matéria condensada (por exemplo, líquidos de Luttinger, isolantes topológicos) e teorias de gauge em reticulado onde a duplicação de férmions deve ser evitada.

O significado reside em demonstrar que operadores densos e não locais — anteriormente considerados difíceis de lidar eficientemente em computadores quânticos devido à sua densidade matricial — podem ser codificados em blocos de forma eficiente. Os autores enfatizam que sua abordagem evita o sobrecusto exponencial de amplitude frequentemente associado a codificações em blocos genéricas de operadores pseudo-diferenciais, explorando a estrutura analítica específica das derivadas SLAC e usando preparação de estado baseada em testes de desigualdade.

Além disso, o trabalho preenche a lacuna entre os benefícios teóricos das derivadas SLAC (física contínua exata) e os requisitos práticos dos algoritmos quânticos (codificação em blocos eficiente e baixos números de condição). Ao integrar a QSWT e o pré-condicionamento diagonal, os autores mostram que a má condição inerente a operadores diferenciais discretizados pode ser mitigada para uma constante, preservando a aceleração exponencial prometida por solucionadores lineares quânticos. As construções são apresentadas como amigáveis à tolerância a falhas, com contagens explícitas de Toffoli fornecidas, tornando-as adequadas para estimativa de recursos em futuras arquiteturas quânticas.