Symmetry-guided quantum state preparation: Branched-Subspaces Adiabatic Preparation (B-SAP)

Résumé Technique : Préparation Adiabatique par Sous-espaces Ramifiés (B-SAP)

Énoncé du Problème
La préparation d'états quantiques est une tâche fondamentale pour la simulation de systèmes complexes à corps multiples en physique, en chimie et en science des données. Les approches actuelles font face à des compromis importants. Les algorithmes quantiques variationnels (VQA) reposent sur des ansatz paramétrés et une optimisation classique, qui souffrent souvent de « plateaux stériles » (gradients s'amont de façon exponentielle), d'une faible extensibilité due au besoin d'espaces de paramètres exponentiellement grands, et de problèmes de convergence liés aux minima locaux. Inversement, la préparation adiabatique (AP) fait évoluer un état initial sous l'influence d'un hamiltonien dépendant du temps pour atteindre un état cible. Bien que conceptuellement distincte, l'AP est très sensible aux écarts spectraux d'énergie ; si les niveaux d'énergie se croisent ou si les écarts se referment exponentiellement avec la taille du système, le théorème adiabatique échoue à garantir la préparation d'un état propre spécifique, entraînant souvent une superposition inconnue d'états dégénérés. Ces limitations entravent la préparation simultanée des états fondamentaux et excités, ce qui est crucial pour la compréhension des spectres optiques, des processus de transfert de charge et des intermédiaires de réaction.

Méthodologie : Préparation Adiabatique par Sous-espaces Ramifiés (B-SAP)
Les auteurs proposent le B-SAP, un algorithme hybride qui intègre les forces structurelles des VQA et de l'AP tout en atténuant leurs faiblesses respectives grâce à des structures de théorie des groupes et un post-traitement classique.

1. Initialisation Guidée par la Symétrie : Contrairement à l'AP conventionnelle, qui sélectionne généralement un hamiltonien initial ($H_0$) avec un groupe de symétrie minimal pour éviter la dégénérescence, le B-SAP sélectionne un $H_0$ dont le groupe de symétrie contient le groupe de symétrie de l'hamiltonien cible ($H_T$) comme un sous-groupe. Par conséquent, le spectre d'énergie initial présente une haute dégénérescence.
2. Exploration de Sous-espaces via des Portes Paramétrées : L'algorithme identifie un sous-espace dégénéré $S(B_n)$ de $H_0$ correspondant à un niveau d'énergie $E_n$. Au lieu d'utiliser un ansatz générique, le B-SAP construit un ensemble spécifique de portes unitaires paramétrées $\{G_j(\alpha_j)\}$ qui agissent comme des endomorphismes sur ce sous-espace. Ces portes sont dérivées des générateurs du groupe unitaire $U(d_n)$ (où $d_n$ est la dimension du sous-espace), réduits à l'espace coset $U(d_n)/U(d_{n-1})$ pour garantir que la mise en correspondance soit bijective pour un état de référence fixé.
   • Crucialement, le nombre de paramètres entraînables évolue selon $2d_n - 1$ (ou $2(d_n-1)$ si la phase globale est non pertinente). Puisque $d_n$ évolue de manière polynomiale avec la taille du système $L$ pour les modèles considérés, le nombre de paramètres évite l'évolution exponentielle qui mène aux plateaux stériles.
