Dynamical structure factor with a pumping approach on a trapped-ion quantum computer

Résumé Technique : Facteur de Structure Dynamique avec une Approche par Pompage sur un Ordinateur Quantique à Ions Piégés

Énoncé du Problème
Le facteur de structure dynamique (FSD), noté $S(q, \omega)$, est une quantité critique en physique de la matière condensée qui décrit la réponse d'un matériau aux neutrons. Il est défini comme la transformée de Fourier espace-temps des corrélations dynamiques $\langle O(x, t)O(0, 0) \rangle$. Le calcul du FSD est difficile pour les ordinateurs classiques car il nécessite à la fois de trouver une représentation de l'état fondamental ou thermique et de simuler la dynamique hamiltonienne au cours du temps, une tâche généralement considérée comme BQP-difficile. Bien que les ordinateurs quantiques offrent une voie pour préparer des états à température finie et effectuer des dynamiques hamiltoniennes efficaces, les algorithmes quantiques existants pour calculer $S(q, \omega)$ reposent généralement sur le calcul de corrélations temporelles à de nombreux points de temps discrets, suivis d'une transformée de Fourier. Cette approche induit un surcoût important en termes de mesures (shots), particulièrement lorsque seuls quelques points de fréquence spécifiques sont d'intérêt.

Méthodologie : L'Approche par Pompage
Les auteurs introduisent une « approche par pompage » qui cible directement des valeurs de fréquence $\omega$ spécifiques sans nécessiter le calcul des corrélations dans le domaine temporel et la réalisation d'une transformée de Fourier. Le cœur de la méthode consiste à conduire le système hors de l'équilibre en utilisant un Hamiltonien dépendant du temps perturbé par un terme de source oscillant à la fréquence cible.

1. Configuration de l'Hamiltonien : Étant donné un Hamiltonien de système $H$ et un observable $O$, les auteurs définissent un Hamiltonien dépendant du temps :
   $$H(t) = H - \epsilon O(0) \sin(\omega(T - t))$$
   où $\epsilon$ est une constante de couplage et $T$ est le temps d'évolution final.
2. Protocole de Mesure : Le système est évolué à partir d'un état d'équilibre initial (par exemple, l'état fondamental) en utilisant cet Hamiltonien dépendant du temps. La valeur attendue d'un opérateur d'impulsion $\hat{O}_q$ est mesurée au temps $T$.
3. Extraction du FSD : En mesurant le système avec un couplage $+\epsilon$ et $-\epsilon$, le FDT est extrait via la relation :
   $$S(q, \omega) \approx \frac{\langle \hat{O}_q(T) \rangle_\epsilon - \langle \hat{O}_q(T) \rangle_{-\epsilon}}{\epsilon(1 - e^{-\beta\omega})}$$
   Cette formule est dérivée d'un développement au premier ordre en $\epsilon$, en exploitant la relation de l'équilibre détaillé pour les états thermiques.
4. Détails de Mise en Œuvre :
   • Préparation de l'État : Pour le modèle de Heisenberg 1D, l'état fondamental est préparé à l'aide d'un ansatz inspiré de l'évolution adiabatique. Le système est initialisé dans un état de Néel et évolue via une interpolation linéaire entre un Hamiltonien à champ de structure et le Hamiltonien de Heisenberg. Les paramètres sont optimisés classiquement à l'aide de simulations de l'état (statevector) sur une chaîne de 20 sites, puis appliqués au matériel.
   • Évolution Temporelle : L'opérateur d'évolution temporelle est approximé par une décomposition de Trotter. Le circuit consiste en une séquence de portes à un qubit et à deux qubits (spécifiquement des portes TK2 sur le matériel Quantinuum).
   • Régularisation : Puisque $T \to \infty$ n'est pas réalisable, un $T$ fini est utilisé, ce qui revient à convoluer le FDT réel avec une fonction sinc, lissant le résultat mais introduisant potentiellement des caractéristiques parasites.

Contributions Clés

• Ciblage Direct de la Fréquence : La contribution principale est une méthode permettant aux ordinateurs quantiques de cibler directement des fréquences $\omega$ arbitraires, évitant ainsi le coût computationnel de l'échantillonnage de l'ensemble du domaine temporel.
• Efficacité des Mesures (Shot Efficiency) : Les auteurs démontrent que pour les scénarios nécessitant seulement quelques valeurs de fréquence, cette approche offre une réduction significative du surcoût de mesures par rapport aux méthodes nécessitant des mesures de corrélations dans le domaine temporel.
• Implémentation Hybride : Ce travail présente une implémentation complète sur l'ordinateur quantique à ions piégés Quantinuum Reimei, intégrée dans un environnement hybride quantique-HPC Fugaku–Reimei. L'optimisation classique des paramètres de préparation d'état est effectuée sur le supercalculateur Fugaku, tandis que les circuits quantiques sont exécutés sur Reimei.

Résultats
Les auteurs ont validé la méthode à l'aide de deux cas :

1. Modèle de Heisenberg 1D (20 sites) : Le FDT a été calculé pour une chaîne de 20 sites et comparé à des simulations sans bruit. Les résultats matériels ont montré un excellent accord avec les simulations sans bruit, obtenu sans aucune technique d'atténuation du bruit.
2. Cristal de Sulfate de Cuivre ($CuSO_4$) : Le modèle de Heisenberg 1D est connu pour décrire les excitations magnétiques des cristaux de sulfate de cuivre. Le FDT calculé a été comparé aux données expérimentales de diffusion neutronique. Les résultats du matériel quantique concordent très bien avec les données expérimentales.

Les paramètres spécifiques utilisés comprenaient une taille de système $L=20$, un pas de Trotter $dt=0.24$, un temps d'évolution total $T=2.4$, et un couplage $\epsilon=2$. La préparation de l'état fondamental utilisait un seul tour ($R=1$) de l'ansatz optimisé. L'étude a noté que l'utilisation d'une seule couche de l'ansatz donnait en fait un écart plus faible par rapport aux résultats exacts de l'état (statevector) qu'en utilisant deux couches, probablement en raison de la réduction du nombre de portes.

Signification et Revendications
L'article affirme que cette approche par pompage offre une voie pratique et efficace pour calculer les facteurs de structure dynamique sur le matériel quantique de l'ère NISQ (Near-term). En simulant directement une « version jouet d'une expérience de diffusion neutronique », la méthode évite de recourir à des calculs complexes de corrélations en temps réel.

Les auteurs soulignent la pertinence potentielle de ce protocole au-delà de la diffusion neutronique, spécifiquement pour les simulations de Résonance Magnétique Nucléaire (RMN). En RMN, les spectres sont souvent localisés de manière très nette à des fréquences spécifiques, et cette méthode pourrait offrir des économies de temps d'exécution significatives en échantillonnant des fréquences fixes plutôt que de longs déclinements dans le domaine temporel.

Concernant la scalabilité, les auteurs se montrent modestes quant à l'avantage quantique actuel. Ils reconnaissent que la préparation d'état dans ce travail reposait sur des ansatz optimisés classiquement à partir de simulations exactes de l'état (statevector), un régime où aucun avantage quantique n'existe. Cependant, ils proposent que la combinaison de leur algorithme de FDT avec des méthodes de préparation d'état classiques scalables (telles que les Sparse Pauli Dynamics ou les Matrix Product States via DMRG) au sein d'un cadre hybride puisse éventuellement mener à un régime où la tâche conjointe de préparation d'état et de calcul du facteur de structure devienne difficile pour les ordinateurs classiques seuls.