⚛️ quantum physics

Preparation of the single-spinon wave function on a quantum computer

Cet article présente des algorithmes quantiques basés sur une combinaison linéaire d'unitaires pour préparer et analyser l'énergie des fonctions d'onde de spinons uniques pour les modèles de Heisenberg et de Haldane-Shastry, en comparant trois stratégies distinctes en termes de ressources computationnelles et de performance.

Auteurs originaux : D. Faílde, A. Gómez, J. Fernández-Rossier

Publié 2026-09-02

📖 7 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : D. Faílde, A. Gómez, J. Fernández-Rossier

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 microscopique des matériaux quantiques, les atomes s'organisent souvent en longues chaînes unidimensionnelles, se comportant comme de minus-aimants qui interagissent avec leurs voisins. Pendant des décennies, les physiciens ont étudié ces chaînes pour comprendre comment le magnétisme fonctionne au niveau le plus fondamental. Un puzzle clé dans ce domaine concerne un type spécifique de chaîne magnétique où les interactions sont parfaitement équilibrées, créant un état de désordre quantique. Dans un tel système, les excitations collectives — la façon dont la chaîne réagit lorsqu'elle est perturbée — ne se comportent pas comme les ondes sonores ou lumineuses familières. Au lieu de cela, la théorie suggère que ces excitations se fragmentent en morceaux plus petits et indépendants. Ces morceaux sont appelés spinons. Contrairement à une onde magnétique complète, un spinon ne porte que la moitié de la force magnétique d'un seul atome, une propriété qui semble impossible dans le monde classique mais qui est parfaitement normale dans le domaine quantique. Comprendre comment ces spinons se déplacent et transportent l'énergie est crucial pour saisir le comportement de matériaux exotiques, pourtant les observer directement dans une expérience physique est incroyablement difficile car ils n'existent pas en tant qu'objets isolés dans le monde réel.

Pour résoudre ce problème, une équipe de chercheurs s'est tournée vers la puissance émergente de l'informatique quantique. Plut lieu de construire une chaîne physique d'atomes, ils ont utilisé un processeur quantique numérique pour simuler le comportement de ces chaînes magnétiques et construire le schéma mathématique d'un spinon unique. Les chercheurs se sont concentrés sur deux modèles spécifiques de ces chaînes magnétiques : une où les atomes n'interagissent qu'avec leurs voisins immédiats, et une autre où les atomes interagissent avec chaque autre atome de la chaîne, quelle que soit la distance. Leur objectif était de comprendre comment préparer un état quantique qui représente un spinon unique, une tâche qui était auparavant limitée aux calculs théoriques ou aux simulations sur des ordinateurs classiques. En réussissant à créer cet état sur un ordinateur quantique, ils ont démontré une nouvelle façon d'explorer les propriétés de ces particules évasives, telles que leur énergie et leur impulsion, sans avoir besoin de les isoler dans un laboratoire physique.

Le processus qu'ils ont développé est une procédure en plusieurs étapes qui commence par la préparation de l'état fondamental de la chaîne magnétique. Cet état fondamental est la configuration d'énergie la plus basse que le système puisse contenir, un arrangement stable où les spins sont équilibrés. Une fois cet état stable établi sur l'ordinateur quantique, les chercheurs introduisent une unité magnétique supplémentaire, ajoutant de fait un atome de plus à la chaîne. Cet ajout perturbe l'équilibre parfait, créant les conditions nécessaires à l'existence d'un spinon. Cependant, ajouter simplement un atome ne suffit pas ; les chercheurs avaient besoin de créer un motif spécifique où cet excédent est réparti à travers la chaîne de manière à lui donner une impulsion définie. Pour y parvenir, ils ont utilisé une technique appelée combinaison linéaire d'unitaires. En termes simples, cela signifie qu'ils n'ont pas essayé de forcer l'ordinateur quantique à créer l'état final en une seule étape. Au lieu de cela, ils ont créé une superposition de nombreuses possibilités différentes, où l'excès de spin est situé à différents endroits le long de la chaîne, puis ont combiné ces possibilités d'une manière mathématique précise pour former l'état final du spinon.

