Fast design and scaling of multi-qubit gates in large-scale trapped-ion quantum computers

Technische Samenvatting: Snelle Ontwerp en Schaling van Multiqubit-poorten in Grootschalige Trapped-Ion Kwantumcomputers

Probleemstelling
Trapped-ion kwantumcomputers bieden langafstandsconnectiviteit en lange coherentietijden, wat hen een leidende technologie maakt voor kwantumberekeningen. Het schalen van deze systemen naar honderden ionen vormt echter een aanzienlijke conceptuele en technologische uitdaging. Hoewel de ionen all-to-all connectiviteit bezitten via hun collectieve bewegingsmodi (fononbus), is het ontwerpen van hoogwaardige, programmeerbare multiqubit verstrengelingspoorten voor grote kristallen computationeel prohibitief.

De kern van het poortontwerpprobleem ligt in het vinden van een set controle-drives die simultaan $N(N-1)/2$ kwadratische beperkingen (doelverstrengelingsfasen tussen alle qubitparen) voldoen, terwijl de drive-kracht wordt geminimaliseerd en wordt gewaarborgd dat de bewegingsmodi terugkeren naar de grondtoestand. Dit is geformuleerd als een kwadratisch beperkt optimalisatieprobleem dat kwadratisch schaalt met het aantal ionen $N$, waardoor het een NP-hard probleem is. Bestaande benaderingen beperken het probleem vaak tot kleine subsets van ionen of vertrouwen op gradiëntgebaseerde minimalisatie-algoritmen die moeite hebben met de complexiteit en de beperkingen van grootschalige systemen.

Methodologie: De Large-Scale Fast (LSF) Methode
De auteurs introduceren een methodologisch kader genaamd "Large-Scale Fast" (LSF) om de NP-harde aard van het poortontwerpprobleem te verzachten, wat het polynomiale tijdontwerp van robuuste, programmeerbare poorten voor kristallen van tot honderden ionen mogelijk maakt.

