Deel 1 van de reeks Van quarks tot quantumvoordeel

Een hadron is geen star bolletje met een quark en een antiquark die rustig op hun plaats blijven. Het is een dynamische gebonden toestand. De materie en het veld dat de materie bijeenhoudt, kunnen intern bewegen, energie uitwisselen en deeltjesparen vormen of vernietigen. Precies die real-time dynamica is moeilijk te berekenen. De paper Observation of Robust and Coherent Non-Abelian Hadron Dynamics on Noisy Quantum Processors onderzoekt of een grote, nog foutgevoelige quantumprocessor toch toegang kan geven tot dat gedrag.

De auteurs simuleren geen volledig proton en ook niet de volledige theorie van de sterke kernkracht. Zij kiezen een eenvoudiger, maar nog steeds niet-Abels model: een SU(2)-lattice-gaugetheorie in één ruimtedimensie en één tijdsdimensie. Daarin laten zij een mesonachtige toestand evolueren op een rooster van 60 sites. Voor iedere site worden twee dynamische binaire variabelen op qubits gezet. De berekening gebruikt daardoor 120 logische qubits op de 156-qubit IBM Boston-processor.

Het resultaat is een ruimtetijdbeeld waarin het meson zich uitbreidt, maar zijn bestanddelen niet onbeperkt van elkaar wegvliegen. Binnen die begrensde structuur verschijnt bovendien een vroege interne oscillatie: een breathing mode. Het hadron wordt afwisselend breder en smaller. De auteurs behandelen de frequentie van die oscillatie als een spectroscopische aanwijzing voor de laagste excitatie-energie van het meson.

Dat klinkt als één quantumexperiment, maar de overtuigingskracht komt juist uit vier onderling verbonden berekeningen. Om te begrijpen waarom, moeten we eerst zien waarom deze dynamica klassiek zo lastig is.

Waarom real-time sterke-interactiefysica moeilijk is

De sterke interactie wordt in het standaardmodel beschreven door quantumchromodynamica, of QCD. In die theorie dragen quarks een kleurlading en reageren zij via niet-Abelse gaugevelden. Niet-Abels betekent onder andere dat de lokale symmetrieoperaties niet simpelweg met elkaar verwisselbaar zijn. De volgorde van bewerkingen kan ertoe doen, en de lokale Gauss-wetten koppelen materie en veldconfiguraties sterk aan elkaar.

Voor statische eigenschappen is Euclidische Monte Carlo bijzonder krachtig. Real-time evolutie veroorzaakt echter een fase- of tekenprobleem, waardoor die methode haar gebruikelijke probabilistische interpretatie verliest. Exacte toestandssimulatie loopt tegen een tweede barrière aan: de Hilbertruimte groeit exponentieel met het aantal lokale vrijheidsgraden.

Tensor networks kunnen dat probleem in één dimensie vaak lang uitstellen. Zij comprimeren een quantumtoestand door vooral de relevante verstrengeling te bewaren. Na een quench of andere niet-evenwichtsevolutie kan die verstrengeling echter snel groeien. Dan moet ook de interne bonddimensie van het tensor network groeien om dezelfde nauwkeurigheid te behouden. Dat is de zogeheten entanglement wall.

Een quantumprocessor hoeft de volledige toestand niet als een klassieke lijst van amplitudes op te slaan. De 120 qubits dragen de toestand rechtstreeks. Dat neemt de andere problemen niet weg: het model moet efficiënt worden gecodeerd, de tijdsevolutie moet in gates worden ontbonden en de hardware voegt ruis toe. De vraag is daarom niet alleen of het circuit kan worden uitgevoerd, maar of de gemeten bits nog aantoonbaar iets zeggen over de bedoelde hadronfysica.

Een vereenvoudigd maar niet-triviaal hadronmodel

