⚛️ quantum physics

qiskit-qudits: A Qiskit Extension for Simulating Qudit Circuits

Cet article présente **qiskit-qudits**, une extension de Qiskit qui simule des qudits arbitraires à $d$ niveaux en les encodant dans des qubits avec un remplissage d'identité pour gérer les dimensions qui ne sont pas des puissances de deux, offrant à la fois des décompositions de portes standards et une synthèse unitaire dense exacte pour des dimensions allant jusqu'à 16, tout en validant son implémentation par la vérification numérique de la transformée de Fourier quantique et des décompositions de portes.

Auteurs originaux : Francesco Pudda

Publié 2026-08-25

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Francesco Pudda

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Dans le monde de l'informatique quantique, la brique élémentaire standard est le qubit, un système minuscule qui peut exister dans deux états à la fois, souvent décrits comme zéro et un. Mais la nature ne se limite pas à seulement deux options. De nombreux systèmes physiques, tels que certains atomes ou ions, possèdent naturellement trois, quatre ou même plus de niveaux d'énergie distincts. Les scientifiques appellent ces systèmes multi-niveaux des qudits. Si les qubits sont le fondement des ordinateurs quantiques actuels, les qudits offrent un paysage plus riche pour le calcul, permettant potentiellement des algorithmes plus efficaces et une meilleure utilisation du matériel physique. Cependant, un obstacle important est resté : les outils logiciels les plus populaires pour concevoir des programmes quantiques ont été construits exclusivement pour les qubits à deux niveaux. Cela signifiait que les chercheurs souhaitant expérimenter avec les qudits devaient soit écrire leur propre logiciel complexe à partir de zéro, soit abandonner les puissants outils développés par la communauté élargie.

Francesco Pudda a comblé cette lacune avec une nouvelle extension logicielle appelée qiskit-qudits. Cet outil permet aux scientifiques de concevoir des circuits quantiques utilisant des qudits multi-niveaux tout en utilisant toujours la suite logicielle standard Qiskit, qui est largement utilisée pour construire et tester des programmes quantiques. Le cœur de cette réalisation est une méthode de traduction ingénieuse. Puisque le matériel sous-jacent et le logiciel standard pensent en termes de qubits à deux niveaux, le nouveau système encode chaque qudit multi-niveaux dans un petit groupe de qubits. Par exemple, un système à trois niveaux est représenté par deux qubits, et un système à cinq niveaux par trois. Cette approche crée une légère inefficacité car le groupe de qubits peut représenter plus d'états que le qudit n'en utilise réellement ; mais le logiciel est conçu pour ignorer ces états supplémentaires inutilisés. En traitant ces états inutilisés comme s'ils ne faisaient rien, le système garantit que la logique complexe du qudit fonctionne parfaitement dans le langage familier des qubits.

Les chercheurs n'ont pas seulement construit une idée théorique ; ils ont créé une bibliothèque entièrement fonctionnelle qui s'intègre de manière transparente aux flux de travail quantiques existants. Le logiciel définit un nouvel ensemble d'instructions, ou portes, qui agissent sur ces systèmes multi-niveaux. Ces portes peuvent effectuer des opérations standards comme le décalage de niveaux ou la création de superpositions complexes. Une caractéristique clé de la conception est sa capacité à gérer simultanément différentes tailles de qudits. Un seul circuit peut contenir un système à trois niveaux, un système à quatre niveaux et un système à cinq niveaux travaillant tous ensemble, et le logiciel gère la traduction pour chacun d'eux automatiquement. Lorsque le nombre de niveaux est une puissance de deux, le logiciel décompose l'opération en une série d'étapes de qubits standards que les outils d'optimisation de l'ordinateur peuvent facilement traiter. Lorsque le nombre de niveaux n'est pas une puissance de deux, le logiciel utilise une méthode différente, calculant la transformation mathématique exacte pour l'ensemble du groupe de qubits à la fois. Cela garantit que le système fonctionne correctement pour n'importe quelle taille jusqu'à seize niveaux, une plage qui couvre de nombreux systèmes physiques pratiques.

