Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

Résumé technique : Sondage de l'état fondamental du modèle de Heisenberg antiferromagnétique sur le réseau Kagomé à l'aide d'un algorithme de l'Eigensolver Quantique Variationnel géométriquement informé

Énoncé du problème
L'article traite du défi que représente la recherche de l'état fondamental du modèle de Heisenberg antiferromagnétique sur un réseau de Kagomé, un système connu pour sa frustration géométrique et ses potentiels états de liquide de spin quantique. La détermination de l'état fondamental de tels Hamiltoniens $k$-locaux est un problème QMA-difficile, rendant la simulation classique intraitable pour les grands systèmes en raison de la croissance exponentielle de l'espace de Hilbert. Bien que les algorithmes de l'Eigensolver Quantique Variationnel (VQE) offrent une approche hybride quantique-classique adaptée aux dispositifs NISQ (Noisy Intermediate-Scale Quantum), ils se heurtent à des obstacles importants : le besoin de circuits profonds, la sensibilité au bruit (erreurs cohérentes et incohérentes), les plateaux stériles (barren plateaux) où les gradients s'annulent, et la difficulté de concevoir des ansatzes qui soient à la fois expressifs et efficaces pour le matériel. Plus précisément, les auteurs visent à caractériser les propriétés de l'état fondamental, telles que les corrélations spin-spin et les facteurs de structure statiques, sur un matériel quantique réel sans dépendre de réductions spécifiques au problème ou d'une atténuation d'erreur intensive pendant la boucle d'optimisation.

Méthodologie
Les auteurs proposent une approche VQE géométriquement informée, adaptée au réseau de Kagomé, en se concentrant sur deux géométries fondamentales : un triangle unique et une étoile de Kagomé (six triangles partageant des sommets).

