A Quantum Spectral Solver for Periodic Incompressible Stokes Flow

Resumo Técnico: Um Solucionador Espectral Quântico para Fluxo de Stokes Incompressível Periódico

Enunciado do Problema
O artigo aborda o desafio de resolver as equações de Stokes incompressíveis, estacionárias e bidimensionais em um domínio periódico utilizando computação quântica. Embora os métodos espectrais sejam altamente eficazis para problemas periódicos — transformando operadores diferenciais em multiplicadores algébricos diagonais no espaço de Fourier — representar a restrição de incompressibilidade (condição de divergência nula) e o acoplamento pressão-velocidade associado dentro de uma estrutura unitária quântica permanece um obstáculo significativo. Os autores visam desenvolver um solucionador espectral quântico que trate explicitamente a natureza vetorial do campo de velocidade, imponha a restrição de divergência nula e recupere tanto os campos de velocidade e pressão quanto observáveis homogeneizados, sem reconstruir o campo completo.

Metodologia
O método proposto utiliza a Transformada de Fourier Quântica (QFT) como uma mudança de base coerente para transitar entre o espaço físico e o espaço espectral. O núcleo do algoritmo baseia-se em uma fatoração espectral do operador de Stokes que separa o problema em operações quânticas gerenciáveis:

1. Fatoração Espectral: No espaço de Fourier, as equações de Stokes se desacoplam modo a modo. A solução envolve uma projeção de Helmholtz para impor a incompressibilidade (removendo a componente longitudinal da força) seguida por um escalonamento escalar de Laplaciano inverso.
2. Projeção Baseada em Rotação: Para o caso 2D, a projeção de Helmholtz é realizada por uma rotação dependente do modo que transforma as componentes cartesianas da velocidade $(u_0, u_1)$ em coordenadas longitudinais ($\parallel$) e transversais ($\perp$) relativas ao vetor de onda $k$. A restrição de incompressibilidade é então imposta zerando a componente longitudinal e escalonando a componente transversal pelo fator do Laplaciano inverso ($1/(\mu|k|^2)$).
3. Arquitetura de Circuito Quântico:
   • Registradores: O circuito utiliza registradores de modo ($k_0, k_1$) para modos de Fourier, um registrador de componentes ($c$) para direções cartesianas, um registrador de alvo ($t$) para codificação de amplitude e um registrador de trabalho ($w$) para aritmética.
   • Preparação de Estado: O campo de força de entrada é codificado em um estado quântico normalizado.
   • QFT: Uma QFT 2D transforma o estado para o espaço de Fourier.
   • Rotação Angular ($U_R$): Uma rotação dependente do modo $R(k)$ é aplicada ao registrador de componentes para alinhar com a base longitudinal-transversal. O ângulo de rotação é aproximado via polinômios de partes contínuas (piecewise polynomials) codificados através de rotações controladas.
   • Escalonamento Escalar ($U_\Gamma, U_\Lambda$): A componente transversal é escalonada pela função de Green $\Gamma(k) = 1/(\mu|k|^2)$, e a componente longitudinal é escalonada por $\Lambda(k) = 1/|k|$ (para recuperação da pressão). Esses escalonamentos são implementados via codificação de amplitude polinomial no registrador alvo.
   • Pós-seleção e Reconstrução: A solução desejada é recuperada pós-selecionando ramos específicos dos registradores ancilares (ex: $c=1$ para velocidade transversal, $c=0$ para pressão). A QFT inversa reconstrói o campo físico.
4. Extração de Observáveis: Os autores demonstram que quantidades homogeneizadas, como a energia cinética, podem ser estimadas diretamente do estado espectral antes da QFT inversa, evitando o custo de uma reconstrução total do campo.

Principais Contribuições

1. Primeiro Solucionador de Stokes Quântico Baseado em QFT: Os autores apresentam o que afirmam ser o primeiro solucionador espectral quântico para as equações de Stokes incompressíveis, estacionárias e periódicas em 2D.
2. Fatoração Espectral Explícita: Eles derivam uma fatoração compatível com sistemas quânticos que separa a projeção de Helmholtz, o escalonamento do Laplaciano inverso e a recuperação da pressão, utilizando uma representação baseada em rotação para o projetor 2D.
3. Design de Circuito Modular: O artigo fornece uma arquitetura de circuito explícita apresentando rotações dependentes do modo e blocos espectrais polinomiais condicionados às componentes. Os parâmetros de aproximação (grau polinomial e paginação/tiling) são tornados explícitos.
4. Recuperação de Observável Homogeneizado: O método é validado em um benchmark inspirado em RVE, mostrando que a energia cinética homogeneizada pode ser recuperada sem reconstruir o campo de velocidade completo.

Resultos Numéricos
O método foi validado usando três benchmarks em um simulador de estado de vetor sem ruído:

• Vórtice de Taylor-Green: Uma solução suave manufaturada usada para verificar a consistência dos circuitos montados (convenções de Fourier, normalização e projeção).
• Dipolo de Força Regularizado: Um benchmark periódico com força localizada e amplo suporte de Fourier. Este teste distinguiu entre erros de truncamento de resolução finita e erros introduzidos pela aproximação polinomial dos operadores espectrais. Os resultados mostraram que, embora o solver discreto exato convirja com o refinamento da grade, a aproximação polinomial introduz um platô de erro para modos de alta frequência se a paginação (tiling) for fixa.
• Flutuação Inspirada em RVE: Um benchmark usando um espectro de flutuação do tipo Kolmogorov ($k^{-5/3}$). Os resultados demonstraram que o observável de energia cinética homogeneizada converge para a solução de referência conforme a grade é refinada, desde que a aproximação polinomial seja ajustada ao espectro ponderado pela energia. O estudo destacou que a contagem de portas é dominada pela codificação polinomial da rotação e dos operadores da função de Green, não pela própria QFT.

Significância e Alegações
O artigo posiciona este trabalho não como um substituto direto para DNS clássico ou um modelo de turbulência, mas como um bloco espectral de Stokes para mecânica computacional quântica. Sua principal significância reside em:

• Compatibilidade Multiescala: O solver é projetado para ser incorporado em arquiteturas maiores de elemento de volume representativo (RVE) clássico-quântico, onde o paralelismo quântico pode atualizar simultaneamente todos os estados de RVE associados a pontos de integração macroscópicos.
• Eficiência na Estimativa de Observáveis: Demonstra a capacidade de extrair quantidades médias (como energia cinética) diretamente do estado espectral, o que é crucial para fluxos multiescala onde a reconstrução total do campo é desnecessária e custosa.
• Escalabilidade: Sob suposições padrão de preparação de estado eficiente e estimativa de observáveis, o circuito exibe uma dependência polilogarítmica da resolução da grade ($O(\text{polylog}(N))$), com o grau polinomial e a contagem de tiles servindo como parâmetros de aproximação explícitos.

Os autores reconhecem que os resultados atuais são limitados a simulações sem ruído e que uma avaliação completa de recursos deve considerar custos de preparação de estado, probabilidades de sucesso de pós-seleção e overheads de tolerância a falhas. Eles observam que estender esta abordagem para 3D ou domínios não periódicos provavelmente exigiria formulações de Combinação de Unitários (LCU) para o projetor vetorial, em vez da abordagem baseada em rotação utilizada aqui.