From barren plateaus through fertile valleys: Conic extensions of parameterised quantum circuits

1. Problématique : Le phénomène des plateaux stériles

L'article s'attaque à l'un des principaux obstacles des algorithmes quantiques variationnels (VQA) sur les dispositifs NISQ (Noisy Intermediate-Scale Quantum) : le phénomène des plateaux stériles (barren plateaus).

• Contexte : Les algorithmes comme QAOA (Quantum Approximate Optimization Algorithm) utilisent des circuits quantiques paramétrés (PQC) pour minimiser une fonction de coût.
• Le problème : Pour des circuits profonds (nombre de couches élevé), les gradients de la fonction de coût deviennent exponentiellement petits par rapport au nombre de qubits. Cela rend l'optimisation classique (par descente de gradient) inefficace, car le paysage d'optimisation devient plat, piégeant l'algorithme dans des minima locaux ou l'empêchant de progresser.
• Limitation actuelle : Les méthodes classiques de VQA sont contraintes à une variété unitaire. Si la direction optimale pour réduire l'énergie est orthogonale à la tangente de cette variété, l'algorithme ne peut pas y accéder.

2. Méthodologie : Extensions coniques et combinaisons linéaires d'unitaires

Les auteurs proposent une approche novatrice pour « sauter » hors des plateaux stériles en étendant la classe d'ansatz au-delà des opérations unitaires pures.

• Extension non-unitaire (Cône variationnel) : Au lieu de se limiter à des portes unitaires, l'algorithme introduit des portes non-unitaires $M_\alpha$ définies comme une combinaison linéaire d'unitaires (LCU) :
  $$M_\alpha = \sum_{i=1}^{\ell} \alpha_i U_i$$
  où les $U_i$ sont des portes unitaires fixes (extraites du PQC précédent) et $\alpha_i$ sont des paramètres complexes. Cela transforme l'espace de recherche en un cône variationnel, permettant d'explorer l'intérieur de l'espace de Hilbert plutôt que de se limiter à sa surface (comme le font les portes unitaires).

• Optimisation via Problème Généralisé de Valeurs Propres (GEP) :
  La recherche des paramètres optimaux $\alpha$ pour minimiser l'énergie est reformulée comme un problème de minimisation sous contrainte. Cela se réduit à la résolution d'un problème généralisé de valeurs propres (GEP) de faible dimension ($\ell \times \ell$) :
  $$H \alpha = \lambda E \alpha$$
  où $H$ et $E$ sont des matrices de moments estimées via des tests de Hadamard (mesures intermédiaires).

  • $E_{ij} = \langle \phi | U_i^\dagger U_j | \phi \rangle$
  • $H_{ij} = \langle \phi | U_i^\dagger H U_j | \phi \rangle$
    Cette approche permet de trouver une direction de « saut » optimale de manière efficace sur un ordinateur classique, indépendamment du nombre de qubits du système quantique.

• Implémentation NISQ :
  L'opération non-unitaire est réalisée de manière probabiliste sur le processeur quantique :

  1. Utilisation d'un registre d'ancillaire pour encoder les coefficients $\alpha$.
  2. Application d'une porte contrôlée combinant les $U_i$.
  3. Mesure projective sur l'ancillaire et post-sélection.
     La probabilité de succès est donnée par $p_{succ} = \|\alpha\|_1^{-2}$. Bien que non déterministe, cette probabilité reste raisonnable pour un petit nombre de termes dans la combinaison.

3. Contributions Clés

1. Cadre théorique (QCP) : Introduction du « Programmation Conique Quantique » (Quantum Conic Programming), un cadre général permettant d'étendre les ansatzs variationnels pour inclure des opérations non-unitaires.
2. Évasion des plateaux stériles : Démonstration que les sauts à travers le cône variationnel permettent d'accéder à des directions de gradient qui sont inaccessibles aux méthodes purement unitaires, contournant ainsi le problème de la disparition des gradients.
3. Réduction de complexité : Transformation d'un problème d'optimisation quantique complexe en un problème de valeurs propres classique de petite dimension, robuste au bruit de tir (shot noise).
4. Intégration dans QAOA : Développement d'un protocole hybride où des étapes LCU sont insérées itérativement entre des phases d'optimisation QAOA standard lorsque la convergence stagne.

4. Résultats Expérimentaux (Simulations)

Les auteurs ont testé leur méthode sur des instances du problème MAXCUT (de 4 à 22 qubits) sur des graphes 3-réguliers.

• Performance de l'approximation :
  • Le QAOA standard (sans LCU) stagne rapidement (après ~10 itérations) avec un ratio d'approximation d'environ 0,785, bien en dessous de la borne de Goemans-Williamson (GW) de 0,878.
  • L'ajout d'une seule étape LCU améliore le ratio à 0,842.
  • Après trois étapes LCU, le ratio atteint 0,904, surpassant significativement la borne GW et les meilleurs solveurs classiques pour ces instances.
• Robustesse : La méthode montre une capacité à échapper non seulement aux plateaux stériles, mais aussi aux minima locaux où les méthodes purement unitaires échouent.
• Probabilité de succès :
  • Il existe un compromis : chaque étape LCU réduit la probabilité de succès globale (environ 36,5 % après une étape, 14,1 % après trois étapes).
  • Cependant, la distribution des états finaux est fortement biaisée vers les solutions de haute qualité (ratio > 0,9), contrairement au QAOA standard qui produit des états avec un faible recouvrement sur les solutions optimales.
• Évolutivité : Contrairement aux attentes, les probabilités de succès ne diminuent pas avec la taille du problème (jusqu'à 22 qubits), suggérant que la méthode est prometteuse pour des problèmes de taille modérée.

5. Signification et Conclusion

Cet article propose une avancée conceptuelle majeure pour l'ère NISQ. Il démontre que l'introduction contrôlée d'opérations non-unitaires (via des mesures intermédiaires et des ancillas) peut surmonter les limitations fondamentales des algorithmes variationnels actuels.

• Impact : La méthode offre une voie pratique pour obtenir des solutions de haute qualité pour des problèmes d'optimisation combinatoire, même sur des dispositifs bruyants et de petite échelle, là où les algorithmes purement variationnels échouent.
• Perspectives : Bien que la probabilité de succès soit un défi (nécessitant potentiellement de nombreuses répétitions), les auteurs suggèrent que le contrôle de cette probabilité pourrait mener à des versions itératives où l'algorithme VQA de base deviendrait obsolète. Ce travail pose les bases d'un nouveau paradigme, la programmation conique quantique, qui pourrait être généralisé à d'autres classes d'algorithmes.