Markovian Embeddings of Non-Markovian Open System Dynamics

Sintesi Tecnica: Embedding Markoviani di Dinamiche di Sistemi Aperti Non-Markoviani

Definizione del Problema
La simulazione di sistemi quantistici aperti non-markoviani è una sfida centrale nell'ottica quantistica, nella fisica della materia condensata e nella fisica chimica. Sebbene esistano metodi non perturbativi come il Monte Carlo con integrali di cammino e le reti tensoriali, essi spesso si affidano a effetti ambientali impliciti tramite funzionali di influenza. Alternativamente, gli approcci tempo-locali (ad esempio, l'equazione gerarchica di moto [HEOM], il formalismo dei pseudomodi di Lindblad) rappresentano esplicitamente l'ambiente tramite modi ausiliari discreti. Tuttavia, questi metodi esistenti sembrano distinti nella loro costruzione matematica, nell'inizializzazione degli stati ausiliari (vuoto vs termico) e nelle procedure per ottenere l'operatore di densità ridotto (traccia parziale vs proiezione). Inoltre, le implementazioni pratiche di tali metodi soffrono spesso di instabilità numerica quando la base dei modi ausiliari viene troncata, in particolare durante evoluzioni a lungo termine o in regimi di accoppiamento forte. Le connessioni fisiche tra queste formulazioni apparentemente disparate e l'origine delle loro instabilità numeriche richiedono un chiarimento teorico unificato.

Metodologia
Gli autori sviluppano un quadro unificato basato sulla rappresentazione per integrale di cammino del funzionale di influenza, estendendo il concetto precedentemente introdotto di "Dissipazione Quantistica in Spazio degli Stati Minimalmente Esteso" (QD-MESS). La metodologia centrale prevede:

1. Unraveling Gaussiano: Il kernel dell'auto-energia del bagno $\Sigma(t)$, che caratterizza la memoria non-markoviana, viene fattorizzato in un prodotto di matrici costanti e una matrice di funzione Green tempo-dipendente $G(t)$ (ovvero, $\Sigma(t) = U^\dagger G(t) V$).
2. Introduzione di Campi Ausiliari: Utilizzando un'identità gaussiana, la fase di influenza non locale viene trasformata in un'azione tempo-locale che coinvolge campi ausiliari deterministici $\Psi(t)$. La dimensione e la struttura di questi campi sono dettate dalla scelta della fattorizzazione di $G(t)$.
3. Decomposizione della Funzione Green: Gli autori dimostrano che diverse scelte della struttura della funzione Green (Keldysh, a blocchi diagonali o solo ritardata) portano a diversi "unraveling". Queste scelte corrispondono a diversi "unraveling" dell'auto-energia del bagno.
4. Trasformazioni di Bogoliubov: Il lavoro stabilisce che la libertà di gauge nella fattorizzazione (Eq. 14) corrisponde a trasformazioni lineari nello spazio dei cammini dei modi ausiliari. Queste trasformazioni sono identificate come trasformazioni di Bogoliubov, che mappano tra diversi frame di rappresentazione (ad esempio, tra i formalismi HEOM e dei pseudomodi di Lindblad).

Contributi Chiave e Risultati

• Derivazione Unificata dei Metodi Esistenti: Il quadro dimostra che l'equazione gerarchica di moto (HEOM) e il formalismo dei pseudomodi di Lindblad emergono naturalmente da specifiche scelte di strutture della funzione Green e matrici di fattorizzazione.

  • Rappresentazione Keldysh: L'uso di una matrice della funzione Green di Keldysh porta a una formulazione in cui i modi ausiliari sono inizializzati in uno stato termico (occupazione $n_{ref}$), e l'operatore di densità ridotto è ottenuto tramite una traccia parziale. Ciò recupera le equazioni standard dei pseudomodi di Lindblad quando vengono scelti parametri specifici.
  • Rappresentazione Pure-State (a Blocchi Diagonali): L'uso di una funzione Green a blocchi diagonali decoppia i contorni forward e backward. Ciò porta a una formulazione in cui i modi ausiliari sono inizializzati nello stato di vuoto e l'operatore di densità ridotto è ottenuto tramite una proiezione sul vuoto. Questa struttura produce naturalmente la gerarchia HEOM.
  • Rappresentazione dello Spazio di Hilbert: Una formulazione che utilizza solo una funzione Green ritardata permette una rappresentazione in cui il generatore totale agisce come un Liouvilliano nello spazio del sistema e come un Hamiltoniano nello spazio ausiliario.

• Chiarimento di Inizializzazione e Proiezione: Il lavoro chiarisce che l'inizializzazione dei modi ausiliari (vuoto vs termico) e il metodo di estrazione dello stato ridotto (proiezione vs traccia parziale) non sono arbitrari, ma sono determinati dalle condizioni al contorno imposte dalla specifica struttura della funzione Green utilizzata nell'unraveling.

• Risoluzione dell'Instabilità Numerica: Un risultato critico è la dimostrazione che le instabilità numeriche osservate nelle formulazioni HEOM troncate (particolarmente per modi discreti o forte smorzamento) sono dipendenti dalla rappresentazione. Gli autori mostrano che l'applicazione di specifiche trasformazioni di Bogoliubov (trasformazioni di similitudine sul Liouvilliano) può mappare una rappresentazione instabile in una stabile. Ad esempio, trasformare le equazioni HEOM derivate dalla rappresentazione pure-state ripristina la stabilità numerica per evoluzioni a lungo termine, un risultato verificato rispetto alla letteratura precedente.

• Relazione con la Termodinamica di Campo (Thermofield Dynamics): Il lavoro collega esplicitamente la trasformazione di termofield (usata per mappare ambienti a temperatura finita in spazi di Hilbert ampliati) alle trasformazioni di Bogoliubov inerenti alla fattorizzazione dell'auto-energia.

Significatività e Rivendicazioni
Il lavoro sostiene di fornire una "fondazione teorica trasparente" per le tecniche di embedding. Rivelando che distinti metodi di simulazione non-markoviana sono correlati attraverso trasformazioni lineari (trasformazioni di Bogoliubov) nello spazio dei cammini ausiliari, il lavoro offre uno strumento sistematico ai ricercatori per:

1. Selezionare rappresentazioni ottimali basate sull'intuizione fisica o sull'efficienza della simulazione.
2. Comprendere l'origine delle instabilità numeriche nelle simulazioni con base troncata.
3. Sviluppare nuovi metodi numericamente stabili scegliendo frame di trasformazione appropriati che ridistribuiscono gli errori di troncamento.

Gli autori sottolineano che, sebbene la dinamica fisica esatta sia invariante sotto queste trasformazioni (poiché corrispondono a trasformazioni di similitudine del Liouvilliano completo), la dinamica approssimata ottenuta dopo il troncamento della base è altamente sensibile alla rappresentazione scelta. Pertesi, la capacità di navigare tra questi frame è presentata come un grado di libertà cruciale per ottimizzare la simulazione di sistemi quantistici aperti fortemente accoppiati. Il lavoro non propone nuovi setup sperimentali, ma mira a raffinare lo strumento teorico e computazionale per la dinamica quantistica non-perturbativa.