Designing quantum technologies with a quantum computer

Résumé technique : Concevoir des technologies quantiques avec un ordinateur quantique

Énoncé du problème
Les matériaux à l'état solide contenant des défauts actifs en spin, tels que les centres azote-lacune (NV) dans le diamant, sont essentiels pour les capteurs quantiques, les sources de photons uniques et les processeurs. Cependant, la modélisation de ces systèmes est complexe en raison des interactions complexes entre le spin et l'environnement, notamment la levée de dégénérescence en champ nul, les effets Zeeman, les couplages hyperfins, les interactions dipôle-dipôle et la décohérence électron-phonon. Les simulations classiques peinent face à la nature de corps multiples de ces systèmes, particulièrement lors du calcul de la dynamique à long terme (nanosecondes à microsecondes) nécessaire pour déterminer les temps de cohérence et les spectres d'absorption. Bien que les ordinateurs quantiques numériques offrent une alternative prometteuse, le matériel actuel de type NISQ (Noisy Intermediate-Scale Quantum) et les premiers matériels tolérants aux fautes sont contraints par un nombre limité de qubits, un bruit intrinsèque et une profondeur de circuit limitée. Le défi central consiste à développer des algorithmes efficaces capables d'approcher l'évolution temporelle quantique à long terme tout en minimisant les ressources de calcul (nombre de portes et profondeur de circuit) afin de rester compatibles avec les limitations matérielles actuelles.

Méthodologie
Les auteurs proposent un cadre de conception assisté par ordinateur quantique qui simule des systèmes de spins interagissants à l'aide d'un Hamiltonien de résonance paramagnétique électronique (ESR) général. La méthodologie intègre plusieurs techniques d'optimisation :

1. Modélisation de l'Hamiltonien : Le cadre construit un Hamiltonien intégrant la levée de dégénérescence en champ nul, les termes Zeeman, les interactions hyperfines, les termes de spin-spin dipôle-dipôle et un modèle spin-boson pour la décohérence électron-phonon. Le système est transposé en une représentation de spin-1/2 via un codage de Gray pour minimiser les besoins en qubits ($\lceil \log_2(d) \rceil$).
2. Optimisation du circuit : L'Hamiltonien résultant, exprimé comme une somme de chaînes de Pauli, est partitionné en groupes commutant par qubit (QWC). Cela permet une diagonalisation simultanée via des rotations de Clifford à un qubit, permettant l'utilisation de circuits canoniques pour l'exponentiation de chaînes de Pauli diagonales. Cette étape facilite le réordonnancement optimal des exponentielles pour annuler les opérations adjointes adjacentes, réduisant ainsi le nombre de portes et la profondeur du circuit.
3. Fast-Forwarding de Krylov rapide à sélection multi-référence (sQKFF) : Pour accéder aux dynamiques à long terme sans une profondeur de circuit prohibitive, les auteurs emploient l'algorithme sQKFF. Cette approche hybride approxime l'état évolué temporellement $|\psi(t)\rangle$ au sein d'un sous-espace de Krylov construit à partir d'un ensemble d'états de référence.
   • Composante Quantique : Calcule les éléments de matrice du recouvrement ($S$) et des matrices de l'Hamiltonien ($H$) en utilisant des tests de Hadamard avec une Trotterisation de second ordre.
   • Composante Classique : Résout le problème de la valeur propre généralisée (équation de Schrödinger dans le sous-espace de Krylov) pour déterminer les coefficients dépendants du temps.
   • Sélection de Référence : Les états de référence sont sélectionnés sur la base des probabilités de mesure d'un état évolué à court terme, visant à capturer les caractéristiques pertinentes de la dynamique à long terme.
4. Calcul de Propriétés : Le cadre calcule les fonctions d'autocorrélation, les spectres d'absorption micro-ondes (via la théorie de la réponse linéaire et la partie imaginaire de la susceptibilité) et la norme $\ell_1$ de la cohérence.

Contributions Clés

• Cadre Intégré : Une approche unifiée combinant le codage de Gray, l'agrégation QWC et l'algorithme sQKFF spécifiquement adapté aux défauts de spin à l'état solide sous des conditions expérimentales réalistes (incluant les bains de phonons et les impuretés).
• Réduction des Ressources : L'application de l'agrégation QWC et de l'optimisation de l'évolution temporelle par Trotterisation a entraîné des réductions de 18 à 30 % du nombre de portes et de la profondeur de circuit par rapport aux implémentations non optimisées pour toutes les configurations testées.
• Dynamique à Long Terme : Le cadre a réussi à calculer des fonctions d'autocorrélation jusqu'à environ 100 ns, une échelle de temps souvent difficilement accessible pour une Trotterisation standard sur le matériel NISQ.
• Benchmarking : La méthode a été validée par rapport à des simulations classiques utilisant QuTiP sur trois configurations distinctes de centres NV (un seul NV, trois NV, et un NV avec des impuretés $^{14}$N).

Résultats

• Précision : Le cadre a reproduit les résultats de QuTiP avec une haute fidélité pour les spectres d'absorption micro-ondes et les fonctions d'autocorrélation. Pour les configurations 1 et 3, des résultats précis ont été maintenus jusqu'à des dizaines de nanosecondes, même avec des erreurs de Trotterisation modérées.
• Sensibilité des États de Référence : L'étude a identifié que la sélection des états de référence est le facteur principal régissant la précision pour un coût matériel donné. L'augmentation du nombre d'états de référence ($R$) a considérablement amélioré tant la précision de la simulation que l'étendue temporelle atteignable, souvent plus que la réduction de l'erreur de Trotterisation.
• Complexité des Configurations : Bien que la configuration 2 (taille de système plus grande, norme 1 de l'Hamiltonien plus élevée) ait présenté des défis pour identifier des états de référence efficaces en raison d'une distribution d'information plus éparse, le cadre est resté applicable. La configuration 3, malgré les impuretés additionnelles et les interactions dipôle-dipôle, a atteint une précision comparable à la configuration 1 plus simple lorsque $R=15$ et que l'erreur de Trotter est minimisée.
• Métriques de Cohérence : La norme $\ell_1$ de la cohérence a été calculée efficacement au sein du sous-espace de Krylov, permettant l'estimation des temps de cohérence effectifs $T_1$ et $T_2$ sans surcoût classique exponentiel.

Signification et Revendications
L'article affirme que cette méthodologie fournit un schéma flexible pour utiliser les ordinateurs quantiques afin de concevoir, comparer et optimiser les technologies de qubits à spin à l'état solide. Les auteurs soutiennent que leur approche hybride démontre que les ordinateurs quantiques, lorsqu'ils sont combinés à des techniques de simulation optimisées, peuvent servir de plateformes efficaces pour l'étude de la dynamique des spins à l'état solide, même sous les contraintes de ressources de l'ère NISQ.

Le travail souligne que si les limitations matérielles (bruit, nombre de qubits) sont significatives, les optimisations algorithmiques telles que le sQKFF et l'agrégation QWC peuvent étendre les échelles de temps opérationnelles des simulations. Les auteurs concluent modestement que, bien que la simulation de dynamiques complexes avec des interactions non triviales reste un objectif à long terme, leurs résultats prouvent la faisabilité de telles démarches. Ils soulignent que la capacité de simuler des configurations classiquement insolubles est cruciale pour l'avancement de l'ingénierie des dispositifs quantiques. L'article ne revendique pas un déploiement immédiat sur un matériel réel, mais établit un cadre numérique validé qui traite les surcoûts de mesures (shot overheads) et la sélection des états de référence comme des domaines critiques pour une mise en œuvre pratique future.