Graph automorphisms to obtain Clifford symmetries in open and closed qudit models

Resumo Técnico: Automorfismos de Grafos para Obter Simetrias de Clifford em Modelos de Qudits Abertos e Fechados

Definição do Problema
A identificação de simetrias em sistemas quânticos é crucial para compreender propriedades do estado fundamental, dinâmica e termodinâmica, bem como para reduzir custos computacionais em simulações numéricas através da diagonalização em blocos. Embora as simetrias sejam frequentemente descobertas por inspeção, elas podem ser não óbvias em modelos complexos. As abordagens algorítmicas existentes são amplamente limitadas a simetrias de Pauli ou tornam-se computacionalmente proibitivas à medida que o tamanho do sistema aumenta. Especificamente, há uma carência de algoritmos eficientes para encontrar simetrias de Clifford — uma classe de simetrias eficientemente implementáveis em computadores clássicos — para sistemas de qudits gerais, incluindo sistemas quânticos tanto fechados (Hamiltonianos) quanto abertos (Liouvillianos).

Metodologia
Os autores propõem um algoritmo que mapeia o problema de encontrar simetrias de Clifford para um problema de Automorfismo de Grafo (AG). A metodologia baseia-se na representação simplética de operadores de qudits e procede através das seguintes etapas:

1. Representação Simplética: O sistema é representado usando uma formalidade de tableau. Para um sistema fechado, o Hamiltoniano $\hat{H}$ é codificado como uma soma ponderada de strings de Pauli, definida por um vetor de coeficientes $\vec{c}$, uma matriz de expoentes $p$ (representando strings de Pauli como vetores em $\mathbb{Z}_d^{2n}$) e um vetor de fase $\vec{\eta}$. Para sistemas abertos, isso é estendido ao espaço de Liouville, dobrando o número de qudits para representar o gerador da equação mestre GKSL como um tableau de superoperador.
2. Construção do Grafo: O Hamiltoniano (ou Liouvilliano) é mapeado para um grafo direcionado colorido $G = (V, E, X_V, X_E)$:
   • Vértices ($V$): Cada string de Pauli no tableau corresponde a um vértice.
   • Cores dos Vértices ($X_V$): Os vértices são coloridos por seus coeficientes $c_i$. Como as operações de Clifford preservam a magnitude dos coeficientes, uma simetria válida deve mapear vértices para outros com cores idênticas.
   • Rótulos de Arestas ($X_E$): Arestas direcionadas entre vértices representam o produto simplético $\langle \vec{p}_i, \vec{p}_j \rangle \pmod d$, que codifica as relações de comutação entre as strings de Pauli. Uma simetria válida deve preservar esses rótulos de aresta.
3. Tratamento de Dependência e Fase:
   • Dependências Lineares: Para garantir que a permutação respeite a estrutura de dependência linear das strings de Pauli, os autores empregam duas estratégias: (a) aumentar o grafo com vértices auxiliares representando circuitos mínimos (subconjuntos dependentes), ou (b) realizar uma verificação de "automorfismo de código" post-hoc em permutações candidatas para verificar se preservam o espaço linha do gerador.
   • Correção de Fase: Como as portas de Clifford incluem vetores de fase que não são estritamente invariantes, o algoritmo primeiro encontra uma permutação candidata (AG) que satisfaz os invariantes estruturais. Em seguida, resolve um sistema linear para determinar se existe um vetor de fase $\vec{\phi}$ válido para corrigir o desajuste de fase residual, verificando efetivamente se o candidato é uma verdadeira simetria até uma porta de Pauli.
4. Algoritmos de Busca: O problema de AG é resolvido usando duas abordagens:
   • Bliss/igraph: Utilizando bibliotecas existentes para encontrar os geradores do grupo de automorfismo.
   • Heurística MRV: Uma busca de backtracking personalizada usando uma heurística de Valor Mínimo Restante (Minimum Remaining Values) para podar o espaço de busca de forma eficiente, frequentemente acelerada pela refinação de cores de Weisfeiler-Leman (WL) para distinguir vértices com base em seus vizinhos de comutação local.

Principais Contribuições

