Tomographie d'ensembles de spins : comment transformer un spectre en état quantique
Contexte
Ce travail fait un peu bande à part : dans les archives, ses plus proches voisins traitent d'algorithmes quantiques, de correction d'erreurs et d'information quantique, mais jamais de reconstruction pratique d'états à partir de spectres. Il existe déjà des outils pour modéliser la dynamique des spins (PULSEE) et des bibliothèques générales (QuTiP, nmrglue), mais il manquait le chaînon entre les données brutes du spectromètre et une matrice densité validée. Le paquet Tomography-NMR comble ce vide, avec pour objectif un traitement transparent et reproductible des données RMN pour la tomographie quantique.
Méthodes
Le paquet implémente une chaîne d'analyse en Python 3.11+ : les signaux bruts dans le domaine temporel subissent une transformée de Fourier avec correction de phase automatique ou manuelle, puis les intensités des pics des doublets sont extraites de trois façons. La première est la mesure directe de la hauteur du pic en un point ; la deuxième, une intégration numérique par la méthode de Simpson dans une fenêtre fixe ; la troisième, un balayage systématique des paramètres (décalages de fréquence, largeur de fenêtre, correction de la séparation) afin de maximiser la précision par rapport à un état théorique connu. La somme et la différence des pics gauche et droit du doublet sont associées aux coefficients des opérateurs produits via une table de correspondance, puis la matrice densité est reconstruite en imposant des contraintes physiques (hermiticité, trace normalisée, semi-définie positive).
Résultats
Sur 20 états expérimentaux préparés sur deux noyaux de phosphore 31 couplés dans l'ADP, la précision moyenne de reconstruction selon la métrique de Fortunato s'échelonne de 0,975 à 0,995 selon la catégorie d'état. Les états de base atteignent 0,995, les sorties de la porte de Hadamard 0,992, celles de CNOT 0,983, les états de Bell 0,980, et la combinaison H+T 0,975. L'intégration optimisée donne une précision supérieure à 99 % pour certains états (jusqu'à 0,998 pour |00⟩), tandis que les méthodes sans connaissance de l'état cible restent autour de 98 %. Les erreurs sur les éléments de la matrice densité sont estimées à 10⁻³–10⁻⁴.
Implications
Contrairement aux travaux voisins, centrés sur les algorithmes ou la modélisation, cette étude résout le problème inverse : reconstruire un état quantique à partir de spectres expérimentaux réels. Elle fournit la première chaîne ouverte et entièrement documentée pour la tomographie RMN, ce qui affranchit les chercheurs des logiciels commerciaux fermés et permet une vérification indépendante des résultats publiés. Cela change la pratique de l'informatique quantique expérimentale, en la rendant plus transparente et plus reproductible.
Développements futurs
L'extension à des systèmes de trois qubits ou plus exigera de résoudre le problème du recouvrement des multiplets et de la passage à l'échelle de l'optimisation. Un prototype de module pour trois spins est déjà inclus, mais la prise en charge complète reste un problème ouvert. Sont également prévues une optimisation bayésienne en remplacement de la recherche sur grille, une analyse en temps réel pour les protocoles adaptatifs et de nouvelles métriques de validation, comme les témoins d'intrication.
Impact
Ce travail aura un impact sur l'informatique quantique expérimentale, en particulier sur le calcul quantique par RMN, en fournissant un outil standardisé pour la caractérisation des états. Il accélérera aussi la formation des nouveaux spécialistes et facilitera la comparaison des résultats entre groupes.
Prochaines étapes
Les prochaines étapes comprennent la mise en œuvre d'une optimisation bayésienne pour accélérer la recherche des paramètres et le développement de méthodes de déconvolution spectrale pour les systèmes à trois qubits. Il importe aussi d'élargir la prise en charge des formats de données d'autres spectromètres.
Problèmes ouverts clés
En informatique quantique, la tomographie d'états à grande échelle reste un problème non résolu : le nombre de mesures croît exponentiellement avec le nombre de qubits. Ce travail ne le résout pas, mais il établit une référence de précision et de transparence pour le cas à deux qubits, ce qui est indispensable pour tester les futures méthodes.
Le paquet s'appelle Tomography-NMR et est distribué sous licence MIT, ce qui permet de l'utiliser librement dans des projets académiques et commerciaux. La validation a été menée sur la molécule d'ADP, où deux noyaux de phosphore 31 sont couplés par une interaction J d'environ 22 Hz.
Chiffres clés
- précision de reconstruction (optimisée): jusqu'à 0,998
- précision de reconstruction (sans connaissance de la cible): environ 0,980
- nombre d'états validés: 20
- nombre de coefficients non identité pour deux qubits: 15
- temps d'optimisation pour un état: 2–5 minutes