Variational non-gaussian approach to interacting spin-boson models

Sintesi Tecnica: Approccio Variazionale Non-Gaussiano per Modelli Spin-Bosone Interagenti

Definizione del Problema
I modelli spin-bosone interagenti, in cui gradi di libertà quantistici discreti (spin) si accoppiano a uno o più modi bosonici (fotoni o fononi), sono centrali nella scienza quantistica contemporanea, inclusi la QED a cavità, la QED a circuito e i sistemi spin-fonone. Una sfida primaria nella simulazione di questi sistemi deriva dalla natura a dimensioni infinite dello spazio di Hilbert bosonico combinata con forti correlazioni many-body nel settore di spin. I metodi numerici convenzionali affrontano un compromesso: la troncatura della base di occupazione bosonica introduce un bias, mentre il rilassamento del settore di spin verso descrizioni a campo medio fallisce nel catturare l'essenziale entanglement spin-bosone. Inoltre, gli approcci diretti basati su reti tensoriali (come DMRG) applicati al sistema composito completo soffrono spesso di una rapida crescita della dimensione del legame (bond dimension) a causa della necessità di rappresentare esplicitamente il campo bosonico insieme agli spin.

Metodologia
Gli autori propongono un framework variazionale ibrido che decupla l'ottimizzazione del settore bosonico da quello di spin. Lo stato fondamentale è rappresentato all'interno di un manifold variazionale non-gaussiano compatto definito dall'ansatz:
$$|\psi_{\text{NGS}}\rangle = U_\lambda (U_{\text{GS}} |0_b\rangle \otimes |\phi_s\rangle)$$
dove:

1. Strato Bosonico Gaussiano ($U_{\text{GS}}$): Una trasformazione gaussiana che comprende operatori di spostamento ($\Delta_R$) e di squeezing ($\xi$) che agiscono sul vuoto bosonico $|0_b\rangle$. Questo cattura le fluttuazioni bosoniche e gli spostamenti a campo medio.
2. Trasformazione di Dressing ($U_\lambda$): Un operatore unitario che introduce correlazioni spin-bosone oltre il livello di campo medio. Esso agisce come uno spostamento dipendente dallo spin (trasformazione tipo polarone) controllato da un parametro variazionale $\lambda$. Questo strato mappa l'accoppiamento lineare spin-bosone in interazioni spin-spin efficaci.
3. Settore di Spin Effettivo ($|\phi_s\rangle$): Lo stato di spin non è un semplice stato prodotto, ma lo stato fondamentale esatto di un Hamiltoniano di spin efficace ($H_{\text{eff}}$). Questo Hamiltoniano è ottenuto mediando l'Hamiltoniano completo sul campo gaussiano bosonico e sulla trasformazione di dressing.

Il metodo impiega un ciclo di ottimizzazione auto-coerente:

• Per parametri bosonici fissi ($\Delta_R, \xi, \lambda$), l'Hamiltoniano di spin efficace viene risolto utilizzando DMRG (Density Matrix Renormalization Group) per ottenere $|\phi_s\rangle$ e gli osservabili di spin.
• Questi osservabili di spin vengono utilizzati per aggiornare i parametri variazionali bosonici minimizzando il funzionale dell'energia variazionale totale.
• Il ciclo itera fino alla convergenza dell'energia e di determinati osservabili.

Questo approccio evita la troncatura esplicita del campo fotonico. Invece, i gradi di libertà bosonici vengono integrati variazionalmente, mappando il problema su un modello di spin auto-coerente.

Contributi Chiave e Risultati
Il documento confronta questo framework NGS-DMRG con i modelli Dicke e Dicke-Ising, comparando i risultati con simulazioni DMRG di riferimento (che utilizzano espliciti cutoff fotonici).

1. Accuratezza e Miglioramento Variazionale:

   • L'approccio non-gaussiano produce costantemente energie dello stato fondamentale inferiori rispetto al limite gaussiano (campo medio), confermando il miglioramento variazionale.
   • Per il modello di Dicke, il metodo riproduce accuratamente la transizione di fase quantistica (QPT) superradiante e lo scaling delle dimensioni finite, corrispondendo ai risultati DMRG di riferimento con risorse computazionali significativamente ridotte.
   • Nel modello Dicke-Ising, il framework mappa con successo un diagramma di fase complesso caratterizzato da:
     • Transizioni continue (secondo ordine) tra fasi normali e superradianti.
     • Transizioni discontinue (primo ordine) e punti tricritici nel regime ferromagnetico.
     • Una stretta fase intermedia di coesistenza di ordine antiferromagnetico e superradianza (AFM-SP) nel regime antiferromagnetico, una caratteristica predetta dalla teoria di campo medio ma difficile da risolvere numericamente.
   • Il metodo cattura le fluttuazioni di dimensione finita e converge al limite termodinamico, dove l'entanglement spin-fotone svanisce, in accordo con i limiti analitici stabiliti.

2. Efficienza Computazionale:

   • Riduzione della Dimensione del Legame (Bond Dimension): Trattando i bosoni in modo variazionale, il solver DMRG opera su un Hamiltoniano efficace composto solo di spin. Ciò risulta in una dimensione del legame MPS ($\chi$) significativamente inferiore rispetto al DMRG di riferimento, che deve gestire esplicitamente la catena composita spin-bosone.
   • Tempo di Esecuzione: Nonostante l'overhead del ciclo auto-coerente, la ridotta complessità della rete tensoriale porta a una sostanziale diminuzione del tempo totale di esecuzione. Il ciclo di ottimizzazione esterno dimostra una convergenza esponenziale.
   • Il vantaggio è più pronunciato nei regimi con forti correlazioni e confini di fase complessi (ad esempio, la fase di coesistenza), dove il DMRG di riferimento richiede dimensioni del legame e cutoff fotonici molto elevati.

Significatività
Il documento sostiene che questo framework ibrido offre un percorso algoritmico distinto rispetto ai metodi di rete tensoriale a ottimizzazione della base locale (LBO). Catturando i modi fotonici attraverso parametri variazionali (spostamento, squeezing, dressing) anziché tramite una riconfigurazione iterativa della base, il metodo permette al solver DMRG di concentrarsi esclusivamente sulle correlazioni di spin.

La significatività risiede nel fornire un metodo numerico robusto e non perturbativo capace di accedere alla fisica spin-bosone fortemente correlata senza eliminare o troncare arbitrariamente il campo bosonico. Gli autori dimostrano che questo approccio può riprodurre accuratamente diagrammi di fase complessi e comportamenti critici in modelli come la catena Dicke-Ising, raggiungendo un'alta fedeltà con costi computazionali ridotti. Gli autori forniscono un'implementazione open-source del codice disponibile all'indirizzo: https://github.com/Jpedromend/NGS, che supporta accoppiamenti spin-spin a intervallo arbitrario e disordine, con potenziali estensioni a sistemi multi-modo e geometrie ad alta dimensionalità tramite altri solver (es. PEPS, QMC).