Green's Functions from Sample-based Krylov Quantum Diagonalization: An Impurity Solver for Dynamical Mean-Field Theory

Materie die zich vreemd en koppig gedraagt, zoals de supergeleiders die elektriciteit zonder weerstand geleiden of de magnetische materialen die onze harde schijven aandrijven, tart vaak de eenvoudige regels die we gebruiken om gewone vaste stoffen te verklaren. In deze materialen bewegen elektronen niet onafhankelijk; ze duwen en trekken aan elkaar met een kracht die zo sterk is dat het gedrag van één elektron onlosmakelijk verbonden is met het gedrag van zijn buren. Om te begrijpen hoe deze materialen werken, moeten wetenschappers de "Groen's functie" berekenen, een wiskundig hulpmiddel dat fungeert als een kaart van hoe een elektron door deze drukke, chaotische omgeving beweegt. Deze kaart onthult de energieniveaus van het materiaal en hoe het reageert op licht of elektriciteit, maar het berekenen ervan voor grote systemen is een monumentale taak. Het aantal mogelijke arrangementen voor de elektronen groeit zo explosief met de grootte van het systeem dat zelfs de krachtigste supercomputers snel het geheugen tekortkomen, waardoor onderzoekers gedwongen worden om slechts minuscule, vereenvoudigde versies van de materialen te bestuderen die ze eigenlijk belangrijk vinden.

Een team van onderzoekers heeft nu een nieuwe manier voorgesteld om dit probleem aan te pakken met behof van kwantumcomputers, niet door te proberen het gehele systeem in één keer te simuleren, maar door een reeks slimme, kleine snapshots te nemen. In een studie gericht op een modelsysteem dat het gedrag van elektronen in een metaal nabootst, hebben de onderzoekers een techniek aangepast genaamd sample-gebaseerde Krylov-kwantumdengeneratie (Krylov quantum diagonalization). In plaats van een kwantumcomputer te vragen het volledige, enorme stadium van het systeem in zijn geheugen te houden, vroegen ze de computer om het systeem voor zeer korte perioden van tijd te laten evolueren en vervolgens de resultaten te meten. Door dit proces vele malen te herhalen en de meest voorkomende uitkomsten te verzamelen, waren zij in staat de volledige kaart van de elektronbeweging te reconstrueren zonder ooit de onmogelijke hoeveelheid gegevens te hoeven opslaan die nodig is door traditionele methoden. De aanpak is ontworpen om te werken op de ruisige, imperfecte kwantumcomputers die vandaag de dag beschikbaar zijn, die nog niet in staat zijn om de lange, complexe circuits uit te voeren die vereist zijn door oudere methoden.

De onderzoekers testten hun methode op een specif kind model bekend als het single-impurity Anderson-model, dat een enkel atoom beschrijft dat interageert met een zee van elektronen. Ze zetten een simulatie op waarbij dit centrale atoom verbonden was met een keten van twaalf locaties, die de omgeving vertegenwoordigen. In deze opstelling draagt het centrale atoom een sterke afstotende kracht die voorkomt dat twee elektronen dezelfde plek innemen, een cruciaal kenmerk dat de overgang drijft tussen een metaal dat elektriciteit geleidt en een isolator die dat niet doet. Het team gebruikte een klassieke computer om het kwantumapparaat te simuleren, waardoor ze hun nieuwe methode konden vergelijken met de exacte, perfecte oplossing voor dit systeem. Ze ontdekten dat ze, door de kwantumevolutie te samplen, de energiemap van het elektron met hoge nauwkeurigheid konden reconstrueren met slechts een fractie van de totale mogelijke toestanden. Voor een systeem met twaalf locaties verminderde de methode de vereiste computationele ruimte met drie tot vier ordes van grootte, wat betekent dat het slechts enkele duizenden toestanden moest verwerken in plaats van de honderden miljoenen die nodig zouden zijn voor een volledige berekening.

De studie onthulde dat de methode het beste werkt wanneer het systeem zich in een staat bevindt die noch volledig eenvoudig, noch volledig chaotisch is. Wanneer de interactie tussen elektronen zeer zwak is, gedraagt het systeem zich als een kalme zee van onafhankelijke deeltjes, en wanneer de interactie zeer sterk is, vergrendelen de elektronen zich in een rigide patroon. In beide van deze extreme gevallen vond de methode het antwoord gemakkelijk. Echter, in het middengebied, waar het materiaal op de drempel staat om van een metaal naar een isolator te veranderen, bevinden de elektronen zich in een complexe, fluctuerende staat. Hier ontdekten de onderzoekers dat de methode een iets groter aantal samples nodig had om het volledige beeld te vangen, maar het slaagde er nog steeds in waar traditionele methoden zouden hebben gefaald. De resultaten toonden aan dat de spectrale kenmerken, die de duidelijke pieken en dalen in de energiemap vormen die de eigenschappen van het materiaal definiëren, getrouw werden gereproduceerd over een breed scala aan interactiestatussen.