1. Ansatz géométriquement informé :
   L'innovation centrale est un ansatz personnalisé, peu profond et efficace pour le matériel. Contrairement à l'ansatz de l'Hamiltonien Variationnel (HVA), qui peut être trop profond pour les dispositifs NISQ, cet ansatz est construit à l'aide d'un motif de rotations monoquantiques paramétrées ($R_y$) et de portes d'intrication restreintes à la connectivité native du matériel (évitant ainsi les portes SWAP).

   • Espace de paramètres euclidien : Une caractéristique critique de cette conception est que le tenseur métrique de Fubini-Study est prouvé analytiquement comme étant diagonal et constant. Cela garantit que l'espace de paramètres est naturellement euclidien, exempt de singularités et de sur-paramétrage. Par conséquent, le gradient naturel quantique est équivalent au gradient standard, simplifiant ainsi le paysage d'optimisation.
   • Compatibilité avec le matériel : La profondeur du circuit est minimée (une couche d'intrication), et le nombre de paramètres est réduit en alternant le placement des portes $R_y$.

2. Stratégie d'optimisation (AQNGD/AGD) :
   Les auteurs emploient un optimiseur de descente de gradient adaptative (AGD). Parce que l'ansatz produit un tenseur métrique de Fubini-Study diagonal et constant, l'AQNGD (Adaptive Quantum Natural Gradient Descent) se réduit à l'AGD. Le taux d'apprentissage est adapté dynamiquement à l'aide de la règle de retour en arrière d'Armijo. Cela permet de prendre de grandes étapes initiales pour accélérer la convergence tout en réduisant automatiquement la taille du pas à mesure que l'optimiseur approche du minimum, évitant ainsi la stagnation dans des géométries étroites. Cette méthode est comparée à la méthode SPSA (Simultaneous Perturbation Stochastic Approximation).

3. Atténuation d'erreur (Post-traitement) :
   L'atténuation de l'erreur est appliquée après l'optimisation VQE pour affiner les estimations d'énergie finales et les propriétés structurelles, plutôt que durant l'entraînement itératif.

   • Atténuation de l'erreur de lecture (REM) : Une approche par matrice de réponse partitionnée est utilisée pour inverser les erreurs SPAM (State Preparation and Measurement), ce qui est plus efficace que l'inversion de matrice complète pour les mesures creuses.
   • Extrapolation du niveau de bruit (ZNE) : Les auteurs appliquent la ZNE en utilisant une extrapolation quadratique (et une régression polynomiale bayésienne) pour extrapoler les résultats vers la limite de bruit nul. Ils notent que bien que la ZNE ne respecte pas strictement le principe variationnel de Rayleigh-Ritz, elle est utilisée ici pour améliorer la précision.

Résultats clés
L'étude a été menée sur du matériel quantique réel (IBM Oslo, IBM Kyoto, IBM Torino et IBM Algiers) et validée via des simulations d'état de vecteur.

• Convergence et optimisation :

  • Sur un système de triangle à 3 qubits, le VQE optimisé par AGD a convergé de manière fluide vers l'énergie de l'état fondamental ($E \approx -3$) avec une haute fidélité (99,89 % sur simulation) en seulement 19 itérations. Les paramètres finaux correspondaient étroitement aux valeurs théoriques exactes.
  • Sur une étoile de Kagomé à 12 qubits, l'AGD a démontré une convergence supérieure à la SPSA, nécessitant moins d'itérations pour atteindre l'énergie de l'état fondamental ($E \approx -18$). La fidélité finale sur simulation était de 96,37 %.
  • L'adaptation de la taille du pas via la règle d'Armijo a montré qu'elle permettait de naviguer efficacement dans le paysage de coût, en commenissant par de grands pas et en affinant la convergence à l'approche du minimum.

• Résilience du bruit des propriétés structurelles :

  • Une découverte significative est que les corrélations spin-spin et le facteur de structure statique ($S(q)$) sont résilients au bruit. Même sans atténuation d'erreur, les cartes expérimentales de $S(q)$ conservaient la signature qualitative de l'état fondamental (par exemple, la suppression à moment zéro indiquant l'absence d'ordre ferromagnétique).
  • Le coefficient de corrélation de Pearson entre le $S(q)$ expérimental bruité et le $S(q)$ théorique exact est resté extrêmement élevé (> 99,9 %).
  • Les expériences de simulation ont révélé une dépendance "en escalier" du facteur de structure statique par rapport aux niveaux de bruit : la structure de corrélation reste robuste jusqu'à un seuil de bruit critique, après quoi elle se dégrade rapidement, ressemblant à une transition de phase. Cela contraste avec l'énergie, qui décroît de manière monotone avec le bruit.

• Efficacité de l'atténuation d'erreur :

  • L'application de la REM et de la ZNE a considérablement amélioré les estimations d'énergie. Pour l'étoile de Kagomé, l'énergie brute était d'environ $-14,9$, tandis que la REM seule l'a améliorée à environ $-17,3$, et la ZNE (quadratique) l'a rapprochée de la valeur exacte de $-18$.
  • Les auteurs ont observé que la combinaison de la REM et de la ZNE entraînait parfois un dépassement de l'énergie fondamentale (violant la borne variationnelle), un comportement cohérent avec la nature non variationnelle de la ZNE. La ZNE seule offrait souvent la meilleure précision parmi les configurations testées.

Signification et affirmations
L'article affirme qu'un design d'ansatz géométriquement informé, qui garantit un espace de paramètres naturellement euclidien, permet un entraînement efficace du VQE sur le matériel NISQ actuel sans nécessiter de réductions complexes et spécifiques au problème.

• Robustesse : Ce travail démontre que des propriétés physiques significatives, spécifiquement les corrélations spin-spin et le facteur de structure statique, peuvent être extraites de matériel quantique bruité avec une précision suffisante pour caractériser l'état fondamental, même sans atténuation d'erreur durant la boucle d'optimisation.
• Efficacité de l'optimisation : L'utilisation de l'AGD avec des tailles de pas adaptatives s'est avérée plus efficace que la SPSA pour cet ansatz spécifique, atteignant une convergence plus rapide en termes d'itérations.
• Scalabilité : En évitant les circuits profonds et une atténuation d'erreur lourde pendant l'entraînement, cette approche offre une voie évolutive pour sonder les aimants quantiques frustrés. Les auteurs suggèrent que ce cadre peut être étendu à des réseaux de Kagomé plus grands et à d'autres systèmes quantiques frustrés.
• Modestie : Les auteurs reconnaissent que si leur méthode caractérise avec succès de petites cellules fondamentales (triangle et étoile), l'état fondamental général du réseau de Kagomé infini reste une question ouverte dans la littérature. Ils ne prétendent pas avoir résolu le problème du réseau infini, mais ont démontré une méthodologie viable pour sonder ces systèmes sur les dispositifs actuels. Ils notent également que si la ZNE améliore les estimations d'énergie, elle ne préserve pas le principe variationnel, et que la combinaison de la REM et de la ZNE nécessite une interprétation prudente.