3. Évolution Adiabatique avec Ramification : Le système subit une évolution adiabatique de $H_0$ vers $H_T$. Parce que le groupe de symétrie de $H(s)$ reste constant ou diminue (les symétries sont brisées, pas introduites), les représentations irréductibles (représ) du groupe de symétrie restent inchangées ou se divisent (ramification). Cette structure de « ramification » empêche la fusion de niveaux d'énergie distincts qui cause les croisements de niveaux dans l'AP conventionnelle. Le protocole garantit que la condition d'écart spectral n'est requise qu'entre des sous-espaces d'énergie distincts, plutôt qu'entre chaque état propre individuel, relâchant ainsi les contraintes pour une préparation réussie.
4. Post-traitement Classique (MC-VQE) : Pour optimiser les paramètres variationnels sans souffrir de goulots d'étranglement d'optimisation, les auteurs emploient l'algorithme de l'Eigensolver Quantique Variationnel à Contraction Multi-états (MC-VQE).
   • Au lieu d'optimiser un seul état, l'algorithme prépare des états de recouvrement $|\Psi^p_{ml}\rangle$ (superpositions d'états de base) sur le registre quantique.
   • Il mesure les éléments de matrice de l'hamiltonien cible $H_T$ dans la base évoluée, construisant une matrice $d_n \times d_n$ $[H_T]^{(n)}_{ml}$.
   • Cette matrice est diagonalisée classiquement pour obtenir les valeurs propres et les coefficients des vecteurs propres cibles.
   • Les paramètres variationnels sont ensuite ajustés classiquement pour maximiser le recouvrement entre l'état préparé et le vecteur cible dérivé classiquement. Cela évite le besoin d'un circuit quantique lors de l'étape finale de minimisation et fournit un minimum global connu pour la fonction de perte.

Contributions Clés

• Cadre Hybride : Le papier introduit un nouvel algorithme qui contourne le besoin d'ansatz heuristiques (communs aux VQA) et les exigences strictes de non-croisement des écarts (communes à l'AP).
• Construction par la Théorie des Groupes : Il utilise la théorie des représentations du groupe de symétrie de l'hamiltonien pour construire un circuit paramétré qui explore efficacement un sous-espace dégénéré spécifique, garantissant que le nombre de paramètres évolue de manière polynomiale avec la taille du système.
• Robustesse aux Croisements de Niveaux : En concevant le chemin adiabatique de telle sorte que les symétries soient brisées plutôt qu'introduites, la méthode empêche la formation de dégénérescences accidentelles qui tourmentent typiquement l'évolution adiabatique, permettant la préparation des états excités.
• Contournement de l'Optimisation Classique : L'intégration du MC-VQE permet l'extraction d'informations sur l'état excité via la diagonalisation classique d'une petite matrice, réduisant la dépendance aux boucles d'optimisation quantique-classique difficiles.

Résultats
Les auteurs ont validé le B-SAP en utilisant le modèle Heisenberg XYZ unidimensionnel avec des conditions aux limites périodiques sur un simulateur quantique (Qiskit) pour des tailles de système allant jusqu'à $L=10$ qubits.

• Précision : L'algorithme a préparé avec succès l'état fondamental et les quatre premiers états excités sur une large gamme de constantes de couplage ($J_x, J_y, J_z$).
• Métriques d'Erreur : Pour l'état fondamental, l'erreur (infidélité) était négligeable par rapport au recouvrement d'un état aléatoire ($2^{-L}$) dans la plupart des régions de paramètres, même avec un nombre minimal d'étapes de Trotter ($L/2$). Pour les états excités, les erreurs restaient généralement inférieures à 10 % et souvent de l'ordre du pourcentage.
• Extensibilité : La profondeur de circuit requise évolue de manière polynomiale avec la taille du système $L$.
• Comparaison : Comparé à l'AP standard sur le même modèle et les mêmes paramètres :
  • L'AP standard a performé de manière comparable pour l'état fondamental mais a totalement échoué pour le premier état excité dans la plupart des régions de paramètres en raison des croisements de niveaux d'énergie.
  • Le B-SAP a préparé avec succès les états excités là où l'AP standard a échoué, démontant sa capacité à naviguer à travers les croisements de niveaux via la ramification des sous-espaces dégénérés.

Signification et Revendications
Le papier affirme que le B-SAP offre une voie prometteuse pour la préparation d'états quantiques sur les plateformes d'informatique quantique de l'ère NISQ (Near-term) en équilibrant précision, efficacité des ressources et robustesse.

• Applicabilité à Court Terme : L'extensibilité polynomiale de la profondeur et des paramètres du circuit rend l'algorithme adapté aux dispositifs dotés de temps de cohérence et de connectivité limités.
• Accès aux États Excités : La méthode comble une lacune critique dans les capacités actuelles en permettant la préparation fiable d'états excités d'ordre supérieur, qui sont essentiels pour déterminer les spectres optiques et la réactivité chimique, mais qui sont typiquement inaccessibles via les méthodes adiabatiques ou variationnelles standards en raison des croisements de niveaux et des difficultés d'optimisation.
• Perspectives : Les auteurs suggèrent que l'approche peut être étendue aux simulations quantiques de matériaux complexes, de la chimie quantique et des systèmes de théorie des champs, les travaux futurs se concentrant sur l'intégration de techniques de mitigation d'erreurs pour un déploiement sur matériel réel. Le papier ne prétend pas résoudre le problème général de toute préparation d'état quantique, mais positionne le B-SAP comme une solution spécifique et efficace pour les hamiltoniens à corps multiples possédant des symétries connues.