L'équipe a testé cette méthode sur deux types différents de chaînes magnétiques. Pour la chaîne où les atomes ne communiquent qu'avec leurs voisins, ils ont utilisé une approche flexible et itérative connue sous le nom de solveur variationnel d'énergie propre (VQE) pour trouver l'état fondamental. Pour la chaîne avec des interactions à longue portée, ils ont utilisé une méthode directe et exacte basée sur une projection mathématique spécifique. Dans les deux cas, ils ont généré avec succès l'état de spinon unique et mesuré ses propriétés. Ils ont constaté que l'énergie du spinon change en fonction de son impulsion d'une manière qui correspond parfaitement aux prédictions des formules théoriques vieilles de plusieurs décennies. Cela a confirmé que leur construction numérique était exacte. Ils ont également mesuré une quantité appelée norme, qui nous indique essentiellement à quel point l'état est « réel » ou bien défini. Leurs résultats ont montré que l'état du spinon n'existe que dans une plage spécifique de valeurs d'impulsion, disparaissant en dehors de cette plage, ce qui concorde avec la compréhension théorique selon laquelle les spinons sont confinés à une région spécifique du paysage quantique.

Bien que la méthode directe de création de l'état fonctionne, elle nécessite un nombre important de bits quantiques supplémentaires, appelés qubits ancillaire, pour fonctionner, et le succès de l'opération est probabiliste, ce qui signifie qu'elle ne fonctionne pas à chaque fois. Pour remédier à cela, les chercheurs ont exploré des stratégies alternatives nécessitant moins de ressources. Une approche consistait à calculer l'énergie et les propriétés du spinon en décomposant le problème en morceaux plus petits et gérables, et en les mesurant individuellement, plutôt que d'essayer de créer l'état complet d'un coup. Cette méthode, connue sous le nom de test de Hadamard, utilise un seul qubit supplémentaire pour sonder le système. Bien que cette approche nécessite d'exécuter beaucoup plus de circuits séparés pour recueillir la même quantité de données, elle est beaucoup plus efficace en termes de matériel requis. Les chercheurs ont comparé ces différentes stratégies et ont trouvé que, si la méthode directe est puissante, l'approche de mesure étape par étape est mieux adaptée à la génération actuelle d'ordinateurs quantiques, qui possèdent des connexions limitées entre leurs qubits.

L'étude a été menée entièrement par des simulations informatiques, ce qui signifie que les chercheurs ont exécuté leurs algorithmes sur un ordinateur classique qui imite le comportement d'un processeur quantique. Ils n'ont pas réalisé ces expériences sur un véritable matériel quantique physique, mais les simulations étaient suffisamment rigoureuses pour tester la faisabilité des méthodes. Les résultats indiquent que les techniques proposées sont viables et que la fonction d'onde du spinon unique peut être préparée et analysée avec précision. Les chercheurs ont également étendu leurs travaux au modèle de Haldane-Shastry, un système plus complexe avec des interactions à longue portée, et ont constaté que les mêmes méthodes s'appliquaient avec succès. C'est une avancée significative car les outils mathématiques qui fonctionnent bien pour les chaînes simples à courte portée échouent souvent lorsqu'ils sont appliqués à des systèmes plus complexes à longue portée. En démontrant que ces algorithmes quantiques peuvent gérer les deux types de modèles, l'équipe a ouvert la voie à des études futures de matériaux magnétiques plus complexes.

En fin de compte, ce travail fournit une feuille de route pour étudier les particules quantiques complexes sur un ordinateur quantique. La capacité de préparer et de mesurer l'état d'un spinon unique est une étape fondamentale vers la compréhension de phénomènes plus compliqués, tels que la façon dont ces particules se comportent dans les matériaux bidimensionnels. Dans le monde réel, les techniques mathématiques qui fonctionnent pour les chaînes unidimensionnelles, comme l'ansatz de Bethe ou la renormalisation de la matrice densité, deviennent inefficaces ou échouent complètement lorsqu'elles sont appliquées à des systèmes bidimensionnels. La simulation quantique offre une solution potentielle à ce goulot d'étranglement. En démontrant que ces méthodes fonctionnent dans un environnement simulé, les chercheurs ont montré que les ordinateurs quantiques pourraient éventuellement devenir l'outil principal pour explorer les fonctions d'onde et la dynamique des spinons dans des dimensions supérieures, ouvrant ainsi une nouvelle fenêtre sur le comportement de la matière quantique qui est actuellement hors de portée des ordinateurs classiques.

