Faster matrix product state preparation by exploiting symmetry-induced block-sparsity

Resumo Técnico: Preparação Mais Rápida de Estados Produto de Matriz Explorando Esparsidade em Blocos Induzida por Simetria

Enunciado do Problema
Estados Produto de Matriz (MPS) são uma ferramenta primária para simular sistemas quânticos em química e física da matéria condensada, particularmente para encontrar estados fundamentais de Hamiltonianos locais unidimensionais com gap. A preparação desses estados em computadores quânticos é uma etapa crítica para algoritmos como a Estimativa de Fase Quântica (QPE). Embora existam métodos padrão de preparação tolerante a falhas (por exemplo, a abordagem de profundidade linear assistida por ancila), eles frequentemente tratam os tensores MPS como matrizes densas. No entanto, muitos sistemas físicos de interesse possuem simetrias $U(1)$ (conservação do número de partículas e projeção de spin), que induzem uma estrutura de esparsidade em blocos nos tensores MPS. Algoritmos clássicos como o Grupo de Renormalização de Matriz de Densidade (DMRG) exploram essa esparsidade para eficiência, mas os circuitos padrão de preparação quântica não a aproveitam plenamente, levando a custos de recursos mais altos (especificamente contagens de portas Toffoli) do que o necessário.

Metodologia
Os autores propõem um método para reduzir o custo de preparação de MPS com esparsidade em blocos transformando as operações unitárias subjacentes em uma forma de blocos diagonais. A abordagem baseia-se no esquema padrão de preparação assistida por ancila, onde um MPS é preparado por meio de uma sequência de unitárias $U_i$.

1. Exploração da Esparsidade em Blocos:
   Os tensores MPS $M_i$ para sistemas com simetrias abelianas globais são esparsos em blocos. Os autores demonstram que a matriz $U'_i$ (formada pela pilha das transpostas dos blocos do tensor) herda essa estrutura, onde cada linha e coluna de blocos contém apenas um único bloco não nulo.

   • Estratégia de Permutação: O método emprega matrizes de permutação de linhas e colunas ($P_r$ e $P_c$) para reorganizar os blocos não nulos de $U'_i$.
   • Blocos Diagonalização: Ao preencher as entradas não especificadas (denotadas como $\ast$) da extensão unitária $U_i = [U'_i \mid \ast]$ com blocos retangulares complementares, a matriz é transformada em uma unitária de blocos diagonais $V$. Isso permite que a síntese de $U_i$ seja decomposta na síntese de blocos unitários menores e independentes, em vez de uma única unitária densa grande.

2. Síntese Unitária Otimizada:
   As unitárias de blocos diagonais são sintetizadas usando rotações de Givens.

   • Otimização de Valores Reais: Os autores modificam a abordagem de síntese unitária de Berry et al. [8]. Para unitárias de valores reais (comuns nesses sistemas), eles demonstram que o número de portas Toffoli necessárias pode ser reduzido por um fator de $\sqrt{2}$. Isso é alcançado expressando camadas de Givens usando apenas rotações $R_y$ (removendo portas $R_z$) e utilizando um framework de gradiente de fase.
   • Alinhamento de Camadas: Para sintetizar a matriz de blocos diagonais completa de forma eficiente, os autores alinham as camadas de rotação de Givens dos blocos individuais. Eles lidam com as matrizes de permutação ($W$ e $Q$) necessárias para transformar entre a base original e a base de blocos diagonais usando QROAM (Memória de Acesso Aleatório Quântica) com correções de inversão de sinal, garantindo que o custo permaneça gerenciável.

Principais Contribuições

• Transformação de Blocos Diagonais: Um procedimento novel para converter tensores MPS esparsos em blocos em unitárias de blocos diagonais por meio de permutações de linhas/colunas e preenchimento estratégico da extensão unitária. Isso reduz o custo de síntese de ser determinado pela dimensão total do vínculo para ser determinado pelo tamanho do maior bloco.
• Aceleração de $\sqrt{2}$ para Unitárias Reais: Uma modificação do algoritmo de síntese unitária de Berry et al. que reduz o custo de Toffoli para unitárias de valores reais por um fator de $\sqrt{2}$, eliminando rotações $R_z$ desnecessárias.
• Benchmarking de Recursos: Implementação do algoritmo no framework Qualtran e estimativa rigorosa de recursos para vários sistemas moleculares, incluindo espaços ativos de enzimas P450 e clusters de ferro-enxofre.

Resultados
Benchmarks numéricos em sistemas moleculares fortemente correlacionados (por exemplo, P450, [Fe$_2$S$_2$], [Fe$_4$S$_4$]) com dimensões de vínculo e estados de spin variáveis mostram melhorias significativas:

• Redução de Custo Toffoli: O método proposto alcança fatores de melhoria de custo de Toffoli variando de 10 a 30 em comparação com os métodos de preparação densa mais avançados (tanto com quanto sem a otimização de valores reais).
• Correlação com Esparsidade: A aceleração é mais pronunciada em sistemas com dimensões de vínculo mais altas, onde a densidade de entradas não nulas nos tensores é menor. Os autores observam que a aceleração alcançada é aproximadamente metade do máximo teórico (que seria $\sqrt{2}$ vezes o inverso da densidade do tensor), sugerindo espaço para otimização adicional na etapa de bloqueio.
• Sobrecarga de Qubits: Os circuitos de permutação introduzem um custo comparável a uma única camada da síntese unitária, o que é gerenciável em relação às economias totais.

Significado e Alegações
O artigo afirma que explorar a esparsidade em blocos induzida por simetria é um passo crucial para tornar a estimativa de energia do estado fundamental de ponta a ponta mais eficiente em computadores quânticos tolerantes a falhas.

• Relevância para QPE: Avanços algorítmicos recentes reduziram o custo da parte de estimativa de fase da QPE, estreitando a lacuna com o custo da preparação do estado inicial. Ao reduzir significativamente o custo de preparação do MPS, este trabalho permite o uso de estados iniciais de maior qualidade (com dimensões de vínculo maiores) sem incorrer em sobrecargas de recursos proibitivas.
• Integração de Clássico e Quântico: O trabalho serve como uma ponte entre técnicas de simulação clássica (que já utilizam esparsidade em blocos) e simulação quântica, garantindo que as estimativas de recursos quânticos reflitam as simetrias físicas dos sistemas-alvo.
• Direções Futuras: Os autores sugerem que melhorias adicionais poderiam ser feitas otimizando a própria etapa de bloqueio ou adaptando o esquema para formas $SU(2)$ adaptadas ao spin. Eles também destacam a necessidade de comparar essa preparação direta de MPS contra outros métodos, como preparar somas esparsas de determinantes de Slater ou usar técnicas adiabáticas/dissipativas.

Os autores concluem que sua abordagem é uma parte integral de conectar simulações clássicas e quânticas para química e física da matéria condensada, fornecendo um caminho concreto para reduzir os requisitos totais de recursos para vantagem quântica nesses domínios.