Pour prouver l'exactitude de leur méthode de traduction, l'équipe a mené des tests numériques rigoureux. Ils ont comparé le comportement de leur logiciel par rapport aux règles mathématiques connues pour un type spécifique d'opération quantique appelée transformée de Fourier quantique, qui est un outil fondamental dans de nombreux algorithmes quantiques. Ils ont testé cela à travers diverses tailles, y compris celles qui ne s'insèrent pas proprement dans les puissances de deux. Les résultats ont montré que la sortie du logiciel correspondait aux attentes théoriques avec une précision extrême, avec des erreurs si petites qu'elles étaient effectivement nulles pour des fins pratiques. Ils ont également vérifié que chaque type de porte qu'ils ont créé se comportait exactement comme prévu, qu'il s'agisse d'un simple décalage ou d'une interaction multi-niveaux complexe. Les tests ont confirmé que le logiciel gère correctement les parties « inutilisées » de l'espace encodé, garantissant qu'aucune information ne s'échappe ou ne soit corrompue pendant le processus.

L'article souligne également les compromis impliqués dans cette approche. Parce que le système utilise des qubits supplémentaires pour représenter chaque qudit, il nécessite plus de mémoire pour la simulation sur un ordinateur classique qu'un simulateur de qudit natif. Pour des systèmes très larges, cette surcharge devient significative, et les chercheurs notent que pour la simulation pure de systèmes massifs, un outil construit spécifiquement pour les qudits pourrait être plus efficace. Cependant, pour la grande majorité des cas d'utilisation actuels, le bénéfice d'utiliser l'écosystème Qiskit standard l'emporte sur ce coût. Le logiciel permet aux chercheurs de dessiner des circuits au niveau du qudit, de les faire passer par des processus d'optimisation standards et de les préparer pour l'exécution sur du matériel quantique réel ou des simulateurs sans avoir besoin de réécrire manuellement le code.

L'un des aspects les plus pratiques de ce travail est la gestion de la mesure. Lorsqu'un système quantique est mesuré, le résultat doit être reconverti des données brutes de qubits vers les niveaux significatifs du qudit. Le nouveau logiciel inclut une couche de décodage qui effectue cette traduction automatiquement. Il inclut également un mécanisme de sécurité pour détecter la « fuite » (leakage), une situation où un système se retrouve accidentellement dans un état qui ne correspond à aucun niveau valide du qudit. Dans une simulation parfaite, cela ne devrait jamais arriver, mais sur du matériel réel ou lorsqu'on introduit du bruit, cela peut se produire. Le logiciel permet aux utilisateurs de choisir comment gérer ces événements rares, soit en rejetant les mauvaises données, soit en les signalant pour examen, rendant les résultats de la simulation plus fiables et transparents.

Ce travail représente une étape importante pour rendre l'informatique quantique plus accessible et polyvalente. En transformant les concepts théoriques des systèmes quantiques multi-niveaux en une couche logicielle transposable et fonctionnelle, il lève un obstacle majeur pour les chercheurs. L'équipe a vérifié que leur méthode fonctionne pour des dimensions allant de deux jusqu'à seize, couvrant les besoins de nombreuses plateformes expérimentales actuelles. Bien que le logiciel repose sur une stratégie d'encodage spécifique qui ajoute une couche de complexité, l'auteur démontre que cette complexité est gérée efficacement, permettant aux outils standards de la communauté quantique d'opérer sur une gamme beaucoup plus large de systèmes physiques. Le résultat est un outil qui permet aux scientifiques d'explorer tout le potentiel de la mécanique quantique multi-niveaux sans avoir à abandonner l'infrastructure puissante et établie sur laquelle ils comptent déjà.

Résumé Technique : QISKIT-QUDITS

Énoncé du Problème