• Estrutura Algorítmica: O artigo fornece o primeiro algoritmo detalhado para encontrar simetrias de Clifford para sistemas arbitrários de qudits abertos ou fechados, mapeando o problema para Automorfismos de Grafos.
• Extensão para Sistemas Abertos: A formalidade é generalizada de Hamiltonianos para equações mestras GKSL, permitindo a descoberta de simetrias em sistemas quânticos dissipativos ao tratar como tableaux de Liouvillianos de qudits dobrados.
• Eficiência e Escalabilidade: A abordagem aproveita a representação simplética para manter o problema classicamente eficiente. O uso de coloração de grafos e heurísticas reduz significamente o espaço de busca em comparação com verificações de permutação por força bruta.
• Tratamento de Dependências: O trabalho detalha métodos para incorporar estruturas de dependência linear (circuitos) seja via aumento de grafo ou verificações de automorfismo de código, garantindo que as permutações encontradas correspondam a simetrias físicas válidas.

Resultados
Os autores testaram o algoritmo em vários modelos:

• Hamiltonianos de Pauli Aleatórios: Demonstraram escalonamento quadrático com o tamanho do sistema para encontrar simetrias injetadas, limitado primariamente pelo tempo de construção do grafo ($O(M^2)$).
• Modelos Físicos: Identificaram com sucesso simetrias em modelos de Campo Transverso de Ising (TFI) (geometrias de escada, quadrada, triangular), modelos Fermiônicos (cadeias de Fermi-Hubbard e $t-V$ desordenadas) e modelos Heisenberg XXZ.
• Escalonamento: Para modelos com interações locais (ex: TFI, $t-V$), o tempo para encontrar simetrias escala aproximadamente de forma quadrática com o número de qudits $n$. Para modelos all-to-all (ex: Heisenberg XXZ com acoplamentos todos-com-todos), o escalonamento depende do número de termos de Pauli $M$, que escala como $n^2$, levando a diferentes características de desempenho.
• Código Torico: Em regimes altamente simétricos, o algoritmo encontrou simetrias para até $n=28$ qubits. Para instâncias maiores, os autores usaram a saída algorítmica para extrapolar analiticamente a estrutura de simetria para tamanhos de sistema arbitrários. A estratégia de aumento de grafo provou ser benéfica para o Código Torico para evitar a explosão combinatória ao verificar dependências.
• Sistemas Abertos: O algoritmo lidou com sucesso com sistemas abertos (ex: TFI com desfasagem), observando o aumento do custo computacional devido ao dobro da dimensão do espaço de Hilbert necessária para a representação do Liouvilliano.

Significância e Alegações
O artigo afirma preencher uma lacuna crítica na literatura ao fornecer um algoritmo prático para encontrar simetrias de Clifford, uma tarefa para a qual não existia solução geral anteriormente. Os autores enfatizam que seu método não se limita a modelos específicos, mas aplica-se a qualquer sistema representável como uma soma de Pauli em qudits.

A significância reside na capacidade de:

1. Automatizar a Descoberta de Simetrias: Ir além da inspeção manual para a detecção algorítmica de simetrias de Clifford não triviais.
2. Lidar com Sistemas Abertos: Estender a análise de simetria para a dinâmica dissipativa, identificando setores de espaço de operador invariantes em sistemas quânticos abertos.
3. Permitir Insight Analítico: Usar resultados numéricos de sistemas de pequeno a médio porte para derivar indutivamente formas analíticas de simetrias para tamanhos de sistema arbitrários (como demonstrado com o Código Torico).

Os autores permanecem modestos quanto à complexidade computacional, reconhecendo que, embora o problema de AG seja geralmente difícil, a estrutura específica das simetrias de Clifford (via invariantes como coeficientes e relações de comutação) permite soluções eficientes em muitos regimes físicos. Eles observam que a correção de fase e as verificações de dependência são etapas necessárias que distinguem este problema de um AG padrão, embora, empiricamente, essas verificações não aumentem significativamente a complexidade para os modelos testados. O trabalho fornece uma implementação em Python ("Sympleq") para facilitar pesquisas e aplicações futuras.