Finite size scaling of bitstring probability distributions for Rydberg arrays

Technische Samenvatting: Eindelijke Grootte Schaling van Bitstring-waarschijnlijkheidsverdelingen voor Rydberg-arrays

Probleemstelling
In kwantumsimulaties van deeltjes-, kern- en gecondenseerde materiefysica vereist het schatten van observabelen met statistische nauwkeurigheid een groot aantal metingen (shots). Naarmate de gesimuleerde systeemgrootte ($N_q$, het aantal qubits) toeneemt, groeit het aantal mogelijke bitstring-toestanden exponentieel, waardoor individuele toestandswaarrscheinlichheden ($p_{\{n\}}$) afnemen en dichter worden in het regime van lage waarschijnlijkheid. Dit roept een kritische vraag op: kan de kleinheid van individuele waarschijnlijkheden worden gecompenseerd door hun enorme aantal, en wat is het minimale aantal shots dat vereist is om een gewenste statistische nauwkeurigheid voor grotere systemen te bereiken? Voorgaand werk suggereert dat het karakteriseren van systemen via grootschalige verdelingen structurele gedragingen kan onthullen, maar expliciete schalingswetten voor bitstring-waarschijnlijkheden in Rydberg-atoomarrays zijn nog niet volledig gekwantificeerd.

Methodologie
De auteurs onderzoeken Rydberg-atoomarrays die zijn gerangschikt in een ladderconfiguratie (twee rijen met een verticale-naar-horizontale afstandratio van $\rho=2$). Het systeem wordt beschreven door een Hamiltoniaan met Rabi-frequentie ($\Omega$), detuning ($\Delta$) en van der Waals-interacties ($V_{ij}$).

Om de waarschijnlijkheidsverdelingen te analyseren, maken de auteurs gebruik van Matrix Product States (MPS) berekend via Density Matrix Renormalization Group (DMRG)-methoden. Ze genereren woordenboeken van bitstring-waarschijnlijkheden ($p_{\{n\}}$) voor systemen variërend van 4 tot 22 rungs (tot 44 atomen). De kern van hun analyse betreft de cumulatieve waarschijnlijkheidsverdeling, gedefinieerd als:
$$\Sigma(p_\Lambda, N_q) = \sum_{\{n\}: p_{\{n\}} \leq p_\Lambda} p_{\{n\}}$$
Deze functie vertegenwoordigt de totale waarschijnlijkheid van het observeren van een toestand met een waarschijnlijkheid die kleiner dan of gelijk is aan een drempelwaarde $p_\Lambda$. De auteurs analyseren het schalingsgedrag van $\Sigma$ over verschillende systeemgroottes ($N_q$) en shot-totalen ($N_{sh} = 10^9$).

Belangrijkste Bijdragen en Resultaten

1. Eindelijke Grootte Schaling en Curve Collaps:
   De auteurs demonstreren dat cumulatieve waarschijnlijkheidsverdelingen voor verschillende systeemgroottes een "herhaald gedrag" vertonen, met name in het bereik $0.2 < \Sigma < 0.8$. Door een karakteristieke waarschijnlijkheid $p^*$ te identificeren (waar de verdeling een vaste $\Sigma$-waarde bereikt) en de onafhankelijke variabele te herschalen als $(p_\Lambda/p^*)^B$ (waarbij $B$ de helling van het lineaire gebied is), vallen de verdelingen voor diverse $N_q$ samen op één universele curve.

2. Fermi-Dirac Analogie (Regio II):
   In het centrale gebied van de verdeling (Regio II) lijkt de samengevallen curve nauw overeen te komen met een Fermi-Dirac statistische verdeling. De auteurs tonen aan dat $\Sigma(p_\Lambda)$ benaderd kan worden door:
   $$\Sigma \simeq \frac{1}{1 + e^{\beta(\epsilon - \mu)}}$$
   waarbij $\epsilon = -\ln(p_\Lambda)$ en $\mu = -\ln(p^*)$. Deze analogie maakt het mogelijk om het "inflectiepunt" (waar $\Sigma = 1/2$, aangeduid als $p^*_{1/2}$) te karakteriseren. De schaling van $p^*_{1/2}$ met de systeemgrootte blijkt exponentieel te zijn, wat aangeeft dat naarmate $N_q$ toeneemt, de waarschijnlijkheidsmassa meer uniform verdeeld raakt over een enorm aantal toestanden.

3. Schaling van Maximale Waarschijnlijkheden (Regio I):
   In Regio I (toestanden met hoge waarschijnlijkheid) volgen de auteurs de maximale enkele-toestandswaararschijnlijkheid ($p_{\max}$). Ze bevestigen dat $p_{\max}$ exponentieel schaalt met de systeemgrootte ($p_{\max} \sim 2^{-N_q}$). Dit stelt een fundamentele limiet aan de toegankelijke systeemgrootte voor een gegeven shot-totaal vast, aangezien de minimaal oplosbare waarschijnlijkheid begrensd wordt door $1/N_{sh}$.

4. Lage Waarschijnlijkheid en Single-Occurrence Toestanden (Regio III & IV):

   • Regio III: Bij lagere waarschijnlijkheden volgen de verdelingen een machtswetvorm $\Sigma(p_\Lambda) \simeq A p_\Lambda^B$ wanneer ze op een log-log schaal worden weergegeven. De auteurs merken op dat het gedrag hier sterk wordt beïnvloed door het eindige shot-totaal ($10^9$), wat een harde afkapwaarde creëert bij $p \approx 10^{-9}$.
   • Regio IV: Deze regio vertegenwoordigt de "staart" van de single-occurrence toestanden. De auteurs kwantificeren het aantal toestanden dat exact één keer wordt waargenomen, en tonen aan dat hoewel dit aantal groeit met de systeemgrootte, het langzamer groeit dan de totale Hilbert-ruimte dimensie ($2^{N_q}$).

5. Schatting van Shot-vereisten:
   Door de schaling van $p^*_{1/2}$ en $p_{\max}$ te analyseren, schatten de auteurs de shot-vereisten voor toekomstige simulaties.

   • Op basis van de $p_{\max}$ schaling zou een systeemgrootte van ongeveer 120 qubits gesimuleerd kunnen worden met $10^9$ shots.
   • Op basis van de $p^*_{1/2}$ schaling (waar de cumulatieve waarschijnlijkheid 0,5 bereikt), wordt de limiet voor een specifieke nauwkeurigheid geschat op ongeveer 60 qubits voor hetzelfde shot-totaal.

Betekenis en Claims
Het artikel claimt een methode te bieden voor het voorspellen van noodzakelijke shot-vereisten voor grotere kwantumsystemen door de cumulatieve waarschijnlijkheidsverdelingen van kleinere, computationeel toegankelijke systemen te analyseren. De primaire betekenis ligt in de identificatie van een exponentiële kostenpost die gepaard gaat met het vergroten van de systeemgrootte, welke wordt gekwantificeerd door de schaling van $p^*_{1/2}$ en $p_{\max}$.

De auteurs zijn echter bescheiden in hun conclusies; zij merken op dat hoewel de kosten exponentieel zijn, deze met een fractie van het tempo van de toename van de systeemgrootte zelf toenemen. Dit vormt dus geen "harde stop" voor het gebruik van apparaten, maar motiveert eerder het "intelligente gebruik van middelen". Het werk vestigt dat het gedrag van grote systemen beschreven kan worden door de verdelingen van kleinere systemen via de voorgestelde collaps en de Fermi-Dirac analogie, en biedt hiermee een instrument om statistische fouten en de haalbaarheid van het simuleren van steeds grotere Rydberg-arrays in te schatten.