Les qudits, qui généralisent les qubits à deux niveaux aux qudits à $d$ niveaux de calcul, offrent des avantages potentiels dans le traitement de l'information quantique. Cependant, le framework Qiskit manque de support natif pour les qudits, ce qui constitue un écart de l'extension tierce connue. Les chaînes d'outils existantes sont confrontées à une dichotomie : les simulateurs de qudits natifs (ex. : Cirq, QuTiP) ne peuvent pas exploiter la pile de compilation étendue, les passes de transpileur et les backends matériels de Qiskit, tandis que les outils ciblant les qubits manquent souvent d'une surface de programmation pour les opérations au niveau des qudits. Les utilisateurs doivent actuellement implémenter manuellement les schémas d'encodage, les matrices de portes, les décompositions et le décodage classique, un processus sujet aux erreurs concernant le traitement de l'espace d'état non physique et le décodage de la mesure.

Méthodologie

Le document présente qiskit-qudits, une extension de Qiskit qui simule des qudits à $d$ niveaux en encodant chacun d'eux en $m = \lceil \log_2 d \rceil$ qubits. Le logiciel est construit sur les classes publiques Instruction, Gate et ControlledGate de Qiskit, assurant ainsi la compatibilité avec l'écosystème existant.

Décisions de Conception Fondamentales

Encodage et Remplissage d'Identité (Identity Padding) :
Comme $d$ n'est pas nécessairement une puissance de deux, l'espace de Hilbert encodé ( $D = 2^N$ ) est généralement plus large que le sous-espace logique ( $L = \prod d_i$ ). La bibliothèque adopte la convention de remplissage d'identité : chaque porte agit comme l'identité sur le complément orthogonal du sous-espace logique. Formellement, pour une porte logique $G$ , l'unitaire encodée est $U = E G E^\dagger + (I - \Pi)$ , où $E$ est l'isométrie d'encodage et $\Pi$ est le projecteur sur le sous-espace logique. Ce choix empêche la fuite (leakage) lors de la simulation sans bruit et préserve les propriétés de permutation.
Hiérarchie des Portes et Décomposition :
La bibliothèque définit trois classes de base : QuditGate (qudit simple), QuditControlledGate (multi-contrôle) et QuditMultiGate.
- Voie Puissance de Deux : Lorsque toutes les dimensions des opérandes sont des puissances de deux, les portes se décomposent en cascades fixes de portes de qubits Qiskit standards (ex. : mcx, cp, synth_qft_full). Cela permet au transpileur de Qiskit d'optimiser le circuit.
- Repli sur Synthèse Dense (Dense-Synthesis Fallback) : Pour les dimensions non puissances de deux, les portes basculent sur une synthèse unitaire dense exacte. La bibliothèque génère la matrice unitaire complète $D \times D$ (via le protocole __array__ de NumPy) et l'intègre comme une instruction unitaire opaque unique, que la passe UnitarySynthesis décompose ensuite.
Modèle de Circuit et Mesure :
Le QuditQuantumCircuit maintient deux vues synchronisées : un journal d'instructions au niveau qudit (.data) et un circuit encodé au niveau qubit (.circuit). La mesure est gérée via QuditMeasure, qui mappe les qubits d'encodage vers des bits classiques (ClByte). Une couche de décodage dédiée traduit les comptes du backend vers les niveaux qudit, détectant explicitement les événements de "fuite" (résultats $\ge d$ ) qui peuvent être gérés via des politiques configurables (keep, drop, raise).

Contributions Clés

