Variational optimization of projected entangled-pair states on the triangular lattice

Resumo Técnico: Otimização Variacional de Estados de Pares de Entrelaçamento Projetados em Redes Triangulares

Declaração do Problema
A simulação de sistemas de muitos corpos quânticos frustrados, particularmente em redes triangulares e de kagome, permanece um desafio significativo na física da matéria condensada. Embora os Estados de Pares de Entrelaçamento Projetados (PEPS) tenham se tornado um arcabouço padrão para simular sistemas quânticos bidimensionais, sua otimização variacional em geometrias triangulares tem sido historicamente limitada. Abordagens convencionais tipicamente mapeiam a rede triangular em uma rede quadrada, introduzindo interações de vizinho próximo de segunda ordem diagonais. Esse mapeamento força as correlações ao longo da direção diagonal a serem roteadas através das direções restantes da rede, o que constitui uma representação ineficiente da estrutura de entrelaçamento intrínseca do sistema. Além disso, métodos existentes frequentemente dependem de mapeamentos indiretos que falham em capturar totalmente a simetria de rotação e a conectividade específica da rede de Bravais triangular, potencialmente limitando a precisão das minimizações de energia variacional para uma dimensão de ligação fixa.

Metodologia
Os autores introduzem um algoritmo de Grupo de Renormalização de Matriz de Transferência de Canto (CTMRG) generalizado, especificamente adaptado para PEPS infinitos na rede triangular. Os componentes metodológicos principais incluem:

• Formulação Triangular Nativa: Em vez de mapear para uma rede quadrada, o algoritmo opera diretamente na geometria triangular. O ansatz PEPS é definido usando tensores de sete índices (um índice físico e seis índices virtuais conectando-se a todos os vizinhos próximos), refletindo o número de coordenação da rede triangular.
• Diferenciação Automática (AD): A abordagem integra AD para computar gradientes de energia de forma eficiente, permitindo a minimização direta da energia variacional sem a necessidade de derivação manual de regras de atualização complexas.
• Procedimento CTMRG Radial: Diferente do CTMRG de rede quadrada, que utiliza direções ortogonais, a variante triangular emprega um esquema de absorção radial. O ambiente consiste em seis tensores de canto ($C$) e doze tensores de aresta ($T$). O algoritmo realiza métodos de potência iterativos para computar tensores de ambiente de ponto fixo.
• Absorção e Renormalização: A absorção dos tensores PEPS no ambiente é realizada radialmente. Dev devido aos vetores de rede não ortogonais, as atualizações dos tensores de canto e de aresta são interdependentes. O método utiliza um processo de renormalização de duas etapas envolvendo dois conjuntos de projetores ($P$ e $Q$) para truncar as dimensões de ligação ampliadas de volta a uma dimensão de ligação de ambiente fixa $\chi_E$, garantindo estabilidade numérica.
• Tratamento de Simetria: O algoritmo não assume a simetria de rotação dos tensores de entrada, permitindo acomodar estados arbitrários, incluindo estados de quebra de simetria e espirais. Ele também pode ser estendido para incorporar simetrias físicas para melhorar a eficiência.

Principais Contribuições

1. Desenvolvimento Algorítmico: O artigo apresenta o primeiro método computacionalmente eficiente e preciso para contrair e otimizar variacionalmente PEPS diretamente na rede triangular, contornando a necessidade de mapeamentos para rede quadrada.
2. Representabilidade Aprimorada: Ao utilizar a geometria nativa, o ansatz TL-PEPS captura as correlações quânticas ao longo de três direções de rede independentes de forma mais fiel que os equivalentes SL-PEPS (rede quadrada). A análise teórica mostra que, para uma dada dimensão de ligação $\chi_B$, o limite superior de entropia de entrelaçamento para TL-PEPS é significativamente maior (escalando como $(8L-2)\log\chi_B$) comparado ao SL-PEPS (escalando como $4L\log\chi_B$).
3. Aumento de Parâmetros Variacionais: Em uma dimensão de ligação fixa, o ansatz TL-PEPS possui um número maior de parâmetros variacionais ($p\chi_B^6$) comparado ao SL-PEPS ($p\chi_B^4$), oferecendo maior poder de expressão para os mesmos recursos computacionais em termos de dimensão de ligação.
4. Estrutura de Benchmarking: Os autores fornecem um protocolo de benchmarking rigoroso comparando o TL-PEPS nativo contra o SL-PEPS mapeado, garantindo comparações justas ao igualar parâmetros variacionais ou analisar a convergência em relação às dimensões de ligação de ambiente.

Resultados
O método foi aplicado a dois modelos paradigmáticos de Heisenberg de spin-1/2 antiferromagnético:

• Modelo de Heisenberg de Rede Triangular:

  • A abordagem TL-PEPS produziu energias variacionais mais baixas do que o SL-PEPS para as mesmas dimensões de ligação.
  • Para $\chi_B = 5$, o TL-PEPS alcançou uma energia de estado fundamental de $E = -0,55136$, que é o valor de energia variacional de rede de tensores mais baixo relatado até o momento para este modelo, superando resultados anteriores de DMRG extrapolado e iPEPS.
  • A extrapolação para o limite de dimensão de ligação infinita resultou em $E = -0,55164$, consistente com resultados SL, mas obtido com propriedades de convergência melhoradas.
  • O parâmetro de ordem magnética $m$ mostrou uma redução sistemática no TL-PEPS ($m \approx 0,157$) em comparação ao SL-PEPS ($m \approx 0,182$), atribuído à sensibilidade aprimorada do ansatz à estrutura de entrelaçamento.

• Modelo de Heisenberg de Rede de Kagome:

  • O TL-PEPS consistentemente forneceu energias variacionais mais baixas em todas as dimensões de ligação acessíveis comparado à formulação SL.
  • A energia de estado fundamental extrapolada foi $E = -0,43929$, melhorando ligeiramente os resultados anteriores de DMRG e diagonalização exata.
  • Os dados do TL-PEPS exibiram um comportamento de escala mais limpo em relação ao inverso do comprimento de correlação ($1/\xi$), facilitando extrapolações mais confiáveis mesmo em dimensões de ligação menores, onde os dados do SL-PEPS mostraram imprecisões.
  • Os resultados apoiam um estado fundamental não magnético consistente com um líquido de spin quântico (QSL), embora a natureza com ou sem gap permaneça uma questão aberta na literatura.

Significância e Alegações
O artigo alega que a formulação nativa de rede triangular oferece uma vantagem significativa na precisão variacional para magnetos quânticos frustrados. As melhorias são atribuídas a:

1. Representação Fiel do Entrelaçamento: O ansatz respeita naturalmente a geometria da rede e as simetrias de rotação, levando a uma representação mais precisa da estrutura de entrelaçamento.
2. Eficiência de Parâmetros: O aumento no número de parâmetros variacionais em uma dimensão de ligação fixa permite que o ansatz explore o espaço de Hilbert de forma mais eficaz.
3. Viabilidade: Apesar do maior custo computacional por etapa (escalando como $O(\chi_B^{13})$ comparado à rede quadrada), o método permanece viável e competitivo, particularmente quando confrontado com simulações de rede quadrada com quantidades equivalentes de parâmetros variacionais.

Os autores concluem que este método fornece uma perspectiva imparcial sobre problemas abertos no magnetismo quântico frustrado e pavimenta o caminho para o estudo de outras redes baseadas em triângulos (ex: honeycomb, dice, maple-leaf) e modelos com interações de vizinho próximo de segunda ordem sem roteamento complexo. A implementação está disponível na biblioteca de código aberto variPEPS.