Technische Samenvatting: QISKIT-QUDITS

Probleemstelling

Qudits, die twee-niveau qubits generaliseren naar $d$ computationele niveaus, bieden potentieel voordelen in kwantuminformatieverwerking. De Qiskit-framework mist echter native ondersteuning voor qudits, wat een bekende kloof in de derde partij extensies vormt. Bestaande toolchains kampen met een dichotomie: native qudit-simulatoren (bijv. Cirq, QuTiP) kunnen de uitgebreide compilatie-stack, transpiler-passes en hardware-backends van Qiskit niet benutten, terwijl qubit-gerichte tools vaak een gebrek hebben aan een programmeerinterface voor operaties op qudit-niveau. Gebruikers moeten momenteel handmatig encoding-schema's, gate-matrices, decomposities en klassieke decodering implementeren, een proces dat foutgevoelig is met betrekking tot de behandeling van de onfysische toestandsruimte en meting-decodering.

Methodologie

Het artikel introduceert qiskit-qudits, een Qiskit-extensie die $d$-niveau qudits simuleert door elke qudit te encoderen in $m = \lceil \log_2 d \rceil$ qubits. De software is gebouwd op de publieke Instruction, Gate, en ControlledGate klassen van Qiskit, wat compatibiliteit met het bestaande ecosysteem garandeert.

Kernontwerpbeslissingen

1. Encoding en Identity Padding:
   Omdat $d$ niet noodzakelijkerwijs een macht van twee is, is de geëncodeerde Hilbertruimte ($D = 2^N$) over het algemeen groter dan de logische subruimte ($L = \prod d_i$). De bibliotheek hanteert de identity-padding conventie: elke gate fungeert als de identiteit op het orthogonale complement van de logische subruimte. Formeel, voor een logische gate $G$, is de geëncodeerde unitaire operator $U = E G E^\dagger + (I - \Pi)$, waarbij $E$ de encoding-isometrie is en $\Pi$ de projector op de logische subruimte. Deze keuze voorkomt lekkage tijdens foutloze simulatie en behoudt permutatie-eigenschappen.

2. Gate Hiërarchie en Decompositie:
   De bibliotheek definieert drie basisklassen: QuditGate (enkelvoudige qudit), QuditControlledGate (multi-control), en QuditMultiGate.

   • Macht-van-Twee Pad: Wanneer alle operanden-dimensies machten van twee zijn, decomponeren gates in vaste cascades van standaard Qiskit qubit-gates (bijv. mcx, cp, synth_qft_full). Dit stelt de Qiskit transpiler in staat om de circuit te optimaliseren.
   • Dense-Synthesis Fallback: Voor dimensies die geen macht van twee zijn, valt de bibliotheek terug op exacte dense unitaire synthese. De bibliotheek genereert de volledige $D \times D$ unitaire matrix (via het NumPy __array__ protocol) en embedt deze als een enkele opaak unitaire instructie, die de UnitarySynthesis pass vervolgens decomponeert.

3. Circuit Model en Meting:
   De QuditQuantumCircuit houdt twee gesynchroniseerde weergaven bij: een instructielogboek op qudit-niveau (.data) en een geëncodeerd qubit-niveau circuit (.circuit). Meting wordt afgehandeld via QuditMeasure, die encoding-qubits naar klassieke bits (ClByte) mapt. Een toegewezen decoderinglaag vertaalt backend-counts terug naar qudit-niveaus, waarbij "lekkage" (uitkomsten $\ge d$) expliciet wordt gedetecteerd, wat afgehandeld kan worden via configureerbare policies (keep, drop, raise).

Belangrijkste Bijdragen

• Implementatie van een Transpilbare Qudit Laag: Het werk zet de theoretische gate-definities uit eerder werk [1] om in een werkende, transpilbare Qiskit-extensie. Het ondersteunt dimensies $2 \le d \le 16$ en gemengde dimensie-registers.
• Exacte Simulatie voor Niet-Machten-van-Twee Dimensies: In tegen tegenstelling tot eerdere qubit-simulatiebenaderingen die beperkt zijn tot $d=2^m$, ondersteunt deze bibliotheek exacte synthese voor willekeurige $d$ binnen het ondersteunde bereik met behulp van dense unitaire fallbacks.
• Verificatie van Theorie en Implementatie:
  • Check A: Geverifieerd dat de dense unitaire matrices de Discrete Fourier Transform (DFT) reproduceren op de logische subruimte voor zowel macht-van-twee als niet-macht-van-twee dimensies, waarmee de verificatie van [1] wordt uitgebreid.
  • Check B: Gevalideerd dat de gedecomponeerde circuits (zowel cascade als dense paden) overeenkomen met de dense unitaire definities over 138 gate-instanties, met afwijkingen binnen machineprecisie ($\approx 10^{-15}$).
  • Padding Invariant: Bevestigd dat de identity-padding invariant exact standhoudt, wat zorgt voor geen vermenging tussen geldige en ongeldige subruimtes.
• Lekkage Detectie: Biedt een expliciet mechanisme om lekkage-events te detecteren en af te handelen in simulaties, waarbij onderscheid wordt gemaakt tussen ideaal gedrag en door ruis geïnduceerde fouten.

Resultaten

Numerieke verificatie bevestigt de correctheid van de implementatie:

• QFT vs. DFT: De QFT-implementatie komt overeen met de theoretische DFT-matrix binnen $10^{-15}$ voor dimensies $d \in \{2, 3, 4, 5, 8\}$.
• Gate Decompositie: Alle gates, inclusief SUMX, SUMP, SWAP, en QFT, produceren decomposities die consistent zijn met hun dense unitaire definities.
• Resource Analyse:
  • Macht-van-Twee: Decomposities zijn efficiënt en transpiler-vriendelijk. Bijvoorbeeld, een 3-qutrit QFT ($d=3$) vereist 1.783 CX-gates, terwijl een 3-ququart Qটি ($d=4$) slechts 56 CX-gates vereist.
  • Niet-Macht-van-Twee: De dense fallback brengt aanzienlijke overhead met zich mee. Een 3-qutrit QFT ($d=3$) synthetiseert naar 1.783 CX-gates, terwijl een 3-ququint ($d=5$) QFT 119.383 CX-gates vereist. Dit onderstreept de kosten van dense synthese voor grotere niet-macht-van-twee dimensies.
  • Geheugenbeperkingen: Het dense pad wordt beperkt door geheugen; een gate op acht 16-niveau qudits ($N=32$) zou $\approx 256$ EiB geheugen vereisen, waardoor de praktische limiet $N \approx 13\text{--}14$ is op standaard hardware.

Betekenis en Beperkingen

Het artikel positioneert qiskit-qudits als een brug die gebruikers in staat stelt om qudit-algoritmen op logisch niveau te schrijven terwijl zij volledige toegang behouden tot de Qiskit-compilatie-stack, backends en primitives. Het adresseert expliciet het "onfysische deel" van de geëncodeerde ruimte en meting-decodering, wat vaak foutgevoelig is in handmatige implementaties.

Beperkingen en Trade-offs:

• Dimensie Cap: De bibliotheek dwingt $2 \le d \le 16$ af en een maximum van 8 qudits per gate. Dit zijn nominale limieten; de praktische limiet wordt bepaald door het geheugen dat vereist is voor dense unitaire synthese ($N \lesssim 13\text{--}14$).
• Lekkage Semantiek: De identity-padding conventie bevriest de gelekte amplitude in plaats van fysieke lekkage te modelleren. Dit commuteert niet met tensoring, wat betekent dat circuit-equivalentie gecontroleerd moet worden op de geprojecteerde subruimte ($\Pi U \Pi$).
• Optimalisatie: De huidige implementatie maakt geen gebruik van "don't-care" rijen/kolommen in de transpiler, en de dense fallback voor niet-macht-van-twee dimensies is computationeel duur vergeleken met native qudit-simulatoren of geoptimaliseerde reversibele logica-synthese.

Het werk concludeert dat hoewel de dense fallback exacte simulatie mogelijk maakt voor niet-macht-van-twee dimensies, toekomstige verbeteringen zich moeten richten op dimensie-agnostische reversibele logica-synthese en het blootstellen van de onfysische subruimte aan de transpiler om de resource-overhead te verminderen.