Résumé Technique : Préparation de la fonction d'onde de spinon unique sur un ordinateur quantique

Énoncé du problème
L'article traite du défi que représente la préparation de fonctions d'onde de spinon unique sur des ordinateurs quantiques pour les modèles de spins unidimensionnels à $S=1/2$ . Bien que le modèle de Heisenberg avec des interactions antiferromagnétiques possède un spectre d'excitations de spin sans gap, celles-ci sont souvent décrites comme des particules composites appelées spinons ( $S=1/2$ ) plutôt que comme des magnons isolés. Un ansatz heuristique pour la fonction d'onde de spinon unique, proposé par Kulka et al. [1], construit un état de spin $S=1/2$ dans une chaîne de $L+1$ sites en ajoutant un site avec $S_z=1/2$ aux configurations de l'état fondamental d'une chaîne de $L$ sites. Cependant, l'implémentation de cet ansatz sur le matériel quantique est non triviale car l'état résultant n'est pas normalisé et ne peut être généré par un opérateur unitaire unique ; il nécessite une combinaison linéaire d'unitaires (LCU). De plus, bien que cet ansatz ait été validé pour le modèle XXZ, son applicabilité au modèle de Haldane-Shastry à interactions à longue portée n'avait pas été établie.

Méthodologie
Les auteurs proposent un protocole composé de trois étapes principales, illustrées par la Figure 1 :

Préparation de l'état fondamental : Préparation de l'état fondamental ( $|\Psi_{GS}\rangle$ $∣ Ψ_{GS} ⟩$ ) du modèle de spin cible (Heisenberg ou Haldane-Shastry) sur $L$ $L$ qubits.
- Pour le modèle de Heisenberg, ils utilisent l'algorithme VQE (Variational Quantum Eigensolver). Ils emploient un circuit quantique paramétré peu profond (PQC) initialisé avec un état de fermions libres (obtenu via des rotations de Givens), suivi de couches de portes $R_{ZZ}$ et $R_{XX+YY}$ .
- Pour le modèle de Haldane-Shastry, ils exploitent le fait que l'état fondamental est une mer de fermions libres projetée par Gutzwiller. Ils préparent l'état de fermions libres à l'aide de rotations de Givens et appliquent une projection de Gutzwiller via une séquence de portes CNOT. Cette étape est non déterministe, avec une probabilité de succès décroissant exponentiellement avec la taille du système.
Construction de l'état de spinon unique : Extension du réseau d'un site supplémentaire (ajout d'un qubit ancilla dans l'état $|\uparrow\rangle$ ) et application de l'ansatz de spinon. L'ansatz est défini comme une transformée de Fourier des configurations étendues : $|\Psi(q)\rangle = \frac{1}{\sqrt{L+1}} \sum_{m=0}^{L} e^{iqm} |\Psi(m)\rangle$ , où $|\Psi(m)\rangle$ représente l'état fondamental avec le spin supplémentaire inséré au site $m$ .
Estimation de l'énergie et de la norme : L'article évalue trois stratégies distinctes pour calculer la norme dépendante du moment $N(q)$ $N (q)$ et la dispersion d'énergie $\epsilon(q)$ $ϵ (q)$ :
- Combinaison Linéaire d'Unitaires (LCU) : Un algorithme non déterministe utilisant $L$ qubits ancillas pour préparer directement l'état de superposition $|\Psi(q)\rangle$ . La norme est estimée par la probabilité de mesurer le registre ancilla dans l'état $|0\rangle^{\otimes L}$ .
- Valeurs attendues de chaînes SWAP (SWAP strings) : Une approche déterministe qui calcule directement les recouvrements $\langle \Psi(n) | \Psi(m) \rangle$ sur l'état fondamental sans ancillas. Cependant, cela nécessite la mesure d'un nombre exponentiel de chaînes de Pauli ( $O(4^L)$ ).
- Test de Hadamard : Une méthode déterministe utilisant un seul qubit ancilla pour estimer les recouvrements et les amplitudes de transition terme par terme. Cela nécessite $O(L^2)$ circuits mais évite l'échelle exponentielle des chaînes de Pauli et le nombre élevé d'ancillas de la méthode LCU.