Een van de meest significante bevindingen was dat de methode niet vereiste dat de kwantumcomputer complexe operaties uitvoert die betrokken zijn bij extra "helper"-bits, bekend als ancilla's, die vaak een bron van fouten zijn in huidige apparaten. In plaats daarvan hoefde het kwantumapparaat alleen de toestand voor te bereiden, deze een korte tijd te laten evolueren en vervolgens de posities van de elektronen te meten. Al het zware werk om deze metingen te verbinden met het uiteindelijke antwoord werd gedaan door een klassieke computer. Deze taakverdeling houdt het kwantumcircuit ondiep en eenvoudig, wat het een realistische kandidaat maakt voor de hardware die in de nabije toekomst beschikbaar zal zijn. De onderzoekers merkten op dat hoewel de methode momenteel een simulatie is, de manier waarop het vereiste aantal samples groeit naarmate het systeem groter wordt, bemoedigend is. Naarmate de keten van locaties groeide van acht naar twaalf, nam het fractie van de totale mogelijke toestanden die nodig was om een accuraat antwoord te krijgen, feitelijk af, wat suggereert dat de methode kan schalen naar zelfs grotere systemen die momenteel buiten bereik liggen.

Het werk adresseerde ook een veelvoorkomende vraag over welk startpunt het beste is voor deze berekeningen. Het team testte twee verschillende manieren om de initiële toestand van het systeem op te zetten: één gebaseerd op een eenvoudig, niet-interagerend beeld en een andere gebaseerd op een complexer, zelfconsistent beeld. Ze ontdekten dat het eenvoudigere startpunt net zo effectief was als het complexere voor deze specifieke taak, wat de opstelling voor toekomstige experimenten vereenvoudigt. Bovendien demonstreerden ze dat de methode de wiskundige "vertakking" (branching) kon afhandelen die vereist is om te berekenen hoe elektronen aan of van het systeem worden toegevoegd of verwijderd, een noodzakelijke stap voor het bouwen van de volledige energiemap. Door zorgvuldig te selecteren welke gesamplede toestanden behouden moesten worden en welke weggegooid, zorgden zij ervoor dat het uiteindelijke resultaat de fundamentele wetten van de fysica respecteerde, zoals de behoud van waarschijnlijkheid.

Deze aanpak biedt een veelbelovend pad voorwaarts voor de dynamische veldentheorie (dynamical mean-field theory), een krachtig kader dat wordt gebruikt om sterk gecorreleerde materialen te bestuderen. In dit kader wordt een complex rooster van atomen in kaart gebracht op een enkel onzuiverheidsprobleem (single impurity problem), dat herhaaldelijk moet worden opgelost om de eigenschappen van het hele materiaal te vinden. De flessenhals is altijd de onzuiverheid-oplosser (impurity solver) geweest, de stap die het gedrag van elektronen voor die enkele locatie berekent. Als deze oplosser op een kwantumcomputer met een groter aantal bad-sites kan worden uitgevoerd dan momenteel mogelijk is, zou dit wetenschappers in staat stellen om materialen met veel grotere precisie te modelleren. De onderzoekers suggereren dat hun methode als een nieuwe onzuiverheid-oplosser zou kunnen dienen, wat het bestuderen van grotere en realistischere systemen mogelijk maakt dan ooit tevoren. Hoewel de studie beperkt was tot een specifiek model en nog geen volledige zelfconsistente lus van de theorie bevatte, bieden de resultaten een sterk bewijs van concept dat een dergelijke berekening haalbaar is.

De studie keek ook naar hoe de methode presteert onder verschillende omstandigheden, zoals het variëren van de sterkte van de elektronrepulsie. Ze ontdekten dat het aantal samples dat nodig is om een specifiek nauwkeurigheidsniveau te bereiken geen eenvoudige rechte lijn is; het piekt bij intermediaire interactiestatussen voordat het weer daalt bij zeer hoge sterktes. Dit gedrag is fysiek logisch, aangezien het intermediaire gebied het gebied is waar het systeem het meest complex is en de meeste informatie vereist om te beschrijven. Ondanks deze piek bleef het totale aantal samples dat nodig was een kleine fractie van de totale mogelijke toestanden, zelfs voor de moeilijkste gevallen. De onderzoekers concludeerden dat de methode robuust is en dat de gunstige schaling die zij in hun simulaties observeerden suggereert dat het een praktische tool kan zijn voor kwantumhardware op de korte termijn.

Uiteindelijk demonstreert dit werk dat we niet hoeven te wachten op perfecte, foutvrije kwantumcomputers om enkele van de moeilijkste problemen in de gecondenseerde materiefysica op te lossen. Door de probabilistische aard van kwantummeting te omarmen en slimme samplingtechnieken te gebruiken, kunnen onderzoekers essentiële informatie extraheren over hoe elektronen bewegen in complexe materialen. De methode transformeert een probleem dat voorheen als onhandelbaar werd beschouwd in een beheersbaar probleem, wat de deur opent naar een nieuwe generatie simulaties die kunnen helpen bij het begrijpen en ontwerpen van de geavanceerde materialen van de toekomst. De bevindingen suggereren dat met het juiste algoritme zelfs de ruisige kwantumcomputers van vandaag licht kunnen werpen op de diepste mysteries van sterk gecorreleerde materie.