Spectral functions on a quantum computer through system-environment interaction

Cet article présente un algorithme quantique efficace qui modélise les interactions système-environnement pour mesurer les fonctions spectrales avec une réduction de $O(N)$ de la surcharge d'échantillonnage par rapport aux techniques standard, démontrant son efficacité sur un système de 27 sites utilisant 54 qubits sur un ordinateur quantique à piège à ions Quantinuum.

L'essentiel

Imaginez que vous essayez de comprendre les « notes musicales » (niveaux d'énergie) qu'un matériau complexe peut jouer. Dans le monde réel, les scientifiques utilisent un appareil photo haute technologie appelé ARPES (Spectroscopie de Photoémission à Résolution Angulaire) pour prendre une photo de ces notes. Pour ce faire, ils projettent de la lumière sur le matériau, arrachant des électrons, puis mesurent la vitesse et la direction dans lesquelles ces électrons s'envolent.

Le problème est que simuler ce processus sur un ordinateur est incroyablement difficile. C'est comme essayer de prédire le son d'une symphonie en écoutant chaque instrument individuellement, dans un silence total, puis en essayant de deviner toute la chanson. Sur un ordinateur quantique, l'ancienne méthode consistait à demander à un musicien de jouer une note, de s'arrêter, de réinitialiser, de jouer la note suivante, de s'arrêter et de réinitialiser à nouveau. Si vous avez 1 000 instruments (ou « sites » dans le matériau), vous devez répéter ce processus 1 000 fois juste pour obtenir une image complète. Cela prend une éternité et gaspille une quantité massive de temps.

La Nouvelle Idée : Un Environnement « Faux »

Les auteurs de cet article ont trouvé un raccourci astucieux. Au lieu de demander à l'ordinateur de calculer les notes une par une, ils ont décidé de simuler l'expérience réelle directement sur l'ordinateur quantique.

Pensez-y ainsi :

• Le Système : C'est le matériau que vous voulez étudier (l'orchestre).
• L'Environnement : C'est l'« appareil photo » ou le « vide » qui capture les électrons (le public).

Dans leur nouvelle méthode, ils connectent l'« orchestre » à un « faux public » (un environnement) à l'intérieur de l'ordinateur. Ils laissent l'orchestre interagir avec ce public pendant un court moment. Ensuite, au lieu de mesurer l'orchestre directement, ils regardent simplement le public pour voir qui a attrapé une note.

Comme le public est connecté à tout l'orchestre en même temps, une seule mesure leur indique les « notes » de l'ensemble de l'orchestre simultanément.

Le Grand Gagnant : Vitesse et Efficacité

L'article affirme que c'est un changement de donne pour un type spécifique d'ordinateur quantique appelé ordinateur à piège à ions (qui utilise des atomes piégés comme qubits).

• L'Ancienne Méthode : Pour obtenir une image claire, vous pourriez avoir besoin de prendre 1 000 photos (mesures) car l'appareil photo est lent et flou.
• La Nouvelle Méthode : Vous n'avez besoin que d'une seule photo.

Les auteurs disent que cela économise une quantité massive de temps. Si l'ancienne méthode prenait 100 heures, cette nouvelle méthode pourrait ne prendre qu'une heure. Ils appellent cela une amélioration O(N), ce qui signifie que si vous doublez la taille du matériau que vous étudiez, l'ancienne méthode devient deux fois plus lente, mais cette nouvelle méthode reste tout aussi rapide.

Le Problème : Il Faut Plus de « Qubits »

Il y a un compromis. Pour réussir ce tour, vous devez doubler le nombre de « qubits » (les unités de base de l'ordinateur quantique) car vous devez simuler à la fois le matériau et le faux environnement. C'est comme avoir besoin d'une plus grande pièce pour loger à la fois le groupe et le public. Cependant, les auteurs soutiennent que pour ces ordinateurs spécifiques, gagner du temps sur les mesures est beaucoup plus important que d'avoir quelques qubits supplémentaires.

Le Tour de « Magie » : La Transformée de Fourier Fermionique

Pour faire fonctionner ce « faux public », l'ordinateur doit effectuer une danse mathématique complexe appelée Transformée de Fourier Fermionique (FFT). Imaginez mélanger un jeu de cartes de sorte que tous les cœurs soient ensemble, tous les piques ensemble, etc., mais en le faisant d'une manière qui respecte les règles étranges des particules quantiques (fermions).

Les auteurs n'ont pas simplement utilisé un mélange standard ; ils ont inventé une manière plus efficace de mélanger ces cartes quantiques spécifiques, en particulier pour une configuration où le nombre de cartes n'est pas une puissance de 2 (comme 27 cartes). Ils ont testé ce mélange sur une machine réelle (H2 de Quantinuum) et ont prouvé que cela fonctionne.

Le Test Réel

L'équipe n'a pas seulement écrit de la théorie ; ils ont mené l'expérience sur un véritable ordinateur quantique avec 54 qubits (27 pour le matériau, 27 pour l'environnement). Ils ont mesuré avec succès la « fonction spectrale » (les notes musicales) d'une chaîne de 27 sites de particules.

