Matrix Product Operator Encodings of the Magnus Expansion and Dyson Series

Sintesi Tecnica: Codifiche di Operatori Prodotto di Matrice per lo Sviluppo di Magnus e la Serie di Dyson

Enunciazione del Problema
La simulazione dell'evoluzione temporale reale di modelli reticolari quantistici unidimensionali governati da Hamiltoniane dipendenti dal tempo, $H(t)$, rimane una sfida significativa nella fisica computazionale dei molti corpi. Sebbene l'operatore di evoluzione temporale $U(t, t_0)$ sia formalmente definito tramite un esponenziale ordinato nel tempo, la sua valutazione analitica è generalmente impossibile per sistemi interagenti. Gli approcci numerici standard spesso si basano sulla suddivisione dell'evoluzione temporale in piccoli passi in cui l'Hamiltoniana è trattata come costante (decomposizione di Trotter-Suzuki) o sull'approssimazione dell'operatore di evoluzione tramite metodi di rete tensoriale come il Principio Variazionale Dipendente dal Tempo (TDVP) o approcci basati su Krylov. Tuttavia, per Hamiltoniane che fluttuano rapidamente, questi metodi richiedono passi temporali estremamente piccoli, portando a costi computazionali elevati. Inoltre, i metodi esistenti di rete tensoriale ad alto ordine, come quelli basati sulla serie di Taylor per Hamiltoniane indipendenti dal tempo, non si generalizzano in modo diretto ai casi dipendenti dal tempo senza incorrere in dimensioni di legame proibitive o perdere l'estensività dimensionale nel limite termodinamico.

Metodologia
Gli autori introducono una codifica tramite Operatore Prodotto di Matrice (MPO) sia per lo sviluppo di Magnus che per la serie di Dyson, specificamente adattata per Hamiltoniane dipendenti dal tempo. Il lavoro si fonda su ricerche precedenti (Rif. 22) che hanno codificato con successo la serie di Taylor di Hamiltoniane indipendenti dal tempo come MPO estensivi dimensionalmente.

La metodologia procede lungo due percorsi principali:

