Variational non-gaussian approach to interacting spin-boson models

Technische Samenvatting: Variational Non-Gaussian Benadering van Interagerende Spin-Boson Modellen

Probleemstelling
Interagerende spin-boson modellen, waarbij discrete kwantum vrijheidsgraden (spins) koppelen aan één of meer bosonische modi (fotonen of fononen), staan centraal in de hedendaagse kwantumwetenschap, waaronder cavity QED, circuit QED en spin-fonon systemen. Een primaire uitdaging bij het simuleren van deze systemen komt voort uit de oneindig-dimensionale aard van de bosonische Hilbertruimte in combinatie met sterke veel-deeltjes correlaties in de spinsector. Conventionele numerieke methoden kampen met een afruil: het afkappen van de bosonische bezettingsbasis introduceert bias, terwijl het versoepelen van de spinsector naar mean-field beschrijvingen essentieel spin-boson verstrengeling niet kan vangen. Bovendien lijden directe tensor netwerk benaderingen (zoals DMRG) toegepast op het volledige samengestelde systeem vaak aan een snelle groei van de bond dimensie door de noodzaak om het bosonische veld expliciet naast de spins te representeren.

Methodologie
De auteurs stellen een hybride variationeel framework voor die de optimalisatie van de bosonische sector ontkoppelt van de spinsector. De grondtoestand wordt gerepresenteerd binnen een compact non-Gaussian variationeel manifold gedefinieerd door de ansatz:
$$|\psi_{\text{NGS}}\rangle = U_\lambda (U_{\text{GS}} |0_b\rangle \otimes |\phi_s\rangle)$$
waarbij:

1. Bosonic Gaussian Layer ($U_{\text{GS}}$): Een Gaussische transformatie bestaande uit verschuivings- ($\Delta_R$) en squeezing- ($\xi$) operatoren die werken op de bosonische vacuüm $|0_b\rangle$. Dit vangt de bosonische fluctuaties en mean-field verschuivingen op.
2. Dressing Transformation ($U_\lambda$): Een unitaire operator die spin-boson correlaties introduceert die verder gaan dan het mean-field niveau. Deze werkt als een spin-afhankelijke verschuiving (polaron-achtige transformatie) gecontroleerd door een variationele parameter $\lambda$. Deze laag brengt de lineaire spin-boson koppeling in kaart naar effectieve spin-spin interacties.
3. Effectieve Spin Sector ($|\phi_s\rangle$): De spin toestand is geen eenvoudige producttoestand maar de exacte grondtoestand van een effectieve spin Hamiltonian ($H_{\text{eff}}$). Deze Hamiltonian wordt verkregen door de volledige Hamiltonian te middelen over de bosonische Gaussian toestand en de dressing transformatie.

De methode maakt gebruik van een zelfconsistente optimalisatie-loop:

• Voor vaste bosonische parameters ($\Delta_R, \xi, \lambda$), wordt de effectieve spin Hamiltonian opgelost met behulp van DMRG (Density Matrix Renormalization Group) om $|\phi_s\rangle$ en spin verwachtingswaarden te verkrijgen.
• Deze spin observables worden gebruikt om de bosonische variationele parameters bij te werken door de totale variationele energie functional te minimaliseren.
• De loop iterereert tot convergentie van de energie en geselecteerde observables.

Deze benadering vermijdt expliciete afkapping van het fotonische veld. In plaats daarvan worden de bosonische vrijheidsgraden variationeel weggeïntegreerd, waardoor het probleem wordt gemapt naar een zelfconsistent spin model.

Belangrijkste Bijdragen en Resultaten
De paper benchmarkt dit NGS-DMRG framework tegen de Dicke en Dicke-Ising modellen, waarbij de resultaten worden vergeleken met geconvergeerde referentie DMRG simulaties (die een expliciete foton cutoff gebruiken).

1. Nauwkeurigheid en Variationele Verbetering:

   • De non-Gaussian benadering levert consistent lagere grondtoestandsenergieën op dan de Gaussian limiet (mean-field), wat de variationele verbetering bevestigt.
   • Voor het Dicke model reproduceert de methode accuraat de superradiante kwantumfaseovergang (QPT) en finite-size scaling, waarbij de resultaten overeenkomen met referentie DMRG resultaten met aanzienlijk minder computationele middelen.
   • In het Dicke-Ising model brengt het framework succesvol een complex fasediagram in kaart met:
     • Continue (tweede-orde) transities tussen normale en superradiante fasen.
     • Discontinue (eerste-orde) transities en tricritische punten in het ferromagnetische regime.
     • Een smalle intermediaire fase van coëxisterende antiferromagnetische orde en superradiantie (AFM-SP) in het antiferromagnetische regime, een kenmerk dat voorspeld wordt door mean-field theorie maar moeilijk numeriek te resolveren is.
   • De methode vangt finite-size fluctuaties op en convergeert naar de thermodynamische limiet, waar spin-foton verstrengeling verdwijnt, consistent met gevestigde analytische limieten.

2. Computationele Efficiëntie:

   • Bond Dimensie Reductie: Door bosonen variationeel te behandelen, opereert de DMRG solver op een effectieve spin-only Hamiltonian. Dit resulteert in een aanzienlijk lagere vereiste MPS bond dimensie ($\chi$) vergeleken met referentie DMRG, die het samengestelde spin-boson keten expliciet moet behandelen.
   • Runtime: Ondanks de overhead van de zelfconsistente loop, leidt de verminderde tensor netwerk complexiteit tot een substantiële afname in totale executietijd. De buitenste optimalisatie-loop vertoont exponentiële convergentie.
   • Het voordeel is het meest uitgesproken in regimes met sterke correlaties en complexe fasegrenzen (bijv. de coëxistentie fase), waar referentie DMRG zeer grote bond dimensies en foton cutoffs vereist.

Significantie
De paper claimt dat dit hybride framework een onderscheidend algoritmisch pad biedt vergeleken met local basis optimization (LBO) tensor netwerk methoden. Door fotonische modi te vangen via variationele parameters (displacement, squeezing, dressing) in plaats van iteratieve basis herconfiguratie, stelt de methode de DMRG solver in staat zich uitsluitend te concentreren op spin correlaties.

De significantie ligt in het bieden van een robuuste, non-perturbatieve numerieke methode die in staat is tot het bereiken van sterk gecorreleerde spin-boson fysica zonder het bosonische veld te elimineren of willekeurig af te kappen. De auteurs demonstreren dat deze benadering in staat is om complexe fase diagrammen en kritische gedragingen in modellen zoals de Dicke-Ising keten accuraat te reproduceren, waarbij een hoge getrouwheid wordt bereikt met gereduceerde computationele kosten. De auteurs bieden een open-source implementatie die arbitrary-range spin-spin koppelingen en disorder ondersteunt (beschikbaar op: https://github.com/Jpedromend/NGS), met potentiële extensies naar multi-mode systemen en hogere-dimensionale geometrieën via alternatieve solvers (bijv. PEPS, QMC).