Site Basis Excitation Ansatz for Matrix Product States

Resumo Técnico: Ansatz de Excitação de Base de Sítio (SBEA) para MPS

1. O Problema

O espectro de excitações elementares de sistemas quânticos é fundamental para a compreensão teórica e a conexão com experimentos. No entanto, obter funções espectrais e relações de dispersão completas é computacionalmente desafiador, especialmente em comparação com propriedades estáticas.

• Limitações Atuais: Métodos tradicionais como a diagonalização exata são limitados a sistemas pequenos. O Grupo de Renormalização de Matriz de Densidade (DMRG) convencional obtém apenas um número limitado de estados excitados. Algoritmos dinâmicos (como DMRG dinâmico ou evolução temporal) são alternativos, mas podem ser custosos.
• Ansatz de Excitação (EA): A versão MPS do "aproximação de modo único" (Excitation Ansatz - EA) introduziu uma melhoria drástica, representando excitações como superposições de ondas planas de MPS onde um tensor local difere do estado fundamental.
• Desafios do EA Padrão:
  1. Requer uma diagonalização completa separada para cada momento $k$, o que é ineficiente.
  2. Métodos para obter o estado fundamental uniforme infinito (como VUMPS) podem falhar na convergência dependendo da inicialização.
  3. A imposição de condições de gauge (ortogonalidade esquerda) no EA tradicional pode prejudicar a convergência e exigir bases maiores.

2. Metodologia Proposta

O autor propõe o Ansatz de Excitação de Base de Sítio (SBEA - Site Basis Excitation Ansatz), uma variação simples e eficiente do EA tradicional. A metodologia divide-se em três etapas principais:

A. Obtenção do Estado Fundamental Infinito (iMPS) via DMRG de Sistema Finito

• Em vez de usar métodos complexos como VUMPS ou iDMRG recursivo, o autor utiliza o DMRG de sistema finito em uma cadeia aberta longa ($N \gg \xi$, onde $\xi$ é o comprimento de correlação).
• Técnica de Inserção de Sítio: Após obter o estado fundamental, insere-se um novo sítio no centro da cadeia. Realiza-se uma diagonalização de Lanczos (sem um estado inicial bom) para otimizar o tensor desse novo sítio.
• Resultado: Como o tensor central possui índices esquerdos e direitos idênticos (devido à simetria do estado fundamental obtido), ele pode ser usado para construir diretamente um MPS uniforme infinito na forma $\dots AAA \dots$. Isso simplifica drasticamente a inicialização e a obtenção do estado fundamental.

B. Construção da Base de Excitação (SBEA)

• No EA tradicional, otimiza-se um tensor $B$ para cada momento $k$ individualmente. No SBEA, propõe-se uma base mínima aproximada $\{B_\alpha\}$ que representa simultaneamente todos os $B(k)$.
• Geração da Base: A base é construída a partir de uma única diagonalização de Lanczos em um único sítio, obtendo os $N_\alpha$ autoestados de menor energia do Hamiltoniano efetivo desse sítio.
• Pré-truncamento: Para manter a eficiência, o ambiente do tensor $B$ é pré-truncado (reduzindo a dimensão de ligação de $\chi$ para $\chi_{tr}$) antes da diagonalização. Isso limita o tamanho da base necessária sem perda significativa de precisão na banda de magnons.
• Ausência de Gauge Esquerdo: Um ponto crucial é não impor a condição de gauge esquerdo (ortogonalidade). O autor demonstra que impor essa condição eleva artificialmente as energias dos estados locais e exige uma base muito maior ($N_\alpha$ grande). Manter a base não ortogonal permite um truncamento severo da base mantendo alta precisão.

C. Teoria de Bandas Não Ortogonais

