A Quantum-Classical Surrogate Model for the Collision Operator of the Lattice Boltzmann Method

Resumo Técnico: Um Modelo Substituto Quântico-Clássico para o Operador de Colisão do Método de Lattice Boltzmann

Definição do Problema
O Método de Lattice Boltzmann (LBM) é uma técnica de dinâmica de fluidos computacional (CFD) amplamente utilizada que discretiza espaço, tempo e velocidade. Embora a etapa de streaming do LBM seja linear e passível de algoritmos quânticos, o passo de colisão local introduz dependências não lineares que desafiam soluções puramente quânticas, as quais dependem de operações unitárias. Abordagens quânticas anteriores ao LBM exigiram ou a reformulação do problema como otimização binária, ou o uso de recozimento adiabático, ou o emprego de Circuitos Quânticos Variacionais (VQC) que são limitados a parâmetros de relaxação fixos (por exemplo, $\tau=1$), necessitando de retreinamento para diferentes regimes de escoamento. Além disso, os métodos existentes frequentemente enfrentam dificuldades com a natureza não unitária do operador de colisão de Bhatnagar-Gross-Krook (BGK) e com a alta contagem de portas exigida por técnicas de linearização, como a expansão de Carleman.

Metodologia
Os autores propõem um modelo substituto híbrido quântico-clássico projetado para o hardware NISQ (Noisy Intermediate-Scale Quantum) para aproximar o operador de colisão BGK não linear. O núcleo da metodologia envolve o desacoplamento do processo de colisão em dois componentes distintos:

1. Predição de Equilíbrio via VQC: Um circuito quântico raso e parametrizado prevê a distribuição de equilíbrio pós-colisão. O circuito utiliza um esquema de Codificação por Ângulo (AE) onde variáveis macroscópicas clássicas (densidade $\rho$, velocidade $\mathbf{u}$) são codificadas como ângulos de rotação. A arquitetura emprega um Ansatz Paralelo (PA) com $N=3$ qubits para lidar com os três componentes de velocidade. A estrutura do circuito inclui o re-upload de dados (repetindo camadas de codificação e emaranhamento) para expandir o espectro de frequência acessível, permitindo que o modelo represente não linearidades de ordem superior.
2. Relaxação BGK via Extensão de Circuito: Em vez de aprender toda a dinâmica de relaxação, os autores incorporam a etapa de relaxação BGK diretamente na estrutura do circuito. O parâmetro de relaxação $\omega$ é codificado como uma amplitude de probabilidade quântica que controla uma porta Fredkin (controlled-SWAP). Esta porta mistura o estado pré-colisão com o estado de equilíbrio previsto. Este design permite que o modelo aceite qualquer parâmetro de relaxação $\omega \in [0.5, 2)$ como entrada sem necessidade de retreinamento.
3. Restrições de Pós-Processamento: Para abordar erros de aproximação inerentes ao substituto quântico, o modelo aplica pós-processamento clássico para impor:
   • Conservação de Massa e Momento: Utilizando uma correção algébrica derivada de Corbetta et al.
   • Equivariância Rotacional: Impondo simetria via média de grupo sobre o grupo octaédrico $D_{48}$ (para 3D) ou o grupo diedral $D_8$ (para 2D).
   • Compressão de Órbita Direcional: Reduzindo o número de circuitos treinados ao agrupar direções de velocidade em órbitas de simetria, compartilhando parâmetros entre direções fisicamente equivalentes.

Principais Contribuições

• Parâmetro de Relaxação Variável: O artigo introduz o primeiro substituto híbrido quântico-clássico que codifica explicitamente o parâmetro de relaxação $\omega$ como uma entrada do circuito, permitindo a recuperação do estado pós-colisão BGK completo através de todo o intervalo fisicamente admissível sem retreinamento.
• Análise Arquitetural: Os autores reavaliam criticamente as métricas padrão de VQC (expressividade, emaranhamento, dimensão efetiva) no contexto da dinâmica de fluidos. Eles identificam que o número de repetições de codificação ($L_{SW}$) é o principal impulsionador da capacidade do modelo, enquanto alcançar o emaranhamento de nível Haar não é um pré-requisito para uma precisão competitiva.
• Simulação de Alta Fidelidade 3D: Este trabalho apresenta a primeira simulação substituta de LBM quântico de alta fidelidade em três dimensões, validada pelo benchmark do vórtice de Taylor–Green.
• Compressão de Órbita: A implementação da compressão de órbita direcional reduz significativamente os custos de treinamento e inferência, mantendo a precisão e estabilizando a simulação contra a variância estocástica do treinamento.

Resultados
O substituto foi validado contra o operador BGK clássico usando problemas de referência estabelecidos:

• Vórtice de Taylor–Green (3D): Simulado em $Re \approx 1$ com parâmetros de relaxação $\omega = 0.75$ e $\omega = 1.25$. O substituto reproduziu com precisão o decaimento da energia cinética e a dinâmica de vórtices com erros relativos medianos em torno de $10^{-4}$. Demonstrou-se que a inclusão de correções de conservação e simetria é necessária para a estabilidade de longo prazo.
• Camada de Cisalhamento Dupla (2D): Simulado para números de Reynolds $Re = 100, 400$e$1000$. O modelo capturou com sucesso a instabilidade de Kelvin-Helmholtz, o pareamento de vórtices e a transição para a turbulência. Os erros pontuais permaneceram baixos ($\sim 10^{-4}$ a $10^{-3}$), concentrados em regiões de alta não linearidade (núcleos de vórtices).
• Generalização: O modelo demonstrou generalização robusta através de diferentes números de Reynolds e parâmetros de relaxação sem retreinamento, confirmando a eficácia da codificação de $\omega$ variável.

Significância e Alegações
O artigo afirma que esta estratégia híbrida oferece um caminho prático para realizar o potencial da computação quântica na engenharia de fluidos. Ao delegar a relaxação BGK não unitária para uma extensão de circuito de forma fechada e usar um VQC raso apenas para a predição de equilíbrio não linear, a abordagem mitiga as limitações do hardware NISQ atual (ruído, qubits limitados, profundidade rasa). Os autores enfatizam que sua abordagem evita os custos exponenciais de portas dos métodos de linearização e os requisitos de retreinamento dos substitutos LBM baseados em VQC anteriores. O trabalho sugere que estratégias híbridas quântico-clássicas podem alcançar alta precisão e generalidade, replicando de perto as soluções clássicas para escoamentos de número de Reynolds moderado, validando assim a viabilidade de substitutos de CFD potencializados por computação quântica.