Deel 2 van de reeks Van quarks tot quantumvoordeel

In het eerste deel zagen we het eindresultaat: een mesonachtig object breidt zich uit over een rooster, maar valt niet uiteen in twee vrij wegvliegende delen. In het midden verschijnt bovendien een interne oscillatie. Om te begrijpen waarom de auteurs dat een hadron noemen, moeten we terug naar de fysische bouwstenen: quarks, kleurlading, gaugevelden en confinement.

De woorden lijken vertrouwd uit de deeltjesfysica, maar krijgen in deze paper een precieze en beperkte betekenis. Het experiment simuleert geen volledig proton en geen driedimensionale SU(3)-QCD. Het gebruikt een SU(2)-lattice-gaugetheorie met één ruimtedimensie. Dat model is eenvoudiger, maar behoudt het niet-Abelse karakter en de lokale koppeling tussen materie en gaugeveld die de berekening interessant maken.

Quarks zijn geen losse kralen

In het quarkmodel bestaan mesonen uit een quark en een antiquark, terwijl baryonen in gewone QCD uit drie quarks bestaan. Dat beeld is nuttig, zolang we niet vergeten dat een hadron meer is dan de som van een paar puntdeeltjes. De quantumtoestand bevat ook het sterke veld, virtuele excitatie en superposities van verschillende materie- en veldconfiguraties.

Een quark draagt naast elektrische lading ook kleurlading. Rood, groen en blauw zijn daarvoor namen, geen zichtbare kleuren. In QCD vormt de gaugegroep SU(3) de wiskundige symmetrie van die kleurlading. Fysische hadronen zijn kleurneutraal: hun totale toestand is een singlet onder lokale gauge-transformaties.

De paper vervangt SU(3) door SU(2). De fundamentele representatie heeft dan twee kleurcomponenten in plaats van drie. Ook de singletstructuur verschilt: SU(2) is niet letterlijk een model waarin ieder detail van het QCD-meson of -baryon behouden blijft. De auteurs gebruiken meson, baryon, string en hadron voor de overeenkomstige gauge-invariante LSH-bouwstenen van dit SU(2)-model. Die terminologie beschrijft de rol in de modeltheorie, niet een directe voorspelling voor een fysisch pion of proton.

Waarom het veld onmisbaar is

Twee elektrisch geladen deeltjes kunnen we in eenvoudige gevallen beschrijven met een potentiaal tussen hun posities. Bij een gauge-theorie is het veld zelf een dynamisch onderdeel van de quantumtoestand. Op een ruimtelijk rooster staat de materie op de sites en leeft het gaugeveld op de links tussen naburige sites.

Voor SU(2) hebben de lokale ladingsdichtheid en de elektrische velden drie algebracomponenten, aangeduid met \(a=1,2,3\). Op iedere site geldt een lokale Gauss-wet. In compacte operatorvorm is de fysieke voorwaarde

\[
G^a(r)|\mathrm{phys}\rangle = 0,
\qquad a=1,2,3.
\]

De operator \(G^a(r)\) telt niet alleen de materielading op site \(r\), maar ook de elektrische flux die via de aangrenzende links binnenkomt en vertrekt. Een willekeurig bitpatroon voor materie en flux is daarom meestal onfysisch. Als op een site kleurlading verschijnt, moet de verandering van het omliggende veld daar precies bij passen.

Dit is het lokale karakter van gauge-invariantie. Het is strenger dan alleen het behoud van één totale lading voor het hele rooster. Een toestand kan globaal neutraal zijn en toch lokaal alle Gauss-wetten schenden. Voor een quantumalgoritme is dat een belangrijk risico: ruis of een onhandige encoding kan amplitudes naar een grote onfysische deelruimte lekken.

flowchart LR
    L["inkomende kleurflux"] --> S["site met quarklading"]
    S --> R["uitgaande kleurflux"]
    G["Gauss-wet op de site"] --> L
    G --> S
    G --> R

Van kleurflux naar een string

Stel dat we een quark en antiquark op verschillende plaatsen creëren. De lokale Gauss-wet vereist dan een bijpassende fluxconfiguratie tussen de twee uiteinden. In één ruimtedimensie kan die flux niet om een object heen buigen of zich over een dwarsvlak verspreiden. Zij volgt de lijn van het rooster.

Die aaneengesloten kleurflux noemen we een string. Het is geen mechanisch touw, maar een gauge-invariante beschrijving van hoe de elektrische flux de materieladingen verbindt. De energie van de configuratie hangt af van de flux en de lengte waarover die aanwezig is. Daardoor is het energetisch ongunstig om de uiteinden zonder grens van elkaar te scheiden.

Dat is de intuïtie achter confinement: de relevante excitatie is niet een vrij quark, maar een gebonden, kleurneutrale combinatie van materie en veld. Bij voldoende energie kan een string in een theorie met dynamische materie breken door een extra paar te vormen. De nieuwe quarks schermen de oorspronkelijke uiteinden af en produceren opnieuw gauge-invariante objecten. De afzonderlijke gekleurde constituent verschijnt nog steeds niet als vrij eindproduct.

De paper noemt in de gemeten dynamica processen als paarcreatie, paarannihilatie, string breaking en rejoining. Op de heatmap zijn dat geen afzonderlijke klassieke gebeurtenissen die we per shot volgen. Het gemeten dichtheidsprofiel is een gemiddelde over veel shots van een quantumtoestand die een superpositie van toegestane configuraties bevat.

