Résumé technique : Certification de l'optimisation quantique et de la découpe de circuits par l'utilisation de la dualité moment quantique-classique

1. Énoncé du problème

Le document traite de deux défis critiques de l'informatique quantique à court terme :

1. Manque de certification de performance : Les algorithmes quantiques variationnels (VQA) tels que le l'Algorithme d'Optimisation Quantique Approchée (QAOA) et la Méthode de Puissance Quantique Variationnelle (VQPM) manquent souvent de garanties de performance rigoureuses par état. Leurs ratios d'approximation sont fréquemment dépendants de l'instance ou reposent sur la convergence vers des optima globaux, laissant un vide dans la certification pour les exécutions bruitées ou interrompues prématurément.
2. Surcoût de la découpe de circuits (Circuit Cutting) : La décomposition de grands circuits quantiques en sous-circuits plus petits (découpe de circuit) pour les exécuter sur un matériel limité est un problème NP-complet lorsque l'on recherche des bipartitions optimales. Les heuristiques existantes reposent souvent sur des hypothèses géométriques ou des recherches par force brute, manquant de bornes d'erreur rigoureuses basées sur des données de corrélation mesurables.

Les auteurs proposent un cadre unifié exploitant une « dualité moment quantique-classique » pour résoudre ces deux problèmes en utilisant uniquement des données de corrélation Pauli-Z à deux points.

2. Méthodologie : Le dictionnaire de moments quantique-classique

La méthodologie centrale établit une correspondance directe entre les propriétés physiques de tout état quantique et l'espace de relaxation classique de Somme de Carrés (SoS).

2.1 La correspondance

Soit $\rho$ un état quantique arbitraire (matrice de densité) sur $n$ qubits. Les auteurs définissent une application $\Phi$ qui extrait les valeurs d'espérance Pauli-Z à un qubit et à deux qubits pour construire une matrice $Y_{\text{quantum}}$ :
$$Y_{\text{quantum}} = \Phi(\rho) = \begin{pmatrix} 1 & \langle Z_1 \rangle_\rho & \cdots & \langle Z_n \rangle_\rho \\ \langle Z_1 \rangle_\rho & 1 & \cdots & \langle Z_1 Z_n \rangle_\rho \\ \vdots & \vdots & \ddots & \vdots \\ \langle Z_n \rangle_\rho & \langle Z_n Z_1 \rangle_\rho & \cdots & 1 \end{pmatrix}$$
où $\langle O \rangle_\rho = \text{Tr}(\rho O)$.

2.2 Fondement théorique

Théorème 2.1 (Rétention de faisabilité rigoureuse) : Le document prouve que pour tout état quantique valide $\rho$, la matrice $Y_{\text{quantum}}$ résultante est toujours un point faisable du cône SoS de degré 2 classique (la relaxation SDP de Goemans–Williamson).

• Logique de preuve : Les éléments diagonaux sont égaux à 1 car $Z_i^2 = I$. La matrice est semi-définie positive (PSD) car pour tout vecteur $v$, $v^T Y_{\text{quantum}} v = \text{Tr}(\rho O_v^2)$, où $O_v$ est un observable hermitien. Puisque $\rho$ et $O_v^2$ sont PSD, leur produit de trace est non négatif.

Cette observation crée un « filet de sécurité universel » : l'état quantique lui-même agit comme un proxy d'optimisation actif au sein de la géométrie continue du cône semi-défini, ne nécessitant aucune résolution active de SDP.

3. Contributions clés

3.1 Certification universelle pour Max-Cut

En établissant que $Y_{\text{quantum}}$ est faisable pour le SDP SoS de degré 2, les auteurs appliquent directement l'arrondi par hyperplan aléatoire de Goemans–Williamson (GW) à la matrice dérivée du quantique.

• Résultat : Pour tout état $\rho$ produit par un algorithme variationnel, le poids de coupe attendu $E[\text{Cut}]$ obtenu via l'arrondi GW satisfait :
  $$E[\text{Cut}] \geq \alpha_{GW} \langle H \rangle_\rho$$
  où $\alpha_{GW} \approx 0,87856$ et $\langle H \rangle_\rho$ est la valeur d'espérance quantique de l'Hamiltonien Max-Cut.
• Signification : Cela fournit un taux d'approximation inférieur certifié pour chaque état de la trajectoire d'optimisation, découplant la certification de la performance de la convergence algorithmique. Cela reste valide même si l'algorithme est bruité, piégé dans des minima locaux ou arrêté prématurément.

3.2 Découpe de circuit basée sur la corrélation

La même matrice de moments capture la structure d'intrication de l'état, permettant une procédure en temps polynomial pour détecter les décompositions produit tensoriel de l'unitaire $U$.