Même si l'ordinateur réel présente un certain « bruit » (comme des parasites à la radio), les résultats étaient assez clairs pour voir les caractéristiques principales du matériau. Le « bruit » a rendu le signal un peu plus faible, mais il n'a pas déformé la forme des notes, ce qui signifie que la physique qu'ils recherchaient est restée précise.

Résumé

En bref, cet article présente une nouvelle façon de simuler comment les matériaux interagissent avec la lumière. En simulant l'expérience entière (système + environnement) à la fois, plutôt que de calculer des parties séparément, ils peuvent obtenir la réponse N fois plus vite (où N est la taille du système). Cela rend beaucoup plus pratique l'étude de matériaux grands et complexes sur les ordinateurs quantiques d'aujourd'hui, spécifiquement ceux de type piège à ions.

Résumé technique

Énoncé du problème
Les fonctions spectrales, mesurées par spectroscopie de photoémission résolue en angle (ARPES), sont essentielles pour élucider la structure de bande des matériaux. Théoriquement, ces fonctions correspondent à des fonctions à un point dynamiques dans un état d'équilibre. Bien que leur calcul soit classiquement difficile, leur mesure sur des ordinateurs quantiques fait face à des obstacles majeurs. Les approches standard, qui consistent à calculer les corrélations dynamiques $\langle c^\dagger_j(t)c_0(0)\rangle$ et à effectuer des transformées de Fourier, souffrent d'une surcharge d'échantillonnage importante. Plus précisément, comme les opérateurs $c_j + c^\dagger_j$ pour différents sites $j$ ne commutent pas, une seule corrélation peut être mesurée par tir. Pour reconstruire la fonction spectrale pour tous les $N$ moments, cela nécessite une surcharge d'échantillonnage de $O(N)$. Ceci est particulièrement problématique pour les ordinateurs quantiques à pièges à ions, qui possèdent une haute fidélité de porte et une connectivité élevée, mais des temps de tir relativement lents, rendant le nombre de tirs $O(N)$ un goulot d'étranglement pour les applications à l'échelle de l'utilité.

Méthodologie
Les auteurs proposent une méthode pour mesurer directement les fonctions spectrales $A(k, \omega)$ en modélisant l'interaction physique entre le système et l'environnement telle qu'elle se produit dans les expériences ARPES. L'approche comprend :

1. Couplage Système-Environnement : Le système d'intérêt (Hamiltonien $H_{sys}$) est couplé à un environnement de $N$ sites (Hamiltonien $H_{env}$) via un terme d'interaction $H_{int}$. L'Hamiltonien total est $H = H_{sys} + \epsilon H_{int} + \omega H_{env}$.
2. Préparation et Évolution de l'État : Le système est initialisé dans un état propre $|E\rangle$ de $H_{sys}$, tandis que l'environnement est initialisé dans un état vide $|0\rangle$ (pour $A_+$) ou un état entièrement rempli $|all\rangle$ (pour $A_-$). Le système total évolue sous $H$ pendant un temps fixe $t$.
3. Mesure : Au lieu de mesurer des corrélations dynamiques, le protocole mesure le nombre d'occupation de moment $n(k)$ sur l'environnement. Cela nécessite d'effectuer une transformée de Fourier fermionique (FFT) sur les qubits de l'environnement à la fin de l'évolution, suivie de la mesure de matrices de Pauli locales $Z$.
4. Résultat Théorique : Dans la limite d'un couplage faible ($\epsilon \to 0$) et d'un temps long ($t \to \infty$), la valeur moyenne $\langle n(k) \rangle$ est proportionnelle à la fonction spectrale $A(k, \omega)$ (convoluée avec un noyau qui tend vers une fonction delta lorsque $t \to \infty$).
5. Implémentation de la FFT Fermionique : Pour gérer des tailles de systèmes arbitraires (spécifiquement $N=27$), les auteurs développent une décomposition de porte efficace pour une FFT fermionique générique de radix-$n$. Ils utilisent une structure récursive impliquant des opérations d'« entrelacement » pour rendre les modes fermioniques non adjacents adjacents pour la transformée, optimisée par des techniques de décimation de graphe pour les architectures à couplage tout-à-tout.