1. Sviluppo di Magnus come MPO: L'operatore di evoluzione temporale è espresso come $U(t, t_0) = e^{\Omega(t, t_0)}$, dove $\Omega$ è sviluppato in commutatori dell'Hamiltoniana a istanti temporali diversi. Gli autori dimostrano che i commutatori di MPO di primo grado (che rappresentano Hamiltoniane locali) possono essere essi stessi costruiti come MPO di primo grado. Ciò permette di costruire l'operatore di Magnus $\Omega$ come combinazione lineare di MPO di primo grado. Successivamente, l'esponenziale $e^\Omega$ è calcolato utilizzando la costruzione MPO della serie di Taylor consolidata.
2. Serie di Dyson come MPO: Gli autori propongono un approccio più diretto codificando la serie di Dyson, $U(t, t_0) = \mathcal{T} \exp(-i \int H(t') dt')$, direttamente come MPO. Decompongono l'Hamiltoniana dipendente dal tempo in una somma di operatori locali indipendenti dal tempo $H^{(a)}$ ponderati da funzioni di guida $f_a(t)$.
   • Ricollegamento (Rewiring): Un'innovazione chiave è il "ricollegamento" della rappresentazione MPO di $H(t)$. Invece di una singola struttura triangolare superiore, gli autori introducono livelli virtuali distinti (etichettati $3a$) per ciascun canale di guida $a$. Ciò permette all'MPO di distinguere tra diverse permutazioni delle funzioni di guida.
   • Trasformazione in MPO Estensivo: Analogamente alla costruzione della serie di Taylor, gli autori trasformano l'MPO di primo grado (che scala linearmente con la dimensione del sistema) in un MPO estensivo dimensionalmente. Ciò comporta il reindirizzamento delle frecce dai nuovi livelli $3a$ al livello base (1) e la moltiplicazione per gli opportuni integrali ordinati nel tempo delle funzioni di guida, indicati come $[f_{a_1} \dots f_{a_k}]$.
   • Costruzione di Ordine Superiore: Viene fornito una procedura algoritmica (Algoritmo 2) per costruire MPO di Dyson di ordine $N$ identificando e sommando i contributi dai prodotti non disgiunti di termini dell'Hamiltoniana, ponderati dagli opportuni integrali ordinati nel tempo.

Tecniche di Compressione
Per gestire la dimensione di legame in rapida crescita intrinseca negli sviluppi ad alto ordine, il documento dettaglia due strategie di compressione:

• Compressione Esatta delle Colonne: I livelli nell'MPO che condividono la stessa "storia" di applicazione degli operatori (cioè, le loro etichette sono identiche dopo aver rimosso gli indici '1' che rappresentano termini non avviati) sono identificati come equivalenti. Le loro righe corrispondenti vengono sommate e le colonne ridondanti vengono rimosse.
• Compressione Approssimata delle Righe: Questa tecnica identifica i livelli che sono seguiti dalle stesse stringhe di operatori ma con coefficienti diversi. Costruendo una base per lo spazio degli operatori della metà destra ed eseguendo una decomposizione QR rivelatrice del rango, i livelli che non introducono operatori genuinamente nuovi (fino all'ordine $N$) sono espressi come combinazioni lineari dei livelli mantenuti. Questo passaggio è approssimato ma introduce errori solo a ordini superiori a $N$, preservando l'accuratezza dello sviluppo di ordine $N$ mentre riduce significativamente la dimensione di legame.

Risultati e Benchmark
Gli autori sottopongono a benchmark la costruzione MPO di Dyson su catene di spin-1/2 con Hamiltoniane modulate nel tempo:

• Sistemi Finiti: I confronti con simulazioni quasi-esatte (utilizzando l'integrazione di Verner ad alto ordine) su un sistema di 8 siti mostrano che l'errore scala come $O(dt^N)$ per un MPO di Dyson di ordine $N$, confermando l'ordine di convergenza teorico.
• Sistemi Infiniti: Le simulazioni su una catena XXZ di Heisenberg infinita con anisotropia modulata nel tempo sono state confrontate con simulazioni TDVP standard con passi temporali estremamente piccoli ($dt=10^{-7}$). L'approccio MPO di Dyson ha raggiunto un'accuratezza comparabile con passi temporali significativamente più grandi.
• Efficienza di Esecuzione: Un'analisi del compromesso tempo di esecuzione-accuratezza dimostra che l'accuratezza migliorata per passo temporale si traduce in un costo computazionale più favorevole. Per un'accuratezza target fissa, il metodo MPO di Dyson richiede meno tempo di esecuzione rispetto agli approcci TDVP standard.

Significato e Affermazioni
Il documento afferma di fornire un quadro sistematico per la simulazione della dinamica quantistica dipendente dal tempo che è:

• Arbitrariamente Accurato: L'ordine di troncamento può essere aumentato per ridurre la scala dell'errore con il passo temporale.
• Estensivo Dimensionalmente: La costruzione è valida direttamente nel limite termodinamico, a differenza di aggiunte naive di potenze dell'Hamiltoniana.
• Versatile: Gestisce sia interazioni a corto che a lungo raggio e si applica sia a sistemi finiti che infiniti.
• Applicabile alla Simulazione Quantistica: Gli autori suggeriscono che il metodo offre un quadro naturale per il benchmarking e l'ottimizzazione di circuiti quantistici per Hamiltoniane dipendenti dal tempo, potenzialmente guidando la compilazione di circuiti efficienti.

Il lavoro posiziona la serie di Dyson MPO come una formulazione più naturale ed efficiente per sistemi dipendenti dal tempo rispetto all'approccio di Magnus in due passaggi, sebbene quest'ultimo rimanga uno strumento potente per lo studio della dinamica di Floquet. Gli autori notano che lavori futuri potrebbero esplorare generalizzazioni a dimensioni superiori tramite Operatori di Coppia di Entanglement Proiettati (PEPO) e applicazioni a sistemi guidati periodicamente.