⚛️ quantum physics

Efficient classical simulation of large-scale unitary cluster Jastrow circuits

Dit artikel introduceert een klassiek algoritme met polynomiale tijd dat in staat is om grootschalige single-layer unitary cluster Jastrow-circuits efficiënt te simuleren, waardoor een laptop in minder dan een minuut een recent 77-qubit kwantumexperiment kan reproduceren en overtreffen.

Oorspronkelijke auteurs: Hrishikesh Belagali, Thomas Van Camp, R. Pradeep, Sourin Das, Namit Anand, Ryan LaRose

Gepubliceerd 2026-07-24

📖 6 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Hrishikesh Belagali, Thomas Van Camp, R. Pradeep, Sourin Das, Namit Anand, Ryan LaRose

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je voor dat je probeert een enorme, verwarde knoop van een touw op te lossen die een molecuul voorstelt. In de wereld van de chemie is het uitzoeken hoe deze knopen hun meest stabiele vorm (hun "grondtoestand") vinden, als het voorspellen van het weer tijdens een orkaan: het is ongelooflijk moeilijk omdat elk stukje touw aan elk ander stukje trekt. Decennialang hebben wetenschappers gehoopt dat quantumcomputers — machines die de vreemde regels van minuscule deeltjes gebruiken om wiskunde te bedrijven — deze knopen sneller kunnen ontwarren dan een normale computer. Het idee is dat terwijl een gewone computer elke mogelijke draai één voor één moet controleren, een quantumcomputer veel draaien tegelijk kan controleren. Onlangs hebben onderzoekers deze quantummachines gebouwd om complexe moleculen te simuleren, zoals ijzer-zwavelclusters die in de natuur voorkomen, in de hoop nieuwe medicijnen of materialen te vinden. De grote vraag is geweest: doen deze quantummachines eigenlijk iets wat een gewone computer niet kan, of gebruiken we gewoon een sloophamer om een noot te kraken die een gewone computer met een beetje meer slimheid ook had kunnen afhandelen?

Dit artikel gaat over een team van onderzoekers die besloten om die sloophamer te testen. Ze keken naar een specifiek type quantumexperiment genaamd de "Unitary Cluster Jastrow" (UCJ) circuit, dat is gebruikt in sommige van de grootste en meest indrukwekkende quantumchemische experimenten tot nu toe. Deze experimenten omvatten quantumcomputers met tot 77 qubits (de quantumversie van bits) en duizenden gates, draaiend op supercomputers met duizenden nodes, enkel om de resultaten te verwerken. De onderzoekers vroegen zich af: "Kunnen we deze zelfde berekening op een gewone laptop uitvoeren?" Het antwoord is, verrassend genoeg, ja. Ze ontwikkelden een nieuwe, supersnelle wiskundige truc die een standaardcomputer in staat stelt om de energie van deze specifieke quantumcircuits in polynomiale tijd te berekenen — wat betekent dat de tijd die nodig is om de grootte van het probleem groeit, redelijk toeneemt, in plaats van exponentieel naar oneindig te exploderen.

Het team kwam erachter dat ze de resultaten van het grootste experiment ooit gedaan op een ijzer-zwavelcluster (dat 77 qubits en 10.570 gates gebruikte) in minder dan een minuut konden reproduceren op een laptop. Sterker nog, omdat hun methode zo snel was, konden ze de parameters van het circuit aanpassen om een nog lagere energietoestand te vinden dan de toestand die het quantumexperiment bereikte, wat de nodige tijd had gekost op een supercomputer met 6.400 nodes en uren aan verwerkingstijd. Er is echter een addertje onder het gras. Het artikel stelt expliciet dat deze "overwinning" alleen werkt voor single-layer circuits. Als je meer lagen aan het circuit toevoegt (het dieper en complexer maakt), stort de wiskunde in en wordt het probleem waarschijnlijk weer te moeilijk voor gewone computers. De auteurs concluderen dat single-layer UCJ-circuits niet voldoende zijn om quantumvoordeel te bewijzen (het punt waarop quantumcomputers klassieke computers verslaan), maar dat we veel diepere, multi-layered circuits zullen moeten bouwen om echt te zien wat quantumcomputers kunnen doen wat klassieke computers niet kunnen.

