Matrix Product Operator Encodings of the Magnus Expansion and Dyson Series

Resumo Técnico: Codificações de Operador Produto Matricial para a Expansão de Magnus e a Série de Dyson

Declaração do Problema
Simular a evolução temporal real de modelos de rede quântica unidimensionais governados por Hamiltonianos dependentes do tempo, $H(t)$, permanece um desafio significativo na física computacional de muitos corpos. Embora o operador de evolução temporal $U(t, t_0)$ seja formalmente definido por meio de um exponencial ordenado no tempo, sua avaliação analítica é geralmente impossível para sistemas interagentes. Abordagens numéricas padrão frequentemente dependem de dividir a evolução temporal em pequenos passos onde o Hamiltoniano é tratado como constante (decomposição de Trotter-Suzuki) ou de aproximar o operador de evolução via métodos de rede tensorial, como o Princípio Variacional Dependente do Tempo (TDVP) ou abordagens baseadas em Krylov. No entanto, para Hamiltonianos que flutuam rapidamente, esses métodos exigem passos de tempo extremamente pequenos, levando a altos custos computacionais. Além disso, métodos existentes de rede tensorial de alta ordem, como aqueles baseados na série de Taylor para Hamiltonianos independentes do tempo, não se generalizam de forma direta para casos dependentes do tempo sem incorrer em dimensões de ligação proibitivas ou perder a extensividade de tamanho no limite termodinâmico.

Metodologia
Os autores introduzem uma codificação de Operador Produto Matricial (MPO) tanto para a expansão de Magnus quanto para a série de Dyson, especificamente adaptada para Hamiltonianos dependentes do tempo. O trabalho baseia-se em pesquisas anteriores (Ref. 22) que codificaram com sucesso a série de Taylor de Hamiltonianos independentes do tempo como MPOs extensivos de tamanho.

A metodologia prossegue em duas rotas principais:

1. Expansão de Magnus como um MPO: O operador de evolução temporal é expresso como $U(t, t_0) = e^{\Omega(t, t_0)}$, onde $\Omega$ é expandido em comutadores do Hamiltoniano em diferentes instantes. Os autores demonstram que comutadores de MPOs de primeiro grau (representando Hamiltonianos locais) podem ser construídos como MPOs de primeiro grau. Isso permite que o operador de Magnus $\Omega$ seja construído como uma combinação linear de MPOs de primeiro grau. Subsequentemente, a exponencial $e^\Omega$ é computada usando a construção estabelecida de MPO de série de Taylor.
2. Série de Dyson como um MPO: Os autores propõem uma abordagem mais direta codificando a série de Dyson, $U(t, t_0) = \mathcal{T} \exp(-i \int H(t') dt')$, diretamente como um MPO. Eles decompõem o Hamiltoniano dependente do tempo em uma soma de operadores locais independentes do tempo $H^{(a)}$ ponderados por funções de acionamento $f_a(t)$.
   • Reconexão (Rewiring): Uma inovação chave é a "reconexão" da representação MPO de $H(t)$. Em vez de uma única estrutura triangular superior, os autores introduzem níveis virtuais distintos (rotulados $3a$) para cada canal de acionamento $a$. Isso permite que o MPO distinga entre diferentes permutações das funções de acionamento.
   • Transformação em MPO Extensivo: Similar à construção da série de Taylor, os autores transformam o MPO de primeiro grau (que escala linearmente com o tamanho do sistema) em um MPO extensivo de tamanho. Isso envolve redirecionar setas dos novos níveis $3a$ de volta ao nível base (1) e multiplicar pelos integrais ordenados no tempo apropriados das funções de acionamento, denotados como $[f_{a_1} \dots f_{a_k}]$.
   • Construção de Ordem Superior: Um procedimento algorítmico (Algoritmo 2) é fornecido para construir MPOs de Dyson de ordem $N$ identificando e somando contribuições de produtos não-disjuntos de termos do Hamiltoniano, ponderados pelos integrais ordenados no tempo corretos.

Técnicas de Compressão
Para gerenciar a dimensão de ligação que cresce rapidamente inerente a expansões de alta ordem, o artigo detalha duas estratégias de compressão:

• Compressão Exata de Colunas: Níveis no MPO que compartilham a mesma "história" de aplicação de operadores (ou seja, seus rótulos são idênticos após remover os índices '1' representando termos não iniciados) são identificados como equivalentes. Suas linhas correspondentes são somadas e colunas redundantes são removidas.
• Compressão Aproximada de Linhas: Esta técnica identifica níveis que são seguidos pelas mesmas cadeias de operadores, mas com coeficientes diferentes. Ao construir uma base para o espaço de operadores da metade direita e realizar uma decomposição QR reveladora de posto, níveis que não introduzem operadores genuinamente novos (até a ordem $N$) são expressos como combinações lineares dos níveis mantidos. Esta etapa é aproximada, mas introduz erros apenas em ordens superiores a $N$, preservando a precisão da expansão de ordem $N$ enquanto reduz significativamente a dimensão de ligação.

Resultados e Benchmarks
Os autores avaliam a construção de MPO de Dyson em cadeias de spin-1/2 com Hamiltonianos modulados no tempo:

• Sistemas Finitos: Comparações contra simulações quase-exatas (usando integração de Verner de alta ordem) em um sistema de 8 sítios mostram que o erro escala como $O(dt^N)$ para um MPO de Dyson de ordem $N$, confirmando a ordem de convergência teórica.
• Sistemas Infinitos: Simulações em uma cadeia XXZ de Heisenberg infinita com anisotropia modulada no tempo foram comparadas contra simulações TDVP padrão com passos de tempo extremamente pequenos ($dt=10^{-7}$). A abordagem MPO de Dyson alcançou precisão comparável com passos de tempo significativamente maiores.
• Eficiência de Tempo de Execução: Uma análise de compromisso entre tempo de execução e precisão demonstra que a precisão aprimorada por passo de tempo se traduz em um custo computacional mais favorável. Para uma precisão-alvo fixa, o método MPO de Dyson requer menos tempo de execução do que as abordagens TDVP padrão.

Significado e Alegações
O artigo alega fornecer uma estrutura sistemática para simular dinâmica quântica dependente do tempo que é:

• Arbitrariamente Precisa: A ordem de truncamento pode ser aumentada para reduzir a escala de erro com o passo de tempo.
• Extensiva de Tamanho: A construção é válida diretamente no limite termodinâmico, ao contrário de adições ingênuas de potências do Hamiltoniano.
• Versátil: Lida com interações de curto e longo alcance e aplica-se tanto a sistemas finitos quanto infinitos.
• Aplicável à Simulação Quântica: Os autores sugerem que o método oferece uma estrutura natural para benchmarking e otimização de circuitos quânticos para Hamiltonianos dependentes do tempo, potencialmente guiando a compilação de circuitos eficientes.

O trabalho posiciona a série de Dyson MPO como uma formulação mais natural e eficiente para sistemas dependentes do tempo em comparação com a abordagem de dois passos de Magnus, embora esta última permaneça uma ferramenta poderosa para estudar dinâmica de Floquet. Os autores observam que trabalhos futuros poderiam explorar generalizações para dimensões superiores via Operadores de Pares Entrelaçados Projetados (PEPOs) e aplicações a sistemas periodicamente acionados.