• Détection exacte : Si un unitaire $U$ se factorise en $A_S \otimes B_T$, la corrélation connectée $C_{ij} = \langle Z_i Z_j \rangle - \langle Z_i \rangle \langle Z_j \rangle$ s'annule pour tout $i \in S, j \in T$. Les auteurs proposent de construire un graphe à partir des corrélations seuillées ; les composantes déconnectées indiquent des coupes exactes.
• Découpe approximative : Pour les états qui ne sont pas des états produits exacts, la méthode formule un problème de coupe minimale pondérée en utilisant $|C_{ij}|$ comme poids d'arête pour trouver une bipartition qui minimise les corrélations inter-partitions.
• Bornes d'erreur : Le document fournit des bornes d'erreur rigoureuses quantifiant la distance entre l'état original et la meilleure approximation d'état produit :
  • Sonde unique : La distance est bornée par $\sqrt{I(S:T)}$, où $I(S:T)$ est l'information mutuelle quantique.
  • Certification globale : En utilisant l'état de Choi, les auteurs dérivent une borne sur la distance de la norme d'opérateur $\|U - A_S \otimes B_T\|_\infty$, garantissant que la décomposition est valable pour des entrées arbitraires, et pas seulement pour l'état spécifique préparé.

4. Résultats et validation

Le cadre a été validé numériquement sur deux fronts :

1. Algorithmes variationnels (QAOA et VQPM) :

   • QAOA ($p=1$ sur un cycle de 4) : L'étude a démontré que la borne arrondie par GW $F(\gamma, \beta)$ reste strictement supérieure au plancher analytique $\alpha_{GW} E_\rho$ tout au long de la trajectoire d'optimisation, même lorsque l'énergie brute $E_\rho$ est sous-optimale.
   • VQPM : La méthode a été appliquée à la méthode de puissance quantique variationnelle, incluant son heuristique de « verrouillage de qubit » (qubit locking). Les résultats ont montré que le certificat de la matrice de moments reste valide même après un verrouillage agressif, fournissant un filet de sécurité contre la convergence prématurée. Le plancher inférieur certifié a suivi de manière cohérente la performance réelle, convergeant finalement vers la valeur exacte de Max-Cut lorsque l'état atteint une solution de base computationnelle.

2. Découpe de circuit (États aléatoires) :

   • En utilisant des états de Haar-aléatoires sur $n=6$ qubits, les auteurs ont testé la coupe basée sur la corrélation par rapport à une recherche exhaustive.
   • Précision : L'information mutuelle des coupes trouvées via la matrice de moments ($I_Y$) coïncidait presque parfaitement avec l'information mutuelle globalement optimale ($I_{\text{best}}$) trouvée par recherche exhaustive.
   • Bornes d'erreur : Les erreurs de norme 2 observées entre le véritable état et l'état produit reconstruit étaient systématiquement inférieures à la borne théorique $\sqrt{I(S:T)}$.
   • Taux de succès : Pour les états produits exacts, l'algorithme a identifié la structure de produit tensoriel correcte avec un succès de 100 % à travers diverses tailles de sous-systèmes.

5. Signification et revendications

Le document revendique une intuition structurelle fondamentale : la matrice de corrélation Pauli-Z à deux qubits de tout état quantique est automatiquement un point faisable du cône SoS de degré 2 classique.

Cette découverte a trois implications principales :

1. Certification passive : Contrairement aux travaux antérieurs qui résolvent activement des SDP pour trouver des bornes, cette méthode utilise simplement le dispositif quantique pour générer des données de mesure. La certification est une étape de post-traitement passive qui ne nécessite aucune ressource quantique supplémentaire au-delà de l'estimation de l'énergie standard.
2. Découplage de la convergence et de la certification : Elle fournit un « filet de sécurité » pour les algorithmes variationnels, garantissant un plancher de performance quel que soit l'état de convergence de l'optimiseur classique ou le niveau de bruit du matériel.
3. Découpe de circuit efficace : Elle offre une méthode en temps polynomial, agnostique à la géométrie, pour la découpe de circuits avec des bornes d'erreur rigoureuses dérivées directement de corrélations mesurables à deux points, contournant la complexité NP-complète de la recherche de coupes optimales.

Les auteurs positionnent ce travail comme un « filet de sécurité universel » qui relie l'optimisation quantique et la relaxation convexe classique, fournissant des garanties rigoureuses tant pour la qualité de la solution (Max-Cut) que pour la structure du circuit (découpe), sans nécessiter de mesures de corrélation d'ordre supérieur ou de tomographie complexe.