Efficient time-evolution of matrix product states using average Hamiltonians

Resumo Técnico: Evolução Temporal Eficiente de Estados de Produto de Matrizes Usando Hamiltonianos Médios

Problema
A simulação da dinâmica de sistemas de muitos corpos quânticos (QMBS) regidos por Hamiltonianos dependentes do tempo permanece um desafio significativo devido ao crescimento exponencial do espaço de Hilbert. Embora as redes de tensores, particularmente os Estados de Produto de Matrizes (MPS), ofereçam um arcabouço eficiente para sistemas unidimensionais, os algoritmos padrão de evolução temporal (como o Time-Evolving Block Decimation, TEBD, e MPO $W_{II}$) são otimizados primariamente para Hamiltonianos independentes do tempo. Quando aplicados a cenários dependentes do tempo, as implementações padrão frequentemente dependem do "passo de Riemann" (Riemann stepper), que substitui o Hamiltoniano instantâneo no início de um intervalo de tempo ($H(0)$) no operador de evolução. Esta abordagem degrada a precisão de integradores de ordem superior, reduzindo a taxa de convergência do operador de evolução de terceira ordem (para um único passo) para segunda ordem no caso dependente do tempo, limitando assim a precisidade da simulação para um dado tamanho de passo.

Metodologia
Os autores propõem uma modificação para os steppers de tempo de MPS existentes para restaurar a precisão de ordem superior para Hamiltonianos dependentes do tempo sem alterar a arquitetura fundamental das bibliotecas de redes de tensores. A metodologia central envolve a substituição do stepper de Riemann padrão por um Hamiltoniano cuidadosamente construído através de um "Hamiltoniano médio" derivado da teoria de integração numérica.

Especificamente, os autores substituem o stepper de Riemann padrão por um Hamiltoniano calculado usando a regra de Simpson:
$$H_{\text{Simp}}(\Delta t) \equiv \frac{1}{6} \left( H(0) + 4H(\Delta t/2) + H(\Delta t) \right)$$
Este Hamiltoniano médio é então alimentado em algoritmos TEBD ou MPO $W_{II}$ de segunda ordem padrão. A análise teórica demonstra que, enquanto o stepper de Riemann ($H(0)$) produz um operador de evolução com um limite de erro de $O(\Delta t^2)$, o stepper de Simpson mantém a precisão local nativa de terceira ordem ($O(\Delta t^3)$) dos integradores subjacentes. A prova baseia-se no limite da diferença entre o operador de evolução exato e o operador aproximado usando a norma Hölder-$\infty$, mostrando que a aproximação de Simpson reduz significativamente o termo de erro de integração em comparação com a aproximação de Riemann.

Principais Contribuições e Resultados
O método proposto foi validado através de simulações numéricas de uma cadeia de sítios de spin-1 representando centros de cor nitrogênio-vacância (NV) em diamante, sendo este último impulsionado por pulsos de micro-ondas. O sistema foi modelado usando um Hamiltoniano composto por termos de deriva (drift), interação e controle, com parâmetros derivados de dados experimentais.

1. Melhoria da Precisão: As simulações compararam o "Simpson stepper" contra o "Riemann stepper" padrão através de variados tamanhos de sistema (3 a 8 centros NV) e contagens de passos de integração.

   • Para sistemas pequenos (ex: 3 centros NV), o Simpson stepper reduziu o erro médio por um fator de aproximadamente 1000 para tamanhos de passo moderados em comparação com o Riemann stepper.
   • Para sistemas maiores (até 8 centros NV), o método apresentou consistentemente erros de propagação menores, com melhorias no erro médio variando de fatores de 50 a mais de 4000, dependendo da contagem de passos e do tamanho do sistema.
   • A razão de erro ($E_{\text{Riemann}} / E_{\text{Simpson}}$) permaneceu alta em todas as configurações testadas, confirmando a robustez do método.

2. Custo Computacional: O aumento da precisão vem acompanhado de um aumento modesto no tempo de execução. O Simpson stepper requer três avaliações do Hamiltoniano por passo de tempo (em $t=0$, $t=\Delta t/2$ e $t=\Delta t$) em comparação com uma para o Riemann stepper. Os autores observam que este aumento no tempo de execução é leve e atribuído unicamente às avaliações de função adicionais, tornando o compromisso (trade-off) altamente favorável para a precisão.

3. Escalabilidade: O método demonstrou melhorias de desempenho constantes à medida que o tamanho do sistema aumentava, mantendo sua vantagem sobre a abordagem padrão sem introduzir novos erros de truncamento, desde que a dimensão de ligação (bond dimension) fosse suficientemente grande para manter os erros de truncamento negligenciáveis.

Significância e Limitações
O artigo afirma que esta abordagem fornece um aumento simples, porém altamente eficiente, para os solucionadores de MPS existentes. Ao alavancar ideias clássicas de integração numérica (regra de Simpson), o método restaura as propriedades de convergência de ordem superior de integradores que, de outra forma, seriam perdidas ao simular dinâmicas dependentes do tempo. Isso é particularmente relevante para aplicações em controle quântico, espectroscopia e engenharia de Floquet, onde campos de condução dependentes do tempo são ubíquos.

Os autores explicitamente observam algumas limitações e considerações práticas:

• Requisito de Suavidade: O método depende da existência e limitação das quatro primeiras derivadas do Hamiltoniano $H(t)$ dentro de cada intervalo.
• Custo de Construção do Hamiltoniano: Se o Hamiltoniano dependente do tempo for computacionalmente caro de construir, o requisito de três avaliações por passo pode aumentar significativamente os custos computacionais.
• Limites de Emaranhamento: Embora o método melhore o erro de integração, ele não elimina as limitações impostas pelo crescimento do emaranhamento em sistemas de muitos corpos de grande escala, que são governadas pelos erros de truncamento na representação MPS.

Os autores sugerem que o conceito de média de Simpson é adaptável e poderia ser combinado com regras de quadratura de ordem superior (ex: Gauss–Legendre) ou métodos de composição de múltiplos estágios para construir Hamiltonianos efetivos de ordem ainda mais elevada em trabalhos futuros.