IntegrateUnitary.jl: A Julia package for symbolic integration over Haar measures

Resumo Técnico de IntegrateUnitary.jl: Um pacote Julia para integração simbólica sobre medidas de Haar

Declaração do Problema
A integração simbólica sobre a medida de Haar de grupos compactos é uma técnica fundamental na teoria de matrizes aleatórias (RMT) e na teoria da informação quântica (QIT), essencial para caracterizar propriedades típicas de sistemas quânticos, como a entropia média de emaranhamento, e para analisar o caos quântico e designs unitários $t$. Embora soluções em forma fechada existam por meio do cálculo de Weingarten, sua avaliação prática é dificultada pela complexidade combinatória das somas sobre grupos simétricos e pelas intricâncias da teoria das representações. Ferramentas de software existentes, como o RTNI baseado em Mathematica e a biblioteca Python Haarpy, enfrentam limitações em extensibilidade, desempenho e na capacidade de lidar nativamente com dimensões simbólicas ($d$) em fluxos de trabalho diversos. Além disso, traduzir expressões matemáticas de alto nível e livres de coordenadas (por exemplo, traços de produtos) nas contrações de tensores baseadas em índices exigidas pelo cálculo de Weingarten frequentemente requer um pré-processamento manual propenso a erros por parte do usuário.

Metodologia
O artigo apresenta o IntegrateUnitary.jl, um pacote Julia puro projetado para automatizar e otimizar a integração simbólica sobre medidas de Haar. A metodologia central baseia-se no cálculo de Weingarten, que expressa integrais de funções polinomiais sobre grupos compactos como somas sobre conjuntos combinatórios (permutações para grupos unitários, partições de pares para grupos ortogonais e simpléticos).

Componentes metodológicos-chave incluem:

• Pipeline de Integração Unificado: O pacote normaliza expressões de entrada usando SymbolicUtils.jl, expande-as em monômios e as despacha para engines especializados com base no tipo de medida (Unitária, Ortogonal, Simplética, Gaussiana, etc.).
• Suporte a Dimensões Simbólicas: Aproveitando a natureza polinomial dos caracteres de grupo e das funções de Schur, o pacote calcula muitas integrais como funções racionais exatas da dimensão $d$. Isso permite a derivação de expansões assintóticas exatas e leis de escala universais por meio de séries de Laurent automatizadas, em vez de exigir dimensões inteiras concretas para cada cálculo.
• Lógica de Traço Simbólica: Uma interface de alto nível permite que os usuários insiram expressões livres de coordenadas (por exemplo, tr(U * A * U' * B)). O sistema reconstrói automaticamente os gráficos de Weingarten necessários ou as contrações de Wick a partir dessas entradas sem índices, eliminando a gestão manual de índices.
• Otimizações Algorítmicas:
  • Grupos Unitários: Utiliza a regra de Murnaghan-Nakayama para o cálculo de caracteres e a fórmula Gancho-Conteúdo para dimensões de representação, permitindo o manuseio eficiente de $d$ simbólico.
  • Grupos Ortogonais/Simpléticos: Emprega uma relação de dualidade onde integrais simpléticas são computadas reutilizando a engine ortogonal com um parâmetro de dimensão com sinal ($d \to -d$), evitando inversões de matriz separadas.
  • Ensembles Gaussianos/Ginibre: Implementa o teorema de Wick (teorema de Isserlis) para contração de momentos.
  • Memorização: Armazena em cache caracteres, dimensões e valores de Weingarten para acelerar tarefas repetidas.
• Integração com Redes de Tensores: Uma ponte nativa para ITensors.jl permite a média simbólica de grandes redes de tensores aleatórios analisando topologias de contração e expandindo-as em somas de tensores contraídos deterministicamente ponderados por funções de Weingarten.

Principais Contribuições
O artigo introduz um pacote Julia abrangente que aborda as limitações das ferramentas existentes através de várias características distintas:

1. Suporte a Ensembles Diversos: Integração unificada sobre grupos compactos ($U(d)$, $O(d)$, $Sp(d)$, $SU(d)$), ensembles circulares (CUE, COE, CSE), ensembles gaussianos e Ginibre, grupos de permutação, estados puros aleatórios, variedades de Stiefel e designs unitários $t$.
2. Dimensões Simbólicas Nativas: Ao contrário de muitos predecessores, o IntegrateUnitary.jl trata $d$ como uma variável simbólica nativa, suportando resultados racionais exatos e expansões assintóticas automatizadas (séries $1/d$) para uma ampla gama de observáveis.
3. Abstração de Alto Nível: A lógica de traço simbólica e os tipos SymbolicMatrix permitem que os usuários definam integrais complexas sem contração explícita de índices, reduzindo o potencial de erros e alinhando a estrutura do código com a notação matemática.
4. Desempenho: Ao aproveitar a compilação just-in-time (JIT) do Julia e a memorização eficiente, o pacote lida com momentos de alto grau e computações simbólicas em grande escala de forma mais eficiente do que as ferramentas existentes.
5. Integrais HCIZ: Implementação da integral de Harish-Chandra-Itzykson-Zuber para entradas de autovalores simbólicas e numéricas, com tratamento específico para espectros degenerados.

Resultados e Benchmarks
Os autores fornecem benchmarks extensos comparando o IntegrateUnitary.jl com o Haarpy e o RTNI:

• Desempenho: O pacote demonstra acelerações significativas, particularmente para momentos de alto grau. Para momentos diagonais unitários ($|U_{11}|^{2k}$) com $d$ simbólico, as acelerações sobre o Haarpy variam de 12,0$\times$ a 88,5$\times$. Para grupos ortogonais, a vantagem é ainda mais pronunciada, com acelerações superiores a $2,7 \times 10^4\times$ para momentos de grau 10, atribuídas ao solucionador univariado otimizado do IntegrateUnitary.jl versus a abordagem de inversão de matriz completa do Haarpy.
• Escalabilidade: O pacote lida eficientemente com polinômios de alto grau (até $2k \approx 10\text{--}12$ em hardware padrão) e escala efetivamente com a dimensão $d$ para integrações de redes de tensores.
• Correção: Todos os resultados são verificados contra valores exatos conhecidos, limites assintóticos e expectativas teóricas (por exemplo, o resultado de Page para o emaranhamento médio).
• Limitações: O artigo observa que momentos de traço puro $|tr(U)|^{2k}$ dependem de $d$ como uma função degrau e requerem dimensões inteiras concretas. Da mesma forma, a integração $SU(d)$ está atualmente limitada a polinômios balanceados, e integrais HCIZ com autovalores simbólicos degenerados exigem tratamento manual.

Significado
O artigo posiciona o IntegrateUnitary.jl como uma ferramenta crítica para físicos teóricos e cientistas da informação quântica. Ao preencher a lacuna entre o cálculo de Weingarten abstrato e a computação prática, ele permite a derivação rigorosa de propriedades para canais quânticos aleatórios, estatísticas de emaranhamento e medidas de embaralhamento sem depender de amostragem numérica frágil. A capacidade do pacote de lidar com dimensões simbólicas e lógica de traço de alto nível reduz significativamente a carga cognitiva e o potencial de erro em derivações complexas. Os autores afirmam que a ferramenta facilita o estudo de fenômenos anteriormente intratáveis devido à complexidade combinatória, servindo como um ambiente versátil tanto para derivações teóricas quanto para experimentos numéricos. Trabalhos futuros são delineados para estender o suporte a grupos de Lie excepcionais e refinar correções de $d$ finito para $SU(d)$.