De paper gebruikt SU(2) in plaats van de SU(3)-symmetrie van volledige QCD en beperkt de ruimte tot een lijn. Dat maakt het model niet realistisch genoeg om direct het gedrag van een proton of een botsingsexperiment te voorspellen. Het behoudt wel twee eigenschappen die voor de computationele uitdaging essentieel zijn: de gaugegroep is niet-Abels en materie is dynamisch gekoppeld aan het gaugeveld.

Op het staggered rooster vertegenwoordigen even sites quarkachtige materie en oneven sites antiquarkachtige materie. Een meson bestaat in deze beschrijving uit een gebonden quark-antiquarkconfiguratie met elektrische kleurflux ertussen. De lokale Gauss-wet staat niet iedere willekeurige combinatie van materie en flux toe. Alleen gauge-invariante toestanden zijn fysisch.

Om die beperking efficiënt te behandelen, gebruikt de paper de Loop-String-Hadron- of LSH-basis. In plaats van eerst redundante gaugevelden te coderen en daarna de onfysische toestanden weg te filteren, beschrijft LSH de lokale toestand met gauge-invariante bouwstenen: doorgaande flux, inkomende en uitgaande stringeinden en lokale hadronconfiguraties. In het gebruikte zwakke-koppelingsregime kunnen de twee dynamische stringeindvariabelen per site rechtstreeks in twee qubits worden opgeslagen.

Het gevolg is een fysica-eigen encoding: 60 sites worden 120 qubits. In deel 3 van deze reeks bouwen we deze mapping precies op. Voor nu is het belangrijke punt dat de qubits geen abstracte optimalisatievariabelen voorstellen. Hun bezetting is gekoppeld aan lokale materie- en fluxinformatie van het rooster.

Twee experimenten voor ieder tijdstip

Het gewenste signaal is klein ten opzichte van de achtergrond en de hardware is ruisachtig. Daarom voeren de auteurs voor ieder Trottertijdstip twee evoluties uit.

De eerste begint in het strong-coupling vacuum, afgekort SCV. In de gebruikte deeltjesconventie heeft deze toestand lokaal geen fermionexcitaties. De tweede begint met hetzelfde vacuüm, maar met een meson op de twee middelste sites. Beide toestanden worden met hetzelfde dynamische circuit geëvolueerd en daarna in de Z-basis gemeten.

Uit de bitkansen wordt een lokale fermiondichtheid n_f(r,t) gereconstrueerd. Vervolgens trekken de onderzoekers het geëvolueerde vacuümsignaal af van het geëvolueerde mesonsignaal:

\[
\Delta n_f(r,t) = n_f^{\mathrm{meson}}(r,t)
– n_f^{\mathrm{SCV}}(r,t).
\]

Deze differentiële meting is een centraal onderdeel van het experiment. Ruis die beide circuits ongeveer op dezelfde manier beïnvloedt, en dynamica die al in de vacuümachtergrond aanwezig is, wordt gedeeltelijk onderdrukt. Wat overblijft is beter gericht op de extra dynamica die door het centrale meson is veroorzaakt.

flowchart LR
    A["Strong-coupling vacuum"] --> C["Dezelfde tijdsevolutie U(t)"]
    B["Vacuum plus centraal meson"] --> C
    C --> D["Z-metingen op 120 qubits"]
    D --> E["Reconstructie van n_f(r,t)"]
    E --> F["Meson minus SCV"]
    F --> G["Ruimtetijdprofiel en breathing mode"]

Het is geen algemene foutcorrectie. Fouten die de twee circuits verschillend beïnvloeden verdwijnen niet automatisch, en een verschil van twee ruizige metingen kan ook extra statistische onzekerheid bevatten. De methode is hier bruikbaar omdat zij aansluit bij de fysische vraag: wat verandert er wanneer we één meson boven op dezelfde achtergrond plaatsen?

Wat verschijnt in het ruimtetijdbeeld?