Om te begrijpen hoe ze dit deden, denk aan het quantumcircuit als een complexe machine die een begintoestand (zoals een plat vel papier) transformeert in een uiteindelijke, gekreukte vorm. Meestal moet je de machine laten draaien en vervolgens proberen de kreukeling terug te ontwerpen, wat een nachtmerrie is voor gewone computers omdat het aantal mogelijkheden astronomisch is. De truc van de auteurs was om achteruit te werken. In plaats van te proberen de uiteindelijke gekreukte vorm te voorspellen, namen ze de "regels" van het molecuul (de Hamiltonian) en lieten ze deze achteruit door de machine lopen.

Stel je voor dat je een recept hebt voor een cake, maar in plaats van de cake te bakken en hem daarna te proeven, begin je met de afgewerkte cake en werk je achteruit door het recept om precies te zien welke ingrediënten zijn gebruikt. De auteurs toonden aan dat je voor deze specifieke single-layer circuits het recept achteruit kunt draaien zonder dat het aantal ingrediënten tot oneindig explodeert. Ze gebruikten een wiskundig instrument genaamd de formule van Löwdin, wat een soort speciale rekenmachine is die snel de waarde van een complexe vorm kan bepalen zonder dat elk afzonderlijk punt gemeten hoeft te worden. Door deze achterwaartse methode te combineren met een slimme manier om de "fasen" (de timing en hoeken van de quantumbewegingen) te hanteren, hielden ze de berekening efficiënt.

De resultaten waren opvallend. Toen ze hun methode toepasten op het ijzer-zwavelclusterexperiment, kregen ze een energiewaarde van -326,796 Hartrees, wat lager (beter) is dan de -326,645 Hartrees die het quantumexperiment bereikte met de Sample-based Quantum Diagonalization (SQD) methode. Het quantumexperiment moest een enorme supercomputer (Fugaku) met 6.400 nodes gebruiken om zijn resultaat te krijgen, terwijl de auteurs een beter resultaat op een laptop kregen in minder dan een minuut. Ze testten dit ook op waterstofketens met tot 160 qubits, waarbij ze lieten zien dat hun methode goed schaalt, hoewel de tijd die het kost inder je met het systeem groter wordt (ongeveer volgens een macht van 4,448, wat veel beter is dan de exponentiële explosie die zou optreden met oudere methoden).

De paper is echter heel duidelijk over wat dit niet betekent. Het betekent niet dat quantumcomputers nutteloos zijn. De auteurs wijzen erop dat hun methode alleen werkt voor "weak simulation", wat betekent het berekenen van de energie, maar het kan niet "strong simuleren" hoe het proces van het samplen van willekeurige bitstrings (de ruwe data-output) dat de quantumcomputer produceert, verloopt. Het oorspronkelijke experiment gebruikte die willekeurige bitstrings om fouten op te schonen en het resultaat te verbeteren. Omdat de methode van de auteurs die willekeurige strings niet genereert, kan het dat specifieke type foutcorrectie niet uitvoeren. Maar, door zo snel te zijn, konden ze de instellingen van het circuit simpelweg beter optimaliseren dan de experimentatoren, waardoor ze het eindresultaat toch versloegen.

De belangrijkste les is een grenslijn. De auteurs beargumenteren dat single-layer circuits niet de "heilige graal" van quantumvoordeel zijn. Als een gewone computer ze in een minuut kan simuleren, doen ze niets bijzonders. Om echt klassieke computers te verslaan, moeten we circuits bouwen met twee of meer lagen (L ≥ 2). De wiskunde wordt te complex voor hun truc zodra je die extra lagen toevoegt, en dat is waarschijnlijk waar de echte quantumkracht ligt. Dus, terwijl dit artikel laat zien dat we de huidige "state-of-the-art" quantumchemische experimenten op een laptop kunnen simuleren, vertelt het ons ook dat de echte race pas net is begonnen: we moeten diepere, complexere circuits bouwen om de problemen te vinden die klassieke computers werkelijk niet kunnen oplossen.