Contributions Clés

• Échantillonnage Efficace : La contribution principale est une réduction de la surcharge d'échantillonnage. En mesurant l'occupation de moment de l'environnement, la méthode obtient des données spectrales pour tous les $N$ moments $k$ en un seul tir, réduisant le nombre de tirs requis d'un facteur $O(N)$ par rapport aux approches standard de corrélation dynamique.
• Amélioration du Temps d'Exécution : Cette réduction du nombre de tirs se traduit par une amélioration de $O(N)$ du temps d'exécution global, ce qui est crucial pour des plateformes comme les ordinateurs quantiques à pièges à ions où le temps de tir est le facteur limitant.
• Robustesse à l'Erreur de Trotter : Les auteurs soutiennent que leur approche est plus robuste aux erreurs d'évolution temporelle (Trotterisation) que les méthodes standard. Alors que les méthodes standard peuvent produire des valeurs négatives non physiques pour les fonctions spectrales en raison d'erreurs, la méthode proposée mesure un nombre d'occupation, qui reste une quantité positive même avec une évolution imparfaite.
• Étalonnage Matériel : Les auteurs ont implémenté une FFT fermionique de radix-3 sur 27 qubits et ont évalué diverses stratégies de compilation (CZ locaux, échelles CX, décimation de graphe) sur les systèmes à pièges à ions H2-1 et H2-2 de Quantinuum. Ils ont identifié que la compilation avec une échelle $\log_2(N)$ offrait la plus faible énergie et la plus faible infidélité.
• Démonstration de l'Algorithme Complet : L'algorithme complet a été démontré sur un système de 54 qubits (27 sites système + 27 sites environnement) sur le matériel Quantinuum H2-2, calculant avec succès la fonction spectrale d'un Hamiltonien de fermions libres 1D.

Résultats

• Performance Matérielle : Sur le matériel H2-2, la méthode a reproduit avec succès la forme de la bande spectrale pour un système de 27 sites. Bien que le bruit matériel ait causé une atténuation générale du signal (un « plancher de bruit » proche de 0,5), les caractéristiques physiques, telles que les positions et largeurs de bande, sont restées non biaisées.
• Analyse des Erreurs : Les comparaisons avec des simulations exactes ont montré que la méthode proposée atteint la même précision avec significativement moins d'étapes de Trotter que l'approche standard de corrélation dynamique, en particulier pour un petit couplage $\epsilon$.
• Analyse des Coûts : Pour l'expérience spécifique de 27 sites, la méthode a entraîné une surcharge de portes à deux qubits (due à l'environnement et à la FFT) mais a économisé un facteur d'environ 27 dans le nombre de tirs. Les auteurs estiment que sans cette méthode, le coût total du temps d'exécution aurait été environ 11 fois plus élevé pour la même précision.
• Effets de Couplage Fini : L'analyse théorique et les simulations de fermions libres indiquent qu'un couplage fini $\epsilon$ introduit un élargissement de bande et des « bandes fantômes ». Cependant, pour $\epsilon t \lesssim \pi$, ces effets sont négligeables et la méthode approxime fidèlement la fonction spectrale.

Signification
L'article affirme que cette approche améliore fondamentalement l'évolutivité du calcul des fonctions spectrales sur les ordinateurs quantiques. En échangeant une surcharge de qubits et de portes contre une réduction massive des exigences d'échantillonnage, la méthode s'attaque à un goulot d'étranglement critique pour les ordinateurs quantiques à pièges à ions. Les auteurs postulent que, à mesure que les tailles de systèmes augmentent vers l'échelle de l'utilité (centaines ou milliers de sites), la réduction $O(N)$ du nombre de tirs deviendra un facteur décisif dans la faisabilité pratique de la simulation des propriétés des matériaux sur du matériel quantique. Ce travail fournit également une implémentation concrète et efficace des transformées de Fourier fermioniques pour des tailles de systèmes non puissances de deux, un composant nécessaire pour simuler des modèles de matériaux réalistes.