De gemeten Delta n_f(r,t) wordt als heatmap weergegeven, met de roosterpositie op de ene as en de tijdstap op de andere. Vanuit het midden breidt het signaal zich naar links en rechts uit. De randen vormen een begrensd gebied dat op een light cone lijkt, maar de quarkachtige bestanddelen vliegen niet als vrije deeltjes uit elkaar. De paper interpreteert dit als confined meson propagation.

Binnen dat gebied wisselen positieve en negatieve dichtheidsverschillen elkaar af. Vooral rond de centrale sites ontstaat een vroegtijdige oscillatie. De mesonconfiguratie ademt als het ware. De frequentie wordt op verschillende manieren uit de centrale dichtheid, de ruimtelijke spreiding en de totale dichtheidsfluctuatie gehaald. Binnen het onderzochte tijdvenster leveren die methoden onderling consistente frequenties.

De auteurs koppelen deze breathing frequency aan het energieverschil tussen de initiële mesontoestand en de laagste aangeslagen toestand die door de evolutie wordt aangesproken. Daarmee wordt de tijdreeks een vorm van hadronspectroscopie: niet door een statisch energiespectrum direct te diagonaliseren, maar door een karakteristieke oscillatie in real time te meten.

Waarom één heatmap nog geen fysisch bewijs is

Een ruisachtig quantumcircuit kan patronen produceren die overtuigend ogen maar niet uit de bedoelde Hamiltoniaan voortkomen. Daarom gebruikt de paper een gelaagde validatie.

De eerste klassieke methode is een tensor-networkberekening met MPS, MPO en 2-site TDVP. Deze route evolueert de volledige LSH-Hamiltoniaan, met een bosonische cutoff en een maximale bonddimensie van 200. Zij bevat niet de zwakke-koppelings- en Trotterbenaderingen van het quantumcircuit. De vergelijking met deze route vraagt dus: volgt het benaderde algoritme in het toegankelijke tijdvenster nog de fysica van het volledige LSH-model?

De tweede methode is Pauli-propagatie. Deze berekent in het Heisenbergbeeld hoe de gemeten Pauli-observabelen achterwaarts door het ideale circuit evolueren. Zij simuleert daarmee het ruisvrije quantumcircuit, niet rechtstreeks de volledige LSH-Hamiltoniaan. De vergelijking tussen tensor network en Pauli-propagatie controleert vooral de modelbenadering en digitale implementatie. De vergelijking tussen Pauli-propagatie en QPU laat vervolgens zien hoeveel de hardware van het ideale circuit afwijkt.

Daarnaast volgen alle methoden de globale ladingen. Voor de hoofdinstantie hoort de totale lading in de sector Q=60 te blijven en de netto lattice flux bij q=0. Behoud van die grootheden is een krachtige sanity check: grote afwijkingen zouden betekenen dat de evolutie de bedoelde fysieke sector verlaat. Maar ook hier geldt een grens. Een correcte totale lading bewijst niet dat het volledige ruimtelijke mesonprofiel correct is. Veel verschillende lokale toestanden kunnen dezelfde globale lading hebben.

De vier zichtbare routes in de paper hebben daarom elk een andere rol:

Route	Wat wordt berekend?	Wat wordt ermee gecontroleerd?
TN/TDVP	volledige LSH-Hamiltoniaan	fysieke referentiedynamica
PP-CPU	ideaal circuit met termlimiet	digitale circuitdynamica
PP-GPU	ideaal circuit met coëfficiëntcutoff	onafhankelijke PP-controle
QPU	werkelijk circuit met hardwarefouten	experimenteel signaal

Juist de keten van vergelijkingen maakt de interpretatie sterker dan een losse hardwaremeting.

De schaal van het experiment

De paper behandelt een rooster van 60 sites in het zwakke-koppelingsregime, met x=100 als hoofdinstantie. De Trotterstap is delta tau = 0.0015. Het circuit heeft twee logische qubits per site en gebruikt een lokale gateconstructie waarvan de diepte per Trotterstap constant is ten opzichte van de roosterlengte.