• Uma vez construída a base $\{B_\alpha\}$, calculam-se os elementos de sobreposição ($O_{\alpha\alpha'}$) e do Hamiltoniano ($H_{\alpha\alpha'}$) entre os estados de base em diferentes sítios.
• Devido ao gap no espectro, esses kernels decaem exponencialmente, permitindo truncamento espacial.
• Para qualquer momento $k$, a relação de dispersão $E(k)$ é obtida resolvendo um problema de autovalor generalizado de tamanho pequeno ($N_\alpha \times N_\alpha$):
  $$\tilde{H}(k)c = E(k) \tilde{O}(k)c$$
  Isso é análogo à teoria de bandas não ortogonais, onde a diagonalização é extremamente rápida comparada à otimização completa por $k$.

D. Funções de Wannier

• O autor demonstra como construir excitações de Wannier (estados localizados no espaço real) a partir das excitações de momento.
• Essas funções são ortonormais e localizadas, permitindo a reconstrução exata dos modos de magnon para todos os momentos. Elas podem servir como base para hamiltonianos efetivos renormalizados de múltiplas excitações.

3. Resultados Principais

O método foi testado na cadeia de Heisenberg antiferromagnética com spin $S=1$ (Cadeia de Haldane), um sistema com gap e comprimento de correlação finito ($\xi \approx 6$).

• Precisão da Dispersão: O SBEA reproduz a dispersão de um magnon com alta precisão. Com apenas $N_\alpha = 7$ estados na base (obtidos com pré-truncamento $\chi_{tr}=8$), o método cobre toda a curva de dispersão, concordando com resultados de DMRG temporal (tDMRG) de alta precisão.
• Gap de Haldane: O valor calculado para o gap de Haldane é $\Delta/J = 0.4107$, muito próximo do valor exato conhecido de $0.410479$.
• Eficiência Computacional:
  • A construção da base e dos kernels leva menos de um minuto.
  • A diagonalização para 500 pontos de $k$ leva apenas alguns segundos em um laptop comum.
  • A precisão é mantida mesmo com uma base pequena, desde que a condição de gauge esquerdo seja omitida.
• Ineficiência do Gauge Esquerdo: O artigo mostra que impor a condição de gauge esquerdo (como feito em trabalhos anteriores) torna os estados locais menos úteis para o truncamento, exigindo $N_\alpha$ muito maiores e prejudicando a convergência.

4. Contribuições Chave

1. Simplicidade na Obtenção de iMPS: Demonstra que um MPS infinito uniforme pode ser obtido de forma trivial e robusta a partir de um DMRG de sistema finito, sem necessidade de algoritmos de inicialização complexos ou recursivos.
2. Eficiência via Base de Sítio (SBEA): Substitui a otimização individual por $k$ por uma diagonalização única de Lanczos local, seguida de uma diagonalização de matriz pequena para cada $k$.
3. Importância do Gauge Não Ortogonal: Identifica que a não imposição da condição de gauge esquerdo é crítica para a eficiência do método, permitindo bases compactas e evitando estados de alta energia artificialmente induzidos.
4. Construção de Wannier: Introduz a construção de excitações de Wannier localizadas a partir do espaço de excitação do MPS, facilitando o estudo de interações de múltiplas partículas.

5. Significado e Impacto

O SBEA representa um avanço significativo na aplicação de MPS a sistemas excitados. Ao reduzir a complexidade algorítmica e aumentar a robustez (evitando falhas de convergência comuns no VUMPS), o método torna-se acessível e eficiente.

• Aplicabilidade: Embora testado em cadeias unidimensionais com gap, a estrutura do método é generalizável para células unitárias maiores e, potencialmente, para sistemas bidimensionais (onde a obtenção de espectros é atualmente um grande desafio).
• Facilitação de Estudos Avançados: A capacidade de gerar rapidamente espectros precisos e construir funções de Wannier localizadas abre caminho para o desenvolvimento de modelos efetivos renormalizados para excitações de múltiplas partículas (ex: espalhamento de dois magnons), algo difícil de tratar com métodos anteriores.

Em resumo, o trabalho oferece uma ferramenta computacional poderosa, simples e altamente eficiente para o cálculo de espectros de excitação em sistemas quânticos unidimensionais, superando limitações técnicas de abordagens anteriores.