Symplectic H2 Model Reduction for High-Dimensional Linear Quantum Systems

Resumo Técnico: Redução de Modelo H2 Simpética para Sistemas Quânticos Lineares de Alta Dimensão

Enunciado do Problema
O artigo aborda o problema de redução de modelo $H_2$ para sistemas quânticos lineares de alta dimensão. Um desafio primário neste domínio é a restrição de Realizabilidade Física (PR). Diferentemente dos sistemas clássicos, os sistemas quânticos lineares devem satisfazer identidades algébricas específicas que acoplam as matrizes do sistema $(A, B, C, D)$ para preservar as relações de comutação canônicas e a estrutura de entrada-saída quântica. Métodos de redução baseados em projeção padrão frequentemente falham em preservar essas restrições não lineares de PR, levando a modelos de ordem reduzida que são fisicamente inválidos. A imposição direta dessas restrições como condições de otimização resulta em problemas algébricos não lineares difíceis que não são escaláveis para sistemas grandes.

Metodologia
Os autores propõem um framework de Petrov-Galerkin Simpético que garante que a PR seja satisfeita por construção, e não por imposição post-hoc.

1. Projeção Simpética: O método baseia-se na construção de um modelo de ordem reduzida via projeção sobre um subespaço simpético. Se a base de tentativa $V \in \mathbb{R}^{2n \times 2r}$ satisfaz a condição simpética $V^\top J_n V = J_r$ (onde $J_k$ é a matriz simpética canônica), e a base de teste é escolhida como $U = J_r^{-1} V^\top J_n$, as matrizes reduzidas resultantes $(A_r, B_r, C_r, D_r)$ satisfazem automaticamente as identidades de PR.
2. Q-IRKA Quântico (Q-IRKA): Os autores desenvolvem uma variante simpética do Algoritmo Iterativo de Krylov Racional (IRKA). O algoritmo opera como uma iteração de ponto fixo:
   • Atualização de Deslocamentos: Os pontos de interpolação (deslocamentos) são atualizados recursivamente com base nos polos do modelo de ordem reduzida atual. Especificamente, os deslocamentos são definidos como o negativo dos polos selecionados (espelhamento de polos), respeitando a simetria espectral do sistema reduzido real.
   • Interpolação Tangencial: Em cada iteração, um conjunto candidato é gerado resolvendo sistemas lineares deslocados $(A - \sigma_i I)^{-1} B t_{i,\ell}$ para deslocamentos e direções tangenciais prescritos.
   • Extração de Base Simpética: Do conjunto candidato, uma base simpética é extraída usando um procedimento do tipo Gram-Schmidt que emparelha vetores com seus conjugados simpéticos ($J_n^\top v$). Uma etapa subsequente de normalização garante que a base satisfaça estritamente $V^\top J_n V = J_r$.
   • Preservação de Estrutura: As matrizes reduzidas são obtidas exclusivamente através da projeção $A_r = UAV$, $B_r = UB$, $C_r = CV$ e $D_r = D$. Nenhuma modificação independente das matrizes reduzidas ocorre, garantindo que a PR seja preservada até a precisão da máquina ao longo da iteração.

Principais Contribuições

• Algoritmo Q-IRKA: A introdução de um algoritmo recursivo escalável e que preserva estrutura para redução de modelo $H_2$ ótima de sistemas quânticos lineares. Ele combina a eficiência do IRKA (usando resoluções lineares deslocadas) com restrições simpéticas estritas.
• PR por Construção: O framework elimina a necessidade de resolver problemas de otimização não linear com restrições. A realizabilidade física é uma propriedade inerente da geometria de projeção, garantindo que o modelo reduzido permaneça um sistema quântico válido.
• Análise Numérica Sistemática: O artigo fornece um estudo abrangente de como a qualidade da redução depende de características estruturais físicas, especificamente:
  • Geometria de Dissipação: Amortecimento homogêneo versus heterogêneo.
  • Posicionamento de Canais: Configurações de baixo canal onde múltiplos osciladores compartilham caminhos dissipativos.
  • Heterogeneidade do Sistema: Variações nas frequências locais e nas intensidades de amortecimento.

Resultados
Experimentos numéricos foram conduzidos em duas famílias de referência:

1. Sistemas Hamiltonianos de Porta de Baixo Canal: Cadeias de osciladores com canais externos agregados.
2. Sistemas Inspirados na Cadeia de Kitaev Bosônica: Sistemas estruturados dominados por transporte com acoplamento entre vizinhos mais próximos.

Principais Descobertas:

• Precisão e Estabilidade: O Q-IRKA produz modelos de ordem reduzida com aproximação $H_2$ precisa. Os resíduos de simplicidade e PR ($\Delta_{symp}$ e $\Delta_{PR}$) permanecem na precisão da máquina ($\approx 10^{-14}$ a $10^{-16}$) ao longo das iterações.
• Dependência da Heterogeneidade: A qualidade da redução é fortemente influenciada pela geometria de dissipação do sistema. Sistemas homogêneos com decaimento rápido dos valores singulares de Hankel produzem aproximações de baixa ordem com alta precisão (erro $H_2 \approx 10^{-5}$). Sistemas heterogêneos exibem decaimento espectral mais lento, resultando em ranks dinâmicos efetivos mais altos e erros maiores (erro $H_2 \approx 10^{-3}$), mesmo com a mesma ordem reduzida.
• Custo Computacional: O custo computacional escala principalmente com a dimensão do sistema $n$ e a ordem reduzida $r$. A dependência do número de canais $m$ é suave. Casos heterogêneos geralmente requerem mais iterações para convergência e incorrem em custos computacionais mais altos devido ao decaimento espectral mais lento.
• Comparação: Em uma comparação em pequena escala, o Q-IRKA alcançou um erro $H_2$ menor (1,2130) em comparação com a truncagem quase balanceada (1,86) e um método de otimização não convexa (1,25).

Significado
O artigo afirma que a redução de modelo $H_2$ escalável para sistemas quânticos lineares de grande escala pode ser alcançada sem sacrificar a estrutura física subjacente. Ao incorporar as restrições de PR diretamente na geometria de projeção, o método evita a intratabilidade computacional da otimização com restrições. Os resultados demonstram que, embora o desempenho da redução seja sensível a parâmetros físicos como geometria de dissipação e heterogeneidade, o framework proposto mantém robustamente a realizabilidade física e fornece modelos de ordem reduzida precisos, adequados para ciclos de projeto, síntese de controladores e estudos de escalabilidade em engenharia quântica. O framework é notado por se estender naturalmente a sistemas quânticos lineares não quadrados, embora isso seja reservado para trabalhos futuros.