Principales contributions

Extension au modèle de Haldane-Shastry : Les auteurs démontrent que l'ansatz de spinon unique de Kulka et al., initialement développé pour le modèle XXZ, décrit avec précision les spinons dans le modèle de Haldane-Shastry. Ils vérifient que l'ansatz reproduit correctement la dispersion d'énergie analytique $\epsilon(q) = \frac{J}{2}(q_0^2 - q^2)$ et le domaine de moment $(-\pi/2 < q < \pi/2)$ où la norme est non nulle.
Stratégies d'implémentation quantique : Le papier fournit une feuille de route détaillée pour l'implémentation de l'ansatz de spinon unique sur le matériel quantique. Il compare les coûts de ressources (qubits, portes, profondeur de circuit, nombre de circuits) des méthodes LCU, SWAP-string et du test de Hadamard.
Validation In-Silico : Les méthodes proposées sont testées via des simulations numériques (« in-silico ») pour des tailles de système allant jusqu'à $L=16$ (Heisenberg) et $L=8$ (Haldane-Shastry). Les simulations confirment que les approches du test de Hadamard et de la LCU donnent des résultats cohérents avec la diagonalisation exacte.

Résultats

Préparation de l'état fondamental : L'approche VQE pour le modèle de Heisenberg atteint une infidélité faible pour des longueurs de chaîne moyennes ( $L \leq 8$ ), avec une performance qui se dégrade légèrement à mesure que $L$ augmente mais reste en dessous des seuils de bruit d'échantillonnage. Pour le modèle de Haldane-Shastry, la méthode de projection de Gutzwiller produit une fidélité élevée mais souffre d'une probabilité de succès décroissant exponentiellement ( $p \sim 2^{-\sqrt{2}L/2}$ ), limitant son application pratique à des longueurs de chaîne modérées.
Dispersion de spinon : La dispersion d'énergie calculée $\epsilon(q)$ et la norme $N(q)$ pour les deux modèles correspondent aux prédictions analytiques. La norme $N(q)$ s'annule correctement pour $|q| > \pi/2$ à mesure que $L$ augmente, confirmant le domaine de moment du spinon.
Comparaison des méthodes :
- La méthode LCU permet la préparation directe de l'état mais nécessite $L$ qubits ancillas et une étape de post-sélection non déterministe. Sa probabilité de succès suit une échelle linéaire avec $N(q)/(L+1)$ , ce qui la rend plus scalable que la projection de Gutzwiller pour l'état fondamental.
- Le test de Hadamard nécessite un seul qubit ancilla et est déterministe. Bien qu'il exige $O(L^2)$ exécutions de circuits, il évite l'échelle exponentielle des termes de Pauli trouvée dans la méthode directe de SWAP-string et le surcoût de qubits de la LCU. Les auteurs notent que pour le modèle de Haldane-Shastry, la décroissance exponentielle de la probabilité de succès de la préparation de l'état fondamental domine le budget d'erreur, nécessant une augmentation exponentielle du nombre de tirages (shots) pour maintenir la précision.

Signification et Revendications
Le papier affirme fournir un cadre pratique pour réaliser des fonctions d'onde de spinon unique dans des systèmes physiques quantiques, une étape nécessaire pour étudier les quasi-particules dans des régimes où les méthodes classiques (comme le Bethe Ansatz ou la DMRG) deviennent inefficaces, particulièrement dans des dimensions supérieures. Les auteurs précisent explicitement que leurs méthodes sont actuellement testées « in-silico » et ne revendiquent pas de réalisation expérimentale sur un matériel physique. La principale contribution est la proposition algorithmique et l'analyse des ressources pour préparer ces fonctions d'onde de spin spécifiques non triviales, servant de précurseur pour de futures implémentations dans des modèles bidimensionnels où les solutions exactes sont indisponibles. Le travail souligne les compromis entre l'estimation déterministe terme par terme (test de Hadamard) et la préparation d'état probabiliste (LCU) dans le contexte des contraintes du matériel quantique de l'ère NISQ.

Articles similaires