1. Formulering in het Spectrale Domein: Het poortontwerp wordt geformuleerd in het spectrale domein, waarbij de drive-velden worden gerepresenteerd als een som van Fourier-componenten die gedetuned zijn rond de qubit-transitiefrequentie. Het probleem wordt gereduceerd tot het vinden van een vector $R$ (amplitudes van drive-tonen) die de norm $\|R\|$ minimaliseert onder lineaire beperkingen (verdwijnende bewegingsverplaatsingen) en kwadratische beperkingen (doelverstrengelingsfasen $\phi_{nn'}$).
2. Zero-Phase Solution (ZPS) Seed: De LSF-methode vertrouwt op het genereren van een "seed" oplossing bekend als de Zero-Phase Solution (ZPS). Een ZPS is een niet-triviale oplossing van het homogene probleem waarbij alle doelverstrengelingsfasen nul zijn ($\phi^{(0)}_s = 0$). Deze oplossingen worden geaccumuleerd voorafgaand aan het oplossen van specifieke doelpoorten.
3. Conversie en Optimalisatie: Het algoritme verloopt in twee hoofdfasen:
   • Conversie: Een ZPS wordt geschaald met een grote factor $\lambda$ en verstoord door een kleine vector $D$ om aan de kwadratische doelbeperkingen te voldoen. Door de beperkingen te lineariseren rond de geschaalde ZPS, leiden de auteurs een lineair systeem af om voor $D$ op te lossen, wat de ZPS converteert naar een benaderde oplossing $R_{conv}$ die de doelfasen voldoet.
   • Normreductie: Een iteratief proces minimaliseert de norm van de oplossingsvector terwijl strikt aan de beperkingen wordt voldaan. Dit omvat het afwisselen tussen "norm-reducerende" stappen (bewegen langs de oplossingsmanifold om vermogen te reduceren) en "fout-reducerende" stappen (het corrigeren van afwijkingen van de doelfasen). De stapgroottes worden analytisch berekend met behulp van polynomiale wortelzoekmethoden (Cardano's en Ferrari's methoden) om te garanderen dat de fout binnen een gespecificeerde tolerantie blijft (doorgaans $10^{-4}$).
4. Aggregatie: De methode aggregeert meerdere ZPS'en om een diverse zoekruimte te waarborgen, waardoor het algoritme efficiënte oplossingen kan vinden voor diverse verstrengelingskaarten.

Belangrijkste Bijdragen en Resultaten

• Polynomiale Schaling: De LSF-methode transformeert het NP-harde ontwerpprobleem in een taak van het genereren van ZPS'en (wat efficiënt kan worden gedaan voor specifieke symmetrieën), gevolgd door efficiënte polynomiale tijd optimalisatie. Benchmarktests laten zien dat LSF een tijd-tot-oplossing bereikt die ongeveer een orde van grootte sneller is dan standaard trust-region constrained en conjugate-gradient optimizers, met name voor systemen met $N > 20$ ionen.
• Drive-kracht Schatting: De auteurs leiden een heuristische estimator af voor de totale Rabi-frequentie ($\Omega_{nuc}$) die nodig is om een poort aan te drijven. Deze estimator hangt af van de nucleaire norm van de doelverstrengelingsmatrix $\|\phi\|_{nuc}$ en schaalt lineair met het aantal ionen $N$. Opmerkelijk genoeg laten simulaties zien dat de vereiste drive-kracht grotendeels onafhankelijk is van de specifieke modusstructuur of frequentieafstand, en convergeert naar een universele curve wanneer deze wordt genormaliseerd door de gemiddelde modusdensiteit tijdsschaal $T_\nu$.
• Minimale Poortduur ($T_{min}$): De studie identificeert een fundamentele ondergrens op de poortduur. Onder een kritieke tijd $T_{min} \approx 2.9 \cdot T_\nu$ (waarbij $T_\nu$ proportioneel is aan de inverse van de gemiddelde modusafstand), kunnen geen hoogwaardige oplossingen worden gevonden. Cruciaal is dat $T_{min}$ lineair schaalt met het aantal ionen $N$. Dit impliceert dat hoewel het aantal gelijktijdige verstrengelingsoperaties schaalt als $N^2$, de poortduur slechts schaalt als $N$, wat een schaalvoordeel biedt ten opzichte van sequentiële poortimplementaties.
• Foutschalingsanalyse: Het artikel kwantificeert de impact van drie dominante foutbronnen op de poortgetrouwheid (fidelity):
  • Drift van bewegingsmodusfrequenties: Fouten schalen kwadratisch met de frequentiedrift en het aantal koppelingen ($N^2$).
  • Fluctuaties in drive-amplitude: Fouten schalen met de variantie van de ruis en het aantal koppelingen.
  • Bewegingsverhitting (Motional heating): De schaling van door verhitting veroorzaakte fouten hangt af van de ruimtelijke correlatie van de ruis. Voor ongecorreleerde ruis neemt de getrouwheid kwadratisch af met $N$; voor sterk gecorreleerde ruis neemt deze lineair af.
    Het LSF-kader staat toe om lineaire beperkingen op te nemen om poorten te robuuster te maken tegen deze fouten, wat de gevoeligheid voor drift en verhitting aanzienlijk vermindert.
• Surface Code Voorbeeld: De methode wordt gedemonstreerd op een 49-ion kristal om een surface-code stabilizer poort te implementeren (het verstrengelen van 33 ionen in een enkele puls). De oplossing vertoont een efficiënte machtsverdeling, waarbij hogere drive-amplitudes worden toegewezen aan "node" ionen en lagere aan "edge" ionen, en onderdrukt effectief ongewenste cross-verstrengeling tussen niet-verbonden paren.

Betekenis en Claims
Het artikel beweert een pad te hebben uitgestippeld naar het schalen van trapped-ion kwantumcomputers naar honderden qubits door de conceptuele barrière van het ontwerpen van multiqubit-poorten te overwinnen. Door een snelle, robuuste algoritme (LSF) te bieden, stellen de auteurs het offline compileren van circuits met tientallen ionen en complexe multiqubit-poorten mogelijk, een taak die voorheen als onoverkomelijk werd beschouwd vanwege de NP-harde aard van de optimalisatie.

Het werk stelt vast dat hoogwaardige multiqubit-poorten fysiek haalbaar zijn voor grote kristallen, mits de poortduur de lineair schalende drempel $T_{min}$ overschrijdt. Bovendien bieden de afgeleide relaties tussen drive-kracht, poortduur en foutbronnen praktische ontwerprichtlijnen en toleranties voor experimentatoren. De auteurs merken op dat hoewel hun methode de dominante foutmechanismen binnen de de Lamb-Dicke benadering behandelt, interacties met qubit-dragers en hogere-orde termen in de Magnus expansie voor hoge-amplitude oplossingen verdere verificatie vereisen. Het kader wordt gepresenteerd als een fundamenteel instrument voor het realiseren van het potentieel van langafstandsconnectiviteit in trapped-ion architecturen voor kwantumfoutcorrectie en optimalisatie-toepassingen.