Staggered fermions op een lijn

Een continue fermionveldtheorie rechtstreeks op een eindig rooster zetten kan ongewenste extra fermionsoorten opleveren. De Kogut-Susskind- of staggered-fermionformulering verdeelt de componenten daarom over afwisselende sites. In de conventie van de paper dragen even sites de quarkinterpretatie en oneven sites de antiquarkinterpretatie.

De lokale LSH-variabelen \(n_i(r)\) en \(n_o(r)\) geven aan of op een site een inkomend of uitgaand stringeinde aanwezig is. Hun som is nog niet zonder meer de fysische fermiondichtheid, omdat de gevulde oneven sites het referentievacuüm vormen. De paper gebruikt

\[
n_f(r)=
\begin{cases}
n_i(r)+n_o(r), & r\ \text{even},\\
2-[n_i(r)+n_o(r)], & r\ \text{oneven}.
\end{cases}
\]

Voor de strong-coupling-vacuümtoestand zijn op even sites beide stringeinden afwezig en op oneven sites beide aanwezig. Door de pariteitsafhankelijke definitie is \(n_f(r)=0\) op iedere site. Een lokaal bitpatroon dat er op oneven sites gevuld uitziet, betekent in de deeltjesinterpretatie dus juist vacuüm.

Dit is een bron van veel verwarring bij het uitlezen van het QASM. Zonder de sitepariteit zou de totale deeltjesdichtheid verkeerd worden gereconstrueerd, ook wanneer alle gemeten bits technisch correct zijn.

Wat is het centrale meson?

De tweede initial state verschilt van het vacuüm op de twee middelste sites. Daar wordt \(n_f=1\) geplaatst op een even-onevenpaar. In de LSH-beschrijving zijn dat twee stringeinden die samen een centraal meson vormen. De rest van het rooster blijft in dezelfde achtergrondsector als het SCV.

Tijdens de evolutie kan het mesonprofiel over meer sites worden uitgesmeerd. Dat betekent niet dat in ieder shot een klassiek object steeds iets groter wordt. De unitaire evolutie bouwt een superpositie op waarin verschillende toegestane materie- en fluxconfiguraties verschillende amplitudes hebben. De gemeten \(n_f(r,t)\) is de verwachting van de lokale bezetting in die toestand.

Door de SCV-evolutie af te trekken, wordt gevraagd welke dichtheidsverandering specifiek door het centrale meson is veroorzaakt. De resulterende

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

kan zowel positief als negatief zijn. Een negatieve waarde is geen negatieve kans of negatief aantal deeltjes. Zij betekent dat de lokale dichtheid in de mesonrun lager is dan in de afzonderlijk geëvolueerde vacuümrun.

Confinement in het gemeten profiel

Bij vrije deeltjes zou een lokaal pakket zich verspreiden volgens hun dispersierelatie. In de paper blijven de relevante fluctuaties binnen een begrensd, gekromd ruimtetijdgebied. De auteurs interpreteren de vorm van die rand samen met de klassieke referenties als confined propagation: de constituenten bewegen, maar blijven deel van een gebonden mesonprofiel.

Binnen de begrensde regio verschijnen afwisselende rode en blauwe banden in de differentieheatmap. Die structuur wordt gebruikt om een breathing frequency te fitten. Fysisch is dat vergelijkbaar met een gebonden object dat periodiek energie uitwisselt tussen verschillende interne configuraties. De frequentie vormt een proxy voor het verschil \(E_1-E_0\) tussen de aangesproken lage excitatie en de initiële mesontoestand.

Het woord proxy is belangrijk. De quantumprocessor diagonaliseert de volledige Hamiltoniaan niet en levert niet rechtstreeks een exact energiespectrum. De energie-informatie wordt indirect uit een beperkte, gedempte tijdreeks gehaald. De betrouwbaarheid hangt dus af van het tijdvenster, het fitmodel, de ruis en de overeenkomst tussen de verschillende methoden.

Wat is hier aangetoond?

De gebruikte SU(2)-theorie bevat lokale niet-Abelse Gauss-wetten, dynamische fermionische materie en fluxstrings. Het centrale object is een gauge-invariante, gebonden mesonconfiguratie binnen die theorie. De gemeten ruimtetijddynamica is in het vroege tijdvenster consistent met begrensde voortplanting en een interne oscillatiemodus.

Wat is nog niet aangetoond?

De simulatie reproduceert niet alle kleuren, dimensies en interacties van fysische QCD. De visuele begrenzing van één heatmap is op zichzelf ook geen bewijs voor confinement. De interpretatie steunt op de Hamiltoniaan, de gauge-invariante encoding en de vergelijking met de klassieke referenties.

In deel 3 maken we die onderliggende berekening concreet. We beginnen bij de Kogut-Susskind-Hamiltoniaan en laten zien hoe de LSH-basis de lokale niet-Abelse Gauss-wetten omzet in direct bruikbare fysieke bouwstenen.

Bronnen

1. Fran Ilčić et al., paper v3, 2026.
2. I. Raychowdhury en J. R. Stryker, Loop, string, and hadron dynamics in SU(2) Hamiltonian lattice gauge theories, Phys. Rev. D 101, 114502 (2020).
3. C. W. Bauer et al., Quantum simulation of fundamental particles and forces, Nature Reviews Physics 5, 420-432 (2023).
4. J. Kogut en L. Susskind, Hamiltonian Formulation of Wilson's Lattice Gauge Theories, Phys. Rev. D 11, 395 (1975).