Hadamard regularization of open quantum systems coupled to unstructured environments in the Schwinger-Keldysh formalism

Dit paper introduceert een tijdstap-algoritme voor de Kadanoff-Baym-vergelijkingen in het Schwinger-Keldysh-formalisme, gebaseerd op Hadamard-regularisatie, dat de numerieke schaalproblemen bij gescheiden tijdschalen oplost door zowel niet-Markovse effecten bij lage temperaturen als renormalisatie door een snelle omgeving nauwkeurig te beschrijven.

De kern

Hoe je een trillende veer in een storm kunt simuleren zonder je computer te laten ontploffen

Stel je voor dat je een heel klein, kwetsbaar object hebt: een quantum-massa aan een veer (een harmonische oscillator). Dit object zit in een kamer die vol zit met onzichtbare, razendsnelle deeltjes (het "milieu" of de "omgeving"). Deze deeltjes botsen voortdurend tegen je veer, waardoor deze gaat trillen, afkoelen en uiteindelijk tot rust komt.

In de wereld van de quantumfysica willen wetenschappers precies berekenen hoe deze veer zich gedraagt. Maar hier zit een groot probleem: de deeltjes in de kamer bewegen zo snel dat ze een "ruis" veroorzaken die veel sneller is dan de beweging van je veer zelf.

Het probleem: De "Stijve" Rekenmachine
Om dit te simuleren, gebruiken wetenschappers een geavanceerde wiskundige methode (de Schwinger-Keldysh theorie). Het probleem is dat deze methode normaal gesproken vereist dat je de tijd in heel, heel kleine stukjes verdeelt. Je moet elke botsing van die snelle deeltjes vastleggen.

Stel je voor dat je een video wilt maken van een vlinder die vliegt, maar de camera moet ook elke microscopische trilling van het stofje in de lucht vastleggen. Omdat de trillingen zo snel zijn, moet je de video met een extreem hoge snelheid opnemen (bijvoorbeeld 1 biljoen beelden per seconde). Zelfs de snelste supercomputer zou hierdoor vastlopen; de rekentijd zou exponentieel toenemen. Dit noemen we een "stijf" probleem.

De oplossing: Hadamard-regulering (De "Scheermes-methode")
Jakob Dolgner, de auteur van dit paper, heeft een slimme truc bedacht. In plaats van elke snelle botsing één voor één te tellen, kijkt hij naar het gemiddelde effect van die botsingen.

Hij gebruikt een wiskundig concept dat Hadamard-regulering heet. Je kunt dit vergelijken met het gebruik van een scheermes om een oneindig scherpe punt af te snijden.

1. De oneindige piek: De wiskunde zegt dat op het exacte moment van botsing de krachten oneindig groot worden (een "divergentie").
2. Het afsnijden: In plaats van te proberen die oneindige piek op te lossen, snijdt hij het "oneindige" deel er wiskundig af. Hij behandelt dit deel als een bekende, snelle "kick" die hij apart berekent.
3. De rest: Wat overblijft is een rustig, langzaam veranderend effect dat je met grote stappen kunt simuleren.

De Analogie: De Trillende Veer en de Storm
Laten we de analogie verder uitwerken:

• De Veer: Je quantum-systeem (bijvoorbeeld een elektron of een atoom).
• De Storm: De omgeving met de snelle deeltjes.
• De Oude Methode: Je probeert elke regendruppel in de storm te tellen om te zien hoe de veer beweegt. Je moet je computer laten rennen met de snelheid van de druppels.
• De Nieuwe Methode: Je zegt: "Oké, de storm is zo snel dat ik niet elke druppel kan zien. Maar ik weet wiskundig precies hoe de gemiddelde windkracht werkt."
  • Hij berekent de "stoot" van de wind (de snelle, lokale effecten) apart.
  • Hij berekent de "duw" van de wind (de langzame, historische effecten) met grote stappen.

Waarom is dit belangrijk?
Met deze nieuwe methode kunnen wetenschappers nu simuleren hoe quantum-systemen werken in situaties waar de omgeving veel sneller is dan het systeem zelf, zonder dat hun computers vastlopen.