Technische Samenvatting: Efficiënte Klassieke Simulatie van Grootschalige Unitary Cluster Jastrow-circuits

Probleemstelling
Recente experimenten in de kwantumchemie hebben gebruikgemaakt van de Unitary Cluster Jastrow (UCJ) ansatz, met name de lokale variant (LUCJ), om de grondtoestanden van sterk gecorreleerde moleculen op kwantunhardware te simuleren. Opmerkelijke experimenten, zoals gerapporteerd in Sci. Adv. 11, 25 (2025), hebben circuits uitgevoerd met tot 77 qubits en meer dan 10.000 gates, wat enorme klassieke supercomputerbronnen vereist (bijv. 6400 nodes op Fugaku) voor post-processing via Sample-based Quantum Diagonalization (SQD). Hoewel de klassieke hardheid van het samplen uit UCJ-circuits is vastgesteld (gelinkt aan IQP-berekeningen), is de specifieke taak van het berekenen van de energie-verwachtingswaarde voor single-layer UCJ-circuits nog niet volledig gekarakteriseerd. Bestaande klassieke simulatiestrategieën voor deze circuits kampen met exponentiële schaling of vertrouwen op afkapmethoden (truncations) die de chemische nauwkeurigheid in gevaar brengen. Het centrale probleem dat wordt aangepakt, is of de energie van single-layer UCJ-circuits exact en efficiënt kan worden berekend op klassieke hardware, onafhankelijk van de lokaliteitsbeperkingen die aan de kwantumhardware worden opgelegd.

Methodologie
De auteurs presenteren een algoritme met een polynomiale looptijd voor de exacte berekening van de energie van elk single-layer ( $L=1$ ) UCJ-circuit. Het algoritme werkt in de tweede-gekwantiseerde formalisme en vermijdt de exponentiële groei van termen die gewoonlijk gepaard gaat met backpropagation in andere simulatiemethoden. De aanpak verloopt in drie afzonderlijke stappen:

Backpropagation door Orbitale Rotaties: Het algoritme propageert eerst de fermionische Hamiltonian achterwaarts door de laatste orbitale rotatie-operatoren ( $e^{K_\mu}$ ). Met behulp van Lemma 1 worden de één-elektron ( $h$ ) en twee-elektron ( $g$ ) integralen exact bijgewerkt via unitaire transformaties. Deze stap wordt exact uitgevoerd in respectievelijk $O(N^3)$ en $O(N^5)$ tijd, waarbij $N$ het aantal orbitalen is, zonder het aantal termen in de Hamiltonian te vergroten.
Propagatie door Jastrow-operatoren: De Hamiltonian wordt vervolgens gepropageerd door de Jastrow-operatoren ( $e^{iJ_\mu}$ ). Lemma 2 demonstreert dat deze operatie ook exact kan worden uitgevoerd in $O(N^5)$ tijd. Cruciaal is dat deze stap het aantal termen niet vergroot; in plaats daarvan accumuleert elke term in de Hamiltonian een scalaire fase en een vectorfase afhankelijk van de aantallen-operatoren (number operators). Dit geldt voor zowel algemene UCJ als ijle (sparse) LUCJ-circuits.
Energie-evaluatie via Löwdin's Formule: De uiteindelijke energie wordt berekend als de verwachtingswaarde van de gewijzigde Hamiltonian tegenover de referentie-Slater-determinant. Omdat de Jastrow-propagatie niet-orthogonale fasen introduceert, worden de matrixelementen geëvalueerd tussen niet-orthogonale Slater-determinanten. De auteurs gebruiken Löwdin's formule (Lemma 3) om deze matrixelementen efficiënt te berekenen. De computationele bottleneck betreft het vormen en factoriseren van een matrix $S$ (gerelateerd aan de overlap van determinanten), wat leidt tot een algemene algoritmische complexiteit van $O(N^7)$ . De auteurs merken op dat in de praktijk het memoïseren van matrixelementen de geobserveerde schaling reduceert tot ongeveer $O(N^{4.45})$ .

