Analytical Series Expansion for Efficient Gradient Evaluation in Multi-Qubit Optimal Control

Technische Samenvatting: Analytische Reeksontwikkeling voor Efficiënte Gradiëntevaluatie in Multi-Qubit Optimale Controle

Probleemstelling
Quantum optimale controle (QOC) is essentieel voor het kalibreren van grootschalige quantumprocessors, met name voor het mitigeren van fouten die voortvloeien uit spuria true interacties en crosstalk in multi-qubit systemen. Ho'ewel gradiëntgebaseerde QOC-methoden (bijv. GRAPE, Krotov) efficiënter zijn dan gradiëntvrije benaderingen, kampen zij met aanzienlijke computationele knelpunten bij het schalen naar grote Hilbert-ruimten of open quantum systemen. De kernuitdaging ligt in het berekenen van de gradiënt van de tijd-geordende propagator, $U(t)$, ten opzichte van willekeurige pulsparameters. Bestaande methoden zoals GOAT (Gradient Optimization of Analytic conTrols) vereisen het co-propageren van een set afgeleiden naast de propagator, wat de keuzes voor parametrisatie beperkt of hoge computationele kosten met zich meebrengt door de herhaalde evaluatie van matrix-exponenten. Piecewise-constant ansätze (gebruikt in GRAPE) omzeilen sommige beperkingen, maar missen de flexibiliteit van continue parametrisaties.

Methodologie
De auteurs presenteren een verenigend kader afgeleid uit eerste principes om de gradiënt van de propagator te berekenen voor willekeurige puls-parametrisaties. De methodologie verloopt in drie hoofdfasen:

1. Formele Oplossing via de Heisenberg-representatie:
   Beginnend vanuit de Schrödinger-vergelijking, leiden de auteurs een exacte integraalexpressie af voor de gradiënt van de propagator, $\nabla_\alpha U(t)$. Deze expressie onthult dat de gradiënt afhankelijk is van de evolutie van de afgeleide van de Hamiltoniaan in de Heisenberg-representatie:
   $$\frac{\partial U(t)}{\partial \alpha_k} = -i U(t) \int_0^t d\tau \frac{\partial u_m(\tau)}{\partial \alpha_k} H_m^{(H)}(\tau)$$
   waarbij $H_m^{(H)}(\tau) = U^\dagger(\tau) H_m U(\tau)$. Deze formulering verbindt de gradiëntberekening met de studie van operator-evolutie, een centraal onderwerp in velden zoals operator spreading en quantum chaos.

2. Reeksontwikkeling:
   Om de dure numerieke integratie van de Heisenberg-operator bij elke tijdstap te vermijden, integreren de auteurs de expressie herhaaldelijk door delen (integration by parts). Dit levert een afgekorte reeksontwikkeling op:
   $$\frac{\partial U(t)}{\partial \alpha_k} \simeq \left[ \sum_{n=0}^N \Theta_n(t) \right] U(t)$$
   Hierbij bestaat $\Theta_n(t)$ uit tijd-onafhankelijke commutatoren van de Hamiltoniaan-termen (bijv. $[H_{m'}, [H_{m''}, H_m]]$) vermenigvuldigd met tijd-afhankelijke coëfficiënten $\beta^{(n)}_{m,m}(t)$. Cruciaal is dat de commutatoren eenmaal aan het begin van de simulatie worden berekend, terwijl alleen de coëfficiënten (die eenvoudige klassieke integralen zijn) tijdens de optimalisatie moeten worden bijgewerkt.

3. Benutten van Localiteit:
   Voor systemen met lokale interacties (bijv. 1D-ketens of 2D-roosters), maken de auteurs gebruik van de grafenstructuur van de Hamiltoniaan. Door Lieb-Robinson-grenzen toe te passen, demonstreren zij dat de invloed van lokale operatoren zich met een eindige snelheid verspreidt. Dit maakt de afkorting van de reeksontwikkeling mogelijk, aangezien commutatoren die betrokken zijn bij verre qubits verwaarloosbaar worden. De methode maakt gebruik van ijle (sparse) representaties (Pauli-strings) en heuristische pruning-strategieën om het aantal termen dat vereist is verder te reduceren, wat de benadering schaalbaar maakt voor grote multi-qubit platformen.

Belangrijkste Bijdragen

• Verenigend Kader: Het artikel biedt een formele oplossing voor de gradiënt van de propagator die van toepassing is op elke puls-parametrisatie (continu of piecewise-constant) en kan worden uitgebreid naar open systemen en hogere-orde afgeleiden.
• Generalisatie van GRAPE: De afgeleide reeksontwikkeling wordt getoond als een generalisatie van de exacte gradiënt gebruikt in GRAPE, waarbij de toepasbaarheid wordt uitgebreid naar continue puls-vormen.
• Computationele Efficiëntie: Door de tijd-onafhankelijke commutatoren te ontkoppelen van de tijd-afhankelijke coëfficiënten, vermindert de methode aanzienlijk het aantal vereiste matrix-exponenten.
• Schaalbaarheid: De integratie van localiteitsbeperkingen en ijle Pauli-dynamica maakt de methode in staat om grote Hilbert-ruimten te verwerken waar traditionele co-propagatie methoden falen.

Resultaten
De auteurs valideren hun methode via numerieke simulaties op multi-qubit systemen:

• Nauwkeurigheid: Vergelijkingen met de GOAT-methode op een twee-qubit systeem met een CRAB-puls ansatz laten zien dat de reeksontwikkeling (zowel dense als sparse implementaties) de exacte gradiënt met hoge nauwkeurigheid reproduceert (relatieve fout $< 10^{-4}$ voor de dense implementatie).
• Prestaties: In simulaties van qubit-ketens en ladder-geometrieën (tot 10 qubits) vertoont de reeksontwikkelingsmethode een versnelling van een orde van grootte in wandtijd (wall time) en een significante reductie in geheugengebruik vergeleken met GOAT.
• Toepassing: De methode werd succesvol toegepast voor het prepareren van een 4-qubit GHZ-toestand op de centrale qubits van een 2D-ladder-geometrie met 8 qubits (inclusief spectatorkubits) en afleidende ZZ-interacties. De optimalisatie, bestaande uit 264 parameters, bereikte een fidelity van 99,95%, wat het vermogen van de methode aantoont om complexe control-landschappen met honderden parameters aan te kunnen.

Betekenis
Het artikel stelt dat dit werk de kloof overbrugt tussen quantum optimale controle en de studie van operator-evolutie. Door de gradiëntberekening in kaart te brengen naar de evolutie van Hamiltoniaan-afgeleiden in de Heisenberg-representatie, stellen de auteurs het gebruik van gevestigde instrumenten uit de theorie van operator spreading (zoals tensor netwerken en sparse Pauli propagatie) mogelijk om quantum-systemen te optimaliseren. De auteurs benadrukken dat deze benadering een praktisch pad biedt voor het schalen van QOC naar grotere quantumprocessors, waarbij het computationele overhead dat momenteel het ontwerp van hoog-frequente pulsen voor complexe, interagerende multi-qubit systemen beperkt, wordt aangepakt. Het werk motiveert tevens toekomstige ontwikkelingen in ijle representaties voor operator- en staat-evolutie om deze optimalisatietechnieken verder te schalen.