• Geen meer "stijve" rekentijd: Je hoeft niet meer te wachten tot je computer de snelle deeltjes heeft geteld. Je kunt de simulatie draaien op de snelheid van de veer zelf.
• Koudere temperaturen: Het werkt zelfs als het heel koud is (waar quantum-effecten belangrijk zijn), iets wat eerdere methoden vaak niet goed konden.
• Nauwkeurigheid: Het houdt rekening met het feit dat de omgeving het systeem iets verandert (zoals een "Lamb-shift", een kleine verandering in de frequentie van de veer), maar doet dit op een slimme manier zonder de rekenkracht te verspillen.

Conclusie
Dit paper is als het vinden van een nieuwe manier om een storm te voorspellen. In plaats van elke windvlaag te meten, begrijp je de natuurwetten van de storm zo goed dat je de gevolgen voor je huis (het quantum-systeem) precies kunt voorspellen, zelfs als je de storm zelf niet in detail kunt zien. Dit opent de deur voor betere quantumcomputers en nieuwe materialen, omdat we nu veel complexere situaties kunnen simuleren dan voorheen mogelijk was.

Technische samenvatting

Probleemstelling

Het artikel adresseert een fundamentele numerieke uitdaging bij het modelleren van open kwantumsystemen die gekoppeld zijn aan een "ongestructureerde" (featureless) omgeving, zoals een Ohmse bad. Hoewel de Schwinger-Keldysh-veldtheorie een krachtig raamwerk biedt voor niet-evenwichtsdynamica, stuiten numerieke oplossingen van de Kadanoff-Baym-vergelijkingen (KBE) op een ernstige beperking: de rekentijd schaalt kubisch met het aantal tijdstappen.

Dit probleem wordt verergerd door het bestaan van een sterk gescheiden tijdschaal:

1. Systeemschaal: De karakteristieke frequentie van het systeem ($\omega_0$) en de demping ($\gamma$).
2. Omgevingschaal: De frequentiecut-off van het bad ($\omega_c$), die vaak veel groter is dan $\omega_0$ ($\omega_c \gg \omega_0$).

In de standaard benadering moet de tijdsstap $\Delta t$ klein genoeg zijn om de snelle dynamica van het bad ($1/\omega_c$) op te lossen. Dit leidt tot een "stijf" (stiff) probleem waarbij de rekentijd onpraktisch groot wordt. Bovendien leiden onbeperkte spectra (zoals in de wide-band limiet) tot divergenties in de impulsvariatie (logaritmische divergentie) en frequentie-renormalisatie, wat een regulatie via een cut-off $\omega_c$ vereist. De uitdaging is om een methode te vinden die de $\omega_c \to \infty$ limiet kan nemen zonder de kubische schaling of de noodzaak van een extreem kleine tijdsstap.

Methodologie

De auteur gebruikt een gedempte kwantumharmonische oscillator (DQHO) als testmodel binnen het Schwinger-Keldysh-formalisme. De kern van de methode is een splitsing van schalen gebaseerd op Hadamard-regulering.