Belangrijkste Bijdragen

Algoritme met Polynomiale Tijd: Het artikel stelt vast dat het berekenen van de energie van een single-layer UCJ-circuit klassiek hanteerbaar is in $O(N^7)$ tijd, onafhankelijk van de connectiviteitsbeperkingen van de kwantumhardware.
Exacte Simulatie: In tegen afstelling tot eerdere methoden die vertrouwen op afkapmethoden of heuristische benaderingen, biedt dit algoritme een exacte berekening van de energie-verwachtingswaarde.
Onderscheid van Sterke Simulatie: De auteurs verduidelijken dat hun methode een "zwakke simulatie" (het berekenen van observabelen) is in plaats van een "sterke simulatie" (het samplen van bitstrings). Bijgevolg kan het algoritme geen bitstring-gebaseerde post-processing zoals SQD uitvoeren, die in experimenten wordt gebruikt om hardware-ruis te mitigeren.

Resultaten
De auteurs valideren hun algoritme via twee primaire casestudy's:

IJzer-Zwavel Cluster (72 Qubits): Het algoritme werd toegepast op het 77-qubit (72 actieve qubits) ijzer-zwavel cluster experiment uit Sci. Adv. 11, 25 (2025). Terwijl het oorspronkelijke experiment 6400 nodes op Fugaku nodig had voor 1,44 uur om een energie van $-326,645$ Ha te bereiken via SQD, simuleerde het klassieke algoritme van de auteurs hetzelfde circuit op een laptop in minder dan één minuut. Bovendien bereikten zij, door de single-layer UCJ-parameters te optimaliseren met de snelle klassieke energie-evaluator (L-BFGS-B), een grondtoestandsenergie van $-326,796$ Ha, waarmee zij de experimentele resultaten verkregen via SQD overtroffen.
Waterstofketens (tot 160 Qubits): Het algoritme werd geschaald naar de simulatie van waterstofketens met tot 160 qubits, wat de huidige grenzen van de meest geavanceerde supergeleidende kwantumcomputers vertegenwoordigt. De resultaten toonden aan dat de single-layer UCJ-ansatz, geïnitialiseerd met Coupled Cluster Singles and Doubles (CCSD) parameters, de Hartree-Fock energieën nauwgezet volgt (gemiddeld verschil van ~11 mHa). Optimalisatie van parameters voor systemen tot 40 qubits leverde energieën op die vergelijkbaar zijn met CCSD.

Betekenis en Claims
Het artikel concludeert dat single-layer UCJ-circuits onvoldoende zijn om kwantumvoordeel in de chemie aan te tonen. Het vermogen om deze circuits exact en efficiënt te simuleren op klassieke hardware impliceert dat het geclaimde "kwantumvoordeel" in recente grootschalige experimenten berust op de specifieke post-processing (SQD) of de potentie van diepere circuits.

De auteurs stellen expliciet dat het bereiken van berekeningen die de klassieke limieten overstijgen waarschijnlijk $L \ge 2$ lagen van de (L)UCJ-ansatz vereist. Zij argumenteren dat hun algoritme faalt voor $L \ge 2$ omdat de orbitale rotaties tussen de Jastrow-lagen de diagonale structuur verstoren die vereist is voor de formule van Löwdin. Gegeven dat alternerende lagen van Givens-rotaties en CZ-gates bekend staan als universeel voor de kwantumchemie, suggereren de auteurs dat hoewel $L \ge 2$ noodzakelijk is voor kwantumvoordeel, het mogelijk niet voldoende is zonder de aanzienlijke gate-overhead op de huidige hardware aan te pakken. Dit werk dient als benchmark en toont aan dat huidige single-layer experimenten de capaciteiten van geoptimaliseerde klassieke algoritmen nog niet overstijgen.

Technische Samenvatting: Efficiënte Klassieke Simulatie van Grootschalige Unitary Cluster Jastrow-circuits

Meer zoals dit