De tweequbitdiepte groeit als 13t-1. Bij 20 stappen is dat 259; bij 25 stappen

1. Het volledige 25-stappencircuit bevat volgens de paper 17.660 tweequbitgates

en meer dan 90.000 enkelqubitgates. Voor ieder gemeten tijdstip worden 10.000 shots gebruikt.

De Quantum Advantage Tracker legt specifiek de 20-stappeninstantie vast. Dat is belangrijk voor latere timingtabellen: de paper bespreekt dynamica tot 25 stappen, terwijl de trackerclaim en de gedownloade stap-20-QASM een afgebakend benchmarkpunt vormen. We zullen die twee scopes in de hele reeks uit elkaar houden.

Wat is hier aangetoond?

Binnen het onderzochte vroege tijdvenster toont de paper een ruisbestendig, differentieel signaal dat overeenkomt met begrensde mesondynamica en een interne oscillatiemodus. De overeenkomst tussen de volledige LSH-tensorberekening, de ideale circuitsimulatie en de hardware ondersteunt drie opeenvolgende stappen: de gebruikte fysische benadering, de circuitconstructie en de experimentele realisatie.

Het experiment laat ook zien waarom een physics-native encoding belangrijk kan zijn. Doordat gauge-invariante bouwstenen direct worden gecodeerd, kan een groot niet-Abels systeem met twee qubits per site en een lokaal circuit worden uitgevoerd.

Wat is nog niet aangetoond?

Dit is geen volledige simulatie van SU(3)-QCD in drie ruimtedimensies. De zwakke-koppelingsbenadering, eindige roosterafstand, open grenzen, eerste-orde-Trotterisatie en beperkte vroege tijd bepalen allemaal het geldigheidsgebied.

Evenmin volgt een algemene quantum-advantageclaim alleen uit de circuitgrootte of de hardwaretijd. Een eerlijke snelheidsclaim moet dezelfde taak, dezelfde observabele, de klassieke convergentie en de gekozen tijdsdefinitie vergelijken. De paper laat zien dat de beschikbare klassieke baselines bij grotere tijd en zwakkere koppeling snel duurder of moeilijker controleerbaar worden. Hoe sterk dat bewijs voor praktisch quantumvoordeel is, behandelen we pas in deel 7.

Onze lokale reproductie voegt daar later een tweede laag aan toe. In dit project zijn de publieke QASM-circuits en data verzameld, zijn de circuits via de zelfstandige Fire Opal-API op ibm_fez uitgevoerd, en zijn controles met Aer, ITensor, MPO en Majorana voorbereid. Die resultaten worden niet gebruikt om de eerste uitleg mooier te maken dan zij is. Ze komen pas aan bod wanneer de observabelen en timinggrenzen volledig zijn gedefinieerd.

Van waarneming naar begrip

We hebben nu het eindbeeld: een centraal meson ontwikkelt een begrensd ruimtetijdprofiel en een interne oscillatie, gemeten met 120 logische qubits en gecontroleerd met twee verschillende klassieke strategieën. Wat nog ontbreekt, is een precieze betekenis van de woorden quark, kleur, gaugeveld, confinement en meson binnen dit model.

Dat is het onderwerp van deel 2. Pas daarna bouwen we de Kogut-Susskind-Hamiltoniaan, de LSH-basis en uiteindelijk het quantumcircuit stap voor stap op.

Bronnen

1. Fran Ilčić et al., Observation of Robust and Coherent Non-Abelian Hadron Dynamics on Noisy Quantum Processors, arXiv:2602.18080v3 (2026).
2. Quantum Advantage Tracker, issue #149.
3. LSH-IBM, publieke circuitimplementatie.
4. lsh_data, publieke klassieke benchmarkdata.
5. Hadron LSH Quantum Advantage Workbench, lokale reproductie en claimgrenzen.