Symmetry-Adapted Eigensolver from Real-Space Sampling of Tensor-Network-Bridged Quantum Circuits

No mundo microscópico da física quântica, átomos e elétrons não se comportam como pequenas bolas de bilhar; eles existem em um estado de conexão profunda conhecido como emaranhamento. Quando as partículas estão emaranhadas, o estado de uma influencia instantaneamente o estado de outra, não importa o quão distantes estejam. Esse fenômeno cria uma complexidade matemática que cresce explosivamente à medida que mais partículas são adicionadas a um sistema. Durante décadas, cientistas lutaram para simular esses sistemas fortemente conectados em computadores clássicos porque a quantidade de informação necessária para descrevê-los rapidamente excede a capacidade até mesmo dos supercomputadores mais poderosos. Para enfrentar isso, pesquisadores frequentemente recorrem a uma técnica chamada rede de tensores, que atua como um algoritmo de compressão, eliminando detalhes desnecessários para manter apenas as conexões mais essenciais. No entanto, mesmo esses métodos comprimidos enfrentam dificuldades quando as conexões se estendem por todo o sistema, exigindo um poder computacional imenso para capturar o quadro completo.

Uma equipe de físicos da Universidade de Zhejiang propôs agora uma nova maneira de navegar por essa complexidade ao separar o ato de amostrar um estado quântico do ato de defini-lo. A abordagem deles, detalhada em um estudo recente, combina os pontos fortes de circuitos quânticos e redes de tensores clássicas para resolver um tipo específico de problema difícil: encontrar o estado de menor energia de um sistema de spins interagentes, que serve como um modelo para materiais magnéticos. Em vez de tentar calcular todo o sistema de uma vez, eles dividem a rede de spins em blocos menores. Dentro de cada bloco, um circuito quântico captura o emaranhamento local, enquanto uma estrutura matemática clássica conecta esses blocos para restaurar as correlações globais. A inovação principal reside em como eles lidam com a simetria do sistema. Em vez de forçar o computador a calcular cada arranjo possível de spins, eles agrupam arranjos que são equivalentes sob operações de simetria, como rotação ou reflexão. Eles então definem o estado quântico não pela soma das amplitudes desses arranjos, mas pela agregação de suas probabilidades. Essa distinção permite que eles gerem amostras do sistema de forma direta e eficiente, sem os processos lentos de caminhada aleatória normalmente exigidos para explorar tais espaços vastos.

Os pesquisadores testaram seu método em um modelo bem conhecido de frustração magnética, onde forças competitivas tornam difícil para o sistema se estabelecer em uma ordem simples. Eles aplicaram sua técnica em cadeias unidimensionais e grades bidimensionais de spins, incorporando simetrias como translação e reflexão para garantir que os resultados respeitassem as leis físicas que regem o sistema. Em suas simulações, eles usaram uma quantidade muito pequena de recursos computacionais, especificamente um parâmetro conhecido como dimensão de ligação, que controla a complexidade das conexões entre os blocos. Mesmo com esse parâmetro definido para um valor tão baixo quanto seis, o método produziu resultados com um erro absoluto de energia por sítio de apenas uma parte em cem mil. Esse nível de precisão é comparável aos métodos tradicionais que exigem dimensões de ligação na casa dos milhares, uma diferença que representa uma redução massiva no custo computacional. A equipe descobriu que, ao usar uma base adaptada à simetria, o emaranhamento de longo alcance que normalmente sobrecarrega o cálculo era efetivamente carregado pela própria estrutura matemática da base, deixando a estrutura computacional para lidar apenas com os detalhes restantes.

Um aspecto crucial deste trabalho é como os pesquisadores construíram os circuitos quânticos usados para gerar as amostras. Em vez de projetar circuitos complexos do zero, o que é uma tarefa difícil e propensa a erros, eles derivaram as conexões necessárias da matriz de densidade reduzida de um sistema menor e exatamente solúvel. Esta matriz descreve as propriedades estatísticas de um pequeno bloco de spins dentro de um sistema maior. Ao extrair os padrões mais significativos desta matriz, eles criaram um mapa isométrico simplificado que captura eficientemente a estrutura essencial do emaranhamento. Esta abordagem evita a necessidade de circuitos profundos e altamente parametrizados que são difíceis de otimizar, incorporando efetivamente a física do problema diretamente no processo de amostragem. O resultado é uma estrutura híbrida onde os circuitos quânticos atuam como um gerador de configurações no espaço real, e a rede de tensores clássica serve como uma ponte que organiza essas configurações em um estado coerente e que respeita a simetria.

O estudo demonstra que essa separação de papéis permite uma simulação altamente eficiente de sistemas fortemente correlacionados. Os pesquisadores mostraram que, para um sistema de sessenta e quatro spins em uma cadeia unidimensional e uma grade de sessenta e quatro spins em duas dimensões, seu método poderia reproduzir a energia do estado fundamental com uma precisão notável. Eles verificaram essas descobertas comparando seus resultados com a diagonalização exata e o grupo de renormalização de matriz de densidade, que são considerados os padrões ouro para tais problemas. A concordância foi tão próxima que as pequenas discrepâncias puderam ser atribuídas ao tamanho finito da dimensão de ligação, que poderia ser futuramente reduzida por extrapolação dos resultados. Isso sugere que o método não é apenas preciso, mas também escalável, já que a base adaptada à simetria lida naturalmente com as correlações de longo alcance que tipicamente exigem recursos computacionais massivos para serem simuladas.

Olhando para o futuro, os autores observam que, embora sua implementação atual utilize uma estrutura de rede de tensores clássica, a estrutura é projetada para ser compatível com hardware quântico distribuído. Os blocos do sistema podem ser processados em paralelo em processadores quânticos separados, com a rede clássica coordenando os resultados. Esta arquitetura evita a necessidade de processadores quânticos compartilharem emaranhamento diretamente, o que é um obstáculo significativo na computação quântica distribuída atual. A etapa generativa do algoritmo pode ser executada sequencialmente em um dispositivo quântico, produzindo amostras exatas do sistema sem a necessidade de estimativa de amplitude complexa. No entanto, a etapa final de calcular o estimador de energia continua sendo um desafio para a implementação direta em hardware, pois requer consultar a amplitude de configurações específicas em vez de medir valores de expectativa. Os pesquisadores sugerem que reformular este estimador para depender de quantidades diretamente mensuráveis é o próximo passo crítico em direção a uma realização totalmente quântico-clássica.

Em última análise, este trabalho oferece uma nova perspectiva sobre como simular a matéria quântica. Ao reconhecer que a definição de um estado quântico e a amostragem desse estado podem ser desacopladas, os pesquisadores encontraram uma maneira de aproveitar o poder da simetria para reduzir a carga computacional. O método prova que, mesmo com recursos muito limitados, é possível capturar a intrincada dança das partículas emaranhadas, desde que se saiba onde olhar e como organizar a informação. As descobertas sugerem que o caminho para simular sistemas quânticos maiores e mais complexos pode não exigir computadores cada vez maiores, mas sim maneiras mais inteligentes de organizar a informação que já está disponível. À medida que o campo avança para dispositivos quânticos mais poderosos, esta abordagem fornece uma ponte robusta e eficiente entre as técnicas de simulação clássica e as capacidades emergentes do hardware quântico.