1. Identificatie van Singulariteiten:
   De symmetrische component van de inbeddings-self-energie ($\Sigma^S_{em}$), die verantwoordelijk is voor fluctuaties, vertoont een singulariteit op de tijddiagonaal ($t_1 = t_2$). In de limiet $\omega_c \to \infty$ gedraagt deze zich als een Hadamard eindig-deel (HFP) distributie (een "nascent" HFP), specifiek gerelateerd aan $1/(t'-t)^2$.

2. Taylor-expansie en Substitutie:
   De convolutie-integrale in de KBE's worden geanalyseerd rond het kritieke punt $t' = t_1$. De correlator $G(t', t_2)$ wordt ontwikkeld in een Taylorreeks rond dit punt:
   $$G(t', t_2) = G(t_1, t_2) + \partial_{t_1}G(t_1, t_2)(t' - t_1) + R_2(t_1, t', t_2)$$
   Hierbij is $R_2$ de restterm die kwadratisch schaal met $(t'-t_1)$.

3. Splitsing in Snelle en Trage Termen:
   De convolutie wordt opgesplitst in:

   • Lokale, snelle termen: De nulde- en eerste-orde termen van de expansie. Deze bevatten de $\omega_c$-afhankelijkheid en de divergenties (massa-renormalisatie en impulsvariatie). Deze kunnen analytisch worden geïntegreerd en leiden tot lokale termen zoals $\gamma \partial_t G$ en een logaritmische "kick" op de diagonaal.
   • Restterm (Trage geschiedenis): De term met $R_2$. Omdat $R_2 \sim (t'-t)^2$, wordt de singulariteit van de self-energie ($\sim 1/(t'-t)^2$) opgeheven. Deze term is goed geconditioneerd, schaalbaar met de systeemtijdschaal en vertoont exponentieel afnemend geheugen.

4. Numerieke Implementatie:

   • Een tijdstap-algoritme wordt ontwikkeld dat de lokale stijfheid (de snelle termen) exact integreert via Filon-type kwadratuur.
   • De trage restterm wordt geïntegreerd met een grovere tijdsstap ($\Delta t \sim 1/\omega_0$).
   • Door de singulariteit te regulariseren via de Hadamard-regulering, kan de limiet $\omega_c \to \infty$ analytisch worden genomen voor alle tijden behalve op de diagonaal, waar de logaritmische divergentie wordt "renormaliseerd" door de discretisatie van het rooster.

Belangrijkste Bijdragen

• Hadamard-Regularisatie in KBE's: De paper introduceert een rigoureuze wiskundige behandeling van de divergenties in de Kadanoff-Baym-vergelijkingen door de self-energie te interpreteren als een Hadamard eindig-deel distributie. Dit lost het probleem van de logaritmische divergentie in de impulsvariatie op zonder kunstmatige cut-offs.
• Efficiënt Tijdstap-algoritme: Er wordt een nieuw algoritme gepresenteerd dat de kubische schaling van de KBE's doorbreekt voor ongestructureerde omgevingen. Door de snelle dynamica van het bad te analytisch te behandelen en alleen de trage geschiedenis numeriek op te lossen, wordt de rekentijd gereduceerd tot een schaling die afhankelijk is van de systeemfrequentie, niet van de bad-cut-off.
• Overbrugging van Schalen: De methode maakt het mogelijk om systemen te simuleren met een extreme scheiding van schalen ($\omega_c / \omega_0 \gg 1$) zonder dat de tijdsstap verkleind hoeft te worden, wat essentieel is voor het modelleren van koude, coherentie-rijke omgevingen.

Resultaten

• Thermalisatie: Het algoritme convergeert correct naar de thermische evenwichtstoestand voor een breed scala aan parameters (temperatuur $T$, demping $\gamma$, en frequentie $\omega_0$).
• Niet-Markoviaanse Dynamica: In het ultra-koude regime ($T \ll \gamma$) toont de simulatie een algebraïsche verval ($t^{-2}$) van de symmetrische correlator, een fenomeen dat door Markoviaanse benaderingen (zoals de Lindblad-vergelijking) niet kan worden vastgelegd.
• Renormalisatie van Divergenties: De numerieke resultaten tonen aan dat de impulsvariatie $\langle \pi^2 \rangle$ logaritmisch schaalt met $\omega_c$ zoals theoretisch voorspeld, maar dat de dynamica op de systeemtijdschaal onafhankelijk is van de exacte waarde van $\omega_c$ zolang deze groot genoeg is. De "coarse-graining" door de tijdsstap fungeert als een natuurlijke renormalisatie.
• Validatie: De resultaten worden vergeleken met exacte analytische oplossingen in het evenwicht en met Lindblad-benaderingen, waarbij de superioriteit van de KBE-benadering voor niet-Markoviaanse en koude regimes wordt bevestigd.

Betekenis en Toekomstperspectief

Deze studie biedt een cruciale doorbraak voor de numerieke behandeling van open kwantumsystemen in de Schwinger-Keldysh-formalisme. Het maakt het mogelijk om:

1. Realistische omgevingen te modelleren: Systemen met een breed bandbreedte-bad (zoals elektronen in een metaal of fononen) kunnen worden gesimuleerd zonder de onpraktische tijdsstap-beperkingen.
2. Koude en Coherente Regimes: Het is nu haalbaar om de overgang van kwantum naar klassiek gedrag in zeer koude, sterk gedempte systemen te bestuderen, waar niet-Markoviaanse effecten dominant zijn.
3. Uitbreiding: De methode is direct toepasbaar op fermionische systemen (waar de wide-band limiet ook leidt tot logaritmische divergenties) en kan worden uitgebreid naar andere bad-spectra (sub-Ohmisch en super-Ohmisch), hoewel voor super-Ohmisch gedrag met $s \geq 2$ de divergentie niet-integreerbaar blijft en een fysieke cut-off nodig blijft.

Kortom, de paper levert een theoretisch en numeriek raamwerk dat de kloof overbrugt tussen de exacte veldtheoretische beschrijvingen en de praktische berekenbaarheid van complexe, open kwantumsystemen.