Implémentation d'une Couche Qudit Transpilable : Ce travail convertit les définitions théoriques de portes issues de travaux antérieurs [1] en une extension Qiskit fonctionnelle et transpilable. Elle supporte les dimensions $2 \le d \le 16$ et les registres à dimensions mixtes.
Simulation Exacte pour les Dimensions Non-Puissances de Deux : Contrairement aux approches de simulation par qubits précédentes limitées à $d=2^m$ , cette bibliothèque supporte la synthèse exacte pour des $d$ arbitraires dans la plage supportée en utilisant des replis unitaires denses.
Vérification de la Théorie et de l'Implémentation :
- Vérification A : Vérifié que les matrices unitaires denses reproduisent la transformée de Fourier discrète (DFT) sur le sous-espace logique pour les dimensions tant puissances de deux que non-puissances de deux, étendant ainsi la vérification de [1].
- Vérification B : Validé que les circuits décomposés (voies cascade et dense) correspondent aux définitions unitaires denses à travers 138 instances de portes, avec des écarts de l'ordre de la précision machine ( $\approx 10^{-15}$ ).
- Invariante de Remplissage : Confirmé que l'invariante de remplissage d'identité est respectée exactement, garantissant l'absence de mélange entre les sous-espaces valides et invalides.
Détection de Fuite (Leakage) : Fournit un mécanisme explicite pour détecter et gérer les événements de fuite dans les simulations, distinguant le comportement idéal des erreurs induites par le bruit.

Résultats

La vérification numérique confirme la correction de l'implémentation :

QFT vs. DFT : L'implémentation QFT correspond à la matrice DFT théorique à $10^{-15}$ près pour les dimensions $d \in \{2, 3, 4, 5, 8\}$ .
Décomposition de Porte : Toutes les portes, incluant SUMX, SUMP, SWAP et QFT, produisent des décompositions cohérentes avec leurs définitions unitaires denses.
Analyse des Ressources :
- Puissance de Deux : Les décompositions sont efficaces et favorables au transpileur. Par exemple, une QFT de 3 qutrits ( $d=3$ ) nécessite 1 783 portes CX, tandis qu'une QFT de 3 ququarts ( $d=4$ ) n'en nécessite que 56.
- Non-Puissance de Deux : Le repli dense induit un surcoût significatif. Une QFT de 3 qutrits ( $d=3$ ) se synthétise en 1 783 portes CX, tandis qu'une QFT de 3 ququint ( $d=5$ ) en nécessite 119 383. Cela souligne le coût de la synthèse dense pour les dimensions non-puissances de deux plus importantes.
- Contraintes de Mémoire : La voie dense est limitée par la mémoire ; une porte sur huit qudits à 16 niveaux ( $N=32$ ) nécessiterait $\approx 256$ EiB de mémoire, ce qui rend la limite pratique $N \approx 13\text{--}14$ sur du matériel courant.

Signification et Limites

Le document positionne qiskit-qudits comme un pont permettant aux utilisateurs d'écrire des algorithmes qudit au niveau logique tout en conservant un accès complet à la pile de compilation, aux backends et aux primitives de Qiskit. Il traite explicitement la "partie non physique" de l'espace encodé et le décodage de la mesure, qui sont souvent sources d'erreurs dans les implémentations manuelles.

Limites et Compromis :

Plafond de Dimension : La bibliothèque impose $2 \le d \le 16$ et un maximum de 8 qudits par porte. Ce sont des plafonds nominaux ; la limite pratique est déterminée par la mémoire requise pour la synthèse unitaire dense ( $N \lesssim 13\text{--}14$ ).
Sémantique de la Fuite : La convention de remplissage d'identité fige l'amplitude fuitée plutôt que de modéliser la fuite physique. Elle ne commute pas avec le produit tensoriel, ce qui signifie que l'équivalence des circuits doit être vérifiée sur le sous-espace projeté ( $\Pi U \Pi$ ).
Optimisation : L'implémentation actuelle n'exploite pas les lignes/colonnes "don't-care" du transpileur, et le repli dense pour les dimensions non-puissances de deux est coûteux en calcul par rapport aux simulateurs de qudits natifs ou à la synthèse de logique réversible optimisée.

Le travail conclut que bien que le repli dense permette une simulation exacte pour les dimensions non-puissances de deux, les améliorations futures devraient se concentrer sur la synthèse de logique réversible agnostique à la dimension et sur l'exposition du sous-espace non physique au transpileur afin de réduire la surcharge de ressources.