Het model van deze pagina is het volledige, spinvolle tweedimensionale Fermi–Hubbardmodel met hopping t, diagonale hopping t′ en lokale afstoting U. De parameters zijn U/t = 8 en t′/t = −0,25. “Volledig” betekent hier dat alle termen van dit gekozen model én de fermionische tekens worden meegenomen. Het is geen model met alle denkbare materiaalinteracties. De historische bare-XY-pilots en de latere volledige fermionroute zijn verschillende berekeningen; resultaten van de eerste bewijzen niet de juistheid van de tweede.
De uitleg hieronder bouwt het model op vanuit roosterplaatsen en bezetting. Daarna volgen de energie- en tijdseenheden, de omzetting naar qubits en een exact oplosbaar voorbeeld. Het onderscheid tussen een fysisch model, zijn qubitrepresentatie en de nauwkeurigheid van een uitvoering loopt door de hele uitleg heen.
1. Wat betekenen 2D, XY en de dimensie van de toestandsruimte?
2D gaat over de geometrie: een roosterplaats heeft twee coördinaten, x en y. Op ons vierkante 6 × 6-rooster zijn er 36 plaatsen. Elektronen kunnen horizontaal en verticaal bewegen; t′ voegt diagonale verbindingen toe. Open randen betekenen dat de linker- en rechterrand, en de boven- en onderrand, niet met elkaar worden verbonden.
i=x+L_x y,\qquad 0\leq x\lt L_x,\quad 0\leq y\lt L_y,\qquad L=L_xL_y=36.
\]
XY gaat over operatoren: een XY-koppeling bevat Pauli-X- en Pauli-Y-operatoren, bijvoorbeeld XX + YY. De letters zijn geen roostercoördinaten. Een XY-spinmodel kan zelf op een tweedimensionaal rooster liggen. “XY of 2D?” vergelijkt dus twee verschillende eigenschappen.
De Hilbertruimtedimensie telt mogelijke quantumtoestanden. Een site heeft vier lokale bezettingen, zodat 36 sites een volledige bezettingsruimte met 436 basisvectoren hebben. Dat is een derde betekenis van dimensie. “Full dimensional” is hier geen duidelijke modelnaam; gebruik “volledig 2D Fermi–Hubbardmodel”.
2. Eén roosterplaats: vier toestanden en twee qubits
We houden per site één ruimtelijke orbitaal bij. Een elektron kan daarin spin omhoog (↑) of spin omlaag (↓) hebben. Een modus is de combinatie van site en spin. Het uitsluitingsprincipe van Pauli laat maximaal één elektron per modus toe. Twee elektronen met tegengestelde spin mogen daardoor wel dezelfde site bezetten.
|0\rangle_i\leftrightarrow|00\rangle,\quad |\uparrow\rangle_i\leftrightarrow|10\rangle,\quad |\downarrow\rangle_i\leftrightarrow|01\rangle,\quad |\uparrow\downarrow\rangle_i\leftrightarrow|11\rangle.
\]
In deze lokale ketnotatie staat de up-bezetting links en de down-bezetting rechts. Een qubitwaarde 1 betekent “deze modus is bezet”; 0 betekent “leeg”. De qubits coderen dus bezetting. Een qubit is niet simpelweg het spinpijltje van één vast elektron. Voor 36 sites zijn er 72 modi en in deze directe codering 72 qubits, ook als er maar 32 elektronen aanwezig zijn.
\dim\mathcal H=4^{36}=2^{72},\qquad \dim\mathcal H_{N_\uparrow=N_\downarrow=16}=\binom{36}{16}^{2}.
\]
De tweede uitdrukking telt alleen toestanden met precies zestien elektronen van elke spin. De dynamica blijft in die sector. Een grote toestandsruimte op zichzelf bewijst geen quantumvoordeel: structuur, verstrengeling, tijd en gewenste observabelen bepalen ook de klassieke kosten.
3. Van bezetting naar operatoren
Een quantumtoestand is een superpositie van bezettingspatronen. De complexe coëfficiënt an heet een amplitude; zijn absolute waarde in het kwadraat is de kans om dat patroon te meten. Een operator werkt op zo’n toestand. De Hamiltoniaan H is de energieoperator en bepaalt tevens hoe de toestand in de tijd verandert.
|\psi\rangle=\sum_{\mathbf n}a_{\mathbf n}|\mathbf n\rangle,\qquad \sum_{\mathbf n}|a_{\mathbf n}|^2=1.
\]
De operator c†iσ voegt een elektron met spin σ op site i toe; ciσ verwijdert er één. Toevoegen aan een al bezette modus geeft nul. De getaloperator niσ telt de bezetting van die modus:
n_{i\sigma}=c^\dagger_{i\sigma}c_{i\sigma},\qquad n_{i\sigma}|\mathbf n\rangle=n_{i\sigma}^{(\mathbf n)}|\mathbf n\rangle,\qquad n_{i\sigma}^{(\mathbf n)}\in\{0,1\}.
\]
Fermionoperatoren voldoen aan de anticommutatieregels hieronder. Hierbij is {A,B} = AB + BA en is δ gelijk aan 1 wanneer de indices overeenkomen, anders 0. Voor verschillende modi geeft het verwisselen van twee creatieoperatoren een minteken. Dat teken verandert interferentie tussen amplitudes en moet dus in de simulatie behouden blijven.
\{c_{i\sigma},c^\dagger_{j\sigma’}\}=\delta_{ij}\delta_{\sigma\sigma’},\qquad \{c_{i\sigma},c_{j\sigma’}\}=0.
\]
4. Hopping: wat doet t precies?
Voor twee verbonden sites i en j beschrijft c†iσcjσ een verplaatsing van j naar i, met behoud van spin. Operatoren werken van rechts naar links: eerst verwijderen op j, dan toevoegen op i. De omgekeerde verplaatsing is de Hermitisch geconjugeerde term. Beide samen maken de energieoperator Hermitisch, zodat zijn eigenwaarden reëel zijn.
H_{ij,\sigma}=-t\left(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma}\right).
\]
t is een energie, geen tijd en geen aantal sprongen. Dat wordt concreet voor één elektron op twee sites. In de basis |L⟩, |R⟩, met het elektron links of rechts, is de Hamiltoniaan een 2 × 2-matrix. De eigenenergieën zijn −t en +t. Vanuit links ontstaat coherent een superpositie:
H_{\rm one} = \begin{pmatrix}0&-t\\-t&0\end{pmatrix},\qquad |\psi(s)\rangle=\cos(ts/\hbar)|L\rangle+\mathrm{i}\sin(ts/\hbar)|R\rangle.
\]
P_R(s)=\sin^2(ts/\hbar).
\]
Hier schrijven we s voor fysieke tijd, zodat t uitsluitend de hoppingenergie blijft. De verhouding t/ℏ bepaalt een frequentieschaal. Het kansverloop is quantummechanische oscillatie; t is geen klassieke sprongfrequentie.
5. De lokale afstoting U en dubbelbezetting
De wisselwerking op één site is U ni↑ni↓. Het product van beide bezettingen is alleen 1 als beide spins aanwezig zijn. Een dubbelbezette site heet ook een doublon. In de volgorde leeg, up, down, dubbelbezet is de interactiematrix:
H_{U,i}=U n_{i\uparrow}n_{i\downarrow}=\operatorname{diag}(0,0,0,U).
\]
U > 0 betekent afstoting: dubbelbezetting voegt energie toe. Voor een bezettingspatroon met d dubbelbezette sites is de interactiebijdrage Ud. Dat is alleen de interactie-energie; de volledige energie bevat ook hopping. In een superpositie kan het gemiddelde doublonaantal niet-geheel zijn.
\langle H_U\rangle=U\sum_i\langle n_{i\uparrow}n_{i\downarrow}\rangle.
\]
Een grote positieve U maakt dubbelbezetting energetisch ongunstig, maar verbiedt haar niet. Alleen een extra projectie of een geschikte oneindige-U-limiet zou dubbelbezette toestanden uitsluiten. Bij U/t = 8 blijven ze in de toestandsruimte en kunnen ze tijdens de dynamica ontstaan.
6. De volledige Hamiltoniaan, term voor term
\begin{aligned}H&=H_t+H_{t’}+H_U,\\ H_t&=-t\sum_{\langle i,j\rangle,\sigma}\left(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma}\right),\\ H_{t’}&=-t’\sum_{\langle\!\langle i,j\rangle\!\rangle,\sigma}\left(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma}\right),\\ H_U&=U\sum_i n_{i\uparrow}n_{i\downarrow}.\end{aligned}
\]
De eerste som loopt over horizontale en verticale buren; de tweede over diagonale buren. Elk ongeordend paar telt één keer en voor elk paar worden beide spins opgeteld. Op een open 6 × 6-rooster zijn er 60 horizontale/verticale verbindingen en 50 diagonale verbindingen. Er zijn daarnaast 36 lokale U-termen. Het getal 8 is dus een coëfficiënt, geen telling van deze termen.
We kiezen t > 0. Omdat t′ = −0,25t, is de coëfficiënt −t′ in de Hamiltoniaan juist +0,25t. Het relatieve teken is relevant voor interferentie. Deze target gebruikt μ = 0. De code ondersteunt ook een uniforme term −μN; binnen een vaste-N-sector geeft die uitsluitend een globale fase aan de tijdsevolutie.
Literatuur bij deze Hamiltoniaan: het oorspronkelijke artikel van Hubbard [1] en het overzicht van Arovas et al. [2].
7. U = 8 betekent U/t = 8: de eenheden volledig uitgeschreven
t, t′ en U hebben allemaal de dimensie energie. Een getal zoals “U = 8” is daarom pas volledig als we de energie-eenheid noemen. Hier kiezen we de positieve hoppingenergie t als eenheid en delen de hele Hamiltoniaan door t. Noem de resulterende dimensieloze operator h:
h=\frac{H}{t},\qquad u=\frac{U}{t}=8,\qquad r=\frac{t’}{t}=-0.25.
\]
\begin{aligned}h={}&-\sum_{\langle i,j\rangle,\sigma}(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma})\\ &+0.25\sum_{\langle\!\langle i,j\rangle\!\rangle,\sigma}(c^\dagger_{i\sigma}c_{j\sigma}+c^\dagger_{j\sigma}c_{i\sigma})\\ &+8\sum_i n_{i\uparrow}n_{i\downarrow}.\end{aligned}
\]
“t = 1, U = 8” is dus een keuze van eenheden: één energie-eenheid is t en de lokale dubbelbezettingsenergie is acht van die eenheden. Als een hypothetische fysieke realisatie t = 0,25 eV had, zou U = 2 eV zijn. Bij t = 0,10 eV zou dezelfde verhouding U = 0,80 eV geven. Deze voorbeelden wijzen ons model niet aan een werkelijk materiaal toe.
Voor tijd beginnen we bij de Schrödingervergelijking. s is fysieke tijd, ℏ is de gereduceerde constante van Planck en τ = ts/ℏ is dimensieloos. Substitutie geeft:
\mathrm{i}\hbar\frac{\partial|\psi\rangle}{\partial s}=H|\psi\rangle,\qquad \tau=\frac{ts}{\hbar}\quad\Longrightarrow\quad \mathrm{i}\frac{\partial|\psi\rangle}{\partial\tau}=h|\psi\rangle.
\]
|\psi(\tau)\rangle=e^{-\mathrm{i}h\tau}|\psi(0)\rangle,\qquad T=0.4\ \Longleftrightarrow\ s=0.4\frac{\hbar}{t}.
\]
De projectnotatie T = 0,4 betekent de eindwaarde τ = 0,4. Bij het hypothetische t = 0,25 eV is dat ongeveer 1,05 femtoseconde aan modeltijd. De quantumcomputer voert daarvoor poorten en herhaalde metingen uit. Die uitvoeringstijd, de wachtrij en de analysetijd zijn andere grootheden. Zonder een fysieke waarde voor t wordt geen tijd in seconden vastgelegd.
8. Waarom XY-termen voorkomen in een fermionmodel
De Pauli-operatoren X, Y en Z werken op de tweedimensionale ruimte van één bezettingsqubit. I is de identiteit. In de basis |0⟩, |1⟩ zijn hun matrices:
X=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad Y=\begin{pmatrix}0&-\mathrm{i}\\\mathrm{i}&0\end{pmatrix},\qquad Z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.
\]
Daarom is n = (I − Z)/2 precies 0 op |0⟩ en 1 op |1⟩. Verder is (X + iY)/2 = |0⟩⟨1| de operator die een bezetting verwijdert. Om ook de fermionalgebra te krijgen, voegt Jordan–Wigner Z-operatoren toe over alle eerdere modi in een gekozen volgorde:
c_p=\left(\prod_{k\lt p}Z_k\right)\frac{X_p+\mathrm{i}Y_p}{2},\qquad c_p^\dagger=\left(\prod_{k\lt p}Z_k\right)\frac{X_p-\mathrm{i}Y_p}{2}.
\]
Voor p < q volgt hieruit de hoppingidentiteit. Operatoren op verschillende qubits worden als tensorproducten gelezen; bijvoorbeeld XpXq betekent X op p én X op q. De Z-string tussen beide uiteinden heet de pariteitsstring:
c_p^\dagger c_q+c_q^\dagger c_p=\frac12\left(X_p Z_{p+1}\cdots Z_{q-1}X_q+Y_p Z_{p+1}\cdots Z_{q-1}Y_q\right).
\]
\prod_{p\lt k\lt q}Z_k\,|\mathbf n\rangle=(-1)^{\sum_{p\lt k\lt q}n_k}|\mathbf n\rangle.
\]
Bij een even aantal bezette tussenliggende modi is de factor +1; bij een oneven aantal −1. Zijn de modi direct opeenvolgend, dan is de string leeg en gelijk aan I. Juist dan reduceert fermionische hopping exact tot een XY-term:
H_{p,p+1}=-\frac{t}{2}(X_pX_{p+1}+Y_pY_{p+1}).
\]
Dat verklaart waarom een correcte Hubbard-simulatie XY-poorten kan bevatten. De omgekeerde redenering werkt niet: een verzameling kale XY-koppelingen bevat niet automatisch de benodigde fermiontekens, de U-termen en de diagonale hopping. Of een string leeg is, hangt af van de volgorde van de modi, niet alleen van de afstand op het rooster. Zelfs een horizontale buurverbinding kan bij spin-afwisselende ordening een tussenliggende modus hebben.
Een stapsgewijze afleiding met uitvoerbare voorbeelden staat in de officiële OpenFermion-documentatie [3].
9. Een voorbeeld waarin kale XY het verkeerde teken geeft
Neem drie geordende modi 0, 1 en 2 en laat een elektron tussen 0 en 2 bewegen. We schrijven hun bezettingen hieronder van links naar rechts als n₀n₁n₂. Met een lege middelste modus geeft hopping de amplitude −t. Met een bezette middelste modus geeft de fermionische pariteit juist +t:
\langle100|H_{02}|001\rangle=-t,\qquad \langle110|H_{02}|011\rangle=+t.
\]
De kale operator −t(X₀X₂ + Y₀Y₂)/2 geeft in beide gevallen −t. Op één geïsoleerde verbinding kan het verschil verborgen blijven in een bezettingskans; bij interfererende paden en superposities is het in het algemeen fysisch relevant. De pariteitsstring weglaten is daarom geen algemene optimalisatie van hetzelfde 2D-fermionmodel.
10. Hoe U en fermionische swaps op qubits werken
Voor de up- en down-modus van dezelfde site volgt de interactie rechtstreeks uit n = (I − Z)/2. Deze term kan als een gecontroleerde fase worden uitgevoerd: alleen |11⟩ krijgt een fase. In dimensieloze tijd Δτ is die fase −uΔτ:
U n_{i\uparrow}n_{i\downarrow}=\frac{U}{4}(I-Z_{i\uparrow}-Z_{i\downarrow}+Z_{i\uparrow}Z_{i\downarrow}).
\]
e^{-\mathrm{i}u\Delta\tau\,n_{i\uparrow}n_{i\downarrow}}=\operatorname{diag}(1,1,1,e^{-\mathrm{i}u\Delta\tau}).
\]
Bij u = 8 en Δτ = 0,2 is de hoek −1,6 radialen. Als deze lokale interactie als enige actief was, zou zij de bezettingskansen niet veranderen: zij verandert alleen de relatieve fase. In combinatie met hopping beïnvloedt die fase de latere interferentie en daarmee de bezettingen.
Een alternatieve correcte implementatie brengt modi naast elkaar met fermionische swaps. Een gewone SWAP verwisselt twee qubits; een fSWAP geeft bovendien |11⟩ een minteken. Zo kan fermionische pariteit in de routing worden verwerkt in plaats van telkens als expliciete lange Z-string:
\mathrm{fSWAP}=\mathrm{SWAP}\,\mathrm{CZ},\qquad |00\rangle\mapsto|00\rangle,\quad |01\rangle\leftrightarrow|10\rangle,\quad |11\rangle\mapsto-|11\rangle.
\]
De projectcode bevat een expliciete Jordan–Wigner-route en fSWAP-routes. Beide kunnen dezelfde volledige Hamiltoniaan vertegenwoordigen. De volgorde van tijdstappen en poorten kan wel verschillen en daarmee ook de benaderingsfout. Modelgelijkheid betekent niet automatisch gelijke circuitfout.
Voor de fermionische swapnetwerken, zie Kivlichan et al. [4].
11. Uitgewerkt: twee sites en twee elektronen bij U/t = 8
Om de wisselwerking tussen t en U exact te zien, nemen we een klein hulpmodel: twee sites, twee elektronen en één hoppingverbinding. Er is hier geen diagonale t′-verbinding. Dit is een didactisch voorbeeld; de uitkomst is geen voorspelling voor ons 6 × 6-rooster of zijn begintoestand.
Kies de modusvolgorde 1↑, 1↓, 2↑, 2↓. Definieer de singlettoestand |S⟩, met één elektron op elke site, en de symmetrische doublontoestand |D₊⟩. De creatieoperatoren hieronder leggen ook de tekenconventie vast:
|S\rangle=\frac{c^\dagger_{1\uparrow}c^\dagger_{2\downarrow}-c^\dagger_{1\downarrow}c^\dagger_{2\uparrow}}{\sqrt2}|\mathrm{vac}\rangle,\qquad |D_+\rangle=\frac{c^\dagger_{1\uparrow}c^\dagger_{1\downarrow}+c^\dagger_{2\uparrow}c^\dagger_{2\downarrow}}{\sqrt2}|\mathrm{vac}\rangle.
\]
De interactie-energie van |S⟩ is nul; die van |D₊⟩ is U. Hopping mengt beide toestanden met matrixelement −2t. In deze tweedimensionale deelruimte geldt:
H_{S,D_+}=\begin{pmatrix}0&-2t\\-2t&U\end{pmatrix},\qquad \det(H-EI)=E^2-UE-4t^2=0.
\]
E_\pm=\frac{U\pm\sqrt{U^2+16t^2}}{2},\qquad \frac{E_-}{t}\bigg|_{U/t=8}=\frac{8-\sqrt{80}}{2}\simeq-0.4721.
\]
De laagste energie is negatief door de menging via hopping, ook al is de lokale afstoting positief. De kans op één doublon in deze grondtoestand is de Hellmann–Feynman-afgeleide ∂E₋/∂U, omdat ∂H/∂U de doublongetaloperator is:
\langle D_{\rm tot}\rangle=\frac{\partial E_-}{\partial U}=\frac12\left(1-\frac{U}{\sqrt{U^2+16t^2}}\right)\simeq0.0528\quad(U/t=8).
\]
Er is dus ongeveer 5,28% kans op dubbelbezetting ergens in deze dimergrondtoestand. Het gemiddelde per site is de helft daarvan. U = 8t onderdrukt dubbelbezetting hier sterk, maar maakt haar niet nul. In de limiet U ≫ t wordt E₋ ≈ −4t²/U en verschijnt de bekende superuitwisselingsschaal J ≈ 4t²/U. Die spinbeschrijving is een lage-energiebenadering; onze volledige Hubbarddynamica houdt ook de ladingsbeweging en doublontoestanden bij.
Voor de sterke-koppelingslimiet en de relatie met spinmodellen, zie Arovas et al., §5.2 [2].
12. Van Hamiltoniaan naar een eindig quantumcircuit
De exacte evolutie is exp(−ihT). De verschillende hopping- en interactietermen commuteren in het algemeen niet: AB ≠ BA. Daarom is het apart na elkaar uitvoeren van alle termen een benadering. Een symmetrische productformule verdeelt een tijdstap over een heen- en teruggaande volgorde. Voor h = a + b is het basisschema:
e^{-\mathrm{i}(a+b)\Delta\tau}=e^{-\mathrm{i}a\Delta\tau/2}e^{-\mathrm{i}b\Delta\tau}e^{-\mathrm{i}a\Delta\tau/2}+O(\Delta\tau^3),\qquad T=m\Delta\tau.
\]
De expliciete referentieroute past elke hoppingterm een halve stap toe, daarna de U-termen een hele stap, en daarna dezelfde hoppinglijst in omgekeerde volgorde een halve stap. Het handboek beschrijft voor deze target T = 0,4 met twee stappen van 0,2. T bepaalt de modeltijd; het aantal stappen bepaalt hoe die evolutie wordt benaderd. Andere runs kunnen een andere stapindeling gebruiken.
Bij vaste T is de gebruikelijke asymptotische foutorde van deze tweede-ordeformule O(Δτ²). De constante hangt onder meer af van commutatoren en modelparameters. Dit levert zonder verdere controle geen numerieke foutgarantie. Hardwarefouten komen daar nog bij. Een volledige Hamiltoniaan in het circuit betekent dus dat de juiste modeltermen aanwezig zijn; de nauwkeurigheid vraagt afzonderlijk om tijdstapverfijning en controle van de uitvoering.
De afhankelijkheid van productformulefouten van commutatoren wordt uitgewerkt door Childs et al. [5].
13. Wat voorspellen we voor het 6 × 6-rooster?
De begintoestand is een bezettingspatroon met afwisselende spins en vier verwijderde elektronen ten opzichte van één elektron per site. Er blijven 32 elektronen over: zestien up en zestien down. “Vier gaten” betekent een tekort ten opzichte van 36 elektronen, niet ten opzichte van de maximale bezetting van 72. De netto gatendoping is 4/36 ≈ 11,1%. Tijdens de evolutie kunnen extra lege sites samen met doublons ontstaan, terwijl N gelijk blijft.
N_\sigma=\sum_i n_{i\sigma},\qquad [H,N_\uparrow]=[H,N_\downarrow]=0,\qquad N_\uparrow=N_\downarrow=16.
\]
Het behoud van N = 32 is vooraf bekend. Het is een nuttige controle op de uitvoering, maar vertelt niet waar de elektronen zitten. De lokale voorspellingen zijn de lading Cᵢ, de spinonbalans Sᵢ en de dubbelbezettingskans Dᵢ. Haakjes ⟨O⟩ betekenen de verwachtingswaarde ⟨ψ(T)|O|ψ(T)⟩:
C_i=\langle n_{i\uparrow}+n_{i\downarrow}\rangle,\qquad S_i=\langle n_{i\uparrow}-n_{i\downarrow}\rangle,\qquad D_i=\langle n_{i\uparrow}n_{i\downarrow}\rangle.
\]
Hier is Sᵢ een dimensieloze spinonbalans. De fysieke z-component van spin is ℏSᵢ/2. De mogelijke bereiken zijn 0 ≤ Cᵢ ≤ 2, −1 ≤ Sᵢ ≤ 1 en 0 ≤ Dᵢ ≤ 1. Dubbelbezetting is een correlatie: Dᵢ is in het algemeen niet gelijk aan het product van de twee gemiddelde bezettingen.
\begin{aligned}C_i&=1-\frac{\langle Z_{i\uparrow}\rangle+\langle Z_{i\downarrow}\rangle}{2},\\ S_i&=\frac{\langle Z_{i\downarrow}\rangle-\langle Z_{i\uparrow}\rangle}{2},\\ D_i&=\frac{1-\langle Z_{i\uparrow}\rangle-\langle Z_{i\downarrow}\rangle+\langle Z_{i\uparrow}Z_{i\downarrow}\rangle}{4}.\end{aligned}
\]
Deze formules verbinden het model rechtstreeks met gemeten qubitbezettingen. De vraag is of het circuit deze lokale tijdsafhankelijke grootheden voldoende nauwkeurig voorspelt. Het model is geïnspireerd door gecorreleerde-elektronenfysica. Deze korte-tijdsimulatie bewijst geen supergeleiding of nauwkeurige beschrijving van een werkelijk cupratenmateriaal.
Een verwant 72-qubitexperiment is de 2D Hubbardstudie van Alam et al. op Google Willow [6]. Dat biedt experimentele context; modelparameters, voorbereiding en validatie moeten per experiment worden vergeleken.
Bronnen en controle van de modelkeuze
De projectconventies zijn gecontroleerd in theory/HANDBOOK.md, cuprate_2d_hubbard_pipeline.py, cuprate_2d_fermionic_dynamics.py en cuprate_2d_fermionic_mapping.py van fermi-hubbard-2d-nighthawk. De expliciete route draagt het modellabel full_2d_fermi_hubbard en vermeldt fermionische pariteit, t′ en U. De kleine-roostercontroles staan in validate_cuprate_2d_fermionic_dynamics.py. Dat identificeert het model en de controleprocedure; het certificeert niet automatisch elke grote hardwarerun.
Literatuur en verdere verdieping
-
J. Hubbard. Electron correlations in narrow energy bands · Vrij toegankelijke tekst. Proceedings of the Royal Society A 276, 238–257 (1963).
Oorspronkelijk artikel over het model van bewegende elektronen met lokale interacties. Historische basis voor de Hamiltoniaan in onderdeel 6; onze specifieke roostermaat en parameters zijn projectkeuzes. -
D. P. Arovas, E. Berg, S. A. Kivelson & S. Raghu. The Hubbard Model · arXiv. Annual Review of Condensed Matter Physics 13, 239–274 (2022).
Overzichtsartikel over de Hubbard-Hamiltoniaan, tweedimensionale roosters en gecontroleerde limieten. Vooral de inleiding en §5.2 over sterke koppeling helpen bij de uitleg van U/t en de superuitwisselingsschaal in onderdeel 11. -
Google Quantum AI / OpenFermion. The Jordan–Wigner and Bravyi–Kitaev Transforms.
Officiële documentatie met de anticommutatieregels, de bezettingscodering en de Jordan–Wigner-pariteitsstring. Bevat uitvoerbare voorbeelden bij de operatoren en XY-relatie in onderdelen 3, 8 en 9. -
I. D. Kivlichan et al.. Quantum Simulation of Electronic Structure with Linear Depth and Connectivity · arXiv. Physical Review Letters 120, 110501 (2018).
Onderzoeksartikel over fermionische swapnetwerken: modi worden naast elkaar gebracht met behoud van de fermionische tekens. Onderbouwt het routingprincipe in onderdeel 10; onze concrete poortenaantallen moeten uit onze eigen circuits komen. -
A. M. Childs, Y. Su, M. C. Tran, N. Wiebe & S. Zhu. Theory of Trotter Error with Commutator Scaling · arXiv. Physical Review X 11, 011020 (2021).
Onderzoeksartikel over foutgrenzen van productformules en hun afhankelijkheid van commutatoren. Verdieping bij de symmetrische tijdstappen in onderdeel 12. De theorie geeft pas een concrete foutgrens voor ons rooster na toepassing op onze Hamiltoniaan en stapgrootte. -
F. Alam et al.. Programmable digital quantum simulation of 2D Fermi-Hubbard dynamics using 72 superconducting qubits. arXiv:2510.26845v2 (2025).
Experimentele preprint over 2D Hubbarddynamica op Google Willow, tot 6 × 6 sites met 72 qubits. Vergelijk vooral vergelijking (1), de tweede-orde-Trotteraanpak en de klassieke controles. De Hamiltoniaan en beginvoorwaarden zijn niet identiek aan onze t–t′–U-target; dit artikel valideert onze IBM-resultaten niet.
De literatuur onderbouwt het model en de methoden. Onze specifieke instellingen U/t = 8, t′/t = −0,25 en T = 0,4, evenals de daadwerkelijk uitgevoerde circuits, worden vastgelegd in de hierboven genoemde projectbronnen. De twee-sitesvoorbeelden zijn didactische afleidingen met hun eigen aannames.