Quantum approximate optimization of finite-state bosonic systems

Probleemstelling

Veel natuurkundige problemen, zoals die in de chemie en materiaalkunde, worden inherent beschreven door systemen met een eindige, maar hoge dimensie $D$ (qudits). Om deze problemen op qubit-gebaseerde quantumhardware uit te voeren, moet de Hilbert-ruimte van de qudits worden gemapt naar een multi-qubit ruimte. Omdat het aantal qubits ($K$) vaak logaritmisch groeit ten opzichte van $D$ (bijv. $K \approx \log_2 D$), is de totale Hilbert-ruimte van de qubits ($2^K$) exponentieel groter dan de fysieke ruimte van het probleem ($D$).

Dit creëert een onhaalbaar deelruimte (infeasible subspace) die niet overeenkomt met fysieke toestanden. Een gebruikelijke aanpak is het toevoegen van straftermen (penalty terms) aan de doelwitfunctie om deze toestanden te onderdrukken. De auteur stelt echter dat deze strategie inefficiënt wordt naarmate de grootte van het illegitieme deelruimte exponentieel toeneemt, waardoor het zoeken naar een oplossing computatiewerkzaam wordt.

Methodologie

Het artikel introduceert een aanpak binnen het kader van het Quantum Approximate Optimization Algorithm (QAOA), specifiek de Hamiltonia-gebaseerde variant (ook wel Quantum Alternating Operator Ansatz genoemd). In plaats van straftermen te gebruiken, wordt de zoekruimte beperkt tot de haalbare deelruimte door het ontwerpen van specifieke mixing Hamiltonians ($\hat{H}_M$).

De kern van de methode is het zorgvuldig construeren van parameteriseerbare quantumcircuits die de dynamica binnen de toegestane toestanden houden. De auteur onderzoekt dit door drie verschillende coderingsmethoden (mappings) te vergelijken:

1. Binary Encoding: Een hardware-efficiënte mapping waarbij qudits worden vertaald naar binaire variabelen.
2. Symmetric Encoding: Een mapping waarbij elke basisstaat $|d\rangle$ wordt gemapt naar een superpositie van toestanden met een Hamming-gewicht van precies $d$.
3. Unary Encoding: Een mapping waarbij een qudit-toestand $d$ wordt gemapt naar een register van $D$ qubits met precies één qubit in de toestand $|1\rangle$.

De prestaties worden geëvalueerd op basis van de CNOT-gate telling (aantal entangling gates) als maatstaf voor implementatiekosten, aangezien deze gates duurder zijn in uitvoering en foutgevoeliger dan single-qubit gates.

Belangrijkste Bijdragen

1. Ontwerp van Constraint-Preserving Mixers: De auteur toont aan dat het mogelijk is om specifieke mixing Hamiltonians te ontwerpen die de haalbare deelruimte invariant houden, waardoor onmogelijke toestanden fysiek worden uitgesloten zonder straftermen.
2. Vergelijking van Coderingsmethoden: Er wordt een gedetailleerde analyse gemaakt van de resourcekosten (qubits en CNOT-gates) voor de drie coderingsmethoden.
3. Optimaliteit van Symmetrische Codering: Het belangrijkste theoretische resultaat is dat voor symmetrische codering de standaard mixing Hamiltonian (de som van Pauli-X operatoren, $\sum X_k$) de optimale keuze is. Deze operator behoudt het symmetrische deelruimte en vereist geen entangling gates (CNOTs) voor de implementatie van de mixer, omdat deze factoriseert in single-qubit operaties.
4. Toepassing op Quantum Thermalisatie en Bose-Hubbard: De framework wordt toegepast op twee specifieke casestudies: het genereren van benaderende thermische toestanden en het vinden van de grondtoestand van het Bose-Hubbard-model.

Resultaten

• CNOT-efficiëntie:
  • Voor symmetrische codering is de standaard mixer optimaal: deze vereist 0 CNOT-gates voor de mixing-stap (alleen single-qubit rotaties).
  • Voor binary en unary codering neemt het aantal benodigde CNOT-gates exponentieel toe met de systeemgrootte $D$. Voor een $p$-laags QAOA neemt het aantal CNOT-gates voor deze methoden $p$-voudig toe, wat leidt tot een aanzienlijk hogere implementatiekosten in vergelijking met de symmetrische methode.
• Quantum Approximate Thermalisatie:
  • Bij het simuleren van twee gekoppelde harmonische oscillatoren bleek de symmetrische codering superieur. De gegenereerde toestanden hadden een hogere fideliteit (0.93) en een lagere quantum relatieve entropie ten opzichte van de exacte Gibbs-toestand, vergeleken met de binary codering (fideliteit 0.89).
  • De symmetrische methode bleek robuuster tegen temperatuursvariaties.
• Bose-Hubbard Model:
  • In het sterk-interagerende regime (waar de grondtoestand een producttoestand is) convergeren beide coderingen snel naar de grondtoestand met hoge fideliteit.
  • In het kinetisch gedomineerde regime (waar de grondtoestand sterk verstrekt is) is een dieper circuit (meer QAOA-lagen) nodig om de complexiteit van de verstrengeling te vangen. Hoewel de binary codering hier iets beter presteerde in termen van fideliteit bij zeer diepe circuits, blijft de symmetrische codering aantrekkelijk vanwege de lagere gate-kosten per laag.

Betekenis en Conclusie

De studie demonstreert dat het kiezen van de juiste codering en mixing Hamiltonian cruciaal is voor de efficiëntie van variational quantum algoritmen op NISQ-hardware.

• Efficiëntie: De symmetrische codering, ondanks het gebruik van meer qubits dan de binary methode, is over het algemeen de meest efficiënte strategie voor systemen met hoge dimensies. Dit komt omdat de standaard mixer geen entangling gates vereist, wat de diepte van het circuit en de gevoeligheid voor ruis drastisch verlaagt.
• Praktische Toepassing: De bevindingen bieden een blauwdruk voor het simuleren van bosonische systemen (zoals in de quantum chemie en materiaalkunde) op huidige quantumcomputers. Het vermijden van straftermen en het minimaliseren van CNOT-gates maakt de methode beter geschikt voor de beperkte coherentietijden van huidige hardware.
• Toekomstperspectief: De auteur suggereert dat de grote Hilbert-ruimte van symmetrische en unary coderingen ook nuttig kan zijn voor foutcorrectie en -mitigatie in de toekomst, en dat de methode een brug slaat tussen theoretische optimalisatie en praktische implementatie op quantumhardware.

Kortom, het artikel pleit voor een verschuiving van "straftermen" naar "geconstrueerde mixers" en identificeert de symmetrische codering als de superieure keuze voor het optimaliseren van eindig-toestand bosonische